"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),
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
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")
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),
############################### # 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(
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.),