## Pythia generator
from Configurables import PythiaInterface
pythia8gen = PythiaInterface(Filename=pythiaConfFile,
OutputLevel=messageLevelPythia)
## Write the HepMC::GenEvent to the data service
pythia8gen.DataOutputs.hepmc.Path = "hepmc"
## Delphes simulator -> define objects to be written out
from Configurables import DelphesSimulation
delphessim = DelphesSimulation(DelphesCard=delphesCard,
ROOTOutputFile=delphesRootOutFile,
OutputLevel=messageLevelDelphes,
outputs=[
"DelphesSaveChargedParticles/muons",
"DelphesSaveChargedParticles/electrons",
"DelphesSaveChargedParticles/charged",
"DelphesSaveNeutralParticles/photons",
"DelphesSaveNeutralParticles/neutral",
"DelphesSaveGenJets/genJets",
"DelphesSaveJets/jets", "DelphesSaveMet/met"
])
delphessim.DataInputs.hepmc.Path = "hepmc"
delphessim.DataOutputs.genParticles.Path = "genParticles"
delphessim.DataOutputs.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"
#out.outputCommands = ["drop *",
# "keep genParticles",
from Gaudi.Configuration import *
from Configurables import ApplicationMgr, FCCDataSvc, AlbersWrite, AlbersOutput
albersevent = FCCDataSvc("EventDataSvc")
from Configurables import DelphesSimulation
delphessimulation = DelphesSimulation(
"Simulation",
filename="/afs/cern.ch/user/d/degrutto/public/forFCC/ZLEP_toAll.hepmc",
detectorcard="Simulation/examples/delphes_card_CMS.tcl",
outputcollections=[
"UniqueObjectFinder/electrons", "UniqueObjectFinder/muons",
"UniqueObjectFinder/jets", "UniqueObjectFinder/photons",
"MissingET/momentum", "ScalarHT/energy"
],
debug_delphes=True)
#delphessimulation = DelphesSimulation("Simulation", filename="example_MyPythia.dat" , detectorcard = "Simulation/examples/delphes_card_FCC_basic.tcl" , outputcollections = [ "ElectronEnergySmearing/electrons", "MuonMomentumSmearing/muons", "JetEnergyScale/jets", "MissingET/momentum", "ScalarHT/energy" ] )
# we shouls add b-jets and track inpact parameters, fixme
delphessimulation.Outputs.genparticles.Path = "genparticles"
delphessimulation.Outputs.particles.Path = "particles"
delphessimulation.Outputs.partons.Path = "partons"
delphessimulation.Outputs.muons.Path = "muons"
delphessimulation.Outputs.electrons.Path = "electrons"
delphessimulation.Outputs.photons.Path = "photons"
delphessimulation.Outputs.jets.Path = "jets"
delphessimulation.Outputs.mets.Path = "met"
delphessimulation.Outputs.hts.Path = "ht"
out = AlbersOutput("out", OutputLevel=DEBUG)
from FCCPileupScenarios import FCCPhase1PileupConf as pileupconf
pythia8gentool = PythiaInterface(Filename=pythiaConfFile, OutputLevel=messageLevelPythia)
## Write the HepMC::GenEvent to the data service
from Configurables import GenAlg
pythia8gen = GenAlg("Pythia8", SignalProvider=pythia8gentool)
pythia8gen.hepmc.Path = "hepmc"
## Delphes simulator -> define objects to be written out
from Configurables import DelphesSimulation
delphessim = DelphesSimulation(DelphesCard=delphesCard,
ROOTOutputFile=delphesRootOutFile,
ApplyGenFilter=True,
OutputLevel=messageLevelDelphes,
outputs=["DelphesSaveChargedParticles/genMuons",
"DelphesSaveChargedParticles/genElectrons",
"DelphesSaveNeutralParticles/genPhotons",
"DelphesSaveGenJets/nonPropagatedGenJets",
"DelphesSaveGenJets/nonPropagatedGenJetsPileUp",
"DelphesSaveGenJets/propagatedGenJets",
"DelphesSaveGenJets/propagatedGenJetsPileUp",
"DelphesSaveMet/genMET"])
delphessim.hepmc.Path = "hepmc"
delphessim.genParticles.Path = "skimmedGenParticles"
delphessim.mcEventWeights.Path = "mcEventWeights"
### Reads an HepMC::GenEvent from the data service and writes a collection of EDM Particles
from Configurables import HepMCToEDMConverter
hepmc_converter = HepMCToEDMConverter("Converter")
hepmc_converter.hepmc.Path="hepmc"
hepmc_converter.genparticles.Path="genParticles"
hepmc_converter.genvertices.Path="genVertices"
PileUpProvider=pythia8gentool,
PileUpTool=pileuptool,
HepMCMergeTool=mergetool)
pythia8gen.hepmc.Path = "hepmc"
## Delphes simulator -> define objects to be written out
from Configurables import DelphesSimulation
delphessim = DelphesSimulation(
DelphesCard=delphesCard,
ROOTOutputFile=delphesRootOutFile,
ApplyGenFilter=True,
OutputLevel=messageLevelDelphes,
outputs=[
"DelphesSaveChargedParticles/muons",
#"DelphesSaveChargedParticles/electrons",
#"DelphesSaveNeutralParticles/photons",
#"DelphesSaveChargedParticles/pfcharged",
#"DelphesSaveNeutralParticles/pfphotons",
#"DelphesSaveNeutralParticles/pfneutrals",
#"DelphesSaveGenJets/genJets",
#"DelphesSaveJets/jets",
#"DelphesSaveJets/fatjets",
#"DelphesSaveMet/met"
])
delphessim.hepmc.Path = "hepmc"
delphessim.genParticles.Path = "skimmedGenParticles"
delphessim.mcEventWeights.Path = "mcEventWeights"
### Reads an HepMC::GenEvent from the data service and writes a collection of EDM Particles
from Configurables import HepMCToEDMConverter
hepmc_converter = HepMCToEDMConverter("Converter")
fatjetSaveTool.subjetsPruningTagged.Path = "subjetsPruningTagged"
fatjetSaveTool.subjetsPruning.Path = "subjetsPruning"
fatjetSaveTool.subjetsSoftDropTagged.Path = "subjetsSoftDropTagged"
fatjetSaveTool.subjetsSoftDrop.Path = "subjetsSoftDrop"
fatjetSaveTool.subjetsTrimming.Path = "subjetsTrimming"
from Configurables import DelphesSaveMet
metSaveTool = DelphesSaveMet("met")
metSaveTool.delphesMETArrayName = "MissingET/momentum"
metSaveTool.delphesSHTArrayName = "ScalarHT/energy"
metSaveTool.missingEt.Path = "met"
## Delphes simulator -> define objects to be written out
from Configurables import DelphesSimulation
delphessim = DelphesSimulation()
## Define Delphes card
delphessim.DelphesCard = os.path.join(os.environ.get("FCCSWSHAREDIR", ""), "Sim/SimDelphesInterface/data/FCChh_DelphesCard_Baseline_v01.tcl")
delphessim.ROOTOutputFile = ""
delphessim.ApplyGenFilter = True
delphessim.outputs = [
"DelphesSaveChargedParticles/muons",
"DelphesSaveChargedParticles/electrons",
"DelphesSaveNeutralParticles/photons",
"DelphesSaveChargedParticles/pfcharged",
"DelphesSaveNeutralParticles/pfphotons",
"DelphesSaveNeutralParticles/pfneutrals",
"DelphesSaveGenJets/genJets",
"DelphesSaveJets/jets",
"DelphesSaveJets/fatjets",
"DelphesSaveMet/met",
# Configure individual modules (algorithms)
#
############################################################
## Pythia generator
from Configurables import PythiaInterface
pythia8gen = PythiaInterface(Filename=pythiaConfFile, OutputLevel=INFO)
## Write the HepMC::GenEvent to the data service
pythia8gen.Outputs.hepmc.Path = "hepmc"
## Delphes simulator
from Configurables import DelphesSimulation
delphessim = DelphesSimulation(DelphesCard=delphesCard,
OutputCollections=delphesOutCol,
HepMCInputFile=delphesHepMCInFile,
ROOTOutputFile=delphesRootOutFile,
OutputLevel=DEBUG)
delphessim.Inputs.hepmc.Path = "hepmc"
delphessim.Outputs.allParticles.Path = "allParticles"
delphessim.Outputs.genPartons.Path = "genPartons"
delphessim.Outputs.genStableParticles.Path = "genStableParticles"
## FCC event-data model output
from Configurables import AlbersWrite, AlbersOutput
out = AlbersOutput("out", OutputLevel=DEBUG)
#out.outputCommands = ["keep *"]
out.outputCommands = [
"drop *", "keep allParticles", "keep genPartons", "keep genStableParticles"
]
## Pythia generator
from Configurables import PythiaInterface
pythia8gen = PythiaInterface(Filename=pythiaConfFile,
OutputLevel=messageLevelPythia)
## Write the HepMC::GenEvent to the data service
pythia8gen.DataOutputs.hepmc.Path = "hepmc"
## Delphes simulator -> define objects to be written out
from Configurables import DelphesSimulation
delphessim = DelphesSimulation(DelphesCard=delphesCard,
HepMCInputFile=delphesHepMCInFile,
ROOTOutputFile=delphesRootOutFile,
MuonsOutArray=delphesMuonsOutArray,
ElectronsOutArray=delphesElectronsOutArray,
ChargedOutArray=delphesChargedOutArray,
NeutralOutArray=delphesNeutralOutArray,
PhotonsOutArray=delphesPhotonsOutArray,
JetsOutArray=delphesJetsOutArray,
METsOutArray=delphesMETsOutArray,
SHTsOutArray=delphesSHTsOutArray,
OutputLevel=messageLevelDelphes)
delphessim.DataInputs.hepmc.Path = "hepmc"
delphessim.DataOutputs.genParticles.Path = "genParticles"
delphessim.DataOutputs.genVertices.Path = "genVertices"
delphessim.DataOutputs.recMuons.Path = "recMuons"
delphessim.DataOutputs.recElectrons.Path = "recElectrons"
delphessim.DataOutputs.recCharged.Path = "recCharged"
delphessim.DataOutputs.recNeutral.Path = "recNeutral"
delphessim.DataOutputs.recPhotons.Path = "recPhotons"
delphessim.DataOutputs.recJets.Path = "recJets"
delphessim.DataOutputs.recMETs.Path = "recMETs"
jetSaveTool.jetsBTagged.Path = "bTags"
jetSaveTool.jetsCTagged.Path = "cTags"
jetSaveTool.jetsTauTagged.Path = "tauTags"
from Configurables import DelphesSaveMet
metSaveTool = DelphesSaveMet("met")
metSaveTool.delphesMETArrayName = "MissingET/momentum"
metSaveTool.delphesSHTArrayName = "ScalarHT/energy"
metSaveTool.missingEt.Path = "met"
## Delphes simulator -> define objects to be written out
from Configurables import DelphesSimulation
delphessim = DelphesSimulation()
## Define Delphes card
delphessim.DelphesCard = os.path.join(os.environ.get("DELPHES_DIR", ""), "cards/delphes_card_IDEA.tcl")
delphessim.ROOTOutputFile = ""
delphessim.ApplyGenFilter = True
delphessim.outputs = [
"DelphesSaveChargedParticles/muons",
"DelphesSaveChargedParticles/electrons",
"DelphesSaveNeutralParticles/photons",
"DelphesSaveChargedParticles/efcharged",
"DelphesSaveNeutralParticles/efphotons",
"DelphesSaveNeutralParticles/efneutrals",
"DelphesSaveGenJets/genJets",
"DelphesSaveJets/jets",
"DelphesSaveMet/met",
]
