RecoPFJets_cff.py

Go to the documentation of this file.

import FWCore.ParameterSet.Config as cms

from RecoJets.JetProducers.sc5PFJets_cfi import sisCone5PFJets
from RecoJets.JetProducers.ic5PFJets_cfi import iterativeCone5PFJets
from RecoJets.JetProducers.ak5PFJets_cfi import ak5PFJets
from RecoJets.JetProducers.ak4PFJets_cfi import ak4PFJets
from RecoJets.JetProducers.ak5PFJetsTrimmed_cfi import ak5PFJetsTrimmed
from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsFiltered, ak5PFJetsMassDropFiltered
from RecoJets.JetProducers.ak5PFJetsPruned_cfi import ak5PFJetsPruned
from RecoJets.JetProducers.ak5PFJetsSoftDrop_cfi import ak5PFJetsSoftDrop
from RecoJets.JetProducers.ak8PFJetsCS_cfi import ak8PFJetsCSConstituents, ak8PFJetsCS, ak8PFJetsCSPruned, ak8PFJetsCSTrimmed, ak8PFJetsCSFiltered
from CommonTools.ParticleFlow.pfNoPileUpJME_cff  import *
from RecoJets.JetProducers.gk5PFJets_cfi import gk5PFJets
from RecoJets.JetProducers.kt4PFJets_cfi import kt4PFJets
from RecoJets.JetProducers.ca4PFJets_cfi import ca4PFJets
from RecoJets.JetProducers.fixedGridRhoProducer_cfi import fixedGridRhoAll
from RecoJets.JetProducers.fixedGridRhoProducerFastjet_cfi import fixedGridRhoFastjetAll
from RecoJets.JetProducers.caTopTaggers_cff import *
from RecoJets.JetProducers.ak8PFJetsCHS_groomingValueMaps_cfi import ak8PFJetsCHSPrunedMass, ak8PFJetsCHSFilteredMass, ak8PFJetsCHSTrimmedMass, ak8PFJetsCHSSoftDropMass
from RecoJets.JetProducers.ca8PFJetsCHS_groomingValueMaps_cfi import ca8PFJetsCHSPrunedMass, ca8PFJetsCHSFilteredMass, ca8PFJetsCHSTrimmedMass, ca8PFJetsCHSSoftDropMass
from CommonTools.PileupAlgos.Puppi_cff import puppi
from CommonTools.PileupAlgos.softKiller_cfi import softKiller
from RecoJets.JetProducers.ak4PFJetsPuppi_cfi import ak4PFJetsPuppi
from RecoJets.JetProducers.ak8PFJetsPuppi_cfi import ak8PFJetsPuppi
from RecoJets.JetProducers.ak4PFJetsSK_cfi import ak4PFJetsSK

sisCone7PFJets = sisCone5PFJets.clone( rParam = 0.7 )
ak7PFJets = ak5PFJets.clone( rParam = 0.7 )
ak8PFJets = ak5PFJets.clone( rParam = 0.8 )
gk7PFJets = gk5PFJets.clone( rParam = 0.7 )
kt6PFJets = kt4PFJets.clone( rParam = 0.6 )
ca8PFJets = ca4PFJets.clone( rParam = 0.8 )

#compute areas for Fastjet PU subtraction  
kt6PFJets.doRhoFastjet = True
kt6PFJets.doAreaFastjet = True
#use active areas and not Voronoi tessellation for the moment
kt6PFJets.voronoiRfact = 0.9
ak4PFJets.doAreaFastjet = True
ak5PFJets.doAreaFastjet = True
ak5PFJetsTrimmed.doAreaFastjet = True
ak7PFJets.doAreaFastjet = True
ak8PFJets.doAreaFastjet = True
ak8PFJetsPuppi.doAreaFastjet = True
ak4PFJetsSK.doAreaFastjet = True

kt6PFJetsCentralChargedPileUp = kt6PFJets.clone(
    src = cms.InputTag("pfPileUpAllChargedParticles"),
    Ghost_EtaMax = cms.double(3.1),
    Rho_EtaMax = cms.double(2.5)
    )

kt6PFJetsCentralNeutral = kt6PFJets.clone(
    src = cms.InputTag("pfAllNeutralHadronsAndPhotons"),
    Ghost_EtaMax = cms.double(3.1),
    Rho_EtaMax = cms.double(2.5),
    inputEtMin = cms.double(0.5)
    )


kt6PFJetsCentralNeutralTight = kt6PFJetsCentralNeutral.clone(
    inputEtMin = cms.double(1.0)
    )


fixedGridRhoFastjetCentral = fixedGridRhoFastjetAll.clone(
    maxRapidity = cms.double(2.5)
    )

fixedGridRhoFastjetCentralChargedPileUp = fixedGridRhoFastjetAll.clone(
    pfCandidatesTag = "pfPileUpAllChargedParticles",
    maxRapidity = 2.5
    )

fixedGridRhoFastjetCentralNeutral = fixedGridRhoFastjetAll.clone(
    pfCandidatesTag = "pfAllNeutralHadronsAndPhotons",
    maxRapidity = 2.5
    )



ak8PFJetsCHSConstituents = cms.EDFilter("PFJetConstituentSelector",
                                        src = cms.InputTag("ak8PFJetsCHS"),
                                        cut = cms.string("pt > 100.0 && abs(rapidity()) < 2.4")
                                        )


# Advanced Algorithms for AK4, AK5, AK8 and CA8 :
#   - CHS, ungroomed
#   - CHS, pruned
#   - 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,
    jetPtMin = 5.0
    )    

ak8PFJetsCHS = ak5PFJetsCHS.clone(
    rParam = 0.8,
    jetPtMin = 50.0
    )

ak8PFJetsCHSPruned = ak5PFJetsCHSPruned.clone(
    rParam = 0.8,
    jetPtMin = 100.0
    )

ak8PFJetsCHSFiltered = ak5PFJetsCHSFiltered.clone(
    rParam = 0.8,
    jetPtMin = 15.0
    )

ak8PFJetsCHSTrimmed = ak5PFJetsCHSTrimmed.clone(
    rParam = 0.8,
    jetPtMin = 15.0
    )

ak8PFJetsCHSSoftDrop = ak5PFJetsCHSSoftDrop.clone(
    rParam = 0.8,
    jetPtMin = 100.0,
    R0 = 0.8
    )

ak8PFJetsPuppiSoftDrop = ak8PFJetsCHSSoftDrop.clone(
    src = cms.InputTag("puppi")
    )



ca8PFJetsCHS = ak8PFJetsCHS.clone(
    jetAlgorithm = cms.string("CambridgeAachen")
    )

ca8PFJetsCHSPruned = ak8PFJetsCHSPruned.clone(
    jetAlgorithm = cms.string("CambridgeAachen")
    )

ca8PFJetsCHSFiltered = ak8PFJetsCHSFiltered.clone(
    jetAlgorithm = cms.string("CambridgeAachen")
    )

ca8PFJetsCHSTrimmed = ak8PFJetsCHSTrimmed.clone(
    jetAlgorithm = cms.string("CambridgeAachen")
    )


ca8PFJetsCHSSoftDrop = ak8PFJetsCHSSoftDrop.clone(
    jetAlgorithm = cms.string("CambridgeAachen")
    )



# Higgs taggers
ca15PFJetsCHSMassDropFiltered = ak5PFJetsMassDropFiltered.clone(
    jetAlgorithm = cms.string("CambridgeAachen"),
    src = cms.InputTag("pfNoPileUpJME"),
    rParam = 1.5,
    jetPtMin=100.0
    )

ca15PFJetsCHSFiltered = ak5PFJetsFiltered.clone(
    jetAlgorithm = cms.string("CambridgeAachen"),
    src = cms.InputTag("pfNoPileUpJME"),
    rParam = 1.5,
    jetPtMin=100.0
    )

cmsTopTagPFJetsCHS.src = cms.InputTag("ak8PFJetsCHSConstituents", "constituents")
hepTopTagPFJetsCHS.src = cms.InputTag("ak8PFJetsCHSConstituents", "constituents")



recoPFJets   =cms.Sequence(fixedGridRhoAll+
                           fixedGridRhoFastjetAll+
                           fixedGridRhoFastjetCentral+
                           fixedGridRhoFastjetCentralChargedPileUp+
                           fixedGridRhoFastjetCentralNeutral+
                           ak4PFJets+
                           pfNoPileUpJMESequence+
                           ak4PFJetsCHS+                           
                           ak8PFJetsCHS+
                           ak8PFJetsCHSConstituents+
                           ak8PFJetsCHSSoftDrop+
                           ak8PFJetsCHSSoftDropMass+
                           cmsTopTagPFJetsCHS
    )

recoAllPFJets=cms.Sequence(sisCone5PFJets+sisCone7PFJets+
                           iterativeCone5PFJets+
                           kt4PFJets+kt6PFJets+
                           kt6PFJetsCentralChargedPileUp+
                           kt6PFJetsCentralNeutral+
                           kt6PFJetsCentralNeutralTight+
                           fixedGridRhoAll+
                           fixedGridRhoFastjetAll+
                           fixedGridRhoFastjetCentral+
                           fixedGridRhoFastjetCentralChargedPileUp+
                           fixedGridRhoFastjetCentralNeutral+
                           iterativeCone5PFJets+
                           ak4PFJets+
                           ak5PFJets+ak7PFJets+ak8PFJets+
                           gk5PFJets+gk7PFJets+
                           ca4PFJets+ca8PFJets+
                           pfNoPileUpJMESequence+
                           ak5PFJetsCHS+
                           ak5PFJetsCHSPruned+
                           ak5PFJetsCHSFiltered+
                           ak5PFJetsCHSTrimmed+
                           ak5PFJetsCHSSoftDrop+
                           ak4PFJetsCHS+                                                      
                           ak8PFJetsCHS+
                           ak8PFJetsCHSPruned+
                           ak8PFJetsCHSFiltered+
                           ak8PFJetsCHSTrimmed+
                           ak8PFJetsCHSSoftDrop+
                           ak8PFJetsCHSPrunedMass+
                           ak8PFJetsCHSTrimmedMass+
                           ak8PFJetsCHSSoftDropMass+
                           ak8PFJetsCHSFilteredMass+
                           ca8PFJetsCHS+
                           ca8PFJetsCHSPruned+
                           ca8PFJetsCHSFiltered+
                           ca8PFJetsCHSTrimmed+
                           ca8PFJetsCHSSoftDrop+
                           cmsTopTagPFJetsCHS+
                           hepTopTagPFJetsCHS+
                           ca15PFJetsCHSMassDropFiltered+
                           ca15PFJetsCHSFiltered+
                           ak8PFJetsCS+
                           ak8PFJetsCSConstituents+
                           ak8PFJetsCSPruned+
                           ak8PFJetsCSTrimmed+
                           ak8PFJetsCSFiltered
    )


recoPFJetsWithSubstructure=cms.Sequence(
                           fixedGridRhoAll+
                           fixedGridRhoFastjetAll+
                           fixedGridRhoFastjetCentral+
                           fixedGridRhoFastjetCentralChargedPileUp+
                           fixedGridRhoFastjetCentralNeutral+
                           ak4PFJets+
                           ak5PFJets+ak7PFJets+ak8PFJets+
                           ca4PFJets+ca8PFJets+
                           pfNoPileUpJMESequence+
                           ak5PFJetsCHS+
                           ak5PFJetsCHSPruned+
                           ak5PFJetsCHSFiltered+
                           ak5PFJetsCHSTrimmed+
                           ak5PFJetsCHSSoftDrop+
                           ak4PFJetsCHS+                                                      
                           ak8PFJetsCHS+
                           ak8PFJetsCHSPruned+
                           ak8PFJetsCHSFiltered+
                           ak8PFJetsCHSTrimmed+
                           ak8PFJetsCHSSoftDrop+
                           ak8PFJetsCHSConstituents+
                           ak8PFJetsCHSPrunedMass+
                           ak8PFJetsCHSTrimmedMass+
                           ak8PFJetsCHSSoftDropMass+
                           ak8PFJetsCHSFilteredMass+
                           ca8PFJetsCHS+
                           ca8PFJetsCHSPruned+
                           ca8PFJetsCHSFiltered+
                           ca8PFJetsCHSTrimmed+
                           ca8PFJetsCHSSoftDrop+
                           cmsTopTagPFJetsCHS+
                           hepTopTagPFJetsCHS+
                           ca15PFJetsCHSMassDropFiltered+
                           ca15PFJetsCHSFiltered+
                           ak8PFJetsCS+
                           ak8PFJetsCSConstituents+
                           ak8PFJetsCSPruned+
                           ak8PFJetsCSTrimmed+
                           ak8PFJetsCSFiltered+
                           puppi+
                           ak4PFJetsPuppi+
                           softKiller+
                           ak4PFJetsSK
    )