In This Package:

pmtbox.py

Go to the documentation of this file.

#!/usr/bin/env python

class MyGenToolsConfig:
    def __init__(self,volume):
        import GaudiPython as gm
        app = gm.AppMgr()
        app.TopAlg = []


        ########### KINEMATICS ##########


        # 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.Strategy = "AvoidDaughters"
        poser.Strategy = "FullVolume"
        poser.Mode = "Smeared"
        poser.Spread = 1.0*units.cm
        poser.Position = [0.0,0.0,0.0];
        poser.OutputLevel = 1

        # Set up gun
        gun = app.property("ToolSvc.GtGunGenTool")
        gun.ParticlesPerEvent = 100
        gun.ParticleName = "opticalphoton"
        gun.Momentum = 2.5*units.eV
        #gun.ParticleName = "mu+"
        #gun.Momentum = 10*units.GeV
        gun.MomentumMode = "Fixed"
        gun.DirectionMode = "Fixed"
        gun.Direction = [ 0.243,0.0,-0.970 ]
        #gun.OutputLevel = 1


        ### Now initi algs

        app.TopAlg += [ "GtGenerator/Generator", "GtHepMCDumper/Dumper" ]


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

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

        ########### SIMULATION ##########

        from GaudiKernel import SystemOfUnits as units
        app.ExtSvc += ["GiGa"]

        modularPL = app.property("GiGa.GiGaPhysListModular")
        modularPL.OutputLevel = 1
        modularPL.CutForElectron = 100*units.micrometer
        modularPL.CutForPositron = 100*units.micrometer
        modularPL.CutForGamma = 1*units.millimeter
        modularPL.PhysicsConstructors = [ 
            "DsPhysConsGeneral", 
            "DsPhysConsOptical" 
            ]

        giga = app.service("GiGa")
        giga.OutputLevel = 1
        giga.PhysicsList = "GiGaPhysListModular"

        gggeo = app.service("GiGaGeo")
        gggeo.OutputLevel = 1

        # Make Geant4 sing!
        ggrm = app.property("GiGa.GiGaMgr")
        ggrm.Verbosity = 9
        event_ac_cmds = app.property("GiGa.GiGaEventActionCommand")
        event_ac_cmds.BeginOfEventCommands = [
            "/control/verbose 2",
            "/run/verbose 1",
            "/event/verbose 2",
            "/tracking/verbose 2",
            "/geometry/navigator/verbose 2"
            ]
        giga.EventAction = "GiGaEventActionCommand"


        app.TopAlg += [ "GaudiSequencer/Simulation" ]
        seq = app.algorithm("Simulation")
        seq.Members = []

        #seq.Members = [ "GiGaInputStream/GGInStream" ]
        ggin = app.algorithm("GGInStream")
        ggin.OutputLevel = 1
        ggin.ExecuteOnce = True
        ggin.ConversionSvcName = "GiGaGeo"
        ggin.DataProviderSvcName = "DetectorDataSvc"
        ggin.StreamItems = [
            "/dd/Structure/world",
            "/dd/Geometry/PmtBox/PmtBoxSurfaces",
            ]

        seq.Members = [ "GiGaInputStream/GGInStream", "DsPushKine/PushKine", "DsPullEvent/PullEvent" ]
        push = app.algorithm("PushKine")
        push.Converter = "HepMCtoG4"


        return

if '__main__' == __name__:
    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)
    gtc = MyGenToolsConfig(volume="/dd/Structure/steel-2/water-2")
    
    import gaudimodule as gm
    app = gm.AppMgr()
    app.EvtSel = "NONE"
    app.run(2)
    

---