JetProducers/python/ak8PFJets_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002 from RecoJets.JetProducers.SubJetParameters_cfi import SubJetParameters
0003
0004
0005 from RecoJets.JetProducers.ak4PFJets_cfi import ak4PFJets
0006
0007 ##############################################################################
0008 # Standard AK8 Jets####################################################
0009 ##########################
0010 ak8PFJets = ak4PFJets.clone(
0011     rParam   = 0.8,
0012     jetPtMin = 50.0
0013     )
0014
0015 ##############################################################################
0016 # AK8 jets with various pileup subtraction schemes
0017 ##############################################################################
0018 ak8PFJetsPuppi = ak8PFJets.clone(
0019     src = "particleFlow",
0020     applyWeight = True,
0021     srcWeights  = cms.InputTag("puppi")
0022     )
0023
0024 ak8PFJetsCHS = ak8PFJets.clone(
0025     src = "pfNoPileUpJME"
0026     )
0027
0028 from Configuration.ProcessModifiers.pp_on_AA_cff import pp_on_AA
0029 pp_on_AA.toModify(ak8PFJetsCHS, src = "pfEmptyCollection")
0030 pp_on_AA.toModify(ak8PFJetsPuppi, src = "pfEmptyCollection")
0031
0032 ak8PFJetsCS = ak8PFJets.clone(
0033     useConstituentSubtraction = cms.bool(True),
0034     csRParam = cms.double(0.4),
0035     csRho_EtaMax = ak8PFJets.Rho_EtaMax,   # Just use the same eta for both C.S. and rho by default
0036     useExplicitGhosts = cms.bool(True),
0037     doAreaFastjet = True,
0038     jetPtMin = 100.0
0039     )
0040
0041
0042 ##############################################################################
0043 # Preclustered constituents for substructure, various subtraction schemes
0044 ##############################################################################
0045 ak8PFJetsCSConstituents = cms.EDProducer("PFJetConstituentSelector",
0046                                          src = cms.InputTag("ak8PFJetsCS"),
0047                                          cut = cms.string("pt > 100.0")
0048                                         )
0049
0050 ak8PFJetsCHSConstituents = cms.EDProducer("PFJetConstituentSelector",
0051                                           src = cms.InputTag("ak8PFJetsCHS"),
0052                                           cut = cms.string("pt > 100.0 && abs(rapidity()) < 2.4")
0053                                          )
0054
0055 ak8PFJetsPuppiConstituents = cms.EDProducer("PFJetConstituentSelector",
0056                                           src = cms.InputTag("ak8PFJetsPuppi"),
0057                                           cut = cms.string("pt > 100.0 && abs(rapidity()) < 2.4")
0058                                          )
0059
0060
0061 ##############################################################################
0062 # Substructure algorithms
0063 ##############################################################################
0064 ak8PFJetsCHSFiltered = ak8PFJets.clone(
0065     src = "ak8PFJetsCHSConstituents:constituents",
0066     useFiltering = cms.bool(True),
0067     nFilt = cms.int32(3),
0068     rFilt = cms.double(0.3),
0069     useExplicitGhosts = cms.bool(True),
0070     writeCompound = cms.bool(True),
0071     jetCollInstanceName=cms.string("SubJets"),
0072     jetPtMin = 100.0
0073     )
0074
0075
0076 ak8PFJetsCHSMassDropFiltered = ak8PFJets.clone(
0077     src = "ak8PFJetsCHSConstituents:constituents",
0078     useMassDropTagger = cms.bool(True),
0079     muCut = cms.double(0.667),
0080     yCut = cms.double(0.08),
0081     writeCompound = cms.bool(True),
0082     jetCollInstanceName=cms.string("SubJets"),
0083     jetPtMin = 100.0
0084     )
0085
0086 ak8PFJetsCHSPruned = ak8PFJets.clone(
0087     SubJetParameters,
0088     src = "ak8PFJetsCHSConstituents:constituents",
0089     usePruning = cms.bool(True),
0090     useExplicitGhosts = cms.bool(True),
0091     writeCompound = cms.bool(True),
0092     jetCollInstanceName=cms.string("SubJets"),
0093     jetPtMin = 100.0,
0094     doAreaFastjet = False
0095     )
0096
0097 ak8PFJetsCHSSoftDrop = ak8PFJets.clone(
0098     useSoftDrop = cms.bool(True),
0099     src = "ak8PFJetsCHSConstituents:constituents",
0100     zcut = cms.double(0.1),
0101     beta = cms.double(0.0),
0102     R0   = cms.double(0.8),
0103     useExplicitGhosts = cms.bool(True),
0104     writeCompound = cms.bool(True),
0105     jetCollInstanceName=cms.string("SubJets"),
0106     jetPtMin = 100.0
0107     )
0108
0109
0110 ak8PFJetsCHSTrimmed = ak8PFJets.clone(
0111     useTrimming = cms.bool(True),
0112     src = "ak8PFJetsCHSConstituents:constituents",
0113     rFilt = cms.double(0.2),
0114     trimPtFracMin = cms.double(0.03),
0115     useExplicitGhosts = cms.bool(True),
0116     jetPtMin = 100.0
0117     )
0118
0119 ak8PFJetsPuppiSoftDrop = ak8PFJetsCHSSoftDrop.clone(
0120     src = "ak8PFJetsPuppiConstituents:constituents",
0121     applyWeight = True,
0122     srcWeights = cms.InputTag("puppi")
0123     )