In This Package:

detsim.py

Go to the documentation of this file.

#!/usr/bin/env python

# This is an EXAMPLE of how to setup DetSim

class DetSimConfig:
    def __init__(self):
        import GaudiPython as gm
        from GaudiKernel import SystemOfUnits as units
        app = gm.AppMgr()
        app.ExtSvc += ["GiGa"]

        # Connection to Geant4
        giga = app.service("GiGa")
        giga.OutputLevel = 1

        # Physics List
        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",
            #"DsPhysConsHadron",
            #"DsPhysConsXxx",
            ]
        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/Sites/far-rock"
            "/dd/Geometry/PMT/PmtSurfaces",
            #"/dd/Structure/pmt-box-array",
            ]

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

        


if '__main__' == __name__:
    import xmldetdesc,gentools,os
    loc = os.getenv('XMLDETDESCROOT')
    if not loc:
        raise RuntimeError,"Location of DetDesc XML file not specified"
    loc += "/DDDB/dayabay.xml"
    xddc = xmldetdesc.XmlDetDescConfig(loc=loc)
    tg = gentools.GenToolsConfig(volume="/dd/Geometry/Pool/lvFarPoolIWS")
    dsc = DetSimConfig()
    
    import gaudimodule as gm
    app = gm.AppMgr()
    app.EvtSel = "NONE"
    app.run(10)
    

---