from Configurables import HepMCConverter
# reads an HepMC::GenEvent from the data service and writes
# a collection of EDM Particles
hepmc_converter = HepMCConverter("Converter")
# the input product name matches the output product name of the previous module
hepmc_converter.Inputs.hepmc.Path = "hepmc"
hepmc_converter.Outputs.genparticles.Path = "all_genparticles"
from Configurables import GenParticleFilter
genfilter = GenParticleFilter("StableParticles")
genfilter.Inputs.genparticles.Path = "all_genparticles"
genfilter.Outputs.genparticles.Path = "genparticles"
from Configurables import JetClustering_MCParticleCollection_GenJetCollection_GenJetParticleAssociationCollection_ as JetClustering
genjet_clustering = JetClustering("GenJetClustering", verbose=False)
genjet_clustering.Inputs.particles.Path = 'genparticles'
# giving a meaningful name for the output product
genjet_clustering.Outputs.jets.Path = 'genjets'
genjet_clustering.Outputs.constituents.Path = 'genjets_particles'
from Configurables import AlbersWrite, AlbersOutput
out = AlbersOutput("out", OutputLevel=DEBUG)
out.outputCommands = ["keep *"]
ApplicationMgr(
TopAlg=[pythia8gen, hepmc_converter, genfilter, genjet_clustering, out],
EvtSel='NONE',
EvtMax=100,
ExtSvc=[albersevent],
OutputLevel=DEBUG)
from Gaudi.Configuration import *
from Configurables import ApplicationMgr, HepMCReader, HepMCDumper, FCCDataSvc, AlbersWrite, AlbersOutput
albersevent = FCCDataSvc("EventDataSvc")
reader = HepMCReader("Reader", Filename="example_MyPythia.dat")
reader.Outputs.hepmc.Path = "hepmc"
dumper = HepMCDumper("Dumper")
dumper.Inputs.hepmc.Path="hepmc"
alberswrite = AlbersWrite("AlbersWrite",OutputLevel=DEBUG)
alberswrite.Outputs.albersJets.Path = "albersJets"
out = AlbersOutput("out")
ApplicationMgr( TopAlg = [reader,dumper,alberswrite,out],
EvtSel = 'NONE',
EvtMax = -1,
ExtSvc = [albersevent],
# EventLoop = eventloopmgr,
OutputLevel=INFO
)
