gun2.py

Go to the documentation of this file.

#!/usr/bin/env python

def gunGenerator2():
    import os
    loc = os.getenv('XMLDETDESCROOT')
    if not loc:
        raise RuntimeError,"Location of DetDesc XML file not specified"
    loc += "/DDDB/pmtbox.xml"
    
    import xmldetdesc
    xddc = xmldetdesc.XmlDetDescConfig(loc=loc)
                                


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

    #processor
    gen = app.algorithm("gunGenerator2")
    #gen.OutputLevel = 2
    #gen.GenTools = [ "GtHepEvtGenTool", "GtPositionerTool", "GtTimeratorTool" ]
    gen.GenTools = [ "GtGunGenTool/Gun2Gen", "GtPositionerTool/Gun2Pos", "GtTimeratorTool/Gun2Tim" ]
    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.Gun2Gen').OutputLevel = 2
    app.property('ToolSvc.Gun2Gen').ParticleName = "opticalphoton"
    app.property('ToolSvc.Gun2Gen').MomentumMode = "Fixed"
    app.property('ToolSvc.Gun2Gen').Momentum = 3.5*units.eV
    app.property('ToolSvc.Gun2Gen').DirectionMode = "Fixed"
    app.property('ToolSvc.Gun2Gen').Direction = [ 0.0,0.0,-1.0 ]

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