In This Package:

test.py

Go to the documentation of this file.

#!/usr/bin/env python

import PyCintex
gbl = PyCintex.makeNamespace("")

import xmldetdesc
xmldetdesc.config()

from Gaudi.Configuration import *
conf = ApplicationMgr()

from DetHelpers.DetHelpersConf import PmtGeomInfoSvc
pgisvc = PmtGeomInfoSvc("PmtGeomInfoSvc")
#pgisvc.OutputLevel = 1
pgisvc.StreamItems = [ "/dd/Structure/DayaBay" ]
print pgisvc
conf.ExtSvc += [ pgisvc ]
print "Done with configure"

from GaudiPython.GaudiAlgs import GaudiAlgo
from GaudiPython import *
class MyAlg(GaudiAlgo):
    def __init__(self,name):
        GaudiAlgo.__init__(self,name)
        print "Making MyAlg",name

    def initialize(self):
        status = GaudiAlgo.initialize(self)
        print "Init MyAlg",self.name
        if status.isFailure(): return status
        app = AppMgr()
        app.DLLs += [ "DetHelpersLib", 'DetHelpers', 'DetHelpersDict' ]
        self.pgis = self.svc('IPmtGeomInfoSvc','PmtGeomInfoSvc')
        print self.pgis
        return self.spin()

    def spin(self):
        # check all AD pmts
        import time
        start = time.time()
        count = 0
        for site,nads in [(0x01,2),(0x02,2),(0x04,4)]:
            print "Site:",site,
            for iad in range(0,nads):
                iad += 1
                print "[",iad,"]",
                for icol in range(0,24):
                    icol +=1
                    for iring in range(0,8):
                        iring += 1
                        pmtid = (site<<24)|(iad<<16)|(iring<<8)|icol

                        pgi = self.pgis.get(pmtid)
                        if not pgi:
                            print "PmtGeomInfoSvc returned None"
                            return FAILURE
                        if pmtid != pgi.pmtid():
                            print "Ask for PMT ID %d, got %d"%(pmtid,pgi.pmtid())
                            return FAILURE
                        count += 1
                        continue
                    continue
                continue
            print "\n",
            continue
        stop = time.time()
        print "Found %d PMTs in %f seconds"%(count,stop-start)
        

        return SUCCESS

    def execute(self):
        for r in range(0,10): self.spin()

    pass



from GaudiPython import AppMgr
app = AppMgr()
app.EvtSel = "NONE"
app.ExtSvc += [ "PmtGeomInfoSvc" ]
#app.DLLs += [ "DetHelpersLib", 'DetHelpers', 'DetHelpersDict' ]

myalg = MyAlg("MyAlg")
print 'myalg="%s"'%str(myalg)
app.setAlgorithms( [ myalg ] )
print ' running'
app.run(1)

---