from RecoJets.Configuration.RecoGenJets_cff import ak5GenJets
from PhysicsTools.JetMCAlgos.GenJetParticles_cfi import genParticlesForJetsPlusNoHadron
from PhysicsTools.JetMCAlgos.GenJetParticles_cff import genParticlesForJetsNoNuPlusNoHadron
from PhysicsTools.JetMCAlgos.GenHFHadronMatcher_cfi import matchGenHFHadron
# Supplies PDG ID to real name resolution of MC particles
from SimGeneral.HepPDTESSource.pythiapdt_cfi import *
# Configuration for matching B-hadrons ================================================================
genParticlesForJetsPlusBHadron = genParticlesForJetsPlusNoHadron.clone()
genParticlesForJetsPlusBHadron.injectHadronFlavours = [5]
ak5GenJetsPlusBHadron = ak5GenJets.clone()
ak5GenJetsPlusBHadron.src = "genParticlesForJetsPlusBHadron"
matchGenBHadron = matchGenHFHadron.clone()
matchGenBHadron.flavour = 5
matchGenBHadron.genJets = cms.InputTag("ak5GenJetsPlusBHadron", "", "")
genBHadronMatchingSequence = cms.Sequence( genParticlesForJetsPlusBHadron * ak5GenJetsPlusBHadron * matchGenBHadron )
# Configuration for matching C-hadrons =================================================================
genParticlesForJetsPlusCHadron = genParticlesForJetsPlusNoHadron.clone()
genParticlesForJetsPlusCHadron.injectHadronFlavours = [4]
ak5GenJetsPlusCHadron = ak5GenJets.clone()
ak5GenJetsPlusCHadron.src = "genParticlesForJetsPlusCHadron"
matchGenCHadron = matchGenHFHadron.clone()
matchGenCHadron.flavour = 4
matchGenCHadron.genJets = cms.InputTag("ak5GenJetsPlusCHadron", "", "")
partonicFinalState = cms.bool(False),
excludeResonances = cms.bool(True),
excludeFromResonancePids = cms.vuint32(12, 13, 14, 16),
tausAsJets = cms.bool(False)
)
genParticlesForJetsNoNu = genParticlesForJets.clone()
genParticlesForJetsNoNu.ignoreParticleIDs += cms.vuint32( 12,14,16)
# Flavour byReference
partons = cms.EDProducer("PartonSelector",
withLeptons = cms.bool(False),
src = cms.InputTag("genParticles")
)
AK7GenJets = ak5GenJets.clone(
rParam = cms.double(0.7)
)
# Flavour byValue PhysDef
AK7byValPhys = cms.EDProducer("JetFlavourIdentifier",
srcByReference = cms.InputTag("AK7byRef"),
physicsDefinition = cms.bool(True),
leptonInfo = cms.bool(True)
)
AK7byRef = cms.EDProducer("JetPartonMatcher",
jets = cms.InputTag("AK7PFJets"),
coneSizeToAssociate = cms.double(0.3),
partons = cms.InputTag("partons")
)
partonicFinalState = cms.bool(False),
excludeResonances = cms.bool(True),
excludeFromResonancePids = cms.vuint32(12, 13, 14, 16),
tausAsJets = cms.bool(False)
)
genParticlesForJetsNoNu = genParticlesForJets.clone()
genParticlesForJetsNoNu.ignoreParticleIDs += cms.vuint32( 12,14,16)
# Flavour byReference
partons = cms.EDProducer("PartonSelector",
withLeptons = cms.bool(False),
src = cms.InputTag("genParticles")
)
AK12GenJets = ak5GenJets.clone(
rParam = cms.double(1.2)
)
# Flavour byValue PhysDef
AK12byValPhys = cms.EDProducer("JetFlavourIdentifier",
srcByReference = cms.InputTag("AK12byRef"),
physicsDefinition = cms.bool(True),
leptonInfo = cms.bool(True)
)
AK12byRef = cms.EDProducer("JetPartonMatcher",
jets = cms.InputTag("AK12PFJets"),
coneSizeToAssociate = cms.double(0.3),
partons = cms.InputTag("partons")
)
partonicFinalState = cms.bool(False),
excludeResonances = cms.bool(True),
excludeFromResonancePids = cms.vuint32(12, 13, 14, 16),
tausAsJets = cms.bool(False)
)
genParticlesForJetsNoNu = genParticlesForJets.clone()
genParticlesForJetsNoNu.ignoreParticleIDs += cms.vuint32( 12,14,16)
# Flavour byReference
partons = cms.EDProducer("PartonSelector",
withLeptons = cms.bool(False),
src = cms.InputTag("genParticles")
)
AK4GenJets = ak5GenJets.clone(
rParam = cms.double(0.4)
)
# Flavour byValue PhysDef
AK4byValPhys = cms.EDProducer("JetFlavourIdentifier",
srcByReference = cms.InputTag("AK4byRef"),
physicsDefinition = cms.bool(True),
leptonInfo = cms.bool(True)
)
AK4byRef = cms.EDProducer("JetPartonMatcher",
jets = cms.InputTag("AK4PFJets"),
coneSizeToAssociate = cms.double(0.3),
partons = cms.InputTag("partons")
)
partonicFinalState = cms.bool(False),
excludeResonances = cms.bool(True),
excludeFromResonancePids = cms.vuint32(12, 13, 14, 16),
tausAsJets = cms.bool(False)
)
genParticlesForJetsNoNu = genParticlesForJets.clone()
genParticlesForJetsNoNu.ignoreParticleIDs += cms.vuint32( 12,14,16)
# Flavour byReference
partons = cms.EDProducer("PartonSelector",
withLeptons = cms.bool(False),
src = cms.InputTag("genParticles")
)
AK9GenJets = ak5GenJets.clone(
rParam = cms.double(0.9)
)
# Flavour byValue PhysDef
AK9byValPhys = cms.EDProducer("JetFlavourIdentifier",
srcByReference = cms.InputTag("AK9byRef"),
physicsDefinition = cms.bool(True),
leptonInfo = cms.bool(True)
)
AK9byRef = cms.EDProducer("JetPartonMatcher",
jets = cms.InputTag("AK9PFJets"),
coneSizeToAssociate = cms.double(0.3),
partons = cms.InputTag("partons")
)
partonicFinalState = cms.bool(False),
excludeResonances = cms.bool(True),
excludeFromResonancePids = cms.vuint32(12, 13, 14, 16),
tausAsJets = cms.bool(False)
)
genParticlesForJetsNoNu = genParticlesForJets.clone()
genParticlesForJetsNoNu.ignoreParticleIDs += cms.vuint32( 12,14,16)
# Flavour byReference
partons = cms.EDProducer("PartonSelector",
withLeptons = cms.bool(False),
src = cms.InputTag("genParticles")
)
AK5GenJets = ak5GenJets.clone(
rParam = cms.double(0.5)
)
# Flavour byValue PhysDef
AK5byValPhys = cms.EDProducer("JetFlavourIdentifier",
srcByReference = cms.InputTag("AK5byRef"),
physicsDefinition = cms.bool(True),
leptonInfo = cms.bool(True)
)
AK5byRef = cms.EDProducer("JetPartonMatcher",
jets = cms.InputTag("AK5PFJets"),
coneSizeToAssociate = cms.double(0.3),
partons = cms.InputTag("partons")
)
partonicFinalState = cms.bool(False),
excludeResonances = cms.bool(True),
excludeFromResonancePids = cms.vuint32(12, 13, 14, 16),
tausAsJets = cms.bool(False)
)
genParticlesForJetsNoNu = genParticlesForJets.clone()
genParticlesForJetsNoNu.ignoreParticleIDs += cms.vuint32( 12,14,16)
# Flavour byReference
partons = cms.EDProducer("PartonSelector",
withLeptons = cms.bool(False),
src = cms.InputTag("genParticles")
)
AK8GenJets = ak5GenJets.clone(
rParam = cms.double(0.8)
)
# Flavour byValue PhysDef
AK8byValPhys = cms.EDProducer("JetFlavourIdentifier",
srcByReference = cms.InputTag("AK8byRef"),
physicsDefinition = cms.bool(True),
leptonInfo = cms.bool(True)
)
AK8byRef = cms.EDProducer("JetPartonMatcher",
jets = cms.InputTag("AK8PFJets"),
coneSizeToAssociate = cms.double(0.3),
partons = cms.InputTag("partons")
)
partonicFinalState = cms.bool(False),
excludeResonances = cms.bool(True),
excludeFromResonancePids = cms.vuint32(12, 13, 14, 16),
tausAsJets = cms.bool(False)
)
genParticlesForJetsNoNu = genParticlesForJets.clone()
genParticlesForJetsNoNu.ignoreParticleIDs += cms.vuint32( 12,14,16)
# Flavour byReference
partons = cms.EDProducer("PartonSelector",
withLeptons = cms.bool(False),
src = cms.InputTag("genParticles")
)
AK6GenJets = ak5GenJets.clone(
rParam = cms.double(0.6)
)
# Flavour byValue PhysDef
AK6byValPhys = cms.EDProducer("JetFlavourIdentifier",
srcByReference = cms.InputTag("AK6byRef"),
physicsDefinition = cms.bool(True),
leptonInfo = cms.bool(True)
)
AK6byRef = cms.EDProducer("JetPartonMatcher",
jets = cms.InputTag("AK6PFJets"),
coneSizeToAssociate = cms.double(0.3),
partons = cms.InputTag("partons")
)
from RecoJets.Configuration.RecoGenJets_cff import ak5GenJets
from TopAnalysis.TopUtils.GenJetParticles_cfi import genParticlesForJetsPlusNoHadron
from TopAnalysis.TopUtils.GenJetParticles_cff import genParticlesForJetsNoNuPlusNoHadron
from TopAnalysis.TopUtils.GenHFHadronMatcher_cfi import matchGenHFHadron
# Supplies PDG ID to real name resolution of MC particles
from SimGeneral.HepPDTESSource.pythiapdt_cfi import *
# Configuration for matching B-hadrons ================================================================
genParticlesForJetsPlusBHadron = genParticlesForJetsPlusNoHadron.clone()
genParticlesForJetsPlusBHadron.injectHadronFlavours = [5]
ak5GenJetsPlusBHadron = ak5GenJets.clone()
ak5GenJetsPlusBHadron.src = "genParticlesForJetsPlusBHadron"
matchGenBHadron = matchGenHFHadron.clone()
matchGenBHadron.flavour = 5
matchGenBHadron.genJets = cms.InputTag("ak5GenJetsPlusBHadron", "", "")
genBHadronMatchingSequence = cms.Sequence( genParticlesForJetsPlusBHadron * ak5GenJetsPlusBHadron * matchGenBHadron )
# Configuration for matching C-hadrons =================================================================
genParticlesForJetsPlusCHadron = genParticlesForJetsPlusNoHadron.clone()
genParticlesForJetsPlusCHadron.injectHadronFlavours = [4]
ak5GenJetsPlusCHadron = ak5GenJets.clone()
ak5GenJetsPlusCHadron.src = "genParticlesForJetsPlusCHadron"
matchGenCHadron = matchGenHFHadron.clone()
matchGenCHadron.flavour = 4
matchGenCHadron.genJets = cms.InputTag("ak5GenJetsPlusCHadron", "", "")
import FWCore.ParameterSet.Config as cms from RecoJets.Configuration.RecoGenJets_cff import ak5GenJets from TopAnalysis.TopUtils.GenJetParticles_cff import genParticlesForJetsPlusHadron from TopAnalysis.TopUtils.GenJetParticles_cff import genParticlesForJetsNoNuPlusHadron ak5GenJetsPlusHadron = ak5GenJets.clone() ak5GenJetsPlusHadron.src = "genParticlesForJetsPlusHadron" ak5GenJetsNoNuPlusHadron = ak5GenJets.clone() ak5GenJetsNoNuPlusHadron.src = "genParticlesForJetsNoNuPlusHadron" improvedJetHadronQuarkMatchingSequence = cms.Sequence(genParticlesForJetsPlusHadron * ak5GenJetsPlusHadron) improvedJetNoNuHadronQuarkMatchingSequence = cms.Sequence(genParticlesForJetsNoNuPlusHadron * ak5GenJetsNoNuPlusHadron)