Example #1
0
    "pt > 5. & abs(eta) < 5.0")

# ------------------------------------------------------------------------------
# Jet Substructure (FastJet 3)
# ------------------------------------------------------------------------------
from RecoJets.JetProducers.ca4GenJets_cfi import ca4GenJets
from RecoJets.JetProducers.ca4PFJets_cfi import ca4PFJets

from RecoJets.JetProducers.ak5PFJetsPruned_cfi import ak5PFJetsPruned
ak5PrunedPFlow = ak5PFJetsPruned.clone(doAreaFastjet=cms.bool(True))

from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsFiltered
ak5FilteredPFlow = ak5PFJetsFiltered.clone(doAreaFastjet=cms.bool(True))

from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsMassDropFiltered
ak5MassDropFilteredPFlow = ak5PFJetsMassDropFiltered.clone(
    doAreaFastjet=cms.bool(True))

#process.ca12GenJetsNoNu = ca4GenJets.clone( rParam = cms.double(1.2),src = cms.InputTag("genParticlesForJetsNoNu"))
process.ca12GenJets = ca4GenJets.clone(rParam=cms.double(1.2),
                                       src=cms.InputTag("genParticlesForJets"))
process.ca12PFJetsPFlow = ca4PFJets.clone(rParam=cms.double(1.2),
                                          src=cms.InputTag("pfNoElectron" +
                                                           postfix),
                                          doAreaFastjet=cms.bool(True),
                                          doRhoFastjet=cms.bool(True),
                                          Rho_EtaMax=cms.double(6.0),
                                          Ghost_EtaMax=cms.double(7.0))
## this thing produces subjets by default
process.ca12PFJetsPrunedPFlow = ak5PrunedPFlow.clone(
    src=cms.InputTag("pfNoElectron" + postfix),
    doAreaFastjet=cms.bool(True),
Example #2
0
patJetsCA8CHSfilteredN3.getJetMCFlavour = False
patJetsCA8CHSfilteredN3.jetCorrFactorsSource = cms.VInputTag(cms.InputTag('patJetCorrFactorsCA8CHSfilteredN3'))

#### Adding Nsubjetiness
patJetsCA8CHSfilteredN3withNsub = cms.EDProducer("NjettinessAdder",
  src=cms.InputTag("patJetsCA8CHSfilteredN3"),
  cone=cms.double(0.8)
)


############################################################### MassDropFilter
from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsMassDropFiltered
ca8PFJetsCHSmassDropFiltered = ak5PFJetsMassDropFiltered.clone(
  src = 'pfNoPileUp',
  jetPtMin = cms.double(20.0),
  doAreaFastjet = cms.bool(True),
  rParam = cms.double(0.8),
  jetAlgorithm = cms.string("CambridgeAachen"),
)

jetSourceMassDropFiltered = 'ca8PFJetsCHSmassDropFiltered'

# corrections 
patJetCorrFactorsCA8CHSmassDropFiltered = patJetCorrFactors.clone()
patJetCorrFactorsCA8CHSmassDropFiltered.src = jetSourceMassDropFiltered
patJetCorrFactorsCA8CHSmassDropFiltered.levels == ['L1FastJet', 'L2Relative', 'L3Absolute']
patJetCorrFactorsCA8CHSmassDropFiltered.payload = 'AK7PFchs'
patJetCorrFactorsCA8CHSmassDropFiltered.useRho = True

patJetsCA8CHSmassDropFiltered = patJets.clone()
patJetsCA8CHSmassDropFiltered.jetSource = jetSourceMassDropFiltered
Example #3
0
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")


Example #4
0
process.selectedPatJetsCAVHFilterPF.cut = cms.string("pt > 5. & abs(eta) < 5.0")

# ------------------------------------------------------------------------------
# Jet Substructure (FastJet 3)
# ------------------------------------------------------------------------------
from RecoJets.JetProducers.ca4GenJets_cfi import ca4GenJets
from RecoJets.JetProducers.ca4PFJets_cfi import ca4PFJets

from RecoJets.JetProducers.ak5PFJetsPruned_cfi import ak5PFJetsPruned
ak5PrunedPFlow = ak5PFJetsPruned.clone(doAreaFastjet = cms.bool(True))

from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsFiltered
ak5FilteredPFlow = ak5PFJetsFiltered.clone(doAreaFastjet = cms.bool(True))

from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsMassDropFiltered
ak5MassDropFilteredPFlow = ak5PFJetsMassDropFiltered.clone(doAreaFastjet = cms.bool(True))

#process.ca12GenJetsNoNu = ca4GenJets.clone( rParam = cms.double(1.2),src = cms.InputTag("genParticlesForJetsNoNu"))
process.ca12GenJets = ca4GenJets.clone( rParam = cms.double(1.2),src = cms.InputTag("genParticlesForJets"))
process.ca12PFJetsPFlow = ca4PFJets.clone(
    rParam = cms.double(1.2),
    src = cms.InputTag("pfNoElectron"+postfix),
    doAreaFastjet = cms.bool(True),
    doRhoFastjet = cms.bool(True),
    Rho_EtaMax = cms.double(6.0),
    Ghost_EtaMax = cms.double(7.0)
    )
## this thing produces subjets by default
process.ca12PFJetsPrunedPFlow = ak5PrunedPFlow.clone(
    src = cms.InputTag("pfNoElectron"+postfix),
    doAreaFastjet = cms.bool(True),
Example #5
0
###############################

# Filtered PF Jets
process.caFilteredPFlow = ak5PFJetsFiltered.clone(
    src=cms.InputTag('pfNoElectron' + postfix),
    jetAlgorithm=cms.string("CambridgeAachen"),
    rParam=cms.double(1.2),
    writeCompound=cms.bool(True),
    doAreaFastjet=cms.bool(True),
    jetPtMin=cms.double(100.0))

from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsMassDropFiltered
process.caMassDropFilteredPFlow = ak5PFJetsMassDropFiltered.clone(
    src=cms.InputTag('pfNoElectron' + postfix),
    jetAlgorithm=cms.string("CambridgeAachen"),
    rParam=cms.double(1.2),
    writeCompound=cms.bool(True),
    doAreaFastjet=cms.bool(True),
    jetPtMin=cms.double(100.0))

process.caFilteredGenJetsNoNu = process.ca8GenJetsNoNu.clone(
    nFilt=cms.int32(2),
    rFilt=cms.double(0.3),
    useFiltering=cms.bool(True),
    useExplicitGhosts=cms.bool(True),
    writeCompound=cms.bool(True),
    rParam=cms.double(1.2),
    jetCollInstanceName=cms.string("SubJets"),
    jetPtMin=cms.double(100.0))

###############################
process.caFilteredPFlowCHS = ak5PFJetsFiltered.clone(
    src=cms.InputTag("pfNoElectron"),
    jetAlgorithm=cms.string("CambridgeAachen"),
    rParam=cms.double(1.2),
    writeCompound=cms.bool(True),
    doAreaFastjet=cms.bool(True),
    jetPtMin=cms.double(100.0),
)

from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsMassDropFiltered

process.caMassDropFilteredPFlow = ak5PFJetsMassDropFiltered.clone(
    src=cms.InputTag("particleFlow"),
    jetAlgorithm=cms.string("CambridgeAachen"),
    rParam=cms.double(1.2),
    writeCompound=cms.bool(True),
    doAreaFastjet=cms.bool(True),
    jetPtMin=cms.double(100.0),
)

process.caMassDropFilteredPFlowCHS = ak5PFJetsMassDropFiltered.clone(
    src=cms.InputTag("pfNoElectron"),
    jetAlgorithm=cms.string("CambridgeAachen"),
    rParam=cms.double(1.2),
    writeCompound=cms.bool(True),
    doAreaFastjet=cms.bool(True),
    jetPtMin=cms.double(100.0),
)


process.caFilteredGenJetsNoNu = process.ca8GenJetsNoNu.clone(
Example #7
0
patJetsAK4CHSfilteredN3.getJetMCFlavour = False
patJetsAK4CHSfilteredN3.jetCorrFactorsSource = cms.VInputTag(
    cms.InputTag('patJetCorrFactorsAK4CHSfilteredN3'))

#### Adding Nsubjetiness
patJetsAK4CHSfilteredN3withNsub = cms.EDProducer(
    "NjettinessAdder",
    src=cms.InputTag("patJetsAK4CHSfilteredN3"),
    cone=cms.double(0.4))

############################################################### MassDropFilter
from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsMassDropFiltered
ak4PFJetsCHSmassDropFiltered = ak5PFJetsMassDropFiltered.clone(
    src='pfNoPileUp',
    jetPtMin=cms.double(20.0),
    doAreaFastjet=cms.bool(True),
    rParam=cms.double(0.4),
    jetAlgorithm=cms.string("AntiKt"),
)

jetSourceMassDropFiltered = 'ak4PFJetsCHSmassDropFiltered'

# corrections
patJetCorrFactorsAK4CHSmassDropFiltered = patJetCorrFactors.clone()
patJetCorrFactorsAK4CHSmassDropFiltered.src = jetSourceMassDropFiltered
patJetCorrFactorsAK4CHSmassDropFiltered.levels == [
    'L1FastJet', 'L2Relative', 'L3Absolute'
]
patJetCorrFactorsAK4CHSmassDropFiltered.payload = 'AK5PFchs'
patJetCorrFactorsAK4CHSmassDropFiltered.useRho = True
patJetsAK7CHSfilteredN3.getJetMCFlavour = False
patJetsAK7CHSfilteredN3.jetCorrFactorsSource = cms.VInputTag(cms.InputTag('patJetCorrFactorsAK7CHSfilteredN3'))

#### Adding Nsubjetiness
patJetsAK7CHSfilteredN3withNsub = cms.EDProducer("NjettinessAdder",
  src=cms.InputTag("patJetsAK7CHSfilteredN3"),
  cone=cms.double(0.7)
)


############################################################### MassDropFilter
from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsMassDropFiltered
ak7PFJetsCHSmassDropFiltered = ak5PFJetsMassDropFiltered.clone(
  src = 'pfNoPileUp',
  jetPtMin = cms.double(20.0),
  doAreaFastjet = cms.bool(True),
  rParam = cms.double(0.7),
  jetAlgorithm = cms.string("AntiKt"),
)

jetSourceMassDropFiltered = 'ak7PFJetsCHSmassDropFiltered'

# corrections 
patJetCorrFactorsAK7CHSmassDropFiltered = patJetCorrFactors.clone()
patJetCorrFactorsAK7CHSmassDropFiltered.src = jetSourceMassDropFiltered
patJetCorrFactorsAK7CHSmassDropFiltered.levels == ['L1FastJet', 'L2Relative', 'L3Absolute']
patJetCorrFactorsAK7CHSmassDropFiltered.payload = 'AK7PFchs'
patJetCorrFactorsAK7CHSmassDropFiltered.useRho = True

patJetsAK7CHSmassDropFiltered = patJets.clone()
patJetsAK7CHSmassDropFiltered.jetSource = jetSourceMassDropFiltered
    yCut = cms.double(0.08),
    useFiltering = cms.bool(True),
    useDynamicFiltering = cms.bool(True),
    nFilt = cms.int32(3),
    rFilt = cms.double(0.3),
    rFiltFactor = cms.double(0.5),
    writeCompound = cms.bool(True),
    jetCollInstanceName=cms.string("SubJets")
)
from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsMassDropFiltered
process.ca8PFJetsCHSMDBDRSFiltered = ak5PFJetsMassDropFiltered.clone(
    jetAlgorithm = cms.string("CambridgeAachen"),
    rParam = cms.double(0.8),
    src = process.ca8PFJetsCHS.src,
    srcPVs = process.ca8PFJetsCHS.srcPVs,
    doAreaFastjet = process.ca8PFJetsCHS.doAreaFastjet,
    writeCompound = cms.bool(True),
    jetCollInstanceName=cms.string("SubJets"),
    jetPtMin = cms.double(20.),
    useDynamicFiltering = cms.bool(True),
    rFiltFactor = cms.double(0.5)
)
## CA8 Kt-BDRS filtered jets (Gen and Reco) (each module produces two jet collections, fat jets and subjets)
## Compared to the above filtered jets, here dynamic filtering radius is used (as in arXiv:0802.2470)
## However, here the mass drop is replaced by finding two Kt subjets which then set the size of the filtering radius
from RecoJets.JetProducers.ca4GenJets_cfi import ca4GenJets
process.ca8GenJetsNoNuKtBDRSFiltered = ca4GenJets.clone(
    rParam = cms.double(0.8),
    src = cms.InputTag("genParticlesForJetsNoNu"),
    usePruning = cms.bool(True),
    useKtPruning = cms.bool(True),
    zcut = cms.double(0.),