Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 13:35:06

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     )