# - CHS, filtered
# - CHS, trimmed
ak5PFJetsCHS = ak5PFJets.clone(
src = cms.InputTag("pfNoPileUpJME")
)
ak5PFJetsCHSPruned = ak5PFJetsPruned.clone(
src = cms.InputTag("pfNoPileUpJME")
)
ak5PFJetsCHSFiltered = ak5PFJetsFiltered.clone(
src = cms.InputTag("pfNoPileUpJME")
)
ak5PFJetsCHSTrimmed = ak5PFJetsTrimmed.clone(
src = cms.InputTag("pfNoPileUpJME")
)
ak5PFJetsCHSSoftDrop = ak5PFJetsSoftDrop.clone(
src = cms.InputTag("pfNoPileUpJME")
)
ak4PFJetsCHS = ak5PFJetsCHS.clone(
rParam = 0.4
)
ak8PFJetsCHS = ak5PFJetsCHS.clone(
rParam = 0.8,
jetPtMin = 15.0
)
patJetsCA8CHSpruned.addGenJetMatch = False
patJetsCA8CHSpruned.getJetMCFlavour = False
patJetsCA8CHSpruned.jetCorrFactorsSource = cms.VInputTag(cms.InputTag('patJetCorrFactorsCA8CHSpruned'))
#### Adding Nsubjetiness
patJetsCA8CHSprunedwithNsub = cms.EDProducer("NjettinessAdder",
src=cms.InputTag("patJetsCA8CHSpruned"),
cone=cms.double(0.8)
)
################################################################ TRIMMING
from RecoJets.JetProducers.ak5PFJetsTrimmed_cfi import ak5PFJetsTrimmed
ca8PFJetsCHStrimmed = ak5PFJetsTrimmed.clone(
src = 'pfNoPileUp',
jetPtMin = cms.double(20.0),
doAreaFastjet = cms.bool(True),
rParam = cms.double(0.8),
jetAlgorithm = cms.string("CambridgeAachen"),
)
jetSourceTrimmed = 'ca8PFJetsCHStrimmed'
# corrections
patJetCorrFactorsCA8CHStrimmed = patJetCorrFactors.clone()
patJetCorrFactorsCA8CHStrimmed.src = jetSourceTrimmed
patJetCorrFactorsCA8CHStrimmed.levels = ['L1FastJet', 'L2Relative', 'L3Absolute']
patJetCorrFactorsCA8CHStrimmed.payload = 'AK7PFchs'
patJetCorrFactorsCA8CHStrimmed.useRho = True
### pats
patJetsCA8CHStrimmed = patJets.clone()
patJetsCA8CHStrimmed.jetSource = jetSourceTrimmed
###############################
###### AK 0.7 jets ############
###############################
process.ak7PFlow = process.pfJetsPFlow.clone(rParam=cms.double(0.7))
###############################
###### AK 0.8 jets ############
###############################
process.ak8PFlow = process.pfJetsPFlow.clone(rParam=cms.double(0.8))
###############################
###### AK 0.5 jets groomed ####
###############################
from RecoJets.JetProducers.ak5PFJetsTrimmed_cfi import ak5PFJetsTrimmed
process.ak5TrimmedPFlow = ak5PFJetsTrimmed.clone(src=process.pfJetsPFlow.src,
doAreaFastjet=cms.bool(True))
from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsFiltered
process.ak5FilteredPFlow = ak5PFJetsFiltered.clone(
src=process.pfJetsPFlow.src, doAreaFastjet=cms.bool(True))
from RecoJets.JetProducers.ak5PFJetsPruned_cfi import ak5PFJetsPruned
process.ak5PrunedPFlow = ak5PFJetsPruned.clone(src=process.pfJetsPFlow.src,
doAreaFastjet=cms.bool(True))
###############################
###### AK 0.7 jets groomed ####
###############################
process.ak7TrimmedPFlow = process.ak5TrimmedPFlow.clone(
src=process.pfJetsPFlow.src, rParam=cms.double(0.7))
process.pvCount = cms.EDFilter(
"EDPileupAna",
src=cms.InputTag('goodOfflinePrimaryVertices')
)
###############################
###### AK 0.5 jets groomed ####
###############################
ak5JetSrc = cms.InputTag("selectedPatJetsPFlow","pfCandidates")
from RecoJets.JetProducers.ak5PFJetsTrimmed_cfi import ak5PFJetsTrimmed
process.ak5TrimmedPFlow = ak5PFJetsTrimmed.clone(
src = ak5JetSrc,
doAreaFastjet = cms.bool(True)
)
from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsFiltered
process.ak5FilteredPFlow = ak5PFJetsFiltered.clone(
src = ak5JetSrc,
doAreaFastjet = cms.bool(True)
)
from RecoJets.JetProducers.ak5PFJetsPruned_cfi import ak5PFJetsPruned
process.ak5PrunedPFlow = ak5PFJetsPruned.clone(
src = ak5JetSrc,
doAreaFastjet = cms.bool(True)
)
patJetsAK4CHSpruned.getJetMCFlavour = False
patJetsAK4CHSpruned.jetCorrFactorsSource = cms.VInputTag(
cms.InputTag('patJetCorrFactorsAK4CHSpruned'))
#### Adding Nsubjetiness
patJetsAK4CHSprunedwithNsub = cms.EDProducer(
"NjettinessAdder",
src=cms.InputTag("patJetsAK4CHSpruned"),
cone=cms.double(0.4))
################################################################ TRIMMING
from RecoJets.JetProducers.ak5PFJetsTrimmed_cfi import ak5PFJetsTrimmed
ak4PFJetsCHStrimmed = ak5PFJetsTrimmed.clone(
src='pfNoPileUp',
jetPtMin=cms.double(20.0),
doAreaFastjet=cms.bool(True),
rParam=cms.double(0.4),
jetAlgorithm=cms.string("AntiKt"),
)
jetSourceTrimmed = 'ak4PFJetsCHStrimmed'
# corrections
patJetCorrFactorsAK4CHStrimmed = patJetCorrFactors.clone()
patJetCorrFactorsAK4CHStrimmed.src = jetSourceTrimmed
patJetCorrFactorsAK4CHStrimmed.levels = [
'L1FastJet', 'L2Relative', 'L3Absolute'
]
patJetCorrFactorsAK4CHStrimmed.payload = 'AK5PFchs'
patJetCorrFactorsAK4CHStrimmed.useRho = True
### pats
jetPtMin = cms.double(100.0)
)
ak8PFJetsCSConstituents = cms.EDProducer("PFJetConstituentSelector",
src = cms.InputTag("ak8PFJetsCS"),
cut = cms.string("pt > 100.0")
)
ak8PFJetsCSPruned = ak5PFJetsPruned.clone(
rParam = 0.8,
jetPtMin = 15.0,
src = cms.InputTag("ak8PFJetsCSConstituents", "constituents")
)
ak8PFJetsCSTrimmed = ak5PFJetsTrimmed.clone(
rParam = 0.8,
jetPtMin = 15.0,
src = cms.InputTag("ak8PFJetsCSConstituents", "constituents")
)
ak8PFJetsCSFiltered = ak5PFJetsFiltered.clone(
rParam = 0.8,
jetPtMin = 15.0,
src = cms.InputTag("ak8PFJetsCSConstituents", "constituents")
)
patJetsAK7CHSpruned.addGenJetMatch = False
patJetsAK7CHSpruned.getJetMCFlavour = False
patJetsAK7CHSpruned.jetCorrFactorsSource = cms.VInputTag(cms.InputTag('patJetCorrFactorsAK7CHSpruned'))
#### Adding Nsubjetiness
patJetsAK7CHSprunedwithNsub = cms.EDProducer("NjettinessAdder",
src=cms.InputTag("patJetsAK7CHSpruned"),
cone=cms.double(0.7)
)
################################################################ TRIMMING
from RecoJets.JetProducers.ak5PFJetsTrimmed_cfi import ak5PFJetsTrimmed
ak7PFJetsCHStrimmed = ak5PFJetsTrimmed.clone(
src = 'pfNoPileUp',
jetPtMin = cms.double(20.0),
doAreaFastjet = cms.bool(True),
rParam = cms.double(0.7),
jetAlgorithm = cms.string("AntiKt"),
)
jetSourceTrimmed = 'ak7PFJetsCHStrimmed'
# corrections
patJetCorrFactorsAK7CHStrimmed = patJetCorrFactors.clone()
patJetCorrFactorsAK7CHStrimmed.src = jetSourceTrimmed
patJetCorrFactorsAK7CHStrimmed.levels = ['L1FastJet', 'L2Relative', 'L3Absolute']
patJetCorrFactorsAK7CHStrimmed.payload = 'AK7PFchs'
patJetCorrFactorsAK7CHStrimmed.useRho = True
### pats
patJetsAK7CHStrimmed = patJets.clone()
patJetsAK7CHStrimmed.jetSource = jetSourceTrimmed
selectedPatJetsCA8CHSprunedPre +
btaggingCA8CHSprunedSubjets +
patJetsCA8CHSprunedSubjets +
selectedPatJetsCA8CHSpruned
)
# JETS TRIMMED CA8 ----------------------------
from RecoJets.JetProducers.ak5PFJetsTrimmed_cfi import ak5PFJetsTrimmed
ca8PFJetsCHStrimmed = ak5PFJetsTrimmed.clone(
src = 'pfNoPileUp',
jetPtMin = cms.double(30.0),
doAreaFastjet = cms.bool(True),
rParam = cms.double(0.8),
jetAlgorithm = cms.string("CambridgeAachen"),
rFilt = cms.double(0.1),
trimPtFracMin = cms.double(0.03),
writeCompound = cms.bool(True),
jetCollInstanceName = cms.string("SubJets")
)
jetSource = 'ca8PFJetsCHStrimmed'
# corrections
from PhysicsTools.PatAlgos.recoLayer0.jetCorrFactors_cfi import *
patJetCorrFactorsCA8CHStrimmed = patJetCorrFactors.clone()
patJetCorrFactorsCA8CHStrimmed.src = jetSource
# will need to add L2L3 corrections in the cfg
patJetCorrFactorsCA8CHStrimmed.levels = ['L1FastJet', 'L2Relative', 'L3Absolute']
patJetCorrFactorsCA8CHStrimmed.payload = 'AK7PFchs'
getattr(process,"pfNoTau"+postfix).enable = False
getattr(process,"pfNoJet"+postfix).enable = True
## Remove MC matching when running over data
if options.runOnData:
removeMCMatching( process, ['All'] )
## Define AK6 jets and subjets (GEN and RECO)
process.ak6GenJetsNoNu = process.ak5GenJetsNoNu.clone( rParam = 0.6 )
process.ak6PFJets = process.pfJetsPFlow.clone( rParam = 0.6 )
## AK6 groomed jets
from RecoJets.JetProducers.ak5PFJetsTrimmed_cfi import ak5PFJetsTrimmed
process.ak6PFJetsTrimmed = ak5PFJetsTrimmed.clone(
src = process.pfJetsPFlow.src,
doAreaFastjet = cms.bool(True),
rParam = cms.double(0.6)
)
from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsFiltered
process.ak6PFJetsFiltered = ak5PFJetsFiltered.clone(
src = process.pfJetsPFlow.src,
doAreaFastjet = cms.bool(True),
rParam = cms.double(0.6),
writeCompound = cms.bool(False),
jetCollInstanceName=cms.string("")
)
from RecoJets.JetProducers.ak5PFJetsPruned_cfi import ak5PFJetsPruned
process.ak6PFJetsPruned = ak5PFJetsPruned.clone(
src = process.pfJetsPFlow.src,
rParam=cms.double(0.8),
useConstituentSubtraction=cms.bool(True),
csRParam=cms.double(0.4),
csRho_EtaMax=ak4PFJets.
Rho_EtaMax, # Just use the same eta for both C.S. and rho by default
useExplicitGhosts=cms.bool(True),
doAreaFastjet=cms.bool(True),
jetPtMin=cms.double(100.0))
ak8PFJetsCSConstituents = cms.EDProducer("PFJetConstituentSelector",
src=cms.InputTag("ak8PFJetsCS"),
cut=cms.string("pt > 100.0"))
ak8PFJetsCSPruned = ak5PFJetsPruned.clone(rParam=0.8,
jetPtMin=15.0,
src=cms.InputTag(
"ak8PFJetsCSConstituents",
"constituents"))
ak8PFJetsCSTrimmed = ak5PFJetsTrimmed.clone(rParam=0.8,
jetPtMin=15.0,
src=cms.InputTag(
"ak8PFJetsCSConstituents",
"constituents"))
ak8PFJetsCSFiltered = ak5PFJetsFiltered.clone(rParam=0.8,
jetPtMin=15.0,
src=cms.InputTag(
"ak8PFJetsCSConstituents",
"constituents"))
src = process.ca8PFJetsCHS.src,
srcPVs = process.ca8PFJetsCHS.srcPVs,
doAreaFastjet = process.ca8PFJetsCHS.doAreaFastjet,
writeCompound = cms.bool(True),
jetCollInstanceName=cms.string("SubJets"),
jetPtMin = cms.double(20.),
useKtPruning = cms.bool(True),
zcut = cms.double(0.),
rcut_factor = cms.double(9999.)
)
## CA8 trimmed jets (Reco only)
from RecoJets.JetProducers.ak5PFJetsTrimmed_cfi import ak5PFJetsTrimmed
process.ca8PFJetsCHSTrimmed = ak5PFJetsTrimmed.clone(
jetAlgorithm = cms.string("CambridgeAachen"),
rParam = cms.double(0.8),
src = process.ca8PFJetsCHS.src,
srcPVs = process.ca8PFJetsCHS.srcPVs,
doAreaFastjet = process.ca8PFJetsCHS.doAreaFastjet,
jetPtMin = cms.double(20.)
)
#-------------------------------------
## PATify the above jets
from PhysicsTools.PatAlgos.tools.jetTools import *
## CA8 jets
switchJetCollection(process,
cms.InputTag('ca8PFJetsCHS'),
jetIdLabel='ca8',
rParam = 0.8,
useLegacyFlavour=False,
doJTA=options.doJTA,
doBTagging=options.doBTagging,