In This Package:

adgun.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.Position = [0,0,0]
        poser.OutputLevel = 5

        # Set up gun
        gun = app.property("ToolSvc.GtGunGenTool")
        gun.OutputLevel = 5
        gun.Volume = volume
        gun.ParticlesPerEvent = 1
        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
        
        # Single optical photon which makes a hit:
        #gun.ParticleName = "opticalphoton"
        #gun.Momentum = 2.5*units.eV
                
        gun.ParticleName = "e+"
        gun.Momentum = 1.0*units.MeV
        
        hepevt = app.property("ToolSvc.GtHepEvtGenTool")
        hepevt.HepEvtDataSource = "seeds.dat" #"|InverseBeta.exe -n 1000"
        

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


        gen = app.algorithm("GenAlg")
        gen.OutputLevel = 5
        gen.GenName = "Stick 'em up, bub, I got a gun."
        gen.GenTools = [ "GtGunGenTool", "GtPositionerTool", "GtTimeratorTool" ]
        #gen.GenTools = [ "GtHepEvtGenTool", "GtPositionerTool", "GtTimeratorTool" ]
        #gen.GenName = "Da Invoise Betah"
        gen.Location = "/Event/Gen/HepMCEvents" # this is default anyways

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

        app.ExtSvc += ["GiGa"]

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

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

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

        # Make Geant4 sing!
        ggrm = app.property("GiGa.GiGaMgr")
        ggrm.Verbosity = 2

        # event_ac_cmds = app.property("GiGa.GiGaEventActionCommand")
        # event_ac_cmds.BeginOfEventCommands = [
        #     "/control/verbose 0",
        #     "/run/verbose 0",
        #     "/event/verbose 0",
        #     "/tracking/verbose 0",
        #     "/geometry/navigator/verbose 0"
        #     ]
        #giga.EventAction = "GiGaEventActionCommand"


        # History!



        app.TopAlg += [ "GaudiSequencer/SimSeq" ]
        simseq = app.algorithm("SimSeq")
        simseq.Members = [ "GiGaInputStream/GGInStream", 
                           "DsPushKine/PushKine", 
                           "DsPullEvent/PullEvent",
                           "DrawHistoryAlg/DrawHistory",
                           "DumpUnobservableStatisticsAlg/DumpUnobserved"
                          ]

        ggin = app.algorithm("GGInStream")
        ggin.OutputLevel = 5
        ggin.ExecuteOnce = True
        ggin.ConversionSvcName = "GiGaGeo"
        ggin.DataProviderSvcName = "DetectorDataSvc"
        ggin.StreamItems = [
            "/dd/Structure/Sites/la-rock",
            #"/dd/Structure/AD/far-ade3",
            ]

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

        # Class name to use is set in DetDesc xml's "sensdet" attribute.
        pmtsd = app.property("GiGaGeo.DsPmtSensDet")
        pmtsd.OutputLevel = 5  # change to 1 to get TouchableHistory of the hits.
        pmtsd.TouchableToDetelem = "TouchableToDetectorElementFast"
        
        
      
        giga.SteppingAction = "GiGaStepActionSequence"
        stepseq = app.property("GiGa.GiGaStepActionSequence")        
        stepseq.Members = [ "HistorianStepAction", "UnObserverStepAction" ]
        
        Historian = app.property("GiGa.GiGaStepActionSequence.HistorianStepAction");
        Historian.TrackSelection = "pdg != 20022"
        Historian.VertexSelection = "all"
        
        UnObserver = app.property("GiGa.GiGaStepActionSequence.UnObserverStepAction")
        
        # This is the default search path, which doesn't include PMTs (for speed)
        UnObserver.DetectorElementSearchPath = [
          "/dd/Structure/DayaBay"
         ,"/dd/Structure/Sites"
        # ,"/dd/Structure/Pool"
         ,"/dd/Structure/AD"
        # ,"/dd/Structure/RPC"
        ]
        Historian.DetectorElementSearchPath = [
           "/dd/Structure/DayaBay"
          ,"/dd/Structure/Sites"
         # ,"/dd/Structure/Pool"
          ,"/dd/Structure/AD"
         # ,"/dd/Structure/RPC"
         ]

        # These are the default parameters.
        UnObserver.Stats = [
         [ "MuonTrkLengthInIws", "dx", " (pdg==13 or pdg==-13) and ((Volume == '/dd/Structure/Pool/far-iws' or (Volume == '/dd/Structure/Pool/db-iws' or Volume == '/dd/Structure/Pool/la-iws')) and MaterialName == '/dd/Materials/Water')"]
        ]
        # UnObserver.Stats = [
        #     ["photon_created_energy" , "E" , "StepNumber==1 and pdg==20022" ]
        #    ,["photon_backscatter_r" , "local_r" , "pdg==20022 and dAngle >= 90" ]
        #    ,["photon_forescatter_r" , "local_r" , "pdg==20022 and (dAngle <= 90 and dAngle > 0)" ]
        #    ,["photon_stop_r" ,        "local_r" , "(pdg==20022) and ((IsStopping == 1) and (LogicalVolumeName == '/dd/Geometry/PMT/lvPmtHemi'))" ]
        #    ,["pmt_hit" , "t" , "pdg==20022 and (IsStopping == 1 and (LogicalVolumeName=='/dd/Geometry/PMT/lvPmtHemi'))" ]
        #    ,["edep-water",  "dE", "pdg!=20022 and (MaterialName == '/dd/Materials/Water')"]
        #    ,["edep-oil",    "dE", "pdg!=20022 and (MaterialName == '/dd/Materials/Oil')"]
        #    ,["edep-ls",     "dE", "pdg!=20022 and (MaterialName == '/dd/Materials/LiquidScintillator')"]
        #    ,["edep-gdls",   "dE", "pdg!=20022 and (MaterialName == '/dd/Materials/GdDopedLS')"]
        #    ,["edep-acryilc","dE", "pdg!=20022 and (MaterialName == '/dd/Materials/Acrylic')"]
        #    ,["edep-ad1"  ,"dE"  ,"pdg!=20022 and ((MaterialName == '/dd/Materials/LiquidScintillator' or MaterialName == '/dd/Materials/GdDopedLS') and AD==1)"  ]
        #    ,["edep-ad2"  ,"dE"  ,"pdg!=20022 and ((MaterialName == '/dd/Materials/LiquidScintillator' or MaterialName == '/dd/Materials/GdDopedLS') and AD==2)"  ]
        #    ,["edep-ad3"  ,"dE"  ,"pdg!=20022 and ((MaterialName == '/dd/Materials/LiquidScintillator' or MaterialName == '/dd/Materials/GdDopedLS') and AD==3)"  ]
        #    ,["edep-ad4"  ,"dE"  ,"pdg!=20022 and ((MaterialName == '/dd/Materials/LiquidScintillator' or MaterialName == '/dd/Materials/GdDopedLS') and AD==4)"  ]
        #    ,["qedep-ad1" ,"qdE" ,"pdg!=20022 and ((MaterialName == '/dd/Materials/LiquidScintillator' or MaterialName == '/dd/Materials/GdDopedLS') and AD==1)"   ]
        #    ,["qedep-ad2" ,"qdE" ,"pdg!=20022 and ((MaterialName == '/dd/Materials/LiquidScintillator' or MaterialName == '/dd/Materials/GdDopedLS') and AD==2)"   ]
        #    ,["qedep-ad3" ,"qdE" ,"pdg!=20022 and ((MaterialName == '/dd/Materials/LiquidScintillator' or MaterialName == '/dd/Materials/GdDopedLS') and AD==3)"   ]
        #    ,["qedep-ad4" ,"qdE" ,"pdg!=20022 and ((MaterialName == '/dd/Materials/LiquidScintillator' or MaterialName == '/dd/Materials/GdDopedLS') and AD==4)"   ]
        #    ,["ad" , "AD" ,
        #      "expos" ,"Ex"   ,
        #      "eypos" ,"Ey"   ,
        #      "ezpos" ,"Ez"   ,
        #      "et"    ,"Et"   ,
        #      "qexpos" ,"qEx" ,
        #      "qeypos" ,"qEy" ,
        #      "qezpos" ,"qEz" ,
        #      "qet"    ,"qEt" , "(pdg!=20022) and (MaterialName == '/dd/Materials/LiquidScintillator' or MaterialName == '/dd/Materials/GdDopedLS')"  ]
        #    ,["mu-path-water", "dx", "(pdg==13 || pdg==-13) and MaterialName == '/dd/Materials/Water'"]
        #    ,["mu-path-oil",   "dx", "(pdg==13 || pdg==-13) and MaterialName == '/dd/Materials/Oil'"]
        #    ,["mu-path-ls",    "dx", "(pdg==13 || pdg==-13) and MaterialName == '/dd/Materials/LiquidScintillator'"]
        #    ,["mu-path-gdls",  "dx", "(pdg==13 || pdg==-13) and MaterialName == '/dd/Materials/GdDopedLS'"]
        # ]

        app.algorithm("DrawHistory").do_hits = 0

        return
        

if '__main__' == __name__:
    import os
    print "\tConfiguring geometry\n"
    import xmldetdesc
    xmldetdesc.config()
 
    from GaudiPython import AppMgr
    app = AppMgr()
    app.EvtSel = "NONE"

    #gtc = MyGenToolsConfig(volume="/dd/Structure/AD/la-oil2")
    gtc = MyGenToolsConfig(volume="/dd/Structure/AD/la-gds2")
    
    app.run(10)
    

---