In This Package:

mymod.py

Go to the documentation of this file.

def run(app):
    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):
            status = GaudiAlgo.initialize(self)
            # just example, could do stuff here too.
            return status

        def execute(self):

            print "Init MyAlg",self.name
            
            # Get a volume from the TDS
            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"
            dboil1 = self.getDet(de_name)
            gi = dboil1.geometry()

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

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

            # convert origin of oil into global coordinates
            point = Gaudi.XYZPoint(0.0,0.0,0.0)
            gpoint = gi.toGlobal(point)
            print 'Global point: (',gpoint.x(),gpoint.y(),gpoint.z(),')'

            # Starting from the top, find what the smalled DE 
            # that holds point is
            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

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

---