In This Package:

simvali.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"

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

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


# 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.418*units.meter]

# Set up gun
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|"


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


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


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

app.property("GiGa.GiGaStepActionSequence.UnObserverStepAction").Stats=[
    ["trkLeng_ows","dx","MaterialName == '/dd/Materials/Water' and Volume == '/dd/Pool/lvFarPoolOWS'"],
    ["trkLeng_iws","dx","MaterialName == '/dd/Materials/Water' and Volume == '/dd/Pool/lvFarPoolIWS'"]
    ]

# Make Geant4 sing!
ggrm = app.property("GiGa.GiGaMgr")
ggrm.Verbosity = 9
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='simvali.root' OPT='RECREATE' TYP='ROOT' "]

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

    
app.run(1)
    

---