from Gaudi.Configuration import *
pileupFilenames = ["dummyEventData2.root", "dummyEventData1.root", "dummyEventData3.root",]
# Data service
from Configurables import FCCDataSvc
podioevent = FCCDataSvc("EventDataSvc")
from Configurables import PileupParticlesMergeTool
particlemergetool = PileupParticlesMergeTool("MyPileupParticlesMergeTool")
particlemergetool.genParticlesBranch = "genParticles"
particlemergetool.genVerticesBranch = "genVertices"
particlemergetool.signalGenParticles.Path = "genParticles"
particlemergetool.signalGenVertices.Path = "genVertices"
particlemergetool.mergedGenParticles.Path = "overlaidGenParticles"
particlemergetool.mergedGenVertices.Path = "overlaidGenVertices"
from Configurables import PileupTrackHitMergeTool
trackhitsmergetool = PileupTrackHitMergeTool("MyTrackHitsMergeTool")
trackhitsmergetool.pileupHitsBranch = "trackHits"
trackhitsmergetool.pileupPositionedHitsBranch = "positionedTrackHits"
trackhitsmergetool.signalHits = "trackHits"
trackhitsmergetool.signalPositionedHits = "positionedTrackHits"
trackhitsmergetool.mergedHits = "overlaidTrackHits"
trackhitsmergetool.mergedPositionedHits = "overlaidPositionedTrackHits"
from Configurables import ConstPileUp
pileuptool = ConstPileUp("MyPileupTool", numPileUpEvents=12)
from Configurables import PileupOverlayAlg
overlay = PileupOverlayAlg()
path_to_detector +
'/Detector/DetFCChhHCalTile/compact/FCChh_HCalBarrel_TileCal.xml'
]
from Configurables import GeoSvc
geoservice = GeoSvc("GeoSvc", detectors=detectors_to_use, OutputLevel=WARNING)
import random
seed = int(filter(str.isdigit, output_name))
print 'seed : ', seed
random.seed(seed)
random.shuffle(pileupFilenames)
# edm data from generation: particles and vertices
from Configurables import PileupParticlesMergeTool
particlemergetool = PileupParticlesMergeTool("ParticlesMerge")
# branchnames for the pileup
particlemergetool.genParticlesBranch = "GenParticles"
particlemergetool.genVerticesBranch = "GenVertices"
# branchnames for the signal
particlemergetool.signalGenParticles.Path = "GenParticles"
particlemergetool.signalGenVertices.Path = "GenVertices"
# branchnames for the output
particlemergetool.mergedGenParticles.Path = "pileupGenParticles"
particlemergetool.mergedGenVertices.Path = "pileupGenVertices"
# edm data from simulation: hits and positioned hits
from Configurables import PileupCaloHitMergeTool
ecalbarrelmergetool = PileupCaloHitMergeTool("ECalBarrelHitMerge")
ecalbarrelmergetool.pileupHitsBranch = prefix + "ECalBarrelCells"
ecalbarrelmergetool.signalHits = "ECalBarrelCells"
collectionNames["trajectoryPoints"] = "trajectoryPoints"
if not args.noSignal:
from Configurables import PodioInput
podioinput = PodioInput("PodioReader",
collections=collectionNames.values(),
OutputLevel=DEBUG,
)
podioinput.AuditExecute = True
algList += [podioinput]
print "# pileup: ", args.npileup
if args.npileup > 0:
pileupFilenames = args.pileupfiles
from Configurables import PileupParticlesMergeTool
particlemergetool = PileupParticlesMergeTool("PileupGenParticlesMergeTool")
# branchnames for the pileup
particlemergetool.genParticlesBranch = "GenParticles"
particlemergetool.genVerticesBranch = "GenVertices"
# branchnames for the signal
particlemergetool.signalGenVertices.Path = "GenVertices"
particlemergetool.signalGenParticles.Path = "GenParticles"
# branchnames for the output
particlemergetool.mergedGenParticles.Path = "mergedGenParticles"
particlemergetool.mergedGenVertices.Path = "mergedGenVertices"
simparticlemergetool = PileupParticlesMergeTool("PileupSimParticlesMergeTool")
# branchnames for the pileup
simparticlemergetool.genParticlesBranch = "SimParticles"
simparticlemergetool.genVerticesBranch = "SimVertices"
# branchnames for the signal
from Gaudi.Configuration import *
pileupFilenames = [
"output_dummyEventData2.root",
"output_dummyEventData1.root",
"output_dummyEventData3.root",
]
# Data service
from Configurables import FCCDataSvc
podioevent = FCCDataSvc("EventDataSvc")
from Configurables import PileupParticlesMergeTool
particlemergetool = PileupParticlesMergeTool("MyPileupParticlesMergeTool")
particlemergetool.genParticlesBranch = "genParticles"
particlemergetool.genVerticesBranch = "genVertices"
particlemergetool.signalGenParticles.Path = "genParticles"
particlemergetool.signalGenVertices.Path = "genVertices"
particlemergetool.mergedGenParticles.Path = "overlaidGenParticles"
particlemergetool.mergedGenVertices.Path = "overlaidGenVertices"
from Configurables import PileupTrackHitMergeTool
trackhitsmergetool = PileupTrackHitMergeTool("MyTrackHitsMergeTool")
trackhitsmergetool.pileupHitsBranch = "trackHits"
trackhitsmergetool.pileupPositionedHitsBranch = "positionedTrackHits"
trackhitsmergetool.signalHits = "trackHits"
trackhitsmergetool.signalPositionedHits = "positionedTrackHits"
trackhitsmergetool.mergedHits = "overlaidTrackHits"
trackhitsmergetool.mergedPositionedHits = "overlaidPositionedTrackHits"
from Configurables import ConstPileUp