In This Package:

elecIoTest.py

Go to the documentation of this file.

#!/usr/bin/env python

use_aes = True

import GaudiPython as gp

if use_aes:

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



import xmldetdesc
xmldetdesc.config()

app = gp.AppMgr(outputlevel=3)
#app.OutputLevel = 3

msg = app.service("MessageSvc")
#msg.Format = "\--%-F%40W%S%7W%R%T %23W%t\n \-> %0W%M"
msg.Format = "% F%25W%S%7W%R%T %0W%M"
msg.useColors = True
msg.fatalColorCode=['red','white']
msg.errorColorCode=['red']
msg.warningColorCode=['yellow']
msg.debugColorCode=['blue']
msg.verboseColorCode=['cyan']


app.EvtMax = 1

app.EvtSel = "NONE"

app.ExtSvc += [ "DybStorageSvc" ]
dss = app.service("DybStorageSvc")
dss.OutputLevel = 3
    
app.ExtSvc += [ "RootIOCnvSvc" ]

per = app.service("EventPersistencySvc")
per.CnvServices = [ "RootIOCnvSvc" ];

eds = app.service("EventDataService")
eds.OutputLevel = 3

rio = app.property("RootIOCnvSvc")
rio.OutputLevel = 3

iomap = { "/Event/Gen/GenHeader": "gen.root",
          "/Event/Elec/ElecHeader": "elec.root",
          "default": "default.root" };


rio.OutputStreams = iomap

app.TopAlg = [ ]

# set up gentools

import xmldetdesc
xmldetdesc.config()

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

# time
import GaudiKernel.SystemOfUnits as units
tim = app.property("ToolSvc.GtTimeratorTool")
tim.LifeTime = int(1*units.second)

# position
poser = app.property("ToolSvc.GtPositionerTool")
poser.Volume = volume
poser.Mode = "Fixed"
poser.Position = [0,0,0]

# momentum
gun = app.property("ToolSvc.GtGunGenTool")
gun.ParticlesPerEvent = 1
gun.ParticleName = "e+"
gun.Momentum = 3*units.MeV
gun.MomentumMode = "Fixed"
gun.DirectionMode = "Fixed"
gun.Direction = [ 1,0,0 ]

# translate from local volume to global
trans = app.property("ToolSvc.GtTransformTool")
trans.Volume = volume

app.TopAlg += [ 'GtGenerator/gen' ]
gen = app.algorithm("gen")
gen.GenTools = [ "GtGunGenTool", "GtPositionerTool", "GtTimeratorTool", "GtTransformTool" ]
gen.GenName = "Bang Bang"

app.ExtSvc += ["GiGa"]

modularPL = app.property("GiGa.GiGaPhysListModular")
#modularPL.OutputLevel = 3
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

giga = app.service("GiGa")
giga.OutputLevel = 3
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"]


TH2DE="TH2DE"       # TH2DE is the fastest
#TH2DE="TouchableToDetectorElementFast"
historian = app.property("GiGa.GiGaStepActionSequence.HistorianStepAction")
historian.TouchableToDetelem = TH2DE

params = {
    'track1':"(id == 1 and ProcessType == 1)",
    'track2':"(id == 2 and ProcessType == 1)",
#    'track1':"(id == 1)",
#    'track2':"(id == 2)",
    'GD':    "MaterialName == '/dd/Materials/GdDopedLS'",
    'LS':    "MaterialName == '/dd/Materials/LiquidScintillator'",
    'oil':   "MaterialName == '/dd/Materials/MineralOil'",
    'iAV':   "MaterialName == '/dd/Materials/Acrylic'",
    'oAV':   "MaterialName == '/dd/Materials/Acrylic'"
    }
unobs = app.property("GiGa.GiGaStepActionSequence.UnObserverStepAction")
unobs.TouchableToDetelem = TH2DE
unobs.Stats=[
    ["pdgId_Trk1","pdg","%(track1)s"%params],
    ["t_Trk1",    "t" , "%(track1)s"%params],
    ["x_Trk1",    "lx", "%(track1)s"%params],
    ["y_Trk1",    "ly", "%(track1)s"%params],
    ["z_Trk1",    "lz", "%(track1)s"%params],
    ["e_Trk1",    "E",  "%(track1)s"%params],
    ["p_Trk1",    "p",  "%(track1)s"%params],
    ["ke_Trk1",   "KE", "%(track1)s"%params],
    ["vx_Trk1",   "vx", "%(track1)s"%params],
    ["vy_Trk1",   "vy", "%(track1)s"%params],
    ["vz_Trk1",   "vz", "%(track1)s"%params],
    ["TrkLength_GD_Trk1",  "dx","%(track1)s and %(GD)s"%params],
    ["TrkLength_iAV_Trk1", "dx","%(track1)s and %(iAV)s"%params],
    ["TrkLength_LS_Trk1",  "dx","%(track1)s and %(LS)s"%params],
    ["TrkLength_oAV_Trk1", "dx","%(track1)s and %(oAV)s"%params],
    ["TrkLength_Oil_Trk1", "dx","%(track1)s and %(oil)s"%params]
    ]
                                                

# Make Geant4 sing!
ggrm = app.property("GiGa.GiGaMgr")
ggrm.Verbosity = 3
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 = [
    "/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"


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 = [ "GiGaInputStream/GGInStream", "DsPushKine/PushKine", "DsPullEvent/PullEvent" ]
pull = app.algorithm("PullEvent")
#pull.OutputLevel = 1 
push = app.algorithm("PushKine")
push.Converter = "HepMCtoG4"
#push.Location = "/Event/Gen/GenHeader" # default anyway

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

app.TopAlg +=["ESFEEAlg"]

fee=app.algorithm("ESFEEAlg")
fee.OutputLevel = 1;
fee.ElecHeaderPath = '/Event/Elec/ElecHeader'

feeTool=app.tool("ToolSvc.ESFeeTool")
feeTool.OutputLevel = 1;


app.TopAlg += [ 'RegSeqDumpAlg/rsd' ]
rsd = app.algorithm("rsd")
rsd.OutputLevel = 3

app.TopAlg += [ 'DybStoreAlg/dsa' ]
dsa = app.algorithm("dsa")
dsa.OutputLevel = 3

# Run...
app.initialize()
app.run(app.EvtMax)

---