In This Package:

ADCalibData.py

Go to the documentation of this file.

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

'''

from UserTagging.UserTaggingAlg import UserTaggingAlg
from GaudiPython import SUCCESS, FAILURE
from GaudiPython import gbl

class ADCalibData(UserTaggingAlg):
    '''Add AD Calibration Info to UserData'''
        
    # ----------------------------------------   
    def initTagList(self):  
        self.addTag('Dummy', ''
            ).addData('ADCalibData', '/Event/UserData/AD/ADCalibData'
            ).addInt('NPESum')

    # ----------------------------------------
    def check(self, evt):
        ADTag = evt['/Event/UserTag/Detector/AD']
        if not ADTag:
            return False
        
        calibHdr = evt['/Event/CalibReadout/CalibReadoutHeader']
        if not calibHdr:
            self.warning('cannot find calibHdr')
            return
        self.SaveCalibData(calibHdr)            

        self.tagIt('Dummy')
        
    # ----------------------------------------
    def SaveCalibData(self, calibHdr):
        '''Save calibration info to UserData'''
                
        NPESum = 0.            
        readout = calibHdr.calibReadout()
        for channel in readout.channelReadout():
            sensorId = channel.pmtSensorId()
            pmtId = gbl.DayaBay.AdPmtSensor(sensorId.fullPackedData())
            if pmtId.ring() == 0: continue  # exclude 2'' PMTs
            for i in range(channel.size()):
                charge = channel.charge(i)
                NPESum += charge
        
        myData = self.tagSet.Dummy.dataSet.ADCalibData            
        myData.set("NPESum", NPESum)            
        # myData.dump()
        
#####  Job Configuration for nuwa.py ########################################

def configure(argv=[]):
    pass

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

---