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:
process.load('CommonTools/PileupAlgos/Puppi_cff')
## e.g. to run on miniAOD
process.puppi.candName = cms.InputTag('packedPFCandidates')
process.puppi.vertexName = cms.InputTag('offlineSlimmedPrimaryVertices')
process.puppi.clonePackedCands = cms.bool(True)
process.puppi.useExistingWeights = cms.bool(True)
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"))
from RecoJets.JetProducers.ak4PFJetsSoftDrop_cfi import ak4PFJetsSoftDrop
process.PFJetsCHSSoftDrop = ak4PFJetsSoftDrop.clone(
jetAlgorithm=cms.string(options.jetAlgo),
rParam=cms.double(options.fatJetRadius),
R0=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))
if options.runFatJets:
## PATify the above jets
addJetCollection(
process,
labelName='PFCHS',
jetSource=cms.InputTag('PFJetsCHS'),
algo=algoLabel, # needed for jet flavor clustering
rParam=options.fatJetRadius, # needed for jet flavor clustering
pfCandidates=cms.InputTag(pfCandidates),
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")
)
from RecoJets.JetProducers.ak4PFJetsSoftDrop_cfi import ak4PFJetsSoftDrop
process.PFJetsCHSSoftDrop = ak4PFJetsSoftDrop.clone(
jetAlgorithm = cms.string(options.jetAlgo),
rParam = cms.double(options.fatJetRadius),
R0 = 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)
)
if options.runFatJets:
## PATify the above jets
addJetCollection(
process,
labelName='PFCHS',
jetSource=cms.InputTag('PFJetsCHS'),
algo=algoLabel, # needed for jet flavor clustering
rParam=options.fatJetRadius, # needed for jet flavor clustering
pfCandidates = cms.InputTag(pfCandidates),