In This Package:

beamer.py

Go to the documentation of this file.

#!/usr/bin/env python
"""Test out GtBeamerTool"""

from GaudiKernel import SystemOfUnits as units

def config_geom():
    print "\tConfiguring geometry...",
    import xmldetdesc
    xmldetdesc.config()
    print "done"

def config_kin(app):


    # Set up timerator
    import GaudiKernel.SystemOfUnits as units
    tim = app.property("ToolSvc.GtTimeratorTool")
    tim.LifeTime = int(1*units.second)
    tim.OutputLevel = 1

    beam_ad = True
    beam = app.property("ToolSvc.GtBeamerTool")
    beam.OutputLevel = 1
    beam.Radius = 5*units.cm
    beam.ParticlesPerEvent = 10000
    beam.ParticleName = "opticalphoton"
    beam.Momentum = 2.5*units.eV
    # move from face-RoC to PMT equator
    beam.TargetOffset = [0,0,100.0*units.mm]
    beam.SourceDirection = [0,0,-1.0]
    beam.SourceDistance = 2.0*units.meter

    #beam.Radius=10*units.cm
    beam.TargetOffset = [0, 0, 60*units.cm]
    beam.SourceDistance = 1.0*units.meter 


    if beam_ad:
        beam.TargetElement = "/dd/Structure/AD/far-oil2/far-ad2-ring5-column1"
    else:
        beam.TargetElement = "/dd/Structure/Pool/far-iws/far-InnFin-wall9-id29"
        ang=0*units.degree
        pass



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

    gen = app.algorithm("GenAlg")
    gen.OutputLevel = 1
    gen.GenTools = [ "GtBeamerTool", "GtTimeratorTool" ]
    gen.GenName = "Beam Beam"

    return

def config_sim(app):

    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 = False

    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"]

    # 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"

    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/AdDetails/AdSurfacesAll",
                         "/dd/Geometry/AdDetails/AdSurfacesFar",
                         "/dd/Geometry/PoolDetails/FarPoolSurfaces",
                         "/dd/Geometry/PoolDetails/PoolSurfacesAll",
                         ]

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

    pull = app.algorithm("PullEvent")

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


config_geom()

from  GaudiPython import AppMgr

app = AppMgr()
app.TopAlg = []


config_kin(app)
config_sim(app)


app.EvtSel = "NONE"


app.run(1)

---