Public Member Functions | |
def | __init__ |
def | initialize |
def | execute |
def | finalize |
Public Attributes | |
cableSvc |
Definition at line 28 of file Ge68EnergyStats.py.
def Ge68EnergyStats::EnergyStatsAlg::__init__ | ( | self, | ||
name | ||||
) |
Definition at line 30 of file Ge68EnergyStats.py.
00030 : 00031 DybPythonAlg.__init__(self,name) 00032 return 00033 def initialize(self):
def Ge68EnergyStats::EnergyStatsAlg::initialize | ( | self | ) |
Definition at line 34 of file Ge68EnergyStats.py.
00034 : 00035 status = DybPythonAlg.initialize(self) 00036 if status.isFailure(): return status 00037 self.info("initializing") 00038 00039 self.cableSvc = self.svc('ICableSvc','StaticCableSvc') 00040 if self.cableSvc == None: 00041 self.error("Failed to get StaticCableSvc") 00042 return FAILURE 00043 00044 # Make energy histograms 00045 00046 hist = TH1F("EnergyBias","Energy Bias",500,-2.0,2.5) 00047 hist.GetXaxis().SetTitle("Energy Difference [MeV]") 00048 hist.GetYaxis().SetTitle("Generated Particles") 00049 hist.SetLineColor(4) 00050 self.stats["/file0/energy/EnergyBias"] = hist 00051 00052 hist = TH1F("reconEnergy","Reconstructed Energy",500,0.0,3.0) 00053 hist.GetXaxis().SetTitle("Reconstructed Visible Energy [E_rec]") 00054 hist.GetYaxis().SetTitle("Triggered Readouts") 00055 hist.SetLineColor(4) 00056 self.stats["/file0/energy/reconEnergy"] = hist 00057 return SUCCESS 00058 def execute(self):
def Ge68EnergyStats::EnergyStatsAlg::execute | ( | self | ) |
Definition at line 59 of file Ge68EnergyStats.py.
00059 : 00060 self.info("executing") 00061 00062 evt = self.evtSvc() 00063 00064 00065 # reconstructed energy 00066 recHdr = evt["/Event/Rec/RecHeader"] 00067 00068 00069 roHdr = evt["/Event/Readout/ReadoutHeader"] 00070 readout = roHdr.readout() 00071 if readout != None: 00072 trigTime = readout.triggerTime() 00073 print "seconds=", trigTime.GetSec() 00074 print "nanoseconds=", trigTime.GetNanoSec() 00075 00076 00077 if recHdr == None: 00078 self.error("Failed to get RecHeader") 00079 return FAILURE 00080 recResults = recHdr.recResults() 00081 recTrigger = recResults["AdSimple"] 00082 if recTrigger.energyStatus() == ReconStatus.kGood: 00083 reconEnergy = recTrigger.energy() 00084 self.stats["/file0/energy/reconEnergy"].Fill(reconEnergy) 00085 self.stats["/file0/energy/EnergyBias"].Fill(reconEnergy-1.022) 00086 return SUCCESS 00087 def finalize(self):
def Ge68EnergyStats::EnergyStatsAlg::finalize | ( | self | ) |
Definition at line 88 of file Ge68EnergyStats.py.
00088 : 00089 self.info("finalizing") 00090 00091 status = DybPythonAlg.finalize(self) 00092 return status 00093 00094 ##### Job Configuration for nuwa.py ########################################
Definition at line 39 of file Ge68EnergyStats.py.