gun.py

Go to the documentation of this file.

#!/usr/bin/env python

def gunGenerator(volume):


    import GaudiPython as gm
    app = gm.AppMgr()            
    #msg = self.app.service("MessageSvc")
    #msg.OutputLevel = 1

    #processor
    gen = app.algorithm("gunGenerator")
    #gen.OutputLevel = 2
    #gen.GenTools = [ "GtHepEvtGenTool", "GtPositionerTool", "GtTimeratorTool" ]
    gen.GenTools = [ "GtGunGenTool/GunGen", "GtPositionerTool/GunPos", "GtTimeratorTool/GunTim" ]
    gen.GenName = "Bang Bang"
    gen.Location = "/Event/Gen/HepMCEvents" # this is default anyways
    #
    from GaudiKernel import SystemOfUnits as units


    # Set up gun
    app.property('ToolSvc.GunGen').OutputLevel = 0
    app.property('ToolSvc.GunGen').Volume=volume
    ##app.property('ToolSvc.GunGen').ParticleName = "mu+"
    app.property('ToolSvc.GunGen').ParticleName = "opticalphoton"
    app.property('ToolSvc.GunGen').ParticlesPerEvent = 100

    app.property('ToolSvc.GunGen').MomentumMode = "Fixed"
    app.property('ToolSvc.GunGen').Momentum = 2.5*units.eV
    app.property('ToolSvc.GunGen').DirectionMode = "Fixed"
    from math import sin, cos, pi
    pmt_column_number = 9
    angle = (2*pmt_column_number - 1)*pi/24.0;
    app.property('ToolSvc.GunGen').Direction = [ cos(angle),sin(angle),0 ]

    # Set up positioner
    #print " Creating poser"
    app.property('ToolSvc.GunPos').OutputLevel = 2
    #    app.property('ToolSvc.GunPos').Volume = "/dd/Structure/steel-2/water-2"
    app.property('ToolSvc.GunPos').Volume = volume
    app.property('ToolSvc.GunPos').Strategy = "FullVolume" # also "AvoidDaughters" and "Surface"
    app.property('ToolSvc.GunPos').Mode = "Uniform"
    app.property('ToolSvc.GunPos').Spread = 10*units.cm
    app.property('ToolSvc.GunPos').Position = [0,0,5.0/16.0*units.meter]
    
    # Set up timerator
    app.property('ToolSvc.GunTim').LifeTime = int(5*units.ns)
    
    
    return