In This Package:

Beamer.py

Go to the documentation of this file.

#!/usr/bin/env python

print "\tBuilding detector\n"
import xmldetdesc
xmldetdesc.config()

import os,sys
from math import sin, cos, pi
from  GaudiPython import AppMgr
from GaudiKernel import SystemOfUnits as units

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

eds=app.service("EventDataSvc")
eds.ForceLeaves=True

volume = "/dd/Structure/Pool/far-ows"


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

usebeam=True

if usebeam:
    gen = app.algorithm("GenAlg")
#    gen.OutputLevel = 4
    gen.GenTools = [ "GtBeamerTool", "GtTimeratorTool" ]
    gen.GenName = "Test Bang"
    gen.Location = "/Event/Gen/GenHeader" # this is default anyways


    # Set up beam
    ang=45*units.degree
    beam = app.property("ToolSvc.GtBeamerTool")
    #beam.OutputLevel = 1
    beam.Radius=0.11*units.meter
    beam.ParticlesPerEvent = 50000
    beam.ParticleName = "opticalphoton"
    beam.Momentum = 2.5*units.eV
    beam.TargetElement = "/dd/Structure/Pool/far-iws/far-InnFin-wall9-id29"
    beam.TargetOffset = [0, 0, 0.056*units.meter]
    beam.SourceDirection = [sin(ang),0,-cos(ang)]
    beam.SourceDistance = 0.3*units.meter 

else:
    gen = app.algorithm("GenAlg")
    gen.GenTools = [ "GtHepEvtGenTool", "GtTimeratorTool", "GtTransformTool" ]
    gen.GenName = "Test Bang"
    gen.Location = "/Event/Gen/GenHeader" # this is default anyways


    hep = app.property("ToolSvc.GtHepEvtGenTool")
    #hep.OutputLevel=3
    hep.HepEvtDataSource="/mnt/scratch/02/dayabay/kevin/NuWa-trunk/dybgaudi/Generators/Muon/Linux-i686/Muon.exe -n 1000 -s DYB -music_dir /mnt/scratch/02/dayabay/kevin/NuWa-trunk/dybgaudi/Generators/Muon/data|"

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


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


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.CerenPhotonScaleWeight=3.0
#optical.ScintPhotonScaleWeight=3.0
#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"]

TH2DE="TH2DE"       # TH2DE is the fastest
historian = app.property("GiGa.GiGaStepActionSequence.HistorianStepAction")
historian.TouchableToDetelem = TH2DE
historian.TrackSelection="pdg==13 or pdg==-13"
#historian.VertexSelection="any"

app.property("GiGa.GiGaStepActionSequence.UnObserverStepAction").Stats=[
    ["MuonTrkLengthInOws","dx","(pdg==13 or pdg==-13) and MaterialName == '/dd/Materials/OwsWater'"],
    ["MuonTrkLengthInIws", "dx", " (pdg==13 or pdg==-13) and MaterialName == '/dd/Materials/IwsWater'"],
    ["MuonTrkLengthInLS", "dx", "(pdg==13 or pdg==-13) and MaterialName == '/dd/Materials/LiquidScintillator'"],
    ["MuonTrkLengthInGdLS","dx", "(pdg==13 or pdg==-13) and MaterialName == '/dd/Materials/GdDopedLS'"],
    ["dEInn","dE", "(pdg!=20022) and MaterialName == '/dd/Materials/IwsWater'"],  
    ["dEOut","dE", "(pdg!=20022) and MaterialName == '/dd/Materials/OwsWater'"],
    ["MuonStop", "dx", "(pdg==13 or pdg==-13) and (IsStopping>0)"],    
    ]

# Make Geant4 sing!
ggrm = app.property("GiGa.GiGaMgr")
ggrm.Verbosity = 0
event_ac_cmds = app.property("GiGa.GiGaEventActionCommand")
verbosity_cmds = [
    "/control/verbose 2",
    "/run/verbose 1",
    "/event/verbose 2",
    "/tracking/verbose 2",
    "/geometry/navigator/verbose 2"
    ]
quiet_cmds = [
    "/control/verbose 0",
    "/run/verbose 0",
    "/event/verbose 0",
    "/tracking/verbose 0",
    "/geometry/navigator/verbose 0"
    ]
event_ac_cmds.BeginOfEventCommands = quiet_cmds
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"

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


histsvc = app.service("THistSvc")
histsvc.Output =["file1 DATAFILE='beam.root' OPT='RECREATE' TYP='ROOT' "]

# Add the GenHist algorithm
#app.TopAlg += ["ValiNTuple/vHalg"]
app.TopAlg += ["DetSimVali/vHalg"]
vali = app.algorithm("vHalg")
vali.Volume=volume
#vali.OutPutLevel=1

app.run(1)
    

---