In This Package:

ReadoutDetector.py

Go to the documentation of this file.

#!/usr/bin/env python

# Load DybPython
import math
from DybPython.DybPythonAlg import DybPythonAlg
from GaudiPython import SUCCESS, FAILURE
from GaudiPython import gbl, loaddict
from DybPython.Util import irange

TH1F = gbl.TH1F
TH2F = gbl.TH2F

loaddict("libCLHEPRflx")
loaddict("libHepMCRflx")
Detector = gbl.DayaBay.Detector
AdPmtSensor = gbl.DayaBay.AdPmtSensor
ServiceMode = gbl.ServiceMode
ReconStatus = gbl.ReconStatus


class ReadoutDetectorAlg(DybPythonAlg):
    "Algorithm to make ReadoutDetector file"
    detList = ['DayaBayAD1', 'DayaBayAD2','LingAoAD1','LingAoAD2',
                'FarAD1','FarAD2','FarAD3','FarAD4']   
   
    def findDetector(instance, detname, detlist = detList):
        for i in range(0,len(detlist)):
            if (detlist[i] == detname):
                instance.stats["/file0/readout/DetectorName"].Fill(detlist[i],scale)
    
    def __init__(self,name):
        DybPythonAlg.__init__(self,name)
        return

    def initialize(self):
        status = DybPythonAlg.initialize(self)
        if status.isFailure(): return status
        self.info("initializing")

        self.cableSvc = self.svc('ICableSvc','StaticCableSvc')
        if self.cableSvc == None:
            self.error("Failed to get StaticCableSvc")
            return FAILURE

        hist = TH1F("Events","The Readout Detector" ,16, 0.0, 16.0)
        hist.GetXaxis().SetTitle("Readout Detector")
        hist.GetYaxis().SetTitle("Events/Day")
        for i in range(0,len(self.detList)):
            hist.GetXaxis().SetBinLabel(2 * (i + 1),self.detList[i])
        self.stats["/file0/readout/DetectorName"] = hist
        
        return SUCCESS

    def execute(self):
        self.info("executing")
        evt = self.evtSvc()

        readoutHdr = evt["/Event/Readout/ReadoutHeader"]
        if readoutHdr == None:
            self.error("Failed to get ReadoutHeader")
            return FAILURE

        readout = readoutHdr.readout()
        if readout != None:
            detectorName = readout.detector().detName()
            self.findDetector(detectorName)
        
        return SUCCESS
        
    def finalize(self):
        self.info("finalizing")
        status = DybPythonAlg.finalize(self)
        return status


#####  Job Configuration for nuwa.py ########################################
scale = 1.0

def configure(argv = []):
    import sys,getopt
    opts,args = getopt.getopt(argv,"s:")
    global scale
    for opt,arg in opts:
        if opt == '-s':
            scale = float(arg)    


    from StatisticsSvc.StatisticsSvcConf import StatisticsSvc
    statsSvc = StatisticsSvc()
    statsSvc.Output ={"file0":"ReadoutDetector.root"}
    import DataSvc
    DataSvc.Configure()
    return

def run(app):
    '''
    Configure and add an algorithm to job
    '''
    app.ExtSvc += ["StaticCableSvc", "StatisticsSvc"]
    readoutDetectorAlg = ReadoutDetectorAlg("MyReadoutDetector")
    app.addAlgorithm(readoutDetectorAlg)
    pass

---