In This Package:

sim15.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()

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

import GaudiPython as gp

iapp = gp.iService("ApplicationMgr")
iapp.SvcMapping = [
    'EvtDataSvc/EventDataArchiveSvc',
    'DybDataSvc/EventDataSvc',
    #"EvtDataSvc/EventDataSvc",
    "DetDataSvc/DetectorDataSvc",
    "HistogramSvc/HistogramDataSvc",
    "HbookCnv::PersSvc/HbookHistSvc",
    "RootHistCnv::PersSvc/RootHistSvc",
    "EvtPersistencySvc/EventPersistencySvc",
    "DetPersistencySvc/DetectorPersistencySvc",
    "HistogramPersistencySvc/HistogramPersistencySvc",
    ]
        

app.TopAlg = ["GiGaInputStream/GGInStream",
              'Sim15',
              ]

StageKin=app.property("ToolSvc.Kinematic")
StageKin.Sequencer='GnrtrSeq'
GenSeq=app.property('GnrtrSeq')
GenSeq.Members = ["Gnrtr/Gen1"]

StageDet=app.property("ToolSvc.Detector")
StageDet.Sequencer='DetSimSeq'
DSSeq=app.property('DetSimSeq')
DSSeq.Members = ['DetSimProc']

StageElec=app.property("ToolSvc.Electronic")
StageElec.Sequencer='ElecSimSeq'
ESSeq=app.property('ElecSimSeq')
ESSeq.Members = ['ElecSimProc']

StageTrig=app.property("ToolSvc.Trigger")
StageTrig.Sequencer='TrigSimSeq'
TSSeq=app.property('TrigSimSeq')
TSSeq.Members = ['TrigReadProc']

StageSingle=app.property("ToolSvc.SingleLoader")
StageSingle.Sequencer='SingleSeq'
SingleSeq=app.property('SingleSeq')
SingleSeq.Members = ['SingleLoader']


# Kinematic
gen=app.algorithm('Gen1')
gen.ThisStageName='Kinematic'    
gen.LowerStageName=''            

# Detector
sim=app.algorithm('DetSimProc')
sim.ThisStageName='Detector'     
sim.LowerStageName='Kinematic'   

# Electronic
elec=app.algorithm('ElecSimProc')
elec.ThisStageName='Electronic'  
elec.LowerStageName='Detector'   

# Trigger
trig=app.algorithm('TrigReadProc')
trig.ThisStageName='Trigger'     
trig.LowerStageName='Electronic' 

# SingleLoader
singll=app.algorithm('SingleLoader')
singll.ThisStageName='SingleLoader'
singll.LowerStageName='Trigger'

# Configure generator
#gen.OutputLevel = 1
gen.GenTools = [ "GtGunGenTool", "GtPositionerTool", "GtTimeratorTool", "GtTransformTool" ]
gen.GenName = "ticktack"
#gen.Location = "/Event/Gen/GenHeader" # this is default anyways
gen.OutputLevel=2

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

# Set up timerator
tim = app.property("ToolSvc.GtTimeratorTool")
tim.LifeTime = int(20000*units.nanosecond)
tim.Distribution = "Uniform"
#tim.Mode="Relative"
tim.OutputLevel=2

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

# Set up gun
gun = app.property("ToolSvc.GtGunGenTool")
#gun.OutputLevel = 1
gun.ParticlesPerEvent = 1
gun.ParticleName = "e-"
#gun.Momentum = 6*units.eV   # try to enhance rayleigh scattering
gun.Momentum = 2*units.MeV
gun.MomentumMode = "Smeared"
#gun.MomentumMode = "Fixed"
#gun.MomentumMode = "Uniform"
gun.MomentumSpread = 1.0*units.MeV

gun.DirectionMode = "Uniform"
#from math import sin, cos, pi
#pmt_column_number = 9
#angle = (2*pmt_column_number - 1)*pi/24.0;
#gun.Direction = [ cos(angle),sin(angle),0 ] # aim for PMT
gun.Direction = [ 1, 0, 0 ] # aim for a PMT
gun.DirectionSpread = 3
print 'gun.Direction=',gun.Direction

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

#-------------------------------------------------------------------------
# Configure geant4 simulation
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"
            ]
# in DsPhysConsOptical.cc these properties are declared.
optical = app.property("GiGa.GiGaPhysListModular.DsPhysConsOptical")
#
#optical.UseScintillation = False  # false option will kill all optical photons, should be true all the time
#optical.UseRayleigh = True       # enable rayleigh scattering
optical.CerenPhotonScaleWeight=4.0  # scale photon yield down and scale QE up
optical.ScintPhotonScaleWeight=4.0

giga = app.service("GiGa")
giga.OutputLevel = 4
giga.PhysicsList = "GiGaPhysListModular"

gggeo = app.service("GiGaGeo")
gggeo.OutputLevel = 3
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")
verbosity_cmds = [
        "/control/verbose 0",
        "/run/verbose 0",
        "/event/verbose 2",
        "/tracking/verbose 2",
        "/geometry/navigator/verbose 0"
        ]
silent_cmds = [
        "/particle/property/verbose 0",
        "/control/verbose 0",
        "/run/verbose 0",
        "/event/verbose 0",
        "/tracking/verbose 0",
        "/geometry/navigator/verbose 0"
        ]

event_ac_cmds.BeginOfEventCommands = silent_cmds
giga.EventAction = "GiGaEventActionCommand"


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

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

sim.Converter = "HepMCtoG4"
#push.Location = "/Event/Gen/GenHeader" # default anyway
sim.OutputLevel=2

# other properties for ElecSim
#elec.PmtTool = 'EsPmtEffectPulseTool'
elec.PmtTool = 'EsIdealPulseTool'
elec.RpcTool = 'EsIdealPulseTool'
elec.FeeTool = 'EsIdealFeeTool'
elec.FecTool = 'EsIdealFecTool'
elec.OutputLevel=2

# trigger tools
trig.TrigTools=['TsMultTriggerTool']  # default
trig.RoTools=['ROsFeeReadoutTool',
              #'ROsFecReadoutTool'
              ]
trig.OutputLevel=2
nhits=app.property("ToolSvc.TsMultTriggerTool")
nhits.OutputLevel=2
fee=app.property("ToolSvc.ROsFeePeakOnlyTool")
#fec=app.property("ToolSvc.ROsFecReadoutTool")
#fee.OutputLevel=1
#fec.OutputLevel=1


# Single Loader
singll.OutputLevel=2

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

#
app.EvtMax = 10
app.EvtSel ="NONE"

#
sim15=app.algorithm('Sim15')
#sim15.TopStage='Kinematic'
#sim15.TopStage='Detector'
#sim15.TopStage='Electronic'
#sim15.TopStage='Trigger'
sim15.TopStage='SingleLoader'
sim15.OutputLevel=2

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

msv = app.service('MessageSvc')
msv.OutputLevel = 3
app.initialize()
app.run( app.EvtMax )
app.exit()

---