gun.py

Go to the documentation of this file.

#!/usr/bin/env python
'''
    Note that the volume is set to /dd/Structure/DayaBay,
    So the position is global position which is exactly what is used in Geant4.
    It is good for debugging.
    Zhe Wang Apr 9,2009
'''
class gun:
    '''
       configure gun generator
    '''

    def __init__(self,
                 stage='null',
                 name='gun',
                 volume='/dd/Structure/DayaBay',
                 start_time=0):

        # set up Gnrtr itself
        from Gnrtr.GnrtrConf import Gnrtr
        gnrtr = Gnrtr(name);

        gnrtr.GenTools = [ "GtGunGenTool/"+name+"Gen",
                           "GtPositionerTool/"+name+"Pos",
                           "GtTimeratorTool/"+name+"Tim",
                           "GtTransformTool/"+name+"Tra"]

        gnrtr.ThisStageName = "Kinematic"
        gnrtr.TimeStamp = start_time

        if stage != 'null':
            stage.KinematicSequence.Members.append(gnrtr)

        # set up each tools
        from GaudiKernel import SystemOfUnits as units
        from GenTools.GenToolsConf import GtPositionerTool, GtTransformTool, GtTimeratorTool

        # generator
        from GenTools.GenToolsConf import GtGunGenTool
        gun = GtGunGenTool(name+"Gen")
        gun.ParticlesPerEvent = 1000
        gun.ParticleName = "opticalphoton"
        gun.Momentum = 2*units.eV
        gun.MomentumMode = "Fixed"
        gun.MomentumSpread = 0*units.eV
        gun.DirectionMode = "Fixed"
        gun.Direction = [ 2.52361587010091171,
                          -0.129228715901263058,
                          -3.46788066267345130 ]
        gun.DirectionSpread = 0
        
        # Set up positioner
        pos=GtPositionerTool(name+'Pos',Volume=volume)
        pos.Strategy = "FullVolume" # also "AvoidDaughters" and "Surface"
        #pos.Mode = "Uniform"
        #pos.Mode = "Fixed"
        pos.Mode = "Smeared"        
        pos.Spread = 5*units.mm
        pos.Position = [-411762.29689663457*units.mm,
                        817680.94560857571*units.mm,
                        -1890.6120407181961*units.mm]
        #pos.Position = [-411763.29689663457*units.mm,
        #                817680.94560857571*units.mm,
        #                -1889.6120407181961*units.mm]
        
        # Set up timerator
        tim=GtTimeratorTool(name+'Tim')
        tim.LifeTime = int(30000000*units.ns)
        
        # transform
        tra=GtTransformTool(name+'Tra',Volume=volume)
        
        pass

if __name__ == "__main__":
    obj=gun()