CheckAlg.py

Go to the documentation of this file.

#!/usr/bin/env python

'''Algorithm to test RPC readout simulation

It can be called via nuwa.py:

  nuwa.py -n XXX -m 'ReadoutSim.CheckAlg ["my optional message"]'

'''


from DybPython.DybPythonAlg import DybPythonAlg
from GaudiPython import SUCCESS, FAILURE
#import PyCintex
#ElecFecCrate = PyCintex.makeClass('DayaBay::ElecFecCrate')

class CheckAlg(DybPythonAlg):
    '''
    Algorithm to test RPC readout simulation
    '''

    def __init__(self,myname):
        ''' '''
        DybPythonAlg.__init__(self,myname)
        return

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

        #self.stats.defaultPath = "/file1/examples/"
        ## Make the trigger time histogram
        #self.stats["TrigTime"] = TH1F("TrigTime","Trigger Time [s]",
                                      #100,0.0,10.0)
        ## Make the TDC histogram
        #self.stats["Tdc"] = TH1F("Tdc","TDC Values",300,0,300)

        return SUCCESS

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

        # Elec header info dump
        elhdr = evt["/Event/Elec/ElecHeader"]
        self.info("===== ElecSim =====")
        if elhdr == None:
            self.info("No ElecSim header this cycle")
        else :
            # pulses:
            pulhdr = elhdr.pulseHeader()
            # signals:
            crhdr = elhdr.crateHeader()

            self.info( "Elec time: %s"%(elhdr.timeStamp().GetSeconds()) )
            if pulhdr == None:
                self.info("No pulse header stored this cycle")
            else :
                self.info( "Number of pulse collections: %d"%(pulhdr.pulseCollection().size()) )

            if crhdr == None:
                self.info("No crate header this cycle")
            else :
                self.info( "Number of crates: %d"%(crhdr.crates().size()) )
                det = crhdr.crates().begin().first
                feccr = crhdr.crates()[det]
                self.info("%d FEC boards stored in the crate."%feccr.fecBoards().size())

        # SimTrig header info dump
        trighdr = evt["/Event/SimTrig/SimTrigHeader"]
        self.info("===== TrigSim =====")
        if trighdr == None:
            self.info("No SimElec header this cycle")
        else :
            cmdhdr = trighdr.commandHeader()
            cmdcoll= cmdhdr.collections()[cmdhdr.collections().begin().first]

            if cmdcoll == None:
                self.info("No trigger command collection")
            else :
                self.info("SimTrig: Detector: %s"%(cmdcoll.detector().detName()))
                self.info( "SimTrig: %d commands in TriggerCommandCollection."%(cmdcoll.commands().size()) )
                for cmd in cmdcoll.commands() :
                    self.info("%d input triggers in this command."%cmd.inTriggers().size())
                    self.info("%d output triggers in this command."%cmd.outTriggers().size())

        # SimReadout header info dump
        self.info("===== SimReadout =====")
        simrohdr = evt["/Event/SimReadout/SimReadoutHeader"]
        if simrohdr == None:
            self.info("No SimReadoutHeader header this cycle")
        else :
            simros = simrohdr.readouts()
            self.info("SimReadout: Got %d simreadouts."%simros.size())
            for simro in simros :
                ro = simro.readout()
                if ro == None:
                    self.info("No Readout this cycle")
                else :
                    self.info("Readout: Detector: %s"%(ro.detector().detName()))
                    self.info("Readout: Trigger time: %f"%(ro.triggerTime().GetSeconds()))
                    self.info("Readout: TriggerType: 0x%x"%ro.triggerType())
                    if (ro.detector().isRPC()) :
                        self.info("Readout: RPC readouts: %d"%(ro.fecReadouts().size()))
                        self.info("Readout: Hitmaps: ")
                        i = 0
                        for fec in ro.fecReadouts() :
                            hitmap = fec.second.hitMap()
                            print "%d: 0x%0x"%(i,hitmap)
                            i += 1
                simdaq = ro.daqCrate(1,1,1)
                self.info("DaqCrate: Detector: %s"%(simdaq.detector().detName()))
                self.info("DaqCrate: Trigger time: %f"%(simdaq.triggerTime().GetSeconds()))
                self.info("DaqCrate: TriggerType: 0x%x"%simdaq.triggerType())
                if (simdaq.detector().isRPC()) :
                    self.info( "DaqCrate: RPC readouts: %d"%(simdaq.rpcPanelReadouts().size()) )

        # Readout header info dump
        self.info("===== Readout =====")
        rohdr = evt["/Event/Readout/ReadoutHeader"]
        if rohdr == None:
            self.info("No ReadoutHeader header this cycle")
        else :
            ro = rohdr.readout()
            if ro == None:
                self.info("No Readout this cycle")
            else :
                self.info("Readout: Detector: %s"%(ro.detector().detName()))
                self.info("Readout: Trigger time: %f"%(ro.triggerTime().GetSeconds()))
                self.info("Readout: TriggerType: %x"%ro.triggerType())
                if (ro.detector().isRPC()) :
                    self.info("Readout: RPC readouts: %d"%(ro.fecReadouts().size()))
                    self.info("Readout: Hitmaps: ")
                    i = 0
                    for fec in ro.fecReadouts() :
                        hitmap = fec.second.hitMap()
                        print "%d: 0x%0x"%(i,hitmap)
                        i += 1
            daq = rohdr.daqCrate()
            if daq == None:
                self.info("No daq crate this cycle")
            else :
                self.info("DaqCrate: Detector: %s"%(daq.detector().detName()))
                self.info("DaqCrate: Trigger time: %f"%(daq.triggerTime().GetSeconds()))
                self.info("DaqCrate: TriggerType: 0x%x"%daq.triggerType())
                if (daq.detector().isRPC()) :
                    self.info( "DaqCrate: RPC readouts: %d"%(daq.rpcPanelReadouts().size()) )
                    self.info("DaqCrate: Hitmaps:")
                    for i in range(daq.rpcPanelReadouts().size()) :
                        hitmap = daq.rpcPanelReadouts()[i].hitMap()
                        print "%d: 0x%0x"%(i,hitmap)
                    
        return SUCCESS
    pass



### nuwa.py hooks:

def configure(argv = []):
    pass

def run(app):
    alg = CheckAlg("RpcCheckAlg")
    app.addAlgorithm(alg)
    return