In This Package:

ADReadoutData.py

Go to the documentation of this file.

#!/usr/bin/env python
'''
Usage:
'''

from UserTagging.UserTaggingAlg import UserTaggingAlg
from UserTagging.Models import Tag, Data, Int, Float
from GaudiPython import SUCCESS, FAILURE
from GaudiPython import gbl

class ADReadoutData(UserTaggingAlg):
    '''Add AD Readout (DaqFormat) Info to UserData'''
        
    # ----------------------------------------    
    def initTagList(self):
        
        # Chao recommends method chaning since it's less typing and he is lazy.
        # A tag with no path (default to '') means no tag HeaderObject will be created
        # but this dummy tag is still need in order to manage the UserData 
        #
        self.addTag('Dummy'
            ).addData('ADReadoutData', '/Event/UserData/AD/ADReadoutData'
            ).addInt('nChannels', 'nHits'
            ).addFloat('chargeSum')
        
        # or use Brett's preferred pythonic syntax:
        #
        # self.tagSet.Dummy = Tag()
        # self.tagSet.Dummy.dataSet.ADReadoutData = Data('/Event/UserData/AD/ADReadoutData')
        # paraSet = self.tagSet.Dummy.dataSet.ADReadoutData.paraSet
        # paraSet.nChannels = Int()
        # paraSet.nHits = Int()
        # paraSet.chargeSum = Float()
        
        # you can do other necessary initialization here
        self.cableSvc = self.svc('ICableSvc','StaticCableSvc')
        if self.cableSvc == None:
            self.error("Failed to get StaticCableSvc")
            return FAILURE

    # ----------------------------------------
    def check(self, evt):

        # check if this is an 'AD' event
        ADTag = evt['/Event/UserTag/Detector/AD']
        if not ADTag:
            return False
        
        # save some useful user data
        readoutHdr = evt['/Event/Readout/ReadoutHeader']
        if not readoutHdr:
            self.warning('cannot find readoutHdr')
            return
        self.SaveReadoutData(readoutHdr)            
        
        # tag the event so the user data can be saved
        self.tagIt('Dummy')
        
        # or use Brett's favorate syntax
        # self.tagSet.Dummy.tagIt()
        
    # ----------------------------------------
    def SaveReadoutData(self, readoutHdr):
        '''Save Readout (DaqFormat) Info to UserData'''
                
        nChannels = 0
        nHits = 0
        chargeSum = 0.
        
        # calculate the variables ...    
        svcMode = gbl.ServiceMode(readoutHdr.context(), 0)     
        daqCrate = readoutHdr.daqCrate()
        pmtCrate = daqCrate.asPmtCrate()
       
        for channel in pmtCrate.channelReadouts():
            # channel = channelPair.second
            channelId = channel.channelId()
            board = channelId.board()
            connector = channelId.connector()
            pmtId = self.cableSvc.adPmtSensor(channelId, svcMode)
            ring = pmtId.ring()
            if ring == 0: continue   # exclude 2'' PMT
            
            nChannels += 1
            nHits += channel.hitCount()
            for i in range(channel.hitCount()):
                adc = channel.adc(i)
                preAdc = channel.preAdcAvg(i)

                if channel.isHighGainAdc(i): 
                    adc = adc - preAdc
                else:
                    adc = (adc - preAdc) * 20
                
                chargeSum += adc

        # ...
        
        # save the variables to the User Data
        myData = self.getTag('Dummy').getData('ADReadoutData')
        myData.set("nChannels", nChannels       
            ).set("nHits", nHits
            ).set("chargeSum", chargeSum)
                    
        # or use Brett's preferred pythonic syntax
        #
        # myData = self.tagSet.Dummy.dataSet.ADReadoutData
        # paraSet = myData.paraSet
        # paraSet.nChannels = nChannels
        # paraSet.nHits = nHits
        # paraSet.chargeSum = chargeSum
        
        # you can print all branches of UserData
        # myData.dump()
        # or an individual branch
        # print 'nChannels', myData.get('nChannels')
                       
#####  Job Configuration for nuwa.py ########################################

def configure(argv=[]):
    pass

def run(app):
    app.ExtSvc += ["StaticCableSvc"]
    myAlg = ADReadoutData("UserData::ADReadoutData")
    app.addAlgorithm(myAlg)
    pass                    
            

---