In This Package:

CoincData.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
import GaudiKernel.SystemOfUnits as units

class CoincData(UserTaggingAlg):
    '''Add Coincidence info to UserData'''
        
    # ===========================================    
    def initTagList(self):
        self.addTag('Dummy', '')
        self.getTag('Dummy').addData('CoincData'  , '/Event/UserData/Coinc/CoincData')
        self.getTag('Dummy').addData('PromptData' , '/Event/UserData/Coinc/PromptData')
        self.getTag('Dummy').addData('DelayedData', '/Event/UserData/Coinc/DelayedData')
        
        myData = self.getTag('Dummy').getData('CoincData')
        myData.addInt('multiplicity')
        
        for name in ['PromptData', 'DelayedData']:
            myData = self.getTag('Dummy').getData(name)
            myData.addIntArray('t_sec', 't_nsec', 'triggerType', 'detID'
                 ).addFloatArray('energy', 'x', 'y', 'z')
        
    # ===========================================
    def check(self, evt):
        coincHdr = evt['/Event/UserTag/Coinc/All']
        if not coincHdr:  # not a coincidence event
            return False
        promptHdr  = evt['/Event/UserTag/Coinc/Prompt']
        if not promptHdr: 
            self.warning('missing prompt header')
            return
        delayedHdr = evt['/Event/UserTag/Coinc/Delayed']
        if not delayedHdr: 
            self.warning('missing delayed header')
            return
        self.SaveConicData(coincHdr, promptHdr, delayedHdr)            
        
        self.tagIt('Dummy')
        
    # ===========================================
    def SaveConicData(self, coincHdr, promptHdr, delayedHdr):
        '''Save Coincidence Info to UserData'''
        
        hdr_dispatcher = {
            'PromptData' : promptHdr,
            'DelayedData' : delayedHdr,
        }
        
        for name, tagHeader in hdr_dispatcher.items():
            myData = self.getTag('Dummy').getData(name)
            for recHdr in tagHeader.inputHeaders():
                result = recHdr.recTrigger()
                myData.append("energy", result.energy()/units.MeV
                    ).append("x", result.position().x()/units.mm
                    ).append("y", result.position().y()/units.mm
                    ).append("z", result.position().z()/units.mm
                    ).append("t_sec", recHdr.timeStamp().GetSec()
                    ).append("t_nsec", recHdr.timeStamp().GetNanoSec()
                    ).append("triggerType", result.triggerType()
                    ).append("detID", result.detector().detectorId())
        
        coincData = self.getTag('Dummy').getData('CoincData')
        coincRecHdrs = coincHdr.inputHeaders()
        coincData.set("multiplicity", len(coincRecHdrs))        
            
#####  Job Configuration for nuwa.py ########################################

def configure(argv=[]):
    pass

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

---