TestSimDataSvc.py

Go to the documentation of this file.

#!/usr/bin/env python

from GaudiPython import *
from DybPython.DybPythonAlg import DybPythonAlg
import GaudiKernel.SystemOfUnits as units

from DataSvc.TestCommon import get_sensors, service_mode

GSimPmtSpec = gbl.GSimPmtSpec

class TestAlg(DybPythonAlg):
    def __init__(self,name):
        DybPythonAlg.__init__(self,name)
        print "Making TestAlg",name

    def initialize(self):
        status = DybPythonAlg.initialize(self)
        print "Init TestAlg",self.name()
        if status.isFailure(): return status
        self.cableSvc = self.svc('ICableSvc','CableSvc')
        self.pmtSimSvc = self.svc('ISimDataSvc','SimDataSvc')
        print "===========================================================\n\n"

        svcMode = service_mode()
        adPmtSensors = get_sensors(self.cableSvc,svcMode)

        hrdwTypes = ["kUnknown","kPmt8inch","kPmt2inch","kRpc","kFee","kFec","kRot","kRom","kHighVoltage"]
        print "  number of sensors =",adPmtSensors.size()
        print "  Sensor  Gain FWHM Toff Tspr Eff PreP AfterP Dark"
        for sens in adPmtSensors:
            pmtData = self.pmtSimSvc.pmtSimData(sens, svcMode)
            if not pmtData:
                print "No calibration data for PMT ring=%d,column=%d" % (sens.ring(),sens.column()) 
                continue
            print ' '.join([ str(elem) for elem in [pmtData.m_pmtId.fullPackedData(),
                                                    pmtData.m_gain,
                                                    pmtData.m_sigmaGain,
                                                    pmtData.m_timeOffset/units.ns,
                                                    pmtData.m_timeSpread/units.ns,
                                                    pmtData.m_efficiency,
                                                    pmtData.m_prePulseProb,
                                                    pmtData.m_afterPulseProb,
                                                    pmtData.m_darkRate/units.hertz] ])
        print "\n\n==========================================================="
        return SUCCESS

    def execute(self):
        print "Executing TestAlg",self.name()
        return SUCCESS

    def finalize(self):
        print "Finalizing TestAlg",self.name()
        status = DybPythonAlg.finalize(self)
        return status

def run(app):
    testAlg = TestAlg("TestAlg")
    app.addAlgorithm( testAlg )
    return