In This Package:

poolgun.py

Go to the documentation of this file.

#!/usr/bin/env python

TH2DE="TH2DE"
#TH2DE="TouchableToDetectorElementFast"
spray_photons = True


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*units.meter,0*units.meter,4.35*units.meter]
#        poser.Position = [0.5*units.meter,1.158*units.meter,-3.5*units.meter]
        #poser.OutputLevel = 1

        # Set up gun
        gun = app.property("ToolSvc.GtGunGenTool")
        #gun.OutputLevel = 1
        if spray_photons:
            gun.ParticlesPerEvent = 100
            gun.ParticleName = "opticalphoton"
            gun.Momentum = 2.068*units.eV
            gun.PolarizeMode = "Random"
        else:
            gun.ParticlesPerEvent = 1
            gun.ParticleName = "mu+"
            gun.Momentum = 10*units.GeV
            pass
        gun.PolarizeMode="Random"
        gun.MomentumMode = "Fixed"
        gun.DirectionMode = "Fixed"
        gun.Direction = [1,1,0 ]
        print 'gun.Direction=',gun.Direction

        trans = app.property("ToolSvc.GtTransformTool")
        trans.Volume = volume

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


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

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

        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",
            "DsPhysConsEM",
            #"DsPhysConsElectroNu",
            #"DsPhysConsHadron",
            #"DsPhysConsIon"
            ]
#        optical = app.property("GiGa.GiGaPhysListModular.DsPhysConsOptical")
#        optical.UseScintillation = True


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

        gggeo = app.service("GiGaGeo")
        #gggeo.OutputLevel = 1
        gggeo.XsizeOfWorldVolume = 2.4*units.kilometer
        gggeo.YsizeOfWorldVolume = 2.4*units.kilometer
        gggeo.ZsizeOfWorldVolume = 2.4*units.kilometer

        giga.SteppingAction = "GiGaStepActionSequence"
        stepseq = app.property("GiGa.GiGaStepActionSequence")
        stepseq.Members = ["HistorianStepAction","UnObserverStepAction"]

        historian = app.property("GiGa.GiGaStepActionSequence.HistorianStepAction")
        historian.TouchableToDetelem = TH2DE

        params = {
            'muonid':"(pdg==13 or pdg==-13)",
            'ows':"(MaterialName == '/dd/Materials/Water' and "+\
                "Volume == '/dd/Structure/Pool/far-ows')",
            'iws':"(MaterialName == '/dd/Materials/Water' and "+\
                "Volume == '/dd/Structure/Pool/far-iws')",
            'lso':"MaterialName == '/dd/Materials/LiquidScintillator'",
            'gds':"MaterialName == '/dd/Materials/GdDopedLS'"
            }
        unobs = app.property("GiGa.GiGaStepActionSequence.UnObserverStepAction")
        unobs.TouchableToDetelem = TH2DE
        unobs.Stats=[
            ["MuonTrkLengthInOws","dx","%(muonid)s and %(ows)s"%params],
            ["MuonTrkLengthInIws", "dx", "%(muonid)s and %(iws)s"%params],
            ["MuonTrkLengthInLS", "dx", "%(muonid)s and %(lso)s"%params],
            ["MuonTrkLengthInGdLS","dx", "%(muonid)s and %(gds)s"%params]
            ]


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


        app.TopAlg += [ "GaudiSequencer/SimSeq" ]
        simseq = app.algorithm("SimSeq")
        simseq.Members = [ "GiGaInputStream/GGInStream" ]

        ggin = app.algorithm("GGInStream")
        #ggin.OutputLevel = 1
        ggin.ExecuteOnce = True
        ggin.ConversionSvcName = "GiGaGeo"
        ggin.DataProviderSvcName = "DetectorDataSvc"
        ggin.StreamItems = [ "/dd/Structure/Sites/far-rock",
                             "/dd/Geometry/PoolDetails/FarPoolSurfaces",
                             "/dd/Geometry/PoolDetails/PoolSurfacesAll",
                             ]


        simseq.Members += [ "DsPushKine/PushKine", "DsPullEvent/PullEvent" ]
        push = app.algorithm("PushKine")
        push.Converter = "HepMCtoG4"
        push.Location = "/Event/Gen/GenHeader"

        pull = app.algorithm("PullEvent")
        pull.Location = "/Event/Sim/SimHeader"

        # Class name to use is set in DetDesc xml's "sensdet" attribute.
        pmtsd = app.property("GiGaGeo.DsPmtSensDet")
        #pmtsd.SensorStructures = [ "/dd/Structure/AdPmtStructure",
        #                           "/dd/Structure/PoolPmtStructure" ]
        #pmtsd.SensorStructures = [ "/dd/Structure/AdPmtStructure" ]
        pmtsd.SensorStructures = [ ]
        #pmtsd.OutputLevel = 1
        pmtsd.TouchableToDetelem = "TH2DE"
        #pmtsd.TouchableToDetelem = "TouchableToDetectorElementFast"

        return

if '__main__' == __name__:
    import os

    print "\tConfiguring geometry\n"
    import xmldetdesc
    xmldetdesc.config()

    from GaudiPython import AppMgr
    app = AppMgr()
    app.EvtSel = "NONE"

    eds=app.service("EventDataSvc")
    eds.ForceLeaves=True
    
    #xmldetdesc.dump()

    gtc = MyGenToolsConfig(volume="/dd/Structure/Pool/far-iws")
    
    app.run(1)
    

---