In This Package:

Example.py

Go to the documentation of this file.

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

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

# Make shortcuts to any ROOT classes you want to use
TH1F = gbl.TH1F

# Make your algorithm
class ExampleAlg(DybPythonAlg):
    "Example Python Algorithm"
    def __init__(self,name):
        DybPythonAlg.__init__(self,name)
        # Define Properties and default values here
        self.MyProperty = 1
        return

    def initialize(self):
        status = DybPythonAlg.initialize(self)
        if status.isFailure(): return status
        self.info("initializing")

        # Initialize services
        #  Cable Service: Use for pmt and electronics channel info
        # self.cableSvc = self.svc('ICableSvc','StaticCableSvc')

        # Make a histogram
        self.stats["/file1/path/to/MyHist1"] = TH1F("MyHist1","Test Histogram",
                                                    100,0.0,10.0)
        return SUCCESS

    def execute(self):
        self.info("executing")
        
        # Access current data
        evt = self.evtSvc()
        currentHdr = evt["/Event/Readout/ReadoutHeader"]
        if currentHdr == None:
            self.error("Failed to get current readout header")
            return FAILURE

        # Access recent data in archive
        readoutArchive = self.getAES("/Event/Readout/ReadoutHeader")
        if readoutArchive == None:
            self.error("Failed to get recent readout headers")
            return FAILURE
        # Loop over all readouts in archive, including current readout
        for readoutHdr in readoutArchive:
            readout = readoutHdr.readout()
        # Access current readout (First entry in archive!)
        firstHdr = readoutArchive[0]
        # This is true: firstHdr == currentHdr
        # Loop over previous readouts, not including current readout
        for readoutHdr in readoutArchive[1:]:
            readout = readoutHdr.readout()

        # Do some calculation...
        # ...
        # ...
        myResult = 3.14
        # Add entry to histogram
        self.stats["/file1/path/to/MyHist1"].Fill( myResult )

        # Print messages to the job log
        self.error("an error message")
        self.warning("a warning")
        self.info("some information")
        self.debug("some extra debugging information")
        self.verbose("too much information")
        return SUCCESS
        
    def finalize(self):
        self.info("finalizing")
        status = DybPythonAlg.finalize(self)
        return status


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

def configure():
    from StatisticsSvc.StatisticsSvcConf import StatisticsSvc
    statsSvc = StatisticsSvc()
    statsSvc.Output ={"file1":"readoutResult.root"}
    return

def run(app):
    '''
    Configure and add an algorithm to job
    '''
    app.ExtSvc += ["StaticCableSvc", "StatisticsSvc"]
    example = ExampleAlg("MyExample")
    # Set your algorithm Properties here
    example.MyProperty = 2
    app.addAlgorithm(example)
    pass

---