)
)
###############################################
# CHS JETS
#
# configure additional jet collections, based on chs.
process.chs = cms.EDFilter("CandPtrSelector", src = cms.InputTag("packedPFCandidates"), cut = cms.string("fromPV"))
process.ca8CHSJets = ca8PFJets.clone (src = 'chs', doAreaFastjet = True, jetPtMin = fatjet_ptmin)
process.ak8CHSJets = ak8PFJets.clone (src = 'chs', doAreaFastjet = True, jetPtMin = fatjet_ptmin)
process.ca15CHSJets = process.ca8CHSJets.clone (rParam = 1.5)
from RecoJets.JetProducers.ak4PFJetsPruned_cfi import ak4PFJetsPruned
process.ca8CHSJetsPruned = ak4PFJetsPruned.clone(rParam = 0.8, jetAlgorithm = "CambridgeAachen", doAreaFastjet = True, src = 'chs', jetPtMin = fatjet_ptmin)
from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsFiltered
process.ca15CHSJetsFiltered = ak5PFJetsFiltered.clone(
src = 'chs',
jetAlgorithm = cms.string("CambridgeAachen"),
rParam = cms.double(1.5),
writeCompound = cms.bool(True),
doAreaFastjet = cms.bool(True),
jetPtMin = cms.double(fatjet_ptmin)
)
from RecoJets.JetProducers.AnomalousCellParameters_cfi import *
from RecoJets.JetProducers.PFJetParameters_cfi import *
from RecoJets.JetProducers.CATopJetParameters_cfi import CATopJetParameters
process.cmsTopTagCHS = cms.EDProducer(
)
###############################################
# CHS JETS
#
# configure additional jet collections, based on chs.
process.chs = cms.EDFilter("CandPtrSelector", src = cms.InputTag("packedPFCandidates"), cut = cms.string("fromPV"))
process.ca8CHSJets = ca8PFJets.clone (src = 'chs', doAreaFastjet = True, jetPtMin = fatjet_ptmin)
process.ak8CHSJets = ak8PFJets.clone (src = 'chs', doAreaFastjet = True, jetPtMin = 10.)
process.ak8CHSJetsFat = ak8PFJets.clone (src = 'chs', doAreaFastjet = True, jetPtMin = fatjet_ptmin)
process.ca15CHSJets = process.ca8CHSJets.clone (rParam = 1.5)
from RecoJets.JetProducers.ak4PFJetsPruned_cfi import ak4PFJetsPruned
process.ca8CHSJetsPruned = ak4PFJetsPruned.clone(rParam = 0.8, jetAlgorithm = "CambridgeAachen", doAreaFastjet = True, src = 'chs', jetPtMin = fatjet_ptmin)
from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsFiltered
process.ca15CHSJetsFiltered = ak5PFJetsFiltered.clone(
src = 'chs',
jetAlgorithm = cms.string("CambridgeAachen"),
rParam = cms.double(1.5),
writeCompound = cms.bool(True),
doAreaFastjet = cms.bool(True),
jetPtMin = cms.double(fatjet_ptmin)
)
from RecoJets.JetProducers.AnomalousCellParameters_cfi import *
from RecoJets.JetProducers.PFJetParameters_cfi import *
from RecoJets.JetProducers.CATopJetParameters_cfi import CATopJetParameters
process.cmsTopTagCHS = cms.EDProducer(
process.genJetsNoNuPruned = ak4GenJets.clone(
SubJetParameters,
jetAlgorithm=cms.string(options.jetAlgo),
rParam=cms.double(options.fatJetRadius),
src=(cms.InputTag("packedGenParticlesForJetsNoNu") if options.miniAOD else
cms.InputTag("genParticlesForJetsNoNu" + postfix)),
usePruning=cms.bool(True),
writeCompound=cms.bool(True),
jetCollInstanceName=cms.string("SubJets"))
from RecoJets.JetProducers.ak4PFJetsPruned_cfi import ak4PFJetsPruned
process.PFJetsCHSPruned = ak4PFJetsPruned.clone(
jetAlgorithm=cms.string(options.jetAlgo),
rParam=cms.double(options.fatJetRadius),
src=(getattr(process, "ak4PFJets").src if options.miniAOD else getattr(
process, "pfJetsPFBRECO" + postfix).src),
srcPVs=(getattr(process, "ak4PFJets").srcPVs if options.miniAOD else
getattr(process, "pfJetsPFBRECO" + postfix).srcPVs),
doAreaFastjet=cms.bool(True),
writeCompound=cms.bool(True),
jetCollInstanceName=cms.string("SubJets"),
jetPtMin=cms.double(options.fatJetRawPtMin))
## SoftDrop fat jets (Gen and Reco) (each module produces two jet collections, fat jets and subjets)
process.genJetsNoNuSoftDrop = ak4GenJets.clone(
jetAlgorithm=cms.string(options.jetAlgo),
rParam=cms.double(options.fatJetRadius),
src=(cms.InputTag("packedGenParticlesForJetsNoNu") if options.miniAOD else
cms.InputTag("genParticlesForJetsNoNu" + postfix)),
useSoftDrop=cms.bool(True),
zcut=cms.double(0.1),
beta=cms.double(0.0),
R0=cms.double(options.fatJetRadius),
rParam=0.8)
## Standard AK4 and AK8 jets without CHS
from RecoJets.JetProducers.ak4PFJets_cfi import ak4PFJets
process.ak4PFJets = ak4PFJets.clone(
src='packedPFCandidates',
doAreaFastjet=True) #HERE!!! add jetPtMin = (fatjet_ptmin)
process.ak8PFJets = ak4PFJets.clone(src='packedPFCandidates',
rParam=0.8,
doAreaFastjet=True,
jetPtMin=fatjet_ptmin)
## Pruned AK8 without CHS
from RecoJets.JetProducers.ak4PFJetsPruned_cfi import ak4PFJetsPruned
process.ak8PFJetsPruned = ak4PFJetsPruned.clone(rParam=0.8,
doAreaFastjet=True,
src='packedPFCandidates',
jetPtMin=fatjet_ptmin)
## Softdrop AK8 without CHS
from RecoJets.JetProducers.ak4PFJetsSoftDrop_cfi import ak4PFJetsSoftDrop
process.ak8PFJetsSoftDrop = ak4PFJetsSoftDrop.clone(R0=0.8,
rParam=0.8,
doAreaFastjet=True,
src='packedPFCandidates',
jetPtMin=fatjet_ptmin)
## PUPPI AK8
## Pre-requisites: re-build puppi jets
###Robin's way:
from RecoJets.JetProducers.ca4GenJets_cfi import ca4GenJets
from RecoJets.JetProducers.SubJetParameters_cfi import SubJetParameters
process.ca8GenJetsNoNuPruned = ca4GenJets.clone(
SubJetParameters,
rParam = cms.double(0.8),
src = (cms.InputTag("packedGenParticlesForJetsNoNu") if options.miniAOD else cms.InputTag("genParticlesForJetsNoNu"+postfix)),
usePruning = cms.bool(True),
writeCompound = cms.bool(True),
jetCollInstanceName=cms.string("SubJets")
)
from RecoJets.JetProducers.ak4PFJetsPruned_cfi import ak4PFJetsPruned
process.ca8PFJetsPruned = ak4PFJetsPruned.clone(
jetAlgorithm = cms.string("CambridgeAachen"),
rParam = cms.double(0.8),
src = (getattr(process,"ak4PFJets").src if options.miniAOD else getattr(process,"pfJets"+postfix).src),
srcPVs = (getattr(process,"ak4PFJets").srcPVs if options.miniAOD else getattr(process,"pfJets"+postfix).srcPVs),
doAreaFastjet = cms.bool(True),
writeCompound = cms.bool(True),
jetCollInstanceName=cms.string("SubJets"),
jetPtMin = cms.double(options.fatJetPtMin)
)
if options.runSubJets:
## PATify CA8 jets
addJetCollection(
process,
labelName = 'CA8',
jetSource = cms.InputTag('ca8PFJets'),
trackSource = cms.InputTag(trackSource),
pvSource = cms.InputTag(pvSource),
svSource = svSource,
btagInfos = bTagInfos,
from RecoJets.JetProducers.SubJetParameters_cfi import SubJetParameters
process.genJetsNoNuPruned = ak4GenJets.clone(
SubJetParameters,
jetAlgorithm = cms.string(options.jetAlgo),
rParam = cms.double(options.fatJetRadius),
src = (cms.InputTag("packedGenParticlesForJetsNoNu") if options.miniAOD else cms.InputTag("genParticlesForJetsNoNu"+postfix)),
usePruning = cms.bool(True),
writeCompound = cms.bool(True),
jetCollInstanceName=cms.string("SubJets")
)
from RecoJets.JetProducers.ak4PFJetsPruned_cfi import ak4PFJetsPruned
process.PFJetsCHSPruned = ak4PFJetsPruned.clone(
jetAlgorithm = cms.string(options.jetAlgo),
rParam = cms.double(options.fatJetRadius),
src = (getattr(process,"ak4PFJets").src if options.miniAOD else getattr(process,"pfJetsPFBRECO"+postfix).src),
srcPVs = (getattr(process,"ak4PFJets").srcPVs if options.miniAOD else getattr(process,"pfJetsPFBRECO"+postfix).srcPVs),
doAreaFastjet = cms.bool(True),
writeCompound = cms.bool(True),
jetCollInstanceName=cms.string("SubJets"),
jetPtMin = cms.double(options.fatJetPtMin)
)
## SoftDrop fat jets (Gen and Reco) (each module produces two jet collections, fat jets and subjets)
process.genJetsNoNuSoftDrop = ak4GenJets.clone(
jetAlgorithm = cms.string(options.jetAlgo),
rParam = cms.double(options.fatJetRadius),
src = (cms.InputTag("packedGenParticlesForJetsNoNu") if options.miniAOD else cms.InputTag("genParticlesForJetsNoNu"+postfix)),
useSoftDrop = cms.bool(True),
zcut = cms.double(0.1),
beta = cms.double(0.0),
R0 = cms.double(options.fatJetRadius),
writeCompound = cms.bool(True),
jetCollInstanceName=cms.string("SubJets")