In This Package:

geo Namespace Reference

Functions
def	configure
def	run

---

Function Documentation

def geo::configure

(

)

Definition at line 3 of file geo.py.

00003                :
00004 
00005 #     from Gaudi.Configuration import ApplicationMgr
00006 #     theApp = ApplicationMgr()
00007 #     theApp.TopAlg.append("myalg")
00008     return
00009 
00010 
def run(app):

def geo::run

(

app

)

Definition at line 11 of file geo.py.

            :
    from GaudiPython.GaudiAlgs import GaudiAlgo
    from GaudiPython import SUCCESS, FAILURE

    class MyAlg(GaudiAlgo):
        def __init__(self,name="MyAlg"):
            print "Making MyAlg",name
            GaudiAlgo.__init__(self,name)
            return

        def initialize(self):
            return GaudiAlgo.initialize(self)

        def execute(self):

            print "Init MyAlg",self.name
            
            oil = self.getDet('/dd/Geometry/AD/lvOIL')
            print 'oil=',oil.name()
            print 'material=',oil.materialName()

            #de_name = "/dd/Structure/DayaBay/db-rock/db-ows/db-curtain/db-iws/db-ade1/db-sst1/db-oil1"
            de_name = '/dd/Structure/AD/db-ade1/db-sst1/db-oil1'
            dboil1 = self.getDet(de_name)
            print 'dboil1=',dboil1.name()

            top_name = "/dd/Structure/DayaBay"
            top = self.getDet(top_name)
            top_gi = top.geometry()

            gi = dboil1.geometry()
            print 'dboil1.gi.lvolumename=',gi.lvolumeName()

            print 'types = %s, %s, %s'%(type(oil),type(dboil1),type(gi))

            import PyCintex
            Gaudi = PyCintex.makeNamespace('Gaudi')

            gbl = PyCintex.makeNamespace('')
            #print gbl.std.__dict__
            loc = gi.parentIGeometryInfo().parentIGeometryInfo()

            ret = gbl.ILVolume.ReplicaPath()
            print ret
            stat = gi.location(loc,ret)
            print stat
            print ret,'of size',ret.size()
            for ind in range(0,ret.size()):
                print ret[ind]

            point = Gaudi.XYZPoint(0.0,0.0,0.0)
            gpoint = gi.toGlobal(point)
            path = gbl.std.vector("unsigned int")()

            top_gi.fullGeoInfoForPoint(gpoint,10,loc,path)
            print gi.lvolumeName()
            print loc.lvolumeName()
            print 'path=',path.size()

            s = top_gi.belongsToPath(gpoint,-1)
            print 'belongsToPath:',s
            de = self.getDet(s)
            while de:
                print 's_de.name():',de.name()
                de = de.parentIDetectorElement()
                continue
            

            return SUCCESS
        pass

    #from GaudiPython import AppMgr
    #app = AppMgr()

    ma = MyAlg("myalg")
    app.addAlgorithm(ma)

---