mergelayersHcalFwd.inhits.Path = "HCalFwdHits"
mergelayersHcalFwd.outhits.Path = "mergedHCalFwdHits"
createHcalFwdCells = CreateCaloCells("CreateHcalFwdCaloCells",
doCellCalibration=True,
calibTool=calibHcalFwd,
addCellNoise=False,
filterCellNoise=False,
OutputLevel=INFO)
createHcalFwdCells.hits.Path = "mergedHCalFwdHits"
createHcalFwdCells.cells.Path = "HCalFwdCells"
out = PodioOutput("out", OutputLevel=INFO)
out.outputCommands = [
"drop *", "keep ECalBarrelCells", "keep ECalEndcapCells",
"keep ECalFwdCells", "keep HCalBarrelCells", "keep HCalExtBarrelCells",
"keep HCalEndcapCells", "keep HCalFwdCells", "keep GenParticles",
"keep GenVertices"
]
out.filename = "output_fullCalo_SimAndDigi_e50GeV_" + str(
num_events) + "events.root"
#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
geantsim.AuditExecute = True
createEcalBarrelCellsStep1.AuditExecute = True
positionsEcalBarrel.AuditExecute = True
resegmentEcalBarrel.AuditExecute = True
createEcalBarrelCells.AuditExecute = True
#Fill a collection of CaloHitPositions for detailed Cluster analysis
from Configurables import CreateCaloCellPositions
positionsClusterBarrel =CreateCaloCellPositions("positionsClusterBarrel",
positionsECalBarrelTool = ECalBcells,
positionsHCalBarrelTool = HCalBcells,
positionsHCalExtBarrelTool = HCalExtBcells,
positionsEMECTool = EMECcells,
positionsHECTool = HECcells,
positionsEMFwdTool = ECalFwdcells,
positionsHFwdTool = HCalFwdcells,
hits = "caloClusterBarrelCells",
positionedHits = "caloClusterBarrelCellPositions",
OutputLevel = INFO)
out = PodioOutput("out", filename = "output_BarrelTopo_electrNoise_50GeVe_3ev.root", OutputLevel = DEBUG)
out.outputCommands =["drop *", "keep GenParticles", "keep GenVertices", "keep caloClustersBarrel","keep caloClusterBarrelCells", "keep caloClusterBarrelCellPositions"]
#CPU information
from Configurables import AuditorSvc,ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors =[chra]
podioinput.AuditExecute = True
createemptycells.AuditExecute =True
createEcalBarrelCells.AuditExecute = True
createHcalBarrelCells.AuditExecute =True
createTopoClusters.AuditExecute = True
positionsClusterBarrel.AuditExecute = True
ApplicationMgr(TopAlg =
[podioinput,
positionedHits = "HCalEndcapCellPositions",
OutputLevel = INFO)
positionsEcalFwd = CreateCellPositions("positionsEcalFwd",
positionsTool=ECalFwdcells,
hits = "ECalFwdCells",
positionedHits = "ECalFwdCellPositions",
OutputLevel = INFO)
positionsHcalFwd = CreateCellPositions("positionsHcalFwd",
positionsTool=HCalFwdcells,
hits = "HCalFwdCells",
positionedHits = "HCalFwdCellPositions",
OutputLevel = INFO)
out = PodioOutput("out", OutputLevel=DEBUG)
out.filename = "digi_cellPositions_50GeVelectrons.root"
out.outputCommands = ["keep *","drop ECalBarrelCells","drop ECalEndcapCells","drop ECalFwdCells","drop HCalBarrelCells", "drop HCalExtBarrelCells", "drop HCalEndcapCells", "drop HCalFwdCells"]
#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
podioinput.AuditExecute = True
positionsEcalBarrel.AuditExecute = True
positionsEcalEndcap.AuditExecute = True
positionsEcalFwd.AuditExecute = True
positionsHcalBarrel.AuditExecute = True
positionsHcalExtBarrel.AuditExecute = True
positionsHcalEndcap.AuditExecute = True
positionsHcalFwd.AuditExecute = True
out.AuditExecute = True
savecaltool.DataOutputs.caloHits.Path = "caloHits"
# next, create the G4 algorithm, giving the list of names of tools ("XX/YY")
particle_converter = SimG4PrimariesFromEdmTool("EdmConverter")
particle_converter.DataInputs.genParticles.Path = "allGenParticles"
geantsim = SimG4Alg("SimG4Alg",
outputs = ["SimG4SaveSmearedParticles/saveSmearedParticles", "SimG4SaveCalHits/saveCalHits"],
eventProvider=particle_converter)
from Configurables import SimG4FastSimHistograms
hist = SimG4FastSimHistograms("fastHist")
hist.DataInputs.particlesMCparticles.Path = "particleMCparticleAssociation"
THistSvc().Output = ["rec DATAFILE='histFormula.root' TYP='ROOT' OPT='RECREATE'"]
THistSvc().PrintAll=True
THistSvc().AutoSave=True
THistSvc().AutoFlush=True
THistSvc().OutputLevel=INFO
from Configurables import PodioOutput
## PODIO algorithm
out = PodioOutput("out", filename = "out_fast_tracker_formula_calo_gflash.root")
out.outputCommands = ["keep *"]
# ApplicationMgr
from Configurables import ApplicationMgr
ApplicationMgr( TopAlg = [reader, hepmc_converter, geantsim, hist, out],
EvtSel = 'NONE',
EvtMax = 1,
# order is important, as GeoSvc is needed by SimG4Svc
ExtSvc = [podioevent, geoservice, geantservice],
OutputLevel=INFO)
createClusters.nEtaDuplicates = 7
createClusters.nPhiDuplicates = 13
createClusters.nEtaFinal = 9
createClusters.nPhiFinal = 17
# Minimal energy to create a cluster in GeV (FCC-ee detectors have to reconstruct low energy particles)
createClusters.energyThreshold = 0.1
createClusters.attachCells = True
createClusters.clusters.Path = "CaloClusters"
createClusters.clusterCells.Path = "CaloClusterCells"
ApplicationMgr().TopAlg += [createClusters]
################ Output
from Configurables import PodioOutput
out = PodioOutput("out")
out.outputCommands = ["keep *", "drop ECalBarrelHits", "drop HCal*", "drop ECalBarrelCellsStep*", "drop ECalBarrelPositionedHits", "drop emptyCaloCells", "drop CaloClusterCells"]
import uuid
out.filename = "output_fullCalo_SimAndDigi.root"
ApplicationMgr().TopAlg += [out]
#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
gen.AuditExecute = True
hepmc_converter.AuditExecute = True
geantsim.AuditExecute = True
createEcalBarrelCellsStep1.AuditExecute = True
resegmentEcalBarrel.AuditExecute = True
createEcalBarrelCells.AuditExecute = True
hepmc_converter.genparticles.Path="genParticles"
hepmc_converter.genvertices.Path="genVertices"
## FCC event-data model output -> define objects to be written out
from Configurables import PodioOutput
out = PodioOutput("out",OutputLevel=messageLevelOut)
out.filename = "FCCDelphesOutput.root"
if args.outputfile != '':
out.filename = args.outputfile
#out.outputCommands = ["drop *",
# "keep genParticles",
# "keep genVertices",
# "keep genJets",
# "keep genJetsToMC"]
out.outputCommands = ["keep *", "drop genParticles", "drop genVertices"]
############################################################
#
# Run modules
#
############################################################
# Run Pythia + Delphes
ApplicationMgr( TopAlg = [ pythia8gen, hepmc_converter, delphessim, out ],
EvtSel = 'NONE',
EvtMax = nEvents,
ExtSvc = [podioEvent])
ApplicationMgr.EvtMax=10
genfilter.allGenParticles.Path = "allGenParticles"
genfilter.filteredGenParticles.Path = "GenParticles"
from Configurables import SimG4PrimariesFromEdmTool
particle_converter = SimG4PrimariesFromEdmTool("EdmConverter")
particle_converter.genParticles.Path = "GenParticles"
geantsim.eventProvider = particle_converter
# PODIO algorithm
from Configurables import ApplicationMgr, FCCDataSvc, PodioOutput
podioevent = FCCDataSvc("EventDataSvc")
out = PodioOutput("out")
out.outputCommands = [
"drop *",
"keep GenParticles",
"keep GenVertices",
"keep SimParticles",
"keep SimVertices",
"keep TrackerHits",
"keep TrackerPositionedHits",
"keep TrackerDigiPostPoint",
]
out.filename = output_name
#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
geantsim.AuditExecute = True
out.AuditExecute = True
list_of_algorithms = [genalg, hepmc_converter, genfilter, geantsim, out]
merge = hcalFwdNumberOfLayersToMerge,
OutputLevel = INFO)
mergelayersHcalFwd.inhits.Path = "HCalFwdHits"
mergelayersHcalFwd.outhits.Path = "mergedHCalFwdHits"
createHcalFwdCells = CreateCaloCells("CreateHcalFwdCaloCells",
doCellCalibration=True,
calibTool=calibHcalFwd,
addCellNoise=False, filterCellNoise=False,
OutputLevel=INFO)
createHcalFwdCells.hits.Path="mergedHCalFwdHits"
createHcalFwdCells.cells.Path="HCalFwdCells"
out = PodioOutput("out",
OutputLevel=INFO)
out.outputCommands = ["drop *", "keep ECalBarrelCells", "keep ECalEndcapCells", "keep ECalFwdCells", "keep HCalBarrelCells", "keep HCalBarrelPositions", "keep HCalExtBarrelCells", "keep HCalEndcapCells", "keep HCalFwdCells", "keep GenParticles","keep GenVertices"]
out.filename = "output_fullCalo_SimAndDigi_pi50GeV_"+str(num_events)+"events.root"
#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
geantsim.AuditExecute = True
createEcalBarrelCellsStep1.AuditExecute = True
positionsEcalBarrel.AuditExecute = True
resegmentEcalBarrel.AuditExecute = True
createEcalBarrelCells.AuditExecute = True
createEcalEndcapCells.AuditExecute = True
createEcalFwdCells.AuditExecute = True
createHcalCells.AuditExecute = True
hepmc_converter.genparticles.Path="genParticles"
hepmc_converter.genvertices.Path="genVertices"
## FCC event-data model output -> define objects to be written out
from Configurables import PodioOutput
out = PodioOutput("out",OutputLevel=messageLevelOut)
out.filename = "FCCDelphesOutput.root"
if args.outputfile != '':
out.filename = args.outputfile
#out.outputCommands = ["drop *",
# "keep genParticles",
# "keep genVertices",
# "keep genJets",
# "keep genJetsToMC"]
out.outputCommands = ["keep *", "drop genParticles", "drop genVertices"]
############################################################
#
# Run modules
#
############################################################
# Run Pythia + Delphes
ApplicationMgr( TopAlg = [ pythia8gen, hepmc_converter, delphessim, out ],
EvtSel = 'NONE',
EvtMax = nEvents,
ExtSvc = [podioEvent])
ApplicationMgr.EvtMax=10
]
overlay.PileUpTool = pileuptool
from Configurables import CreateExampleEventData
producer = CreateExampleEventData()
producer.magicNumberOffset = 300
producer.trackhits.Path = "trackHits"
producer.positionedtrackhits.Path = "positionedTrackHits"
producer.calohits.Path = "caloHits"
producer.positionedcalohits.Path = "positionedCaloHits"
producer.genparticles.Path = "genParticles"
producer.genvertices.Path = "genVertices"
# PODIO algorithm
from Configurables import PodioOutput
out = PodioOutput("out", OutputLevel=DEBUG)
out.outputCommands = ["drop *", "keep overlaid*"]
out.filename = "dummyEventDataOverlaid.root"
# ApplicationMgr
from Configurables import ApplicationMgr
ApplicationMgr(
TopAlg=[producer, overlay, out],
EvtSel='NONE',
EvtMax=300,
ExtSvc=[
podioevent,
],
OutputLevel=DEBUG,
)
OutputLevel = INFO)
positionsHcalFwd = CreateCellPositions("positionsHcalFwd",
positionsTool=HCalFwdcells,
hits = prefix+"HCalFwdCells",
positionedHits = "HCalFwdCellPositions",
OutputLevel = INFO)
if addMuons:
positionsTailCatcher = CreateCellPositions("positionsTailCatcher",
positionsTool=TailCatchercells,
hits = "TailCatcherCells",
positionedHits = "TailCatcherCellPositions",
OutputLevel = INFO)
# PODIO algorithm
out = PodioOutput("out", OutputLevel=DEBUG)
out.outputCommands = ["keep *","drop "+prefix+"ECalBarrelCells","drop "+prefix+"ECalEndcapCells","drop "+prefix+"ECalFwdCells","drop "+prefix+"HCalBarrelCells", "drop "+prefix+"HCalExtBarrelCells", "drop "+prefix+"HCalEndcapCells", "drop "+prefix+"HCalFwdCells"]
if addMuons:
out.outputCommands += ["drop TailCatcherCells"]
out.filename = output_name
#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
podioinput.AuditExecute = True
positionsEcalBarrel.AuditExecute = True
positionsEcalEndcap.AuditExecute = True
positionsEcalFwd.AuditExecute = True
positionsHcalBarrel.AuditExecute = True
positionsHcalExtBarrel.AuditExecute = True
delphessim.DataOutputs.recMETs.Path = "recMETs"
delphessim.DataOutputs.recMuonsToMC.Path = "recMuonsToMC"
delphessim.DataOutputs.recElectronsToMC.Path = "recElectronsToMC"
delphessim.DataOutputs.recChargedToMC.Path = "recChargedToMC"
delphessim.DataOutputs.recNeutralToMC.Path = "recNeutralToMC"
delphessim.DataOutputs.recPhotonsToMC.Path = "recPhotonsToMC"
delphessim.DataOutputs.recJetsToMC.Path = "recJetsToMC"
## FCC event-data model output -> define objects to be written out
from Configurables import PodioOutput
out = PodioOutput("out", OutputLevel=messageLevelOut)
out.filename = "FCCDelphesOutput.root"
out.outputCommands = [
"drop *", "keep genParticles", "keep genVertices", "keep recMuons",
"keep recElectrons", "keep recCharged", "keep recPhotons", "keep recJets",
"keep recMETs", "keep recMuonsToMC", "keep recElectronsToMC",
"keep recChargedToMC", "keep recPhotonsToMC", "keep recJetsToMC"
]
#out.outputCommands = ["keep *"]
############################################################
#
# Run modules
#
############################################################
# Run Pythia + Delphes
if delphesHepMCInFile == "":
ApplicationMgr(TopAlg=[pythia8gen, delphessim, out],
EvtSel='NONE',
EvtMax=nEvents,
hepmc_converter2 = HepMCToEDMConverter("Converter2")
hepmc_converter2.hepmc.Path = "hepmc2"
hepmc_converter2.genparticles.Path = "pileupGenParticles"
hepmc_converter2.genvertices.Path = "pileupGenVertices"
from Configurables import GenMerge
mergealg = GenMerge()
mergealg.pileupGenParticles.Path = "pileupGenParticles"
mergealg.pileupGenVertices.Path = "pileupGenVertices"
mergealg.signalGenParticles.Path = "signalGenParticles"
mergealg.signalGenVertices.Path = "signalGenVertices"
mergealg.allGenParticles.Path = "allGenParticles"
mergealg.allGenVertices.Path = "allGenVertices"
from Configurables import PodioOutput
out = PodioOutput("out", filename="output_gen.root")
out.outputCommands = [
"keep *", "drop pileupGenVertices", "drop pileupGenParticles",
"drop signalGenVertices", "drop signalGenParticles"
]
ApplicationMgr(
TopAlg=[gun1, gun2, hepmc_converter1, hepmc_converter2, mergealg, out],
EvtSel='NONE',
ExtSvc=[podioevent],
EvtMax=1,
OutputLevel=VERBOSE,
)
delphessim.DataOutputs.recNeutralToMC.Path = "recNeutralToMC"
delphessim.DataOutputs.recPhotonsToMC.Path = "recPhotonsToMC"
delphessim.DataOutputs.recJetsToMC.Path = "recJetsToMC"
## FCC event-data model output -> define objects to be written out
from Configurables import PodioOutput
out = PodioOutput("out",OutputLevel=messageLevelOut)
out.filename = "FCCDelphesOutput.root"
out.outputCommands = ["drop *",
"keep genParticles",
"keep genVertices",
"keep recMuons",
"keep recElectrons",
"keep recCharged",
"keep recPhotons",
"keep recJets",
"keep recMETs",
"keep recMuonsToMC",
"keep recElectronsToMC",
"keep recChargedToMC",
"keep recPhotonsToMC",
"keep recJetsToMC"]
#out.outputCommands = ["keep *"]
############################################################
#
# Run modules
#
############################################################
# Run Pythia + Delphes
overlay.randomizePileup = False
overlay.mergeTools = ["PileupParticlesMergeTool/MyPileupParticlesMergeTool",
"PileupTrackHitMergeTool/MyTrackHitsMergeTool"]
overlay.PileUpTool = pileuptool
from Configurables import CreateExampleEventData
producer = CreateExampleEventData()
producer.magicNumberOffset = 300
producer.trackhits.Path = "trackHits"
producer.positionedtrackhits.Path = "positionedTrackHits"
producer.calohits.Path = "caloHits"
producer.positionedcalohits.Path = "positionedCaloHits"
producer.genparticles.Path = "genParticles"
producer.genvertices.Path = "genVertices"
# PODIO algorithm
from Configurables import PodioOutput
out = PodioOutput("out",
OutputLevel=DEBUG)
out.outputCommands = ["drop *", "keep overlaid*"]
out.filename = "dummyEventDataOverlaid.root"
# ApplicationMgr
from Configurables import ApplicationMgr
ApplicationMgr( TopAlg=[producer, overlay, out],
EvtSel='NONE',
EvtMax=300,
ExtSvc=[podioevent,],
OutputLevel=DEBUG,
)
mergelayersHcalFwd.inhits.Path = "HCalFwdHits"
mergelayersHcalFwd.outhits.Path = "mergedHCalFwdHits"
createHcalFwdCells = CreateCaloCells("CreateHcalFwdCaloCells",
doCellCalibration=True,
calibTool=calibHcalFwd,
addCellNoise=False,
filterCellNoise=False,
OutputLevel=INFO)
createHcalFwdCells.hits.Path = "mergedHCalFwdHits"
createHcalFwdCells.cells.Path = "HCalFwdCells"
out = PodioOutput("out", OutputLevel=INFO)
out.outputCommands = [
"drop *", "keep ECalBarrelCells", "keep ECalEndcapCells",
"keep ECalFwdCells", "keep newHCalCells", "keep newExtHCalCells",
"keep HCalEndcapCells", "keep HCalFwdCells"
]
out.filename = "output_fullCalo_SimAndDigi_e50GeV_" + str(
num_events) + "events.root"
#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
geantsim.AuditExecute = True
createEcalBarrelCellsStep1.AuditExecute = True
positionsEcalBarrel.AuditExecute = True
resegmentEcalBarrel.AuditExecute = True
createEcalBarrelCells.AuditExecute = True
hepmc_converter2 = HepMCToEDMConverter("Converter2")
hepmc_converter2.hepmc.Path="hepmc2"
hepmc_converter2.genparticles.Path="pileupGenParticles"
hepmc_converter2.genvertices.Path="pileupGenVertices"
from Configurables import GenMerge
mergealg = GenMerge()
mergealg.pileupGenParticles.Path="pileupGenParticles"
mergealg.pileupGenVertices.Path="pileupGenVertices"
mergealg.signalGenParticles.Path="signalGenParticles"
mergealg.signalGenVertices.Path="signalGenVertices"
mergealg.allGenParticles.Path="allGenParticles"
mergealg.allGenVertices.Path="allGenVertices"
from Configurables import PodioOutput
out = PodioOutput("out", filename = "output_gen.root")
out.outputCommands = ["keep *", "drop pileupGenVertices", "drop pileupGenParticles", "drop signalGenVertices", "drop signalGenParticles"]
ApplicationMgr(
TopAlg=[gun1, gun2, hepmc_converter1, hepmc_converter2, mergealg, out],
EvtSel='NONE',
ExtSvc = [podioevent],
EvtMax=1,
OutputLevel=VERBOSE,
)
from Configurables import CreateCaloCells
createcells = CreateCaloCells("CreateCaloCells",
geometryTool=barrelHcalGeometry,
calibTool=calibHcells,
doCellCalibration=True,
addCellNoise=True,
filterCellNoise=False,
noiseTool=noise,
OutputLevel=DEBUG)
createcells.hits.Path = "HCalBarrelCellsStep1"
createcells.cells.Path = "HCalCells"
out = PodioOutput("out",
filename="output_HCalCells_digitisation_allCells.root",
OutputLevel=DEBUG)
out.outputCommands = ["keep *", "drop HCalBarrelCellsStep1"]
#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
podioinput.AuditExecute = True
rewriteHCal.AuditExecute = True
createcells.AuditExecute = True
out.AuditExecute = True
ApplicationMgr(
TopAlg=[podioinput, rewriteHCal, createcells, out],
EvtSel='NONE',
EvtMax=1,
pandoralg.ECalToHadGeVCalibrationBarrel = 1.12 #very small effect
pandoralg.ECalToHadGeVCalibrationEndCap = 1.12
pandoralg.HCalToHadGeVCalibration = 1.07
pandoralg.MuonToMipCalibration = 10.0
pandoralg.DigitalMuonHits = 0
pandoralg.MaxHCalHitHadronicEnergy = 1.0
pandoralg.UseOldTrackStateCalculation = 0
pandoralg.AbsorberRadLengthECal = 0.2854
pandoralg.AbsorberIntLengthECal = 0.0101
pandoralg.AbsorberRadLengthHCal = 0.0569
pandoralg.AbsorberIntLengthHCal = 0.006
pandoralg.AbsorberRadLengthOther = 0.0569
pandoralg.AbsorberIntLengthOther = 0.006
##############################################################################
# write PODIO file
from Configurables import PodioOutput
write = PodioOutput("write")
write.filename = "pan_test.root"
write.outputCommands = ["keep *"]
# ApplicationMgr
from Configurables import ApplicationMgr
ApplicationMgr(TopAlg=[read, pandoralg],
EvtSel='NONE',
EvtMax=10,
ExtSvc=[dsvc, gearSvc],
HistogramPersistency="ROOT",
OutputLevel=INFO)
OutputLevel= DEBUG)
from Configurables import CreateCaloCells
createcells = CreateCaloCells("CreateCaloCells",
geometryTool = barrelHcalGeometry,
calibTool=calibHcells,
doCellCalibration = True,
addCellNoise = True, filterCellNoise = False,
noiseTool = noise,
OutputLevel = DEBUG)
createcells.hits.Path="HCalBarrelCellsStep1"
createcells.cells.Path="HCalCells"
out = PodioOutput("out", filename="output_HCalCells_digitisation_allCells.root",
OutputLevel = DEBUG)
out.outputCommands = ["keep *", "drop HCalBarrelCellsStep1"]
#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
podioinput.AuditExecute = True
rewriteHCal.AuditExecute = True
createcells.AuditExecute = True
out.AuditExecute = True
ApplicationMgr(
TopAlg = [podioinput,
rewriteHCal,
createcells,
positionsHCalBarrelTool=HCalBcells,
positionsHCalExtBarrelTool=HCalExtBcells,
positionsEMECTool=EMECcells,
positionsHECTool=HECcells,
positionsEMFwdTool=ECalFwdcells,
positionsHFwdTool=HCalFwdcells,
hits="caloClusterBarrelCells",
positionedHits="caloClusterBarrelCellPositions",
OutputLevel=INFO)
out = PodioOutput("out",
filename="output_BarrelTopo_50GeVe_3ev.root",
OutputLevel=DEBUG)
out.outputCommands = [
"drop *", "keep GenParticles", "keep GenVertices",
"keep caloClustersBarrel", "keep caloClusterBarrelCells",
"keep caloClusterBarrelCellPositions"
]
#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
podioinput.AuditExecute = True
createemptycells.AuditExecute = True
createTopoClusters.AuditExecute = True
positionsClusterBarrel.AuditExecute = True
out.AuditExecute = True
ApplicationMgr(
createTailCatcherCells = CreateCaloCells("CreateTailCatcherCells",
doCellCalibration=False,
addCellNoise=False,
filterCellNoise=False,
OutputLevel=INFO,
hits="TailCatcherHits",
cells="TailCatcherCells")
# PODIO algorithm
from Configurables import ApplicationMgr, FCCDataSvc, PodioOutput
podioevent = FCCDataSvc("EventDataSvc")
out = PodioOutput("out")
out.outputCommands = [
"drop *", "keep GenParticles", "keep GenVertices", "keep TrackerHits",
"keep TrackerPositionedHits", "keep TrackerDigiPostPoint",
"keep ECalBarrelCells", "keep ECalEndcapCells", "keep ECalFwdCells",
"keep HCalBarrelCells", "keep HCalExtBarrelCells", "keep HCalEndcapCells",
"keep HCalFwdCells", "keep TailCatcherCells"
]
out.filename = output_name
#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
geantsim.AuditExecute = True
createEcalBarrelCellsStep1.AuditExecute = True
positionsEcalBarrel.AuditExecute = True
resegmentEcalBarrel.AuditExecute = True
createEcalBarrelCells.AuditExecute = True
particleName="e-",
energyMin=energy,
energyMax=energy,
etaMin=-0.36,
etaMax=0.36,
OutputLevel=DEBUG)
geantsim = SimG4Alg("SimG4Alg",
outputs=["SimG4SaveCalHits/saveHCalHits"],
eventProvider=pgun,
OutputLevel=DEBUG)
# PODIO algorithm
from Configurables import PodioOutput
out = PodioOutput("out", OutputLevel=DEBUG)
out.outputCommands = ["keep *"]
out.filename = "output_hcalSim_e" + str(int(
energy / 1000)) + "GeV_eta036_1events.root"
#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
geantsim.AuditExecute = True
# ApplicationMgr
from Configurables import ApplicationMgr
ApplicationMgr(
TopAlg=[geantsim, out],
EvtSel='NONE',
filterCellNoise=False)
createHcalEndcapCells.hits.Path = "pileupHCalEndcapCells"
createHcalEndcapCells.cells.Path = "mergedHCalEndcapCells"
createHcalFwdCells = CreateCaloCells("CreateHcalFwdCells",
doCellCalibration=False,
addCellNoise=False,
filterCellNoise=False)
createHcalFwdCells.hits.Path = "pileupHCalFwdCells"
createHcalFwdCells.cells.Path = "mergedHCalFwdCells"
# PODIO algorithm
from Configurables import PodioOutput
out = PodioOutput("out")
out.outputCommands = [
"drop *", "keep mergedECalBarrelCells", "keep mergedECalEndcapCells",
"keep mergedECalFwdCells", "keep mergedHCalBarrelCells",
"keep mergedHCalExtBarrelCells", "keep mergedHCalEndcapCells",
"keep mergedHCalFwdCells"
]
out.filename = output_name
list_of_algorithms += [
overlay, createEcalBarrelCells, createEcalEndcapCells, createEcalFwdCells,
createHcalBarrelCells, createHcalExtBarrelCells, createHcalEndcapCells,
createHcalFwdCells, out
]
# ApplicationMgr
from Configurables import ApplicationMgr
ApplicationMgr(TopAlg=list_of_algorithms,
EvtSel='NONE',
EvtMax=num_events,