In This Package:

ranFarMuon.py

Go to the documentation of this file.

#!/usr/bin/env python

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

import os,sys
from  GaudiPython import AppMgr
from GaudiKernel import SystemOfUnits as units

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


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

gen = app.algorithm("GenAlg")
#gen.OutputLevel = 3

gen.GenTools = [ "GtHepEvtGenTool", "GtTimeratorTool", "GtTransformTool" ]
gen.GenName = "Test Bang"
gen.Location = "/Event/Gen/GenHeader" # this is default anyways

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

trans = app.property("ToolSvc.GtTransformTool")
trans.Volume = volume
trans.Offset=[0., 0., (0.042)*units.meter]

# Set up gun
hep = app.property("ToolSvc.GtHepEvtGenTool")
#hep.OutputLevel=3
hep.HepEvtDataSource="/mnt/scratch/02/dayabay/kevin/NuWa-trunk/dybkevin/Generators/Muon/Linux-i686/Muon.exe -n 1000 -s Far -r Yes -music_dir /mnt/scratch/02/dayabay/kevin/NuWa-trunk/dybkevin/Generators/Muon/data/NewMuonGenerator/data|"

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.StackingAction="DsOpStackAction"

giga.SteppingAction = "GiGaStepActionSequence"
stepseq = app.property("GiGa.GiGaStepActionSequence")
stepseq.Members = ["HistorianStepAction","UnObserverStepAction"]


TH2DE="TH2DE"       # TH2DE is the fastest
#TH2DE="TouchableToDetectorElementFast"
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'"],
    ["MuonStop", "dx", "(pdg==13 or pdg==-13) and (IsStopping>0)"],    
    
    ]


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


# 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 = 3
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 = 3


histsvc = app.service("THistSvc")
p1='/mnt/scratch/04/dayabay/kevin/FastTest/Far'+sys.argv[1]+'.root'
outset="file1 DATAFILE='%s' OPT='RECREATE' TYP='ROOT' "%p1
histsvc.Output =[outset]

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

app.run(200)
    

---