def say(self, det=None):
"""print header/detector information"""
if det is None:
self.sayHeader()
else:
binning = self.detector[det]
say("Bin : ", binning.nb)
say("Title : ", binning.name)
say("Type : ", binning.type)
say("Score : ", binning.score)
say("X : [", binning.xlow, "-", binning.xhigh, "] x",
binning.nx, "dx=", binning.dx)
say("Y : [", binning.ylow, "-", binning.yhigh, "] x",
binning.ny, "dy=", binning.dy)
say("Z : [", binning.zlow, "-", binning.zhigh, "] x",
binning.nz, "dz=", binning.dz)
say("L : ", binning.lntzer)
say("bk : ", binning.bk)
say("b2 : ", binning.b2)
say("tc : ", binning.tc)
def say(self, det=None):
"""print header/detector information"""
if det is None:
self.sayHeader()
else:
binning = self.detector[det]
say("Bin : ", binning.nb)
say("Title : ", binning.name)
say("Type : ", binning.type)
say("Region : ", binning.region)
say("Volume : ", binning.volume)
say("Mhigh : ", binning.mhigh)
say("Zhigh : ", binning.zhigh)
say("NMZmin : ", binning.nmzmin)
def say(self, det=None):
"""print header/detector information"""
if det is None:
self.sayHeader()
else:
det = self.detector[det]
say("BDX : ", det.nb)
say("Title : ", det.name)
say("Type : ", det.type)
say("Dist : ", det.dist)
say("Reg1 : ", det.reg1)
say("Reg2 : ", det.reg2)
say("Area : ", det.area)
say("2way : ", det.twoway)
say("Fluence: ", det.fluence)
say("LowNeu : ", det.lowneu)
say("Energy : [", det.elow, "..", det.ehigh, "] ne=", det.ne,
"de=", det.de)
if det.lowneu:
say("LOWNeut : [", det.egroup[-1], "..", det.egroup[0],
"] ne=", det.ngroup)
say("Angle : [", det.alow, "..", det.ahigh, "] na=", det.na,
"da=", det.da)
say("Total : ", det.total, "+/-", det.totalerror)
def say(self, det=None):
"""print header/detector information"""
if det is None:
self.sayHeader()
else:
bin = self.detector[det]
say("Bin : ", bin.nb)
say("Title : ", bin.name)
say("Type : ", bin.type)
say("Score : ", bin.score)
say("X : [", bin.xlow,"-",bin.xhigh,"] x", bin.nx, "dx=",bin.dx)
say("Y : [", bin.ylow,"-",bin.yhigh,"] x", bin.ny, "dy=",bin.dy)
say("Z : [", bin.zlow,"-",bin.zhigh,"] x", bin.nz, "dz=",bin.dz)
say("L : ", bin.lntzer)
say("bk : ", bin.bk)
say("b2 : ", bin.b2)
say("tc : ", bin.tc)
def sayHeader(self):
say("File : ", self.file)
say("Title : ", self.title)
say("Time : ", self.time)
say("Weight : ", self.weight)
say("NCase : ", self.ncase)
say("NBatch : ", self.nbatch)
def say(self, det=None):
"""print header/detector information"""
if det is None:
self.sayHeader()
else:
det = self.detector[det]
say("BDX : ", det.nb)
say("Title : ", det.name)
say("Type : ", det.type)
say("Dist : ", det.dist)
say("Reg1 : ", det.reg1)
say("Reg2 : ", det.reg2)
say("Area : ", det.area)
say("2way : ", det.twoway)
say("Fluence: ", det.fluence)
say("LowNeu : ", det.lowneu)
say("Energy : [", det.elow,"..",det.ehigh,"] ne=", det.ne, "de=",det.de)
if det.lowneu:
say("LOWNeut : [",det.egroup[-1],"..",det.egroup[0],"] ne=",det.ngroup)
say("Angle : [", det.alow,"..",det.ahigh,"] na=", det.na, "da=",det.da)
say("Total : ", det.total, "+/-", det.totalerror)
def say(self, det=None):
"""print header/detector information"""
if det is None:
self.sayHeader()
else:
bin = self.detector[det]
say("Bin : ", bin.nb)
say("Title : ", bin.name)
say("Type : ", bin.type)
say("Region : ", bin.region)
say("Volume : ", bin.volume)
say("Mhigh : ", bin.mhigh)
say("Zhigh : ", bin.zhigh)
say("NMZmin : ", bin.nmzmin)
def sayHeader(self):
say("File : ",self.file)
say("Title : ",self.title)
say("Time : ",self.time)
say("Weight : ",self.weight)
say("NCase : ",self.ncase)
say("NBatch : ",self.nbatch)
# lastind = ind
# blockname = ""
half = 0
f.close()
# if first:
# self.det.set(first)
#===============================================================================
if __name__ == "__main__":
import sys
# say("="*80)
# mgdraw = Mgdraw("examples/source001_source")
# while mgdraw.readEvent():
# say(mgdraw.data)
say("="*80)
usr = Usrbdx(sys.argv[1])
usr.say()
for i in range(len(usr.detector)):
say("-"*50)
usr.say(i)
data = unpackArray(usr.readData(i))
stat = unpackArray(usr.readStat(i))
#say( len(data), len(stat))
#for j,(d,e) in enumerate(zip(data,stat)):
# say(j,d,e)
#say()
#usr = Resnuclei(sys.argv[1])
#usr.say()
#for i in range(len(usr.detector)):
half = 0
f.close()
# if first:
# self.det.set(first)
#===============================================================================
if __name__ == "__main__":
import sys
# say("="*80)
# mgdraw = Mgdraw("examples/source001_source")
# while mgdraw.readEvent():
# say(mgdraw.data)
say("=" * 80)
usr = Usrbdx(sys.argv[1])
usr.say()
for i in range(len(usr.detector)):
say("-" * 50)
usr.say(i)
data = unpackArray(usr.readData(i))
stat = unpackArray(usr.readStat(i))
#say( len(data), len(stat))
#for j,(d,e) in enumerate(zip(data,stat)):
# say(j,d,e)
#say()
#usr = Resnuclei(sys.argv[1])
#usr.say()
#for i in range(len(usr.detector)):