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 )