# - 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 ) ak8PFJetsCHS = ak5PFJetsCHS.clone( rParam = 0.8, jetPtMin = 15.0 )
patJetsCA8CHSpruned.addGenJetMatch = False patJetsCA8CHSpruned.getJetMCFlavour = False patJetsCA8CHSpruned.jetCorrFactorsSource = cms.VInputTag(cms.InputTag('patJetCorrFactorsCA8CHSpruned')) #### Adding Nsubjetiness patJetsCA8CHSprunedwithNsub = cms.EDProducer("NjettinessAdder", src=cms.InputTag("patJetsCA8CHSpruned"), cone=cms.double(0.8) ) ################################################################ TRIMMING from RecoJets.JetProducers.ak5PFJetsTrimmed_cfi import ak5PFJetsTrimmed ca8PFJetsCHStrimmed = ak5PFJetsTrimmed.clone( src = 'pfNoPileUp', jetPtMin = cms.double(20.0), doAreaFastjet = cms.bool(True), rParam = cms.double(0.8), jetAlgorithm = cms.string("CambridgeAachen"), ) jetSourceTrimmed = 'ca8PFJetsCHStrimmed' # corrections patJetCorrFactorsCA8CHStrimmed = patJetCorrFactors.clone() patJetCorrFactorsCA8CHStrimmed.src = jetSourceTrimmed patJetCorrFactorsCA8CHStrimmed.levels = ['L1FastJet', 'L2Relative', 'L3Absolute'] patJetCorrFactorsCA8CHStrimmed.payload = 'AK7PFchs' patJetCorrFactorsCA8CHStrimmed.useRho = True ### pats patJetsCA8CHStrimmed = patJets.clone() patJetsCA8CHStrimmed.jetSource = jetSourceTrimmed
############################### ###### AK 0.7 jets ############ ############################### process.ak7PFlow = process.pfJetsPFlow.clone(rParam=cms.double(0.7)) ############################### ###### AK 0.8 jets ############ ############################### process.ak8PFlow = process.pfJetsPFlow.clone(rParam=cms.double(0.8)) ############################### ###### AK 0.5 jets groomed #### ############################### from RecoJets.JetProducers.ak5PFJetsTrimmed_cfi import ak5PFJetsTrimmed process.ak5TrimmedPFlow = ak5PFJetsTrimmed.clone(src=process.pfJetsPFlow.src, doAreaFastjet=cms.bool(True)) from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsFiltered process.ak5FilteredPFlow = ak5PFJetsFiltered.clone( src=process.pfJetsPFlow.src, doAreaFastjet=cms.bool(True)) from RecoJets.JetProducers.ak5PFJetsPruned_cfi import ak5PFJetsPruned process.ak5PrunedPFlow = ak5PFJetsPruned.clone(src=process.pfJetsPFlow.src, doAreaFastjet=cms.bool(True)) ############################### ###### AK 0.7 jets groomed #### ############################### process.ak7TrimmedPFlow = process.ak5TrimmedPFlow.clone( src=process.pfJetsPFlow.src, rParam=cms.double(0.7))
process.pvCount = cms.EDFilter( "EDPileupAna", src=cms.InputTag('goodOfflinePrimaryVertices') ) ############################### ###### AK 0.5 jets groomed #### ############################### ak5JetSrc = cms.InputTag("selectedPatJetsPFlow","pfCandidates") from RecoJets.JetProducers.ak5PFJetsTrimmed_cfi import ak5PFJetsTrimmed process.ak5TrimmedPFlow = ak5PFJetsTrimmed.clone( src = ak5JetSrc, doAreaFastjet = cms.bool(True) ) from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsFiltered process.ak5FilteredPFlow = ak5PFJetsFiltered.clone( src = ak5JetSrc, doAreaFastjet = cms.bool(True) ) from RecoJets.JetProducers.ak5PFJetsPruned_cfi import ak5PFJetsPruned process.ak5PrunedPFlow = ak5PFJetsPruned.clone( src = ak5JetSrc, doAreaFastjet = cms.bool(True) )
patJetsAK4CHSpruned.getJetMCFlavour = False patJetsAK4CHSpruned.jetCorrFactorsSource = cms.VInputTag( cms.InputTag('patJetCorrFactorsAK4CHSpruned')) #### Adding Nsubjetiness patJetsAK4CHSprunedwithNsub = cms.EDProducer( "NjettinessAdder", src=cms.InputTag("patJetsAK4CHSpruned"), cone=cms.double(0.4)) ################################################################ TRIMMING from RecoJets.JetProducers.ak5PFJetsTrimmed_cfi import ak5PFJetsTrimmed ak4PFJetsCHStrimmed = ak5PFJetsTrimmed.clone( src='pfNoPileUp', jetPtMin=cms.double(20.0), doAreaFastjet=cms.bool(True), rParam=cms.double(0.4), jetAlgorithm=cms.string("AntiKt"), ) jetSourceTrimmed = 'ak4PFJetsCHStrimmed' # corrections patJetCorrFactorsAK4CHStrimmed = patJetCorrFactors.clone() patJetCorrFactorsAK4CHStrimmed.src = jetSourceTrimmed patJetCorrFactorsAK4CHStrimmed.levels = [ 'L1FastJet', 'L2Relative', 'L3Absolute' ] patJetCorrFactorsAK4CHStrimmed.payload = 'AK5PFchs' patJetCorrFactorsAK4CHStrimmed.useRho = True ### pats
jetPtMin = cms.double(100.0) ) ak8PFJetsCSConstituents = cms.EDProducer("PFJetConstituentSelector", src = cms.InputTag("ak8PFJetsCS"), cut = cms.string("pt > 100.0") ) ak8PFJetsCSPruned = ak5PFJetsPruned.clone( rParam = 0.8, jetPtMin = 15.0, src = cms.InputTag("ak8PFJetsCSConstituents", "constituents") ) ak8PFJetsCSTrimmed = ak5PFJetsTrimmed.clone( rParam = 0.8, jetPtMin = 15.0, src = cms.InputTag("ak8PFJetsCSConstituents", "constituents") ) ak8PFJetsCSFiltered = ak5PFJetsFiltered.clone( rParam = 0.8, jetPtMin = 15.0, src = cms.InputTag("ak8PFJetsCSConstituents", "constituents") )
patJetsAK7CHSpruned.addGenJetMatch = False patJetsAK7CHSpruned.getJetMCFlavour = False patJetsAK7CHSpruned.jetCorrFactorsSource = cms.VInputTag(cms.InputTag('patJetCorrFactorsAK7CHSpruned')) #### Adding Nsubjetiness patJetsAK7CHSprunedwithNsub = cms.EDProducer("NjettinessAdder", src=cms.InputTag("patJetsAK7CHSpruned"), cone=cms.double(0.7) ) ################################################################ TRIMMING from RecoJets.JetProducers.ak5PFJetsTrimmed_cfi import ak5PFJetsTrimmed ak7PFJetsCHStrimmed = ak5PFJetsTrimmed.clone( src = 'pfNoPileUp', jetPtMin = cms.double(20.0), doAreaFastjet = cms.bool(True), rParam = cms.double(0.7), jetAlgorithm = cms.string("AntiKt"), ) jetSourceTrimmed = 'ak7PFJetsCHStrimmed' # corrections patJetCorrFactorsAK7CHStrimmed = patJetCorrFactors.clone() patJetCorrFactorsAK7CHStrimmed.src = jetSourceTrimmed patJetCorrFactorsAK7CHStrimmed.levels = ['L1FastJet', 'L2Relative', 'L3Absolute'] patJetCorrFactorsAK7CHStrimmed.payload = 'AK7PFchs' patJetCorrFactorsAK7CHStrimmed.useRho = True ### pats patJetsAK7CHStrimmed = patJets.clone() patJetsAK7CHStrimmed.jetSource = jetSourceTrimmed
selectedPatJetsCA8CHSprunedPre + btaggingCA8CHSprunedSubjets + patJetsCA8CHSprunedSubjets + selectedPatJetsCA8CHSpruned ) # JETS TRIMMED CA8 ---------------------------- from RecoJets.JetProducers.ak5PFJetsTrimmed_cfi import ak5PFJetsTrimmed ca8PFJetsCHStrimmed = ak5PFJetsTrimmed.clone( src = 'pfNoPileUp', jetPtMin = cms.double(30.0), doAreaFastjet = cms.bool(True), rParam = cms.double(0.8), jetAlgorithm = cms.string("CambridgeAachen"), rFilt = cms.double(0.1), trimPtFracMin = cms.double(0.03), writeCompound = cms.bool(True), jetCollInstanceName = cms.string("SubJets") ) jetSource = 'ca8PFJetsCHStrimmed' # corrections from PhysicsTools.PatAlgos.recoLayer0.jetCorrFactors_cfi import * patJetCorrFactorsCA8CHStrimmed = patJetCorrFactors.clone() patJetCorrFactorsCA8CHStrimmed.src = jetSource # will need to add L2L3 corrections in the cfg patJetCorrFactorsCA8CHStrimmed.levels = ['L1FastJet', 'L2Relative', 'L3Absolute'] patJetCorrFactorsCA8CHStrimmed.payload = 'AK7PFchs'
getattr(process,"pfNoTau"+postfix).enable = False getattr(process,"pfNoJet"+postfix).enable = True ## Remove MC matching when running over data if options.runOnData: removeMCMatching( process, ['All'] ) ## Define AK6 jets and subjets (GEN and RECO) process.ak6GenJetsNoNu = process.ak5GenJetsNoNu.clone( rParam = 0.6 ) process.ak6PFJets = process.pfJetsPFlow.clone( rParam = 0.6 ) ## AK6 groomed jets from RecoJets.JetProducers.ak5PFJetsTrimmed_cfi import ak5PFJetsTrimmed process.ak6PFJetsTrimmed = ak5PFJetsTrimmed.clone( src = process.pfJetsPFlow.src, doAreaFastjet = cms.bool(True), rParam = cms.double(0.6) ) from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsFiltered process.ak6PFJetsFiltered = ak5PFJetsFiltered.clone( src = process.pfJetsPFlow.src, doAreaFastjet = cms.bool(True), rParam = cms.double(0.6), writeCompound = cms.bool(False), jetCollInstanceName=cms.string("") ) from RecoJets.JetProducers.ak5PFJetsPruned_cfi import ak5PFJetsPruned process.ak6PFJetsPruned = ak5PFJetsPruned.clone( src = process.pfJetsPFlow.src,
rParam=cms.double(0.8), useConstituentSubtraction=cms.bool(True), csRParam=cms.double(0.4), csRho_EtaMax=ak4PFJets. Rho_EtaMax, # Just use the same eta for both C.S. and rho by default useExplicitGhosts=cms.bool(True), doAreaFastjet=cms.bool(True), jetPtMin=cms.double(100.0)) ak8PFJetsCSConstituents = cms.EDProducer("PFJetConstituentSelector", src=cms.InputTag("ak8PFJetsCS"), cut=cms.string("pt > 100.0")) ak8PFJetsCSPruned = ak5PFJetsPruned.clone(rParam=0.8, jetPtMin=15.0, src=cms.InputTag( "ak8PFJetsCSConstituents", "constituents")) ak8PFJetsCSTrimmed = ak5PFJetsTrimmed.clone(rParam=0.8, jetPtMin=15.0, src=cms.InputTag( "ak8PFJetsCSConstituents", "constituents")) ak8PFJetsCSFiltered = ak5PFJetsFiltered.clone(rParam=0.8, jetPtMin=15.0, src=cms.InputTag( "ak8PFJetsCSConstituents", "constituents"))
src = process.ca8PFJetsCHS.src, srcPVs = process.ca8PFJetsCHS.srcPVs, doAreaFastjet = process.ca8PFJetsCHS.doAreaFastjet, writeCompound = cms.bool(True), jetCollInstanceName=cms.string("SubJets"), jetPtMin = cms.double(20.), useKtPruning = cms.bool(True), zcut = cms.double(0.), rcut_factor = cms.double(9999.) ) ## CA8 trimmed jets (Reco only) from RecoJets.JetProducers.ak5PFJetsTrimmed_cfi import ak5PFJetsTrimmed process.ca8PFJetsCHSTrimmed = ak5PFJetsTrimmed.clone( jetAlgorithm = cms.string("CambridgeAachen"), rParam = cms.double(0.8), src = process.ca8PFJetsCHS.src, srcPVs = process.ca8PFJetsCHS.srcPVs, doAreaFastjet = process.ca8PFJetsCHS.doAreaFastjet, jetPtMin = cms.double(20.) ) #------------------------------------- ## PATify the above jets from PhysicsTools.PatAlgos.tools.jetTools import * ## CA8 jets switchJetCollection(process, cms.InputTag('ca8PFJetsCHS'), jetIdLabel='ca8', rParam = 0.8, useLegacyFlavour=False, doJTA=options.doJTA, doBTagging=options.doBTagging,