In This Package:

muon.py

Go to the documentation of this file.

#!/usr/bin/env python

from  GaudiPython import AppMgr
from GaudiKernel import SystemOfUnits as units

import os,sys

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

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


volume = "/dd/Structure/AD/far-oil2"


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

# Set up positioner
poser = app.property("ToolSvc.GtPositionerTool")
#poser.OutputLevel = 1
poser.Strategy = "FullVolume"
poser.Volume = volume
poser.Mode = "Fixed"
#poser.Mode = "Smeared"
#poser.Spread = 1*units.meter
poser.Position = [0,10,0.0*units.meter]

# Set up gun
gun = app.property("ToolSvc.GtGunGenTool")
#gun.OutputLevel = 1
gun.Volume = volume
gun.ParticlesPerEvent = 1
gun.ParticleName = "mu+"
gun.Momentum = 10*units.GeV
gun.MomentumMode = "Fixed"
#gun.MomentumSpread = 10*units.GeV
gun.DirectionMode = "Fixed"

gun.Direction = [ -1,0,0 ] # aim for PMT 
print 'gun.Direction=',gun.Direction

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

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


gen = app.algorithm("GenAlg")
#gen.OutputLevel = 1
gen.GenTools = [ "GtGunGenTool", "GtPositionerTool", "GtTimeratorTool", "GtTransformTool" ]
gen.GenName = "Bang Bang"

#print " GtDumper"
gendump = app.algorithm("GenDump")

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

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

app.property("GiGa.GiGaStepActionSequence.UnObserverStepAction").Stats=[
    ["MuonTrkLengthInOws","dx","(pdg==13 or pdg==-13) and (MaterialName == '/dd/Materials/Water' and (Volume == '/dd/Structure/Pool/far-ows' or Volume == '/dd/Structure/Pool/db-ows'))"],
    ["MuonTrkLengthInIws", "dx", " (pdg==13 or pdg==-13) and ((Volume == '/dd/Structure/Pool/far-iws' or Volume == '/dd/Structure/Pool/db-iws') and MaterialName == '/dd/Materials/Water')"],
    ["MuonTrkLengthInLS", "dx", "(pdg==13 or pdg==-13) and MaterialName == '/dd/Materials/LiquidScintillator'"],
    ["MuonTrkLengthInGdLS","dx", "(pdg==13 or pdg==-13) and MaterialName == '/dd/Materials/GdDopedLS'"]
    ]


# 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 1",
    "/run/verbose 1",
    "/event/verbose 1",
    "/tracking/verbose 1",
    "/geometry/navigator/verbose 1"
    ]
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/DayaBay", ]
ggin.StreamItems = [ "/dd/Structure/Sites/far-rock", ]

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='simvali1.root' OPT='RECREATE' TYP='ROOT' "]

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

    
app.run(1)
    

---