In This Package:

ViktorsTest.py

Go to the documentation of this file.

#!/usr/bin/env python
#
# Example analysis script for Viktor
#
#  Usage:
#   nuwa.py -n -1 ViktorsTest simulationData.root

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

# Make shortcuts to ROOT classes
TH1F = gbl.TH1F

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

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

        # Initialize services
        #self.cableSvc = self.svc('ICableSvc','StaticCableSvc')

        self.stats.defaultPath = "/file1/examples/"
        # Make the trigger time histogram
        self.stats["TrigTime"] = TH1F("TrigTime","Trigger Time [s]",
                                      100,0.0,10.0)
        # Make the TDC histogram
        self.stats["Tdc"] = TH1F("Tdc","TDC Values",300,0,300)

        return SUCCESS

    def execute(self):
        self.info("executing")

        evt = self.evtSvc()
        hdr = evt["/Event/Readout/ReadoutHeader"]
        # Exercise 1: Print detector name
        ro = hdr.readout()
        if ro == None:
            self.info("No Readout this cycle")
            return SUCCESS
        print ro.detector().detName()
        # Exercise 2: Histogram the readout trigger time
        self.stats["TrigTime"].Fill( ro.triggerTime().GetSeconds() )        
        # Exercise 3: Histogram the TDC values
        channelIDs = ro.channels() # Get list of channel IDs in the readout
        for channelID in channelIDs:
            channel = ro.channel(channelID) # get the data from the channel
            for tdc in channel.tdc(): # loop over list of TDC values on channel
                self.stats["Tdc"].Fill( tdc )
            for adcClock in channel.adcClocks():
                adc = channel.adcByClock(adcClock)
                print adcClock, adc
        return SUCCESS
        
    def finalize(self):        
        self.info("finalizing")
        status = DybPythonAlg.finalize(self)
        return status

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

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

---