In This Package:

main.py

Go to the documentation of this file.

#!/usr/bin/env python

class MyGenToolsConfig:
    def __init__(self,volume):
        from  GaudiPython import AppMgr
        from GaudiKernel import SystemOfUnits as units

        app = AppMgr()
        app.TopAlg = []

        # Set up timerator
        import GaudiKernel.SystemOfUnits as units
        tim = app.property("ToolSvc.GtTimeratorTool")
        tim.LifeTime = int(1*units.second)

        # Set up positioner
        poser = app.property("ToolSvc.GtPositionerTool")
        poser.Volume = volume
        poser.Mode = "Fixed"
        poser.Position = [0,0,5.0/16.0*units.meter]
        poser.OutputLevel = 7

        # Set up gun
        gun = app.property("ToolSvc.GtGunGenTool")
        gun.OutputLevel = 7
        gun.Volume = volume
        gun.ParticlesPerEvent = 1
        gun.ParticleName = "opticalphoton"
        gun.Momentum = 2.5*units.eV
        #gun.ParticleName = "mu+"
        #gun.Momentum = 10*units.GeV
        gun.MomentumMode = "Fixed"
        gun.DirectionMode = "Fixed"
        from math import sin, cos, pi
        pmt_column_number = 9
        angle = (2*pmt_column_number - 1)*pi/24.0;
        gun.Direction = [ cos(angle),sin(angle),0 ] # aim for PMT 
        print 'gun.Direction=',gun.Direction

        app.TopAlg += [ "GaudiSequencer/GenSeq" ]
        genseq = app.algorithm("GenSeq")
        genseq.Members = [ "GtGenerator/GenAlg", "GtHepMCDumper/GenDump" ]


        gen = app.algorithm("GenAlg")
        gen.OutputLevel = 7
        gen.GenTools = [ "GtGunGenTool", "GtPositionerTool", "GtTimeratorTool" ]
        gen.GenName = "Bang Bang"
        gen.Location = "/Event/Gen/HepMCEvents" # this is default anyways

        #print " GtDumper"
        gendump = app.algorithm("GenDump")
        gendump.Location = "/Event/Gen/HepMCEvents"  # this is default anyways.

        return

if '__main__' == __name__:

    from  GaudiPython import AppMgr

    app = AppMgr()
    app.EvtSel = "NONE"

    gtc = MyGenToolsConfig(volume="/dd/Structure/AD/la-oil2")

    ############### isolate them ##############
    app.TopAlg += [ "GiGaInputStream/GGInStream" ]
    from GGInStream import *
    GGInStream()

    app.TopAlg += [ "DetSimProc/DetSim" ]
    from DetSimProc import *
    DetSim()
    ###########################################

    
    app.run(2)
    

---