In This Package:

ds.py

Go to the documentation of this file.

#!/usr/bin/env python

def DetSim():
    from  GaudiPython import AppMgr
    from GaudiKernel import SystemOfUnits as units
    
    app = AppMgr()

    #######################################################################
    app.ExtSvc += ["GiGa"]
    
    modularPL = app.property("GiGa.GiGaPhysListModular")
    modularPL.OutputLevel = 0
    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 = 0
    giga.PhysicsList = "GiGaPhysListModular"
    
    gggeo = app.service("GiGaGeo")
    gggeo.OutputLevel = 0
    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 = 0
    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"
    

    ds = app.algorithm("DetSim")
    ds.OutputLevel = 1
    ds.Converter = "HepMCtoG4"
    
    # Class name to use is set in DetDesc xml's "sensdet" attribute.
    pmtsd = app.property("GiGaGeo.DsPmtSensDet")
    pmtsd.OutputLevel = 0
    
    return

---