In This Package:

Dump.py

Go to the documentation of this file.

#!/usr/bin/env python
'''
Add an info dumper algorithm
'''

from DybPython.DybPythonAlg import DybPythonAlg, gbl, SUCCESS, FAILURE, irange, units

class SayWhat(DybPythonAlg):
    
    def __init__(self, name = None):
        DybPythonAlg.__init__(self,name)

        self.genPath = '/Event/Gen/GenHeader'
        self.readoutPath = '/Event/Readout/ReadoutHeader'

        return

    def dump_genheader(self,gh):
        if not gh: 
            self.info("GenHeader does not exist")
            return


        evt = gh.event()

        self.info("GenHeader(%s)#%d @ %s" % (gh.generatorName(), evt.event_number(),
                                             gh.context().AsString()))
        vtx = evt.signal_process_vertex()
        for part in irange(vtx.particles_out_const_begin(), vtx.particles_out_const_end()):
            self.info('\t%d (%d) %f MeV' % (part.pdg_id(), part.status(), 
                                            part.generatedMass()/units.MeV) )
            continue
        return

    def dump_roheader(self,rh):
        if not rh: 
            self.info("ReadoutHeader does not exist")
            return

        self.info("ReadoutHeader @ %s" % rh.context().AsString())

        ro = rh.readout()
        if not ro:
            self.info("No Readout object in the ReadoutHeader")
        else:
            self.info("Readout: %s #%d @ %s %d" % (ro.detector().asString(),
                                                   ro.triggerNumber(), 
                                                   ro.triggerTime().AsString(),
                                                   ro.triggerType()) )
            pass
            
        dc = rh.daqCrate()
        if not dc:
            self.info("No DaqCrate object in the ReadoutHeader")
        else:
            self.info("DaqCrate: %s #%d[%d] @ %s %d" % (dc.detector().asString(),
                                                        dc.eventNumber(), 
                                                        dc.localTriggerNumber(),
                                                        dc.triggerTime().AsString(), 
                                                        dc.triggerType() ))
            pass
        return

    def check_roheader(self,roh):
        'Return True if okay'
        dt = roh.latest().GetSeconds() - roh.earliest().GetSeconds()
        if dt <= 1e-9:
            #self.warning('ReadoutHeader at "%s" has too small dt: %f seconds' % \
            #                 (roh.context().AsString(),dt))
            return False
        return True

    def execute(self):
        tes = self.evtSvc()

        self.info("Dumping TES")
        self.dump_genheader(tes[self.genPath])
        roh = tes[self.readoutPath]
        self.dump_roheader(roh)
        if not self.check_roheader(roh):
            pass
            #return FAILURE

        self.info("Dumping AES")
        ghlist = self.getAES(self.genPath) or []
        for gh in ghlist:
            self.dump_genheader(gh)
        rhlist = self.getAES(self.readoutPath) or []
        for rh in rhlist:
            self.dump_roheader(rh)

        return SUCCESS

def configure(argv=None):
    return

def run(app):
    sw = SayWhat()
    app.addAlgorithm(sw)

---