In This Package:

Example_SimEvent.py

Go to the documentation of this file.

#!/usr/bin/env python

# Load GaudiPython
from GaudiPython.GaudiAlgs import GaudiAlgo
from GaudiPython import SUCCESS, FAILURE
from GaudiPython import gbl
from DybPython.Util import irange

# Make shortcuts to ROOT classes
TH1F = gbl.TH1F

# Make your algorithm
class ExampleAlg(GaudiAlgo):
    "Example Python Algorithm"
    def __init__(self,name):
        GaudiAlgo.__init__(self,name)
        return

    def initialize(self):
        status = GaudiAlgo.initialize(self)
        print "Init ExampleAlg",self.name()
        if status.isFailure(): return status
        return SUCCESS

    def execute(self):
        print "Executing ExampleAlg",self.name()
        evt = self.evtSvc()
        hdr = evt["/Event/Sim/SimHeader"]
        shh = hdr.hits()
        for detector in shh.hitDetectors():
            print "Detector: ",detector
            hitCollection = shh.hitsByDetector(detector)
            print "Hits: ",hitCollection
            if hitCollection == None:
                print "No Hit Collection for detector ",detector
            hits = hitCollection.collection()
            for hit in hits:
                print "  PMT",hit.sensDetId(),"hit at time",hit.hitTime() 
        return SUCCESS
        
    def finalize(self):        
        print "Finalizing ExampleAlg",self.name()
        status = GaudiAlgo.finalize(self)
        return status

def configure():
    return

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

---