In This Package:

TestCalibDataSvc.py

Go to the documentation of this file.

#!/usr/bin/env python

from GaudiPython import *
from DybPython.DybPythonAlg import DybPythonAlg

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

CalibPmtSpec = gbl.CalibPmtSpec

class TestAlg(DybPythonAlg):
    def __init__(self,name):
        DybPythonAlg.__init__(self,name)
        self.calibDataSvcName = '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.pmtCalibSvc = self.svc('ICalibDataSvc',self.calibDataSvcName)
        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 Status SpeHigh sigmaSpe SpeLow Toff Tspr Eff PreP AfterP Dark"
        for sens in adPmtSensors:
            pmtData = self.pmtCalibSvc.pmtCalibData(sens, svcMode)
            print ' '.join([ str(elem) for elem in [pmtData.m_pmtId.fullPackedData(),
                                                    pmtData.m_status,
                                                    pmtData.m_speHigh,
                                                    pmtData.m_sigmaSpeHigh,
                                                    pmtData.m_speLow,
                                                    pmtData.m_timeOffset,
                                                    pmtData.m_timeSpread,
                                                    pmtData.m_efficiency,
                                                    pmtData.m_prePulseProb,
                                                    pmtData.m_afterPulseProb,
                                                    pmtData.m_darkRate] ])
        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

calibDataSvcName = "StaticCalibDataSvc"

def configure(argv = []):
    global calibDataSvcName
    import getopt
    opts,args = getopt.getopt(argv,"",["svc-name="])
    for opt,arg in opts:
        if opt == "--svc-name":
            calibDataSvcName = arg
    import DataSvc
    DataSvc.Configure( calibDataSvc = calibDataSvcName )
    if calibDataSvcName == 'DbiCalibDataSvc':
        # Load the reflex dictionary for CalibPmtSpec
        print "Loading DBI class: CalibPmtSpec"
        CalibPmtSpec = gbl.CalibPmtSpec
    return

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

---