ExamplePmtId.py

Go to the documentation of this file.

#!/usr/bin/env python
#
# Template for writing your own algorithm
#
#  Usage:
#   nuwa.py -n 1 DivingIn.ExamplePmtId simulationData.root

# Load DybPython
from DybPython.DybPythonAlg import DybPythonAlg
from GaudiPython import SUCCESS, FAILURE
from GaudiPython import gbl

ServiceMode = gbl.ServiceMode

# Make your algorithm
class ExampleAlg(DybPythonAlg):
    "Example Python Algorithm"
    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','CableSvc')
        self.pmtGeomSvc = self.svc('IPmtGeomInfoSvc','PmtGeomInfoSvc')

        return SUCCESS

    def execute(self):
        self.info("executing")

        # Get the Triggered Readout of the detector        
        evt = self.evtSvc()
        readoutHdr = evt["/Event/Readout/ReadoutHeader"]
        readout = readoutHdr.readout()
        if readout == None:
            self.info("No Triggered Readout for this event")
            return SUCCESS

        # Get the service mode
        svcMode = ServiceMode( readoutHdr.context(), 0 )

        # Look at data from each channel
        for channelPair in readout.channelReadout():
            channel = channelPair.second

            # See #319
            adc = channel.adc()

            channelId = channel.channelId()
            print "==================================="
            print "VME board= ", channelId.board()
            print "connector= ", channelId.connector()
            
            # Find the PMT ID connected this electronics channel
            pmtId = self.cableSvc.adPmtSensor( channelId, svcMode )
            print "pmt ring=", pmtId.ring()
            print "pmt column=", pmtId.column()
            
            # Get the physical PMT number for this PMT,
            # to compare with the PMT Test bench
            pmtHardwareId = self.cableSvc.pmtHardwareId( pmtId, svcMode )
            
            # Get the PMT position
            pmtGeomInfo = self.pmtGeomSvc.get( pmtId.fullPackedData() )
            pmtPosition = pmtGeomInfo.localPosition()
            # position relative to AD center
            print "pmt x= ", pmtPosition.x()
            print "pmt y= ", pmtPosition.y()
            print "pmt z= ", pmtPosition.z()
            print "=================================="

        return SUCCESS
        
    def finalize(self):
        self.info("finalizing")
        status = DybPythonAlg.finalize(self)
        return status


#####  Job Configuration for nuwa.py ########################################

def configure():
    from DetHelpers.DetHelpersConf import PmtGeomInfoSvc
    pgiSvc = PmtGeomInfoSvc()
    pgiSvc.StreamItems = ["/dd/Structure/DayaBay"]
    import DataSvc
    DataSvc.Configure()
    return

def run(app):
    '''
    Configure and add an algorithm to job
    '''
    app.ExtSvc += ["PmtGeomInfoSvc"]
    example = ExampleAlg("MyExample")
    app.addAlgorithm(example)
    pass