In This Package:

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

Detector = gbl.DayaBay.Detector
ServiceMode = gbl.ServiceMode
Context = gbl.Context
Site = gbl.Site
DetectorId = gbl.DetectorId

GSimPmtSpec = gbl.GSimPmtSpec

class TestAlg(DybPythonAlg):
    def __init__(self,name):
        DybPythonAlg.__init__(self,name)
        self.simDataSvcName = 'StaticSimDataSvc'
        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','StaticCableSvc')
        self.pmtSimSvc = self.svc('ISimDataSvc',self.simDataSvcName)
        print "===========================================================\n\n"
        # Select all sites
        context = Context()
        context.SetSite(Site.kAll)
        # Select all detectors
        context.SetDetId(DetectorId.kAll)
        svcMode = ServiceMode(context,0)
        adPmtSensors = self.cableSvc.adPmtSensors(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

simDataSvcName = "StaticSimDataSvc"

def configure(argv = []):
    global simDataSvcName
    import getopt
    opts,args = getopt.getopt(argv,"",["svc-name="])
    for opt,arg in opts:
        if opt == "--svc-name":
            simDataSvcName = arg
    import DataSvc
    DataSvc.Configure( simDataSvc = simDataSvcName )
    return

def run(app):
    global simDataSvcName
    app.ExtSvc += [ "StaticCableSvc", simDataSvcName ]
    testAlg = TestAlg("TestAlg")
    testAlg.simDataSvcName = simDataSvcName
    app.addAlgorithm( testAlg )
    return

---