def jetFlavour(process): process.load("PhysicsTools.JetMCAlgos.CaloJetsMCFlavour_cfi") process.PFAK5byRef = process.AK5byRef.clone( jets = cms.InputTag("ak5PFJetsCorr") ) process.PFAK5byValAlgo = process.AK5byValAlgo.clone( srcByReference = cms.InputTag("PFAK5byRef") ) process.JetSequence += process.myPartons process.JetSequence += process.PFAK5byRef process.JetSequence += process.PFAK5byValAlgo # get gen jets without neutrinos from RecoJets.Configuration.GenJetParticles_cff import genParticlesForJetsNoNu process.genParticlesForJetsNoNu = genParticlesForJetsNoNu from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ak5GenJetsNoNu = ak5GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) process.JetSequence += process.genParticlesForJetsNoNu process.JetSequence += process.ak5GenJetsNoNu
) process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) ) #process.SimpleMemoryCheck = cms.Service("SimpleMemoryCheck", # ignoreTotal = cms.untracked.int32(1) ## default is one #) process.load('Configuration.StandardSequences.Reconstruction_cff') process.load('JetMETCorrections.Configuration.DefaultJEC_cff') process.load('RecoJets.Configuration.GenJetParticles_cff') from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ak5GenJetsNoNu = ak5GenJets.clone(src = cms.InputTag("genParticlesForJetsNoNu")) process.ak5GenJetsSequence = cms.Sequence(process.genParticlesForJetsNoNu*process.ak5GenJetsNoNu) from RecoJets.Configuration.RecoGenJets_cff import * from RecoJets.Configuration.RecoJets_cff import * from RecoJets.Configuration.RecoPFJets_cff import * process.kt6PFJets.doRhoFastjet = True process.ak5PFJets.doAreaFastjet = True process.load("Configuration.StandardSequences.Services_cff") process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff") process.GlobalTag.globaltag = 'START53_V29B::All' #process.GlobalTag.globaltag = 'START53_V7A::All'
maxZ=cms.double(24.0), minNdof=cms.double(4.0) # this is >= 4 ), src=cms.InputTag(pvSrc)) ############################### ########## Gen Setup ########## ############################### process.load("RecoJets.Configuration.GenJetParticles_cff") from RecoJets.JetProducers.ca4GenJets_cfi import ca4GenJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ca8GenJetsNoNu = ca4GenJets.clone( rParam=cms.double(0.8), src=cms.InputTag("genParticlesForJetsNoNu")) process.ak8GenJetsNoNu = ak5GenJets.clone( rParam=cms.double(0.8), src=cms.InputTag("genParticlesForJetsNoNu")) process.load("TopQuarkAnalysis.TopEventProducers.sequences.ttGenEvent_cff") # add the flavor history process.load("PhysicsTools.HepMCCandAlgos.flavorHistoryPaths_cfi") # prune gen particles process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi") process.prunedGenParticles = cms.EDProducer( "GenParticlePruner", src=cms.InputTag("genParticles"), select=cms.vstring( "drop *", "keep status = 3" #keeps particles from the hard matrix element ,
9900012, 9900014, 9900016, 39, ), partonicFinalState=cms.bool(False), excludeResonances=cms.bool(True), excludeFromResonancePids=cms.vuint32(12, 13, 14, 16), tausAsJets=cms.bool(False), ) from RecoJets.JetProducers.ca4GenJets_cfi import ca4GenJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ak5GenJets = ak5GenJets.clone() process.ak5GenJetsNoNu = ak5GenJets.clone(src=cms.InputTag("genParticlesForJetsNoNu")) process.ca8GenJetsNoNu = ca4GenJets.clone(rParam=cms.double(0.8), src=cms.InputTag("genParticlesForJetsNoNu")) process.ak8GenJetsNoNu = ak5GenJets.clone(rParam=cms.double(0.8), src=cms.InputTag("genParticlesForJetsNoNu")) process.load("PhysicsTools.HepMCCandAlgos.flavorHistoryPaths_cfi") process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi") process.prunedGenParticles = cms.EDProducer( "GenParticlePruner", src=cms.InputTag("genParticles"), select=cms.vstring( "drop *", "keep status = 3", # keeps particles from the hard matrix element "keep (abs(pdgId) >= 11 & abs(pdgId) <= 16) & status = 1", # keeps e/mu and nus with status 1 "keep (abs(pdgId) = 15) & status = 3", # keeps taus ),
###### Bare KT 0.6 jets ####### ############################### #from RecoJets.JetProducers.kt4PFJets_cfi import kt4PFJets # For electron (effective area) isolation #process.kt6PFJetsForIsolation = kt4PFJets.clone( rParam = 0.6, doRhoFastjet = True ) #process.kt6PFJetsForIsolation.Rho_EtaMax = cms.double(2.5) ############################### ### Add AK5GenJetsNoMuNoNu #### ############################### from RecoJets.Configuration.GenJetParticles_cff import genParticlesForJetsNoMuNoNu from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ak5GenJetsNoMuNoNu = ak5GenJets.clone( src = cms.InputTag('genParticlesForJetsNoMuNoNu') ) process.ak5GenJetsSeq = cms.Sequence(genParticlesForJetsNoMuNoNu*process.ak5GenJetsNoMuNoNu) ############################### #### Selections Setup ######### ############################### # AK5 Jets # PF process.selectedPatJetsPF2PAT.cut = cms.string("pt > 10") process.selectedPatJetsPF2PATNoLeptonCleaning.cut = cms.string("pt > 10") #process.selectedPatJetsPF2PATNoPFnoPU.cut = cms.string("pt > 10") #process.selectedPatJetsAK5Calo.cut = cms.string("pt > 15") # Flavor history stuff process.load("PhysicsTools.HepMCCandAlgos.flavorHistoryPaths_cfi")
process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(1000)) from RecoJets.JetProducers.ak5PFJets_cfi import ak5PFJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.chs = cms.EDFilter("CandPtrSelector", src=cms.InputTag("packedPFCandidates"), cut=cms.string("fromPV")) process.ak5PFJets = ak5PFJets.clone( src='packedPFCandidates', doAreaFastjet=True ) # no idea while doArea is false by default, but it's True in RECO so we have to set it process.ak5PFJetsCHS = ak5PFJets.clone( src='chs', doAreaFastjet=True ) # no idea while doArea is false by default, but it's True in RECO so we have to set it process.ak5GenJets = ak5GenJets.clone(src='packedGenParticles') process.load("PhysicsTools.PatAlgos.producersLayer1.patCandidates_cff") process.load("Configuration.EventContent.EventContent_cff") process.load('Configuration.StandardSequences.Geometry_cff') process.load('Configuration.StandardSequences.MagneticField_38T_cff') process.load( 'Configuration.StandardSequences.FrontierConditions_GlobalTag_cff') process.GlobalTag.globaltag = 'START70_V6::All' from PhysicsTools.PatAlgos.tools.jetTools import addJetCollection addJetCollection(process, postfix="", labelName='AK5PFCHS', jetSource=cms.InputTag('ak5PFJetsCHS'), trackSource=cms.InputTag('unpackedTracksAndVertices'),
import FWCore.ParameterSet.Config as cms from RecoJets.JetProducers.sc5GenJets_cfi import sisCone5GenJets from RecoJets.JetProducers.ic5GenJets_cfi import iterativeCone5GenJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets from RecoJets.JetProducers.gk5GenJets_cfi import gk5GenJets from RecoJets.JetProducers.kt4GenJets_cfi import kt4GenJets from RecoJets.JetProducers.ca4GenJets_cfi import ca4GenJets from RecoHI.HiJetAlgos.HiGenJets_cff import * sisCone7GenJets = sisCone5GenJets.clone( rParam = 0.7 ) ak7GenJets = ak5GenJets.clone( rParam = 0.7 ) gk7GenJets = gk5GenJets.clone( rParam = 0.7 ) kt6GenJets = kt4GenJets.clone( rParam = 0.6 ) ca6GenJets = ca4GenJets.clone( rParam = 0.6 ) sisCone5GenJetsNoNu = sisCone5GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) sisCone7GenJetsNoNu = sisCone7GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) kt4GenJetsNoNu = kt4GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) kt6GenJetsNoNu = kt6GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) iterativeCone5GenJetsNoNu = iterativeCone5GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) ak5GenJetsNoNu = ak5GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) ak7GenJetsNoNu = ak7GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) gk5GenJetsNoNu = gk5GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) gk7GenJetsNoNu = gk7GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) ca4GenJetsNoNu = ca4GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) ca6GenJetsNoNu = ca6GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") )
import FWCore.ParameterSet.Config as cms from RecoJets.JetProducers.sc5GenJets_cfi import sisCone5GenJets from RecoJets.JetProducers.ic5GenJets_cfi import iterativeCone5GenJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets from RecoJets.JetProducers.gk5GenJets_cfi import gk5GenJets from RecoJets.JetProducers.kt4GenJets_cfi import kt4GenJets from RecoJets.JetProducers.ca4GenJets_cfi import ca4GenJets from RecoHI.HiJetAlgos.HiGenJets_cff import * sisCone8GenJets = sisCone5GenJets.clone(rParam=0.8) ak4GenJets = ak5GenJets.clone(rParam=0.4) ak8GenJets = ak5GenJets.clone(rParam=0.8) gk8GenJets = gk5GenJets.clone(rParam=0.8) kt6GenJets = kt4GenJets.clone(rParam=0.6) ca6GenJets = ca4GenJets.clone(rParam=0.6) sisCone5GenJetsNoNu = sisCone5GenJets.clone( src=cms.InputTag("genParticlesForJetsNoNu")) sisCone8GenJetsNoNu = sisCone8GenJets.clone( src=cms.InputTag("genParticlesForJetsNoNu")) kt4GenJetsNoNu = kt4GenJets.clone(src=cms.InputTag("genParticlesForJetsNoNu")) kt6GenJetsNoNu = kt6GenJets.clone(src=cms.InputTag("genParticlesForJetsNoNu")) iterativeCone5GenJetsNoNu = iterativeCone5GenJets.clone( src=cms.InputTag("genParticlesForJetsNoNu")) ak5GenJetsNoNu = ak5GenJets.clone(src=cms.InputTag("genParticlesForJetsNoNu")) ak8GenJetsNoNu = ak8GenJets.clone(src=cms.InputTag("genParticlesForJetsNoNu")) gk5GenJetsNoNu = gk5GenJets.clone(src=cms.InputTag("genParticlesForJetsNoNu")) gk8GenJetsNoNu = gk8GenJets.clone(src=cms.InputTag("genParticlesForJetsNoNu")) ca4GenJetsNoNu = ca4GenJets.clone(src=cms.InputTag("genParticlesForJetsNoNu"))
# 1000022, # 1000012, 1000014, 1000016, # 2000012, 2000014, 2000016, # 1000039, 5100039, # 4000012, 4000014, 4000016, # 9900012, 9900014, 9900016, # 39), # partonicFinalState = cms.bool(False), # excludeResonances = cms.bool(True), # excludeFromResonancePids = cms.vuint32(12, 13, 14, 16), # tausAsJets = cms.bool(False) #) from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ak5GenJetsptmin2 = ak5GenJets.clone() process.ak5GenJetsptmin2.jetPtMin = cms.double(2.0) process.new_ak5GenJets = cms.Sequence(process.genParticlesForJets* process.ak5GenJetsptmin2) #process.new_ak5GenJets = cms.Sequence(process.ak5GenJetsptmin2) # --- to recover the ak7 GenJets: process.ak7GenJetsptmin2 = ak5GenJets.clone() process.ak7GenJetsptmin2.rParam = cms.double(0.7) process.ak7GenJetsptmin2.jetPtMin = cms.double(2.0) process.new_ak7GenJets = cms.Sequence(process.genParticlesForJets* process.ak7GenJetsptmin2) #process.new_ak7GenJets = cms.Sequence(process.ak7GenJetsptmin2) process.p = cms.Path(process.new_ak5GenJets*process.new_ak7GenJets*process.monster*process.myanalysis)
def reSetJet(process): ## Filter out neutrinos from packed GenParticles process.packedGenParticlesForJetsNoNu = cms.EDFilter( "CandPtrSelector", src=cms.InputTag("packedGenParticles"), cut=cms.string( "abs(pdgId) != 12 && abs(pdgId) != 14 && abs(pdgId) != 16")) ## Define GenJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ak5GenJetsNoNu = ak5GenJets.clone( src='packedGenParticlesForJetsNoNu') ## Select charged hadron subtracted packed PF candidates process.pfCHS = cms.EDFilter("CandPtrSelector", src=cms.InputTag("packedPFCandidates"), cut=cms.string("fromPV")) from RecoJets.JetProducers.ak5PFJets_cfi import ak5PFJets ## Define PFJetsCHS process.ak5PFJetsCHS = ak5PFJets.clone(src='pfCHS', doAreaFastjet=True) ################################################# ## Remake PAT jets ################################################# ## b-tag discriminators bTagDiscriminators = ['pfCombinedInclusiveSecondaryVertexV2BJetTags'] ################################################################# slimmedAddPileupInfo = cms.EDProducer( 'PileupSummaryInfoSlimmer', #src = cms.InputTag('addPileupInfo'), keepDetailedInfoFor=cms.vint32(0)) ####################################################### from PhysicsTools.PatAlgos.tools.jetTools import addJetCollection ## Add PAT jet collection based on the above-defined ak5PFJetsCHS addJetCollection( process, labelName='AK5PFCHS', jetSource=cms.InputTag('ak5PFJetsCHS'), pvSource=cms.InputTag('offlineSlimmedPrimaryVertices'), bsSource=cms.InputTag('offlineBeamSpot'), pfCandidates=cms.InputTag('packedPFCandidates'), svSource=cms.InputTag('slimmedSecondaryVertices'), btagDiscriminators=bTagDiscriminators, #jetCorrections = ('AK5PFchs', ['L1FastJet', 'L2Relative', 'L3Absolute'], 'None'), #commented before adding jet corrections genJetCollection=cms.InputTag('ak5GenJetsNoNu'), genParticles=cms.InputTag('prunedGenParticles'), algo='AK', rParam=0.5) getattr(process, 'selectedPatJetsAK5PFCHS').cut = cms.string('pt > 20.') process.rejet = cms.Path( process.packedGenParticlesForJetsNoNu * process.ak5GenJetsNoNu * process.pfCHS * process.ak5PFJetsCHS * process.patJetCorrFactorsAK5PFCHS * process.patJetPartons * process.patJetFlavourAssociationAK5PFCHS * process.patJetPartonMatchAK5PFCHS * process.patJetGenJetMatchAK5PFCHS * process.pfImpactParameterTagInfosAK5PFCHS * process.pfInclusiveSecondaryVertexFinderTagInfosAK5PFCHS * process.pfCombinedInclusiveSecondaryVertexV2BJetTagsAK5PFCHS * process.patJetsAK5PFCHS * process.selectedPatJetsAK5PFCHS) from PhysicsTools.PatAlgos.tools.pfTools import adaptPVs ## Adapt primary vertex collection and BeamSpot adaptPVs(process, pvCollection=cms.InputTag('offlineSlimmedPrimaryVertices')) adaptBSs(process, bsCollection=cms.InputTag('offlineBeamSpot')) from PhysicsTools.PatAlgos.tools.jetTools import updateJetCollection updateJetCollection( process, jetSource=cms.InputTag('slimmedJets'), labelName='UpdatedJEC', jetCorrections=('AK5PFchs', cms.vstring(['L1FastJet', 'L2Relative', 'L3Absolute']), 'None') # Do not forget 'L2L3Residual' on data ) process.load('Configuration.StandardSequences.Services_cff' ) #Added info from now onwards process.load("JetMETCorrections.Modules.JetResolutionESProducer_cfi") from CondCore.DBCommon.CondDBSetup_cfi import * import os process.jer = cms.ESSource( "PoolDBESSource", CondDBSetup, toGet=cms.VPSet( # Resolution cms.PSet( record=cms.string('JetResolutionRcd'), tag=cms.string('JR_Fall15_25nsV2_MC_PtResolution_AK5PFchs'), label=cms.untracked.string('AK5PFchs_pt')), # Scale factors cms.PSet(record=cms.string('JetResolutionScaleFactorRcd'), tag=cms.string('JR_Fall15_25nsV2_MC_SF_AK4PFchs'), label=cms.untracked.string('AK5PFchs')), ), connect=cms.string('sqlite:Fall15_25nsV2_MC.db')) process.es_prefer_jer = cms.ESPrefer('PoolDBESSource', 'jer')
genParticleCollection = 'genParticles' genJetInputParticleCollection = 'genParticlesForJets' ## producing a subset of genParticles to be used for jet clustering in AOD from RecoJets.Configuration.GenJetParticles_cff import genParticlesForJets process.genParticlesForJets = genParticlesForJets.clone() else: genParticleCollection = 'prunedGenParticles' genJetInputParticleCollection = 'packedGenParticles' # Supplies PDG ID to real name resolution of MC particles process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi") # Producing own jets for testing purposes from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ak5GenJetsCustom = ak5GenJets.clone(src=genJetInputParticleCollection, rParam=cms.double(0.5), jetAlgorithm=cms.string("AntiKt")) # Ghost particle collection used for Hadron-Jet association # MUST use proper input particle collection from PhysicsTools.JetMCAlgos.HadronAndPartonSelector_cfi import selectedHadronsAndPartons process.selectedHadronsAndPartons = selectedHadronsAndPartons.clone( particles=genParticleCollection) # Input particle collection for matching to gen jets (partons + leptons) # MUST use use proper input jet collection: the jets to which hadrons should be associated # rParam and jetAlgorithm MUST match those used for jets to be associated with hadrons # More details on the tool: https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideBTagMCTools#New_jet_flavour_definition from PhysicsTools.JetMCAlgos.sequences.GenHFHadronMatching_cff import genJetFlavourPlusLeptonInfos process.genJetFlavourPlusLeptonInfos = genJetFlavourPlusLeptonInfos.clone( jets=genJetCollection,
# load low-pt charged GenJet configuration here # get charged genParticles process.load('QCDAnalysis.UEAnalysis.UEAnalysisParticles_cfi') from RecoJets.JetProducers.sc5GenJets_cfi import sisCone5GenJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets from RecoJets.JetProducers.FastjetParameters_cfi import * process.load("RecoJets.Configuration.GenJetParticles_cff") process.chargeParticles.cut = cms.string('charge != 0 & pt > 0.5 & status = 1') process.sisCone5ChgGenJets = sisCone5GenJets.clone(rParam = 0.5, jetPtMin=1.0, src = cms.InputTag("chargeParticles"), inputEtMin = cms.double(0.5) ) process.sisCone7ChgGenJets = sisCone5GenJets.clone(rParam = 0.7, jetPtMin=1.0, src = cms.InputTag("chargeParticles"), inputEtMin = cms.double(0.5) ) process.ak5ChgGenJets = ak5GenJets.clone(rParam = 0.5, jetPtMin=1.0, src = cms.InputTag("chargeParticles"), inputEtMin = cms.double(0.5) ) process.ak4ChgGenJets = ak5GenJets.clone(rParam = 0.4, jetPtMin=1.0, src = cms.InputTag("chargeParticles"), inputEtMin = cms.double(0.5) ) process.ak7ChgGenJets = ak5GenJets.clone(rParam = 0.7, jetPtMin=1.0, src = cms.InputTag("chargeParticles"), inputEtMin = cms.double(0.5) ) process.ak10ChgGenJets = ak5GenJets.clone(rParam = 1., jetPtMin=1.0, src = cms.InputTag("chargeParticles"), inputEtMin = cms.double(0.5) ) process.chargedgenjets = cms.Path(process.UEAnalysisParticles*process.sisCone5ChgGenJets* process.sisCone7ChgGenJets*process.ak4ChgGenJets *process.ak5ChgGenJets*process.ak7ChgGenJets*process.ak10ChgGenJets) from CommonTools.RecoAlgos.TrackWithVertexRefSelector_cfi import * from RecoJets.JetProducers.TracksForJets_cff import * process.trackWithVertexRefSelector=trackWithVertexRefSelector.clone() process.trackRefsForJets=trackRefsForJets.clone() process.trackRefsForJets.cut= cms.string('charge != 0 & pt > 0.5 & status = 1 & quality = highPurity & ptErrorCut = 0.05 ')
from RecoJets.JetProducers.AnomalousCellParameters_cfi import * from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets from RecoJets.JetProducers.ak5TrackJets_cfi import ak5TrackJets FastjetWithAreaPU = cms.PSet( Active_Area_Repeats = cms.int32(5), GhostArea = cms.double(0.01), Ghost_EtaMax = cms.double(6.0), UE_Subtraction = cms.string('no') ) ueAk5ChgGenJet500 = ak5GenJets.clone( src = cms.InputTag("chargeParticles"), jetPtMin = cms.double(1.0), inputEtMin = cms.double(0.5) ) ueAk5TracksJet500 = ak5TrackJets.clone( src = cms.InputTag("goodTracks"), jetPtMin = cms.double(1.0), inputEtMin = cms.double(0.5) ) #ueAk5TracksJet.jetType = 'BasicJet' UEAnalysisJetsAkOnlyMC = cms.Sequence(ueAk5ChgGenJet500) UEAnalysisJetsAkOnlyReco = cms.Sequence(ueAk5TracksJet500) UEAnalysisJetsAk = cms.Sequence(UEAnalysisJetsAkOnlyMC*UEAnalysisJetsAkOnlyReco)
) ) ############################################### # GEN JETS from RecoJets.JetProducers.ca4GenJets_cfi import ca4GenJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ca8GenJets = ca4GenJets.clone( rParam = cms.double(0.8), src = cms.InputTag(genjetparticles)) process.ca15GenJets = ca4GenJets.clone( rParam = cms.double(1.5), src = cms.InputTag(genjetparticles)) process.ak3GenJets = ak5GenJets.clone( rParam = cms.double(0.3), src = cms.InputTag(genjetparticles)) process.ca15FilteredGenJets = process.ca8GenJets.clone( nFilt = cms.int32(3), rFilt = cms.double(0.3), src = cms.InputTag(genjetparticles), useFiltering = cms.bool(True), useExplicitGhosts = cms.bool(True), writeCompound = cms.bool(True), rParam = cms.double(1.5), jetCollInstanceName=cms.string("SubJets"), jetPtMin = cms.double(100.0) ) process.ca8GenJetsPruned = ca4GenJets.clone( rParam = cms.double(0.8),
from RecoJets.JetProducers.AnomalousCellParameters_cfi import * from RecoJets.JetProducers.CATopJetParameters_cfi import * from RecoJets.JetProducers.GenJetParameters_cfi import * from RecoJets.JetProducers.caTopTaggers_cff import * # CATopTagInfos, HEPTopTagInfos # process.CATopTagInfos = CATopTagInfos.clone() ############################################### # GEN JETS from RecoJets.JetProducers.ca4GenJets_cfi import ca4GenJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ak3GenJets = ak5GenJets.clone( rParam = cms.double(0.3), src = cms.InputTag(genjetparticles)) process.ca8GenJets = ca4GenJets.clone( rParam = cms.double(0.8), src = cms.InputTag(genjetparticles)) process.ak15GenJets = ak5GenJets.clone( rParam = cms.double(1.5), src = cms.InputTag(genjetparticles)) #ADD NO NU process.ak15GenJetsFiltered = ak5GenJets.clone( rParam = cms.double(1.5), src = cms.InputTag(genjetparticles), useFiltering = cms.bool(True), nFilt = cms.int32(3), rFilt = cms.double(0.3),
import FWCore.ParameterSet.Config as cms ### PP RECO does not include R=3 or R=5 jets. ### re-RECO is only possible for PF, RECO is missing calotowers from RecoJets.JetProducers.ak5PFJets_cfi import ak5PFJets ak5PFJets.doAreaFastjet = True ak3PFJets = ak5PFJets.clone(rParam = 0.3) from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets ak3GenJets = ak5GenJets.clone(rParam = 0.3) from RecoJets.Configuration.GenJetParticles_cff import * from RecoHI.HiJetAlgos.HiGenJets_cff import * from HeavyIonsAnalysis.JetAnalysis.makePartons_cff import * from HeavyIonsAnalysis.JetAnalysis.jets.ak3PFJetSequence_pp_mc_cff import * from HeavyIonsAnalysis.JetAnalysis.jets.ak4PFJetSequence_pp_mc_cff import * from HeavyIonsAnalysis.JetAnalysis.jets.ak5PFJetSequence_pp_mc_cff import * from HeavyIonsAnalysis.JetAnalysis.jets.ak4CaloJetSequence_pp_mc_cff import * highPurityTracks = cms.EDFilter("TrackSelector", src = cms.InputTag("generalTracks"), cut = cms.string('quality("highPurity")') ) # Other radii jets and calo jets need to be reconstructed jetSequences = cms.Sequence( genParticlesForJets + ak3GenJets + ak5GenJets + ak3PFJets + ak5PFJets + makePartons +
process.chargeParticles.cut = cms.string('charge != 0 & pt > 0.5 & status = 1') process.sisCone5ChgGenJets = sisCone5GenJets.clone( rParam=0.5, jetPtMin=1.0, src=cms.InputTag("chargeParticles"), inputEtMin=cms.double(0.5)) process.sisCone7ChgGenJets = sisCone5GenJets.clone( rParam=0.7, jetPtMin=1.0, src=cms.InputTag("chargeParticles"), inputEtMin=cms.double(0.5)) process.ak5ChgGenJets = ak5GenJets.clone(rParam=0.5, jetPtMin=1.0, src=cms.InputTag("chargeParticles"), inputEtMin=cms.double(0.5)) process.ak4ChgGenJets = ak5GenJets.clone(rParam=0.4, jetPtMin=1.0, src=cms.InputTag("chargeParticles"), inputEtMin=cms.double(0.5)) process.ak7ChgGenJets = ak5GenJets.clone(rParam=0.7, jetPtMin=1.0, src=cms.InputTag("chargeParticles"), inputEtMin=cms.double(0.5)) process.ak10ChgGenJets = ak5GenJets.clone(rParam=1., jetPtMin=1.0, src=cms.InputTag("chargeParticles"), inputEtMin=cms.double(0.5)) process.chargedgenjets = cms.Path(
def initialise(runOnMC, decayMode, doOutModule=False, doPAT=True): process = cms.Process("KrAFT") process.load("Configuration.StandardSequences.Services_cff") process.load("Configuration.Geometry.GeometryDB_cff") process.load("Configuration.StandardSequences.MagneticField_cff") process.options = cms.untracked.PSet(wantSummary=cms.untracked.bool(True)) process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(-1)) process.load("FWCore.MessageLogger.MessageLogger_cfi") process.MessageLogger.cerr.FwkReport.reportEvery = 10000 process.source = cms.Source("PoolSource", fileNames=cms.untracked.vstring()) process.source.inputCommands = cms.untracked.vstring( "keep *", "drop *_MEtoEDMConverter_*_*") process.load( "Configuration.StandardSequences.FrontierConditions_GlobalTag_cff") from Configuration.AlCa.autoCond import autoCond if runOnMC: process.GlobalTag.globaltag = autoCond['startup'] else: process.GlobalTag.globaltag = autoCond['com10'] outputModuleForTriggerMatch = "" outputModules = [] if doOutModule: from PhysicsTools.PatAlgos.patEventContent_cff import patEventContentNoCleaning process.out = cms.OutputModule( "PoolOutputModule", fileName=cms.untracked.string("out.root"), outputCommands=cms.untracked.vstring( 'drop *', 'keep recoPFCandidates_particleFlow_*_*', *patEventContentNoCleaning)) process.outPath = cms.EndPath(process.out) outputModuleForTriggerMatch = "out" outputModules.append(process.out) ## Load PAT process.load("PhysicsTools.PatAlgos.patSequences_cff") ## Apply MVA process.load('EgammaAnalysis.ElectronTools.electronIdMVAProducer_cfi') process.eidMVASequence = cms.Sequence(process.mvaTrigV0 + process.mvaNonTrigV0) process.patElectrons.electronIDSources.mvaTrigV0 = cms.InputTag( "mvaTrigV0") process.patElectrons.electronIDSources.mvaNonTrigV0 = cms.InputTag( "mvaNonTrigV0") process.patDefaultSequence.replace( process.patElectrons, process.eidMVASequence * process.patElectrons) ## Load trigger matching process.load("KrAFT.Configuration.hltFilters_cff") #from PhysicsTools.PatAlgos.tools.trigTools import * #switchOnTriggerMatchEmbedding(process, outputModule="") ## Apply PF2PAT from PhysicsTools.PatAlgos.tools.pfTools import usePF2PAT if runOnMC: jecLevels = ['L1FastJet', 'L2Relative', 'L3Absolute'] else: jecLevels = ['L1FastJet', 'L2Relative', 'L3Absolute', 'L2L3Residual'] postfix = "PFlow" #usePFBRECO(process,runPFBRECO=True, usePF2PAT(process, runPF2PAT=True, runOnMC=runOnMC, outputModules=outputModules, postfix=postfix, jetAlgo="AK5", jetCorrections=("AK5PFchs", jecLevels), pvCollection=cms.InputTag('goodOfflinePrimaryVertices'), typeIMetCorrections=True) # process.pfPileUpPFlow.checkClosestZVertex = cms.bool(False) # top projections in PF2PAT: process.pfNoPileUpPFlow.enable = True process.pfNoMuonPFlow.enable = True process.pfNoElectronPFlow.enable = True process.pfNoTauPFlow.enable = False process.pfNoJetPFlow.enable = False # verbose flags for the PF2PAT modules process.pfNoMuonPFlow.verbose = False # Use non-isolated muons and electrons process.patMuonsPFlow.pfMuonSource = "pfMuonsPFlow" process.patElectronsPFlow.pfElectronSource = "pfElectronsPFlow" # And turn on delta-beta corrections while building pfIsolated*PFlow process.pfIsolatedMuonsPFlow.doDeltaBetaCorrection = True process.pfIsolatedElectronsPFlow.doDeltaBetaCorrection = True # Change DR cone size to 0.3 process.pfIsolatedMuonsPFlow.isolationValueMapsCharged = cms.VInputTag( cms.InputTag('muPFIsoValueCharged03PFlow')) process.pfIsolatedMuonsPFlow.deltaBetaIsolationValueMap = cms.InputTag( 'muPFIsoValuePU03PFlow') process.pfIsolatedMuonsPFlow.isolationValueMapsNeutral = cms.VInputTag( cms.InputTag('muPFIsoValueNeutral03PFlow'), cms.InputTag('muPFIsoValueGamma03PFlow'), ) process.pfMuonsPFlow.isolationValueMapsCharged = cms.VInputTag( cms.InputTag('muPFIsoValueCharged03PFlow')) process.pfMuonsPFlow.deltaBetaIsolationValueMap = cms.InputTag( 'muPFIsoValuePU03PFlow') process.pfMuonsPFlow.isolationValueMapsNeutral = cms.VInputTag( cms.InputTag('muPFIsoValueNeutral03PFlow'), cms.InputTag('muPFIsoValueGamma03PFlow'), ) process.patMuonsPFlow.isolationValues.pfNeutralHadrons = cms.InputTag( 'muPFIsoValueNeutral03PFlow') process.patMuonsPFlow.isolationValues.pfChargedAll = cms.InputTag( 'muPFIsoValueChargedAll03PFlow') process.patMuonsPFlow.isolationValues.pfPUChargedHadrons = cms.InputTag( 'muPFIsoValuePU03PFlow') process.patMuonsPFlow.isolationValues.pfPhotons = cms.InputTag( 'muPFIsoValueGamma03PFlow') process.patMuonsPFlow.isolationValues.pfChargedHadrons = cms.InputTag( 'muPFIsoValueCharged03PFlow') process.pfIsolatedElectronsPFlow.isolationValueMapsCharged = cms.VInputTag( cms.InputTag('elPFIsoValueCharged03PFIdPFlow')) process.pfIsolatedElectronsPFlow.deltaBetaIsolationValueMap = cms.InputTag( 'elPFIsoValuePU03PFIdPFlow') process.pfIsolatedElectronsPFlow.isolationValueMapsNeutral = cms.VInputTag( cms.InputTag('elPFIsoValueNeutral03PFIdPFlow'), cms.InputTag('elPFIsoValueGamma03PFIdPFlow')) process.pfElectronsPFlow.isolationValueMapsCharged = cms.VInputTag( cms.InputTag('elPFIsoValueCharged03PFIdPFlow')) process.pfElectronsPFlow.deltaBetaIsolationValueMap = cms.InputTag( 'elPFIsoValuePU03PFIdPFlow') process.pfElectronsPFlow.isolationValueMapsNeutral = cms.VInputTag( cms.InputTag('elPFIsoValueNeutral03PFIdPFlow'), cms.InputTag('elPFIsoValueGamma03PFIdPFlow')) process.patElectronsPFlow.isolationValues.pfNeutralHadrons = cms.InputTag( 'elPFIsoValueNeutral03PFIdPFlow') process.patElectronsPFlow.isolationValues.pfChargedAll = cms.InputTag( 'elPFIsoValueChargedAll03PFIdPFlow') process.patElectronsPFlow.isolationValues.pfPUChargedHadrons = cms.InputTag( 'elPFIsoValuePU03PFIdPFlow') process.patElectronsPFlow.isolationValues.pfPhotons = cms.InputTag( 'elPFIsoValueGamma03PFIdPFlow') process.patElectronsPFlow.isolationValues.pfChargedHadrons = cms.InputTag( 'elPFIsoValueCharged03PFIdPFlow') ## Add common filters process.primaryVertexFilter = cms.EDFilter( "GoodVertexFilter", vertexCollection=cms.InputTag('offlinePrimaryVertices'), minimumNDOF=cms.uint32(4), maxAbsZ=cms.double(24), maxd0=cms.double(2)) from PhysicsTools.SelectorUtils.pvSelector_cfi import pvSelector process.goodOfflinePrimaryVertices = cms.EDFilter( "PrimaryVertexObjectFilter", filterParams=pvSelector.clone(maxZ=cms.double(24.0)), src=cms.InputTag('offlinePrimaryVertices')) #process.load( "TopQuarkAnalysis.Configuration.patRefSel_goodVertex_cfi" ) #process.goodOfflinePrimaryVertices.filter = True process.load('TopQuarkAnalysis.Configuration.patRefSel_eventCleaning_cff') process.trackingFailureFilter.VertexSource = cms.InputTag( 'goodOfflinePrimaryVertices') if runOnMC: process.eventCleaning += process.eventCleaningMC else: process.eventCleaning += process.eventCleaningData #reason : removing many events. process.eventCleaning.remove(process.CSCTightHaloFilter) ## Lepton veto filters for L+J channels process.muonVetoFilter = cms.EDFilter( "PATCandViewCountFilter", src=cms.InputTag("selectedPatMuonsPFlow"), maxNumber=cms.uint32(0), minNumber=cms.uint32(0), ) process.electronVetoFilter = cms.EDFilter( "PATCandViewCountFilter", src=cms.InputTag("selectedPatElectronsPFlow"), maxNumber=cms.uint32(0), minNumber=cms.uint32(0), ) # event counters process.nEventsTotal = cms.EDProducer("EventCountProducer") process.nEventsClean = cms.EDProducer("EventCountProducer") process.nEventsPAT = cms.EDProducer("EventCountProducer") process.nEventsHLTElEl = cms.EDProducer("EventCountProducer") process.nEventsHLTMuMu = cms.EDProducer("EventCountProducer") process.nEventsHLTMuEl = cms.EDProducer("EventCountProducer") process.nEventsHLTMuJets = cms.EDProducer("EventCountProducer") process.nEventsHLTElJets = cms.EDProducer("EventCountProducer") ############################### ### Add AK5GenJetsNoNu #### ############################### #process.load("RecoJets.Configuration.GenJetParticles_cff") from RecoJets.Configuration.GenJetParticles_cff import genParticlesForJetsNoNu from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ak5GenJetsNoNu = ak5GenJets.clone( src=cms.InputTag("genParticlesForJetsNoNu")) process.ak5GenJetsSeq = cms.Sequence(process.genParticlesForJets * process.genParticlesForJetsNoNu * process.ak5GenJetsNoNu) # Flavor history stuff process.load("PhysicsTools.HepMCCandAlgos.flavorHistoryPaths_cfi") process.flavorHistoryFilter.pathToSelect = cms.int32(-1) process.cFlavorHistoryProducer.matchedSrc = cms.InputTag("ak5GenJetsNoNu") process.bFlavorHistoryProducer.matchedSrc = cms.InputTag("ak5GenJetsNoNu") #ptMinParticle = cms.double(0.0) #ptMinShower = cms.double(0.0) #process.load( "RecoJets.Configuration.GenJetParticles_cff") #process.load( "RecoJets.Configuration.RecoGenJets_cff") #process.ReGenJetSequence = cms.Sequence( # process.genParticlesForJets # + process.genParticlesForJetsNoNu # + process.ak5GenJets # + process.ak5GenJetsNoNu #) process.commonFilterSequence = cms.Sequence( process.goodOfflinePrimaryVertices * process.eventCleaning * process.ak5GenJetsSeq * process.primaryVertexFilter * process.nEventsClean * getattr(process, "patPF2PATSequence" + postfix) # main PF2PAT * process.nEventsPAT * process.flavorHistorySeq) # process.patSequenceComplete = cms.Sequence( # + process.patDefaultSequence # + process.patPFBRECOSequencePFlow # process.patPF2PATSequencePFlow # + process.nEventsPAT # ) process.load("KrAFT.GeneratorTools.ttbar2bFilter_cfi") ## Defile paths if decayMode in ("all", "dilepton", "ElEl", "ee"): process.pElEl = cms.Path( process.nEventsTotal * process.commonFilterSequence * process.ttbar2bFilter * process.hltElEl * process.nEventsHLTElEl) # if doPAT: process.pElEl += process.patSequenceComplete if decayMode in ("all", "dilepton", "MuMu", "mumu"): process.pMuMu = cms.Path(process.nEventsTotal * process.commonFilterSequence * process.hltMuMu * process.nEventsHLTMuMu) # if doPAT: process.pMuMu += process.patSequenceComplete if decayMode in ("all", "dilepton", "MuEl", "emu"): process.pMuEl = cms.Path(process.nEventsTotal * process.commonFilterSequence * process.hltMuEl * process.nEventsHLTMuEl) # if doPAT: process.pMuEl += process.patSequenceComplete if decayMode in ("all", "MuJets"): process.selectedPatMuonsPFlow.cut = 'isPFMuon && (isGlobalMuon || isTrackerMuon) && pt > 10 && abs(eta) < 2.5 && (chargedHadronIso + max(0.,neutralHadronIso+photonIso-0.5*puChargedHadronIso))/pt < 0.15' process.selectedPatElectronsPFlow.cut = 'pt > 20 && abs(eta) < 2.5 && electronID("mvaTrigV0") > 0. && (chargedHadronIso + max(0.,neutralHadronIso+photonIso-0.5*puChargedHadronIso))/pt < 0.15' process.muonVetoFilter.maxNumber = 1 process.pMuJets = cms.Path( process.nEventsTotal * process.commonFilterSequence * process.hltMuJets * process.nEventsHLTMuJets) # if doPAT: process.pMuJets += process.patSequenceComplete process.pMuJets *= process.muonVetoFilter process.pMuJets += process.electronVetoFilter if runOnMC: process.pMuJets.remove(process.hltMuJets) if decayMode in ("all", "ElJets"): process.selectedPatMuonsPFlow.cut = 'isPFMuon && (isGlobalMuon || isTrackerMuon) && pt > 10 && abs(eta) < 2.5 && (chargedHadronIso + max(0.,neutralHadronIso+photonIso-0.5*puChargedHadronIso))/pt < 0.15' process.selectedPatElectronsPFlow.cut = 'pt > 20 && abs(eta) < 2.5 && electronID("mvaTrigV0") > 0. && (chargedHadronIso + max(0.,neutralHadronIso+photonIso-0.5*puChargedHadronIso))/pt < 0.15' process.electronVetoFilter.maxNumber = 1 process.pElJets = cms.Path( process.nEventsTotal * process.commonFilterSequence * process.hltElJets * process.nEventsHLTElJets) # if doPAT: process.pElJets += process.patSequenceComplete process.pElJets *= process.muonVetoFilter process.pElJets += process.electronVetoFilter if runOnMC: process.pElJets.remove(process.hltElJets) return process
"RecoLocalTracker.SiStripRecHitConverter.SiStripRecHitConverter_cfi") process.load("MLAnalyzer.RecHitAnalyzer.RHAnalyzer_cfi") process.fevt.mode = cms.string(options.processMode) if options.isTTbar: process.fevt.isTTbar = cms.bool(True) else: process.fevt.isTTbar = cms.bool(False) jet_radius = 0.8 from RecoJets.Configuration.GenJetParticles_cff import genParticlesForJets process.genParticlesForJets = genParticlesForJets.clone() from RecoJets.JetProducers.ak5PFJets_cfi import ak5PFJets process.ak8PFJets = ak5PFJets.clone(rParam=jet_radius) process.ak8PFJets.doAreaFastjet = True from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ak8GenJets = ak5GenJets.clone(rParam=jet_radius) if options.UseAK8: process.fevt.PFJetCollection = cms.InputTag('ak8PFJets') process.fevt.genJetCollection = cms.InputTag('ak8GenJets') else: process.fevt.PFJetCollection = cms.InputTag('ak5PFJets') process.fevt.genJetCollection = cms.InputTag('ak5GenJets') #process.fevt.mode = cms.string("JetLevel") # for when using crab #process.fevt.mode = cms.string("EventLevel") # for when using crab print " >> Processing as:", (process.fevt.mode) process.TFileService = cms.Service( "TFileService", fileName=cms.string(options.outputFile)
fatjet_ptmin = 50. ##Dummy copy for having pf candidates into the process #process.packed = cms.EDFilter("CandPtrSelector", src = cms.InputTag("packedGenParticles"), cut = cms.string("abs(pdgId) != 12 && abs(pdgId) != 14 && abs(pdgId) != 16")) ## Filter out neutrinos from packed GenParticles process.packedGenParticlesForJetsNoNu = cms.EDFilter( "CandPtrSelector", src=cms.InputTag("packedGenParticles"), cut=cms.string("abs(pdgId) != 12 && abs(pdgId) != 14 && abs(pdgId) != 16")) ## Define GenJets from RecoJets.JetProducers.ak4GenJets_cfi import ak4GenJets process.ak4GenJetsNoNu = ak4GenJets.clone(src='packedGenParticlesForJetsNoNu') from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ak8GenJetsNoNu = ak5GenJets.clone(src='packedGenParticlesForJetsNoNu', rParam=0.8) ## Standard AK4 and AK8 jets without CHS from RecoJets.JetProducers.ak4PFJets_cfi import ak4PFJets process.ak4PFJets = ak4PFJets.clone( src='packedPFCandidates', doAreaFastjet=True) #HERE!!! add jetPtMin = (fatjet_ptmin) process.ak8PFJets = ak4PFJets.clone(src='packedPFCandidates', rParam=0.8, doAreaFastjet=True, jetPtMin=fatjet_ptmin) ## Pruned AK8 without CHS from RecoJets.JetProducers.ak4PFJetsPruned_cfi import ak4PFJetsPruned process.ak8PFJetsPruned = ak4PFJetsPruned.clone(rParam=0.8, doAreaFastjet=True,
import FWCore.ParameterSet.Config as cms # Remake AK5 jets from packedPFcandidates # See: https://twiki.cern.ch/twiki/bin/view/CMSPublic/WorkBookMiniAOD#Advanced_topics_re_clustering_ev ## Filter out neutrinos from packed GenParticles packedGenParticlesForJetsNoNu = cms.EDFilter( "CandPtrSelector", src=cms.InputTag("packedGenParticles"), cut=cms.string("abs(pdgId) != 12 && abs(pdgId) != 14 && abs(pdgId) != 16")) ## Define GenJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets ak5GenJetsNoNu = ak5GenJets.clone(src='packedGenParticlesForJetsNoNu') ## Select charged hadron subtracted packed PF candidates pfCHS = cms.EDFilter("CandPtrSelector", src=cms.InputTag("packedPFCandidates"), cut=cms.string("fromPV")) from RecoJets.JetProducers.ak5PFJets_cfi import ak5PFJets ## Define ak5PFJetsCHS ak5PFJetsCHS = ak5PFJets.clone(src='pfCHS', doAreaFastjet=True) ## Define ak5PFJets ak5PFJets = ak5PFJets.clone(src='packedPFCandidates') ak5PFJetsSequence = cms.Sequence(packedGenParticlesForJetsNoNu * ak5GenJetsNoNu * pfCHS * ak5PFJetsCHS * ak5PFJets)
import FWCore.ParameterSet.Config as cms from RecoJets.Configuration.GenJetParticles_cff import genParticlesForJetsNoMuNoNu from RecoJets.Configuration.GenJetParticles_cff import genParticlesForJetsNoNu from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets cms2antikt5GenJets = ak5GenJets.clone(src = cms.InputTag("genParticlesForJetsNoMuNoNu") ) cms2antikt5PFGenJets = ak5GenJets.clone(src = cms.InputTag("genParticlesForJetsNoNu") ) genJetSequence = cms.Sequence( genParticlesForJetsNoMuNoNu * cms2antikt5GenJets * genParticlesForJetsNoNu * cms2antikt5PFGenJets)
import FWCore.ParameterSet.Config as cms from RecoJets.JetProducers.sc5GenJets_cfi import sisCone5GenJets from RecoJets.JetProducers.ic5GenJets_cfi import iterativeCone5GenJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets from RecoJets.JetProducers.ak4GenJets_cfi import ak4GenJets from RecoJets.JetProducers.gk5GenJets_cfi import gk5GenJets from RecoJets.JetProducers.kt4GenJets_cfi import kt4GenJets from RecoJets.JetProducers.ca4GenJets_cfi import ca4GenJets from RecoHI.HiJetAlgos.HiGenJets_cff import * sisCone8GenJets = sisCone5GenJets.clone(rParam=0.8) ak8GenJets = ak5GenJets.clone(rParam=0.8, jetPtMin=15.0) gk8GenJets = gk5GenJets.clone(rParam=0.8) kt6GenJets = kt4GenJets.clone(rParam=0.6) ca8GenJets = ca4GenJets.clone(rParam=0.8) sisCone5GenJetsNoNu = sisCone5GenJets.clone( src=cms.InputTag("genParticlesForJetsNoNu")) sisCone8GenJetsNoNu = sisCone8GenJets.clone( src=cms.InputTag("genParticlesForJetsNoNu")) kt4GenJetsNoNu = kt4GenJets.clone(src=cms.InputTag("genParticlesForJetsNoNu")) kt6GenJetsNoNu = kt6GenJets.clone(src=cms.InputTag("genParticlesForJetsNoNu")) iterativeCone5GenJetsNoNu = iterativeCone5GenJets.clone( src=cms.InputTag("genParticlesForJetsNoNu")) ak5GenJetsNoNu = ak5GenJets.clone(src=cms.InputTag("genParticlesForJetsNoNu")) ak4GenJetsNoNu = ak4GenJets.clone(src=cms.InputTag("genParticlesForJetsNoNu")) ak8GenJetsNoNu = ak8GenJets.clone(src=cms.InputTag("genParticlesForJetsNoNu")) gk5GenJetsNoNu = gk5GenJets.clone(src=cms.InputTag("genParticlesForJetsNoNu")) gk8GenJetsNoNu = gk8GenJets.clone(src=cms.InputTag("genParticlesForJetsNoNu"))
# ignoreParticleIDs = cms.vuint32( # 1000022, # 1000012, 1000014, 1000016, # 2000012, 2000014, 2000016, # 1000039, 5100039, # 4000012, 4000014, 4000016, # 9900012, 9900014, 9900016, # 39), # partonicFinalState = cms.bool(False), # excludeResonances = cms.bool(True), # excludeFromResonancePids = cms.vuint32(12, 13, 14, 16), # tausAsJets = cms.bool(False) #) from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ak5GenJetsptmin2 = ak5GenJets.clone() process.ak5GenJetsptmin2.jetPtMin = cms.double(2.0) process.new_ak5GenJets = cms.Sequence(process.genParticlesForJets * process.ak5GenJetsptmin2) #process.new_ak5GenJets = cms.Sequence(process.ak5GenJetsptmin2) # --- to recover the ak7 GenJets: process.ak7GenJetsptmin2 = ak5GenJets.clone() process.ak7GenJetsptmin2.rParam = cms.double(0.7) process.ak7GenJetsptmin2.jetPtMin = cms.double(2.0) process.new_ak7GenJets = cms.Sequence(process.genParticlesForJets * process.ak7GenJetsptmin2) #process.new_ak7GenJets = cms.Sequence(process.ak7GenJetsptmin2)
from HeavyIonsAnalysis.JetAnalysis.jets.HiReRecoJets_pp_cff import * ak1PFJetAnalyzer.doSubEvent = True ak1CaloJetAnalyzer.doSubEvent = True ak2PFJetAnalyzer.doSubEvent = True ak2CaloJetAnalyzer.doSubEvent = True ak3PFJetAnalyzer.doSubEvent = True ak3CaloJetAnalyzer.doSubEvent = True ak4PFJetAnalyzer.doSubEvent = True ak4CaloJetAnalyzer.doSubEvent = True ak5PFJetAnalyzer.doSubEvent = True ak5CaloJetAnalyzer.doSubEvent = True ak6PFJetAnalyzer.doSubEvent = True ak6CaloJetAnalyzer.doSubEvent = True from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets ak5GenJets = ak5GenJets ak1GenJets = ak5GenJets.clone(rParam=0.1) ak2GenJets = ak5GenJets.clone(rParam=0.2) ak3GenJets = ak5GenJets.clone(rParam=0.3) ak4GenJets = ak5GenJets.clone(rParam=0.4) ak6GenJets = ak5GenJets.clone(rParam=0.6) from RecoJets.Configuration.GenJetParticles_cff import * akGenJets = cms.Sequence(genParticlesForJets + ak1GenJets + ak2GenJets + ak3GenJets + ak4GenJets + ak5GenJets + ak6GenJets) from HeavyIonsAnalysis.JetAnalysis.makePartons_cff import * highPurityTracks = cms.EDFilter("TrackSelector", src=cms.InputTag("generalTracks"), cut=cms.string('quality("highPurity")')) jetSequences = cms.Sequence(akGenJets + ppReRecoPFJets + # ppReRecoCaloJets +
import FWCore.ParameterSet.Config as cms from RecoJets.JetProducers.sc5GenJets_cfi import sisCone5GenJets from RecoJets.JetProducers.ic5GenJets_cfi import iterativeCone5GenJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets from RecoJets.JetProducers.ak4GenJets_cfi import ak4GenJets from RecoJets.JetProducers.gk5GenJets_cfi import gk5GenJets from RecoJets.JetProducers.kt4GenJets_cfi import kt4GenJets from RecoJets.JetProducers.ca4GenJets_cfi import ca4GenJets from RecoHI.HiJetAlgos.HiGenJets_cff import * sisCone8GenJets = sisCone5GenJets.clone( rParam = 0.8 ) ak8GenJets = ak5GenJets.clone( rParam = 0.8, jetPtMin = 15.0 ) gk8GenJets = gk5GenJets.clone( rParam = 0.8 ) kt6GenJets = kt4GenJets.clone( rParam = 0.6 ) ca8GenJets = ca4GenJets.clone( rParam = 0.8 ) sisCone5GenJetsNoNu = sisCone5GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) sisCone8GenJetsNoNu = sisCone8GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) kt4GenJetsNoNu = kt4GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) kt6GenJetsNoNu = kt6GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) iterativeCone5GenJetsNoNu = iterativeCone5GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) ak5GenJetsNoNu = ak5GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) ak4GenJetsNoNu = ak4GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) ak8GenJetsNoNu = ak8GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) gk5GenJetsNoNu = gk5GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) gk8GenJetsNoNu = gk8GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) ca4GenJetsNoNu = ca4GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") )
import FWCore.ParameterSet.Config as cms # Remake AK5 jets from packedPFcandidates # See: https://twiki.cern.ch/twiki/bin/view/CMSPublic/WorkBookMiniAOD#Advanced_topics_re_clustering_ev ## Filter out neutrinos from packed GenParticles packedGenParticlesForJetsNoNu = cms.EDFilter("CandPtrSelector", src = cms.InputTag("packedGenParticles"), cut = cms.string("abs(pdgId) != 12 && abs(pdgId) != 14 && abs(pdgId) != 16")) ## Define GenJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets ak5GenJetsNoNu = ak5GenJets.clone(src = 'packedGenParticlesForJetsNoNu') ## Select charged hadron subtracted packed PF candidates pfCHS = cms.EDFilter("CandPtrSelector", src = cms.InputTag("packedPFCandidates"), cut = cms.string("fromPV")) from RecoJets.JetProducers.ak5PFJets_cfi import ak5PFJets ## Define ak5PFJetsCHS ak5PFJetsCHS = ak5PFJets.clone(src = 'pfCHS', doAreaFastjet = True) ## Define ak5PFJets ak5PFJets = ak5PFJets.clone(src = 'packedPFCandidates') ak5PFJetsSequence = cms.Sequence( packedGenParticlesForJetsNoNu* ak5GenJetsNoNu* pfCHS* ak5PFJetsCHS* ak5PFJets )
def reSetJet(process): ## Filter out neutrinos from packed GenParticles process.packedGenParticlesForJetsNoNu = cms.EDFilter("CandPtrSelector", src = cms.InputTag("packedGenParticles"), cut = cms.string("abs(pdgId) != 12 && abs(pdgId) != 14 && abs(pdgId) != 16")) ## Define GenJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ak5GenJetsNoNu = ak5GenJets.clone(src = 'packedGenParticlesForJetsNoNu') ## Select charged hadron subtracted packed PF candidates process.pfCHS = cms.EDFilter("CandPtrSelector", src = cms.InputTag("packedPFCandidates"), cut = cms.string("fromPV")) from RecoJets.JetProducers.ak5PFJets_cfi import ak5PFJets ## Define PFJetsCHS process.ak5PFJetsCHS = ak5PFJets.clone(src = 'pfCHS', doAreaFastjet = True) ################################################# ## Remake PAT jets ################################################# ## b-tag discriminators bTagDiscriminators = [ 'pfCombinedInclusiveSecondaryVertexV2BJetTags' ] ################################################################# slimmedAddPileupInfo = cms.EDProducer( 'PileupSummaryInfoSlimmer', src = cms.InputTag('addPileupInfo'), keepDetailedInfoFor = cms.vint32(0) ) ####################################################### from PhysicsTools.PatAlgos.tools.jetTools import addJetCollection ## Add PAT jet collection based on the above-defined ak5PFJetsCHS addJetCollection( process, labelName = 'AK5PFCHS', jetSource = cms.InputTag('ak5PFJetsCHS'), pvSource = cms.InputTag('offlineSlimmedPrimaryVertices'), bsSource = cms.InputTag('offlineBeamSpot'), pfCandidates = cms.InputTag('packedPFCandidates'), svSource = cms.InputTag('slimmedSecondaryVertices'), btagDiscriminators = bTagDiscriminators, jetCorrections = ('AK5PFchs', ['L1FastJet', 'L2Relative', 'L3Absolute'], 'None'), genJetCollection = cms.InputTag('ak5GenJetsNoNu'), genParticles = cms.InputTag('prunedGenParticles'), algo = 'AK', rParam = 0.5 ) getattr(process,'selectedPatJetsAK5PFCHS').cut = cms.string('pt > 20.') process.rejet = cms.Path(process.packedGenParticlesForJetsNoNu * process.ak5GenJetsNoNu * process.pfCHS * process.ak5PFJetsCHS * process.patJetCorrFactorsAK5PFCHS * process.patJetPartons * process.patJetFlavourAssociationAK5PFCHS * process.patJetPartonMatchAK5PFCHS * process.patJetGenJetMatchAK5PFCHS * process.pfImpactParameterTagInfosAK5PFCHS * process.pfInclusiveSecondaryVertexFinderTagInfosAK5PFCHS * process.pfCombinedInclusiveSecondaryVertexV2BJetTagsAK5PFCHS * process.patJetsAK5PFCHS * process.selectedPatJetsAK5PFCHS) from PhysicsTools.PatAlgos.tools.pfTools import adaptPVs ## Adapt primary vertex collection and BeamSpot adaptPVs(process, pvCollection=cms.InputTag('offlineSlimmedPrimaryVertices')) adaptBSs(process, bsCollection=cms.InputTag('offlineBeamSpot'))
from RecoJets.JetProducers.FastjetParameters_cfi import * from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets from RecoJets.JetProducers.ak5TrackJets_cfi import ak5TrackJets FastjetWithAreaPU = cms.PSet( Active_Area_Repeats=cms.int32(5), GhostArea=cms.double(0.01), Ghost_EtaMax=cms.double(6.0), UE_Subtraction=cms.string("no"), ) # MC jet ueAntiKt5GenJet = ak5GenJets.clone( src=cms.InputTag("goodParticles"), jetPtMin=cms.double(1.0), inputEtMin=cms.double(0.9) ) ueAntiKt5ChgGenJet = ak5GenJets.clone( src=cms.InputTag("chargeParticles"), jetPtMin=cms.double(1.0), inputEtMin=cms.double(0.9) ) ueAntiKt5GenJet500 = ak5GenJets.clone( src=cms.InputTag("goodParticles"), jetPtMin=cms.double(1.0), inputEtMin=cms.double(0.5) ) ueAntiKt5ChgGenJet500 = ak5GenJets.clone( src=cms.InputTag("chargeParticles"), jetPtMin=cms.double(1.0), inputEtMin=cms.double(0.5) ) ueAntiKt5GenJet1500 = ak5GenJets.clone(
from RecoJets.JetProducers.kt4PFJets_cfi import kt4PFJets from RecoHI.HiJetAlgos.hiFJRhoProducer import hiFJRhoProducer from RecoHI.HiJetAlgos.hiFJGridEmptyAreaCalculator_cff import hiFJGridEmptyAreaCalculator kt4PFJets.src = cms.InputTag('particleFlow') kt4PFJets.doAreaFastjet = True kt4PFJets.jetPtMin = cms.double(0.0) kt4PFJets.GhostArea = cms.double(0.005) kt2PFJets = kt4PFJets.clone(rParam=cms.double(0.2)) ak4PFJetAnalyzer.doSubEvent = True ak4CaloJetAnalyzer.doSubEvent = True from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets ak5GenJets = ak5GenJets ak4GenJets = ak5GenJets.clone(rParam=0.4, src=cms.InputTag("genParticlesForJetsNoNu")) from RecoJets.Configuration.GenJetParticles_cff import * ak3GenJets = ak5GenJets.clone(rParam=0.3, src=cms.InputTag("genParticlesForJetsNoNu")) from RecoJets.Configuration.GenJetParticles_cff import * akGenJets = cms.Sequence(genParticlesForJetsNoNu + ak4GenJets + ak3GenJets) from HeavyIonsAnalysis.JetAnalysis.makePartons_cff import * highPurityTracks = cms.EDFilter("TrackSelector", src=cms.InputTag("generalTracks"), cut=cms.string('quality("highPurity")')) #custom for 80X PFTowers.src = cms.InputTag("particleFlow")
from RecoJets.JetProducers.PFJetParameters_cfi import * from RecoJets.JetProducers.CaloJetParameters_cfi import * from RecoJets.JetProducers.AnomalousCellParameters_cfi import * from RecoJets.JetProducers.CATopJetParameters_cfi import * from RecoJets.JetProducers.GenJetParameters_cfi import * from RecoJets.JetProducers.caTopTaggers_cff import * ############################################### # GEN JETS from RecoJets.JetProducers.ca4GenJets_cfi import ca4GenJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ak3GenJets = ak5GenJets.clone( rParam = cms.double(0.3), src = cms.InputTag(genjetparticles)) process.ak15GenJets = ak5GenJets.clone( rParam = cms.double(1.5), src = cms.InputTag(genjetparticles)) process.ak15GenJetsFiltered = ak5GenJets.clone( rParam = cms.double(1.5), src = cms.InputTag(genjetparticles), useFiltering = cms.bool(True), nFilt = cms.int32(3), rFilt = cms.double(0.3), writeCompound = cms.bool(True), jetCollInstanceName=cms.string("SubJets") )
ak1PFJetAnalyzer.doSubEvent = True ak1CaloJetAnalyzer.doSubEvent = True ak2PFJetAnalyzer.doSubEvent = True ak2CaloJetAnalyzer.doSubEvent = True ak3PFJetAnalyzer.doSubEvent = True ak3CaloJetAnalyzer.doSubEvent = True ak4PFJetAnalyzer.doSubEvent = True ak4CaloJetAnalyzer.doSubEvent = True ak5PFJetAnalyzer.doSubEvent = True ak5CaloJetAnalyzer.doSubEvent = True ak6PFJetAnalyzer.doSubEvent = True ak6CaloJetAnalyzer.doSubEvent = True from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets ak5GenJets = ak5GenJets ak1GenJets = ak5GenJets.clone(rParam=0.1) ak2GenJets = ak5GenJets.clone(rParam=0.2) ak3GenJets = ak5GenJets.clone(rParam=0.3) ak4GenJets = ak5GenJets.clone(rParam=0.4) ak6GenJets = ak5GenJets.clone(rParam=0.6) from RecoJets.Configuration.GenJetParticles_cff import * akGenJets = cms.Sequence( genParticlesForJets + ak1GenJets + ak2GenJets + ak3GenJets + ak4GenJets + ak5GenJets + ak6GenJets ) from HeavyIonsAnalysis.JetAnalysis.makePartons_cff import * highPurityTracks = cms.EDFilter( "TrackSelector", src=cms.InputTag("generalTracks"), cut=cms.string('quality("highPurity")') )
###################### # load low-pt charged GenJet configuration here # get charged genParticles process.load('QCDAnalysis.UEAnalysis.UEAnalysisParticles_cfi') process.chargeParticles.cut = cms.string('charge != 0 & pt > 0.3 & status = 1') # get genjet definitions from RecoJets.JetProducers.ak4GenJets_cfi import ak4GenJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets # clone normal ak4 genjets process.ak4ChgGenJets = ak4GenJets.clone( src = cms.InputTag("chargeParticles"), jetPtMin = cms.double(0.3) ) # clone normal ak5 genjets process.ak5ChgGenJets = ak5GenJets.clone( src = cms.InputTag("chargeParticles"), jetPtMin = cms.double(0.3) ) # add the processes process.chargedgenjets = cms.Path(process.UEAnalysisParticles*process.ak4ChgGenJets*process.ak5ChgGenJets) ###################### # Here starts the CFF specific part import CommonFSQFramework.Core.customizePAT process = CommonFSQFramework.Core.customizePAT.customize(process) # GT customization process = CommonFSQFramework.Core.customizePAT.customizeGT(process) # define treeproducer process.GenLevelTree = cms.EDAnalyzer("CFFTreeProducer") import CommonFSQFramework.Core.GenLevelViewsConfigs
def reSetJet(process): ## Filter out neutrinos from packed GenParticles process.packedGenParticlesForJetsNoNu = cms.EDFilter("CandPtrSelector", src = cms.InputTag("packedGenParticles"), cut = cms.string("abs(pdgId) != 12 && abs(pdgId) != 14 && abs(pdgId) != 16")) ## Define GenJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ak5GenJetsNoNu = ak5GenJets.clone(src = 'packedGenParticlesForJetsNoNu') ## Select charged hadron subtracted packed PF candidates process.pfCHS = cms.EDFilter("CandPtrSelector", src = cms.InputTag("packedPFCandidates"), cut = cms.string("fromPV")) from RecoJets.JetProducers.ak5PFJets_cfi import ak5PFJets ## Define PFJetsCHS process.ak5PFJetsCHS = ak5PFJets.clone(src = 'pfCHS', doAreaFastjet = True) ################################################# ## Remake PAT jets ################################################# ## b-tag discriminators bTagDiscriminators = [ 'pfCombinedInclusiveSecondaryVertexV2BJetTags' ] ################################################################# slimmedAddPileupInfo = cms.EDProducer( 'PileupSummaryInfoSlimmer', #src = cms.InputTag('addPileupInfo'), keepDetailedInfoFor = cms.vint32(0) ) ####################################################### from PhysicsTools.PatAlgos.tools.jetTools import addJetCollection ## Add PAT jet collection based on the above-defined ak5PFJetsCHS addJetCollection( process, labelName = 'AK5PFCHS', jetSource = cms.InputTag('ak5PFJetsCHS'), pvSource = cms.InputTag('offlineSlimmedPrimaryVertices'), bsSource = cms.InputTag('offlineBeamSpot'), pfCandidates = cms.InputTag('packedPFCandidates'), svSource = cms.InputTag('slimmedSecondaryVertices'), btagDiscriminators = bTagDiscriminators, #jetCorrections = ('AK5PFchs', ['L1FastJet', 'L2Relative', 'L3Absolute'], 'None'), #commented before adding jet corrections genJetCollection = cms.InputTag('ak5GenJetsNoNu'), genParticles = cms.InputTag('prunedGenParticles'), algo = 'AK', rParam = 0.5 ) getattr(process,'selectedPatJetsAK5PFCHS').cut = cms.string('pt > 20.') process.rejet = cms.Path(process.packedGenParticlesForJetsNoNu * process.ak5GenJetsNoNu * process.pfCHS * process.ak5PFJetsCHS * process.patJetCorrFactorsAK5PFCHS * process.patJetPartons * process.patJetFlavourAssociationAK5PFCHS * process.patJetPartonMatchAK5PFCHS * process.patJetGenJetMatchAK5PFCHS * process.pfImpactParameterTagInfosAK5PFCHS * process.pfInclusiveSecondaryVertexFinderTagInfosAK5PFCHS * process.pfCombinedInclusiveSecondaryVertexV2BJetTagsAK5PFCHS * process.patJetsAK5PFCHS * process.selectedPatJetsAK5PFCHS) from PhysicsTools.PatAlgos.tools.pfTools import adaptPVs ## Adapt primary vertex collection and BeamSpot adaptPVs(process, pvCollection=cms.InputTag('offlineSlimmedPrimaryVertices')) adaptBSs(process, bsCollection=cms.InputTag('offlineBeamSpot')) from PhysicsTools.PatAlgos.tools.jetTools import updateJetCollection updateJetCollection( process, jetSource = cms.InputTag('slimmedJets'), labelName = 'UpdatedJEC', jetCorrections = ('AK5PFchs', cms.vstring(['L1FastJet', 'L2Relative', 'L3Absolute']), 'None') # Do not forget 'L2L3Residual' on data ) process.load('Configuration.StandardSequences.Services_cff') #Added info from now onwards process.load("JetMETCorrections.Modules.JetResolutionESProducer_cfi") from CondCore.DBCommon.CondDBSetup_cfi import * import os process.jer = cms.ESSource("PoolDBESSource", CondDBSetup, toGet = cms.VPSet( # Resolution cms.PSet( record = cms.string('JetResolutionRcd'), tag = cms.string('JR_Fall15_25nsV2_MC_PtResolution_AK5PFchs'), label = cms.untracked.string('AK5PFchs_pt') ), # Scale factors cms.PSet( record = cms.string('JetResolutionScaleFactorRcd'), tag = cms.string('JR_Fall15_25nsV2_MC_SF_AK4PFchs'), label = cms.untracked.string('AK5PFchs') ), ), connect = cms.string('sqlite:Fall15_25nsV2_MC.db') ) process.es_prefer_jer = cms.ESPrefer('PoolDBESSource', 'jer')
import FWCore.ParameterSet.Config as cms process = cms.Process("S2") process.source = cms.Source("PoolSource", fileNames = cms.untracked.vstring("file:patTuple_mini.root") ) process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(1000) ) from RecoJets.JetProducers.ak5PFJets_cfi import ak5PFJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.chs = cms.EDFilter("CandPtrSelector", src = cms.InputTag("packedPFCandidates"), cut = cms.string("fromPV")) process.ak5PFJets = ak5PFJets.clone(src = 'packedPFCandidates', doAreaFastjet = True) # no idea while doArea is false by default, but it's True in RECO so we have to set it process.ak5PFJetsCHS = ak5PFJets.clone(src = 'chs', doAreaFastjet = True) # no idea while doArea is false by default, but it's True in RECO so we have to set it process.ak5GenJets = ak5GenJets.clone(src = 'packedGenParticles') process.load("PhysicsTools.PatAlgos.producersLayer1.patCandidates_cff") process.load("Configuration.EventContent.EventContent_cff") process.load('Configuration.StandardSequences.Geometry_cff') process.load('Configuration.StandardSequences.MagneticField_38T_cff') process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff') process.GlobalTag.globaltag = 'START70_V6::All' from PhysicsTools.PatAlgos.tools.jetTools import addJetCollection addJetCollection( process, postfix = "", labelName = 'AK5PFCHS', jetSource = cms.InputTag('ak5PFJetsCHS'), trackSource = cms.InputTag('unpackedTracksAndVertices'),
* process.pfMetT0pcT1Txy ) process.JetSequence = cms.Sequence(process.ak5PFJetsCorr * process.pileupJetIdProducer * process.MyAk5PFJetTracksAssociatorAtVertex * process.MyImpactParameterPFTagInfos #* process.MySecondaryVertexTagInfos #* process.MyCombinedSecondaryVertexBJetTags ) if not(isRealData): from RecoJets.Configuration.GenJetParticles_cff import genParticlesForJetsNoNu process.genParticlesForJetsNoNu = genParticlesForJetsNoNu from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ak5GenJetsNoNu = ak5GenJets.clone( src = cms.InputTag("genParticlesForJetsNoNu") ) process.JetSequence += process.genParticlesForJetsNoNu process.JetSequence += process.ak5GenJetsNoNu # JEC uncertainty files #process.ak5PFJets.doAreaFastjet = True ### To get b-tags from ak5PFJets - this works with the 'Reconstruction_cff' above #process.path = cms.Path(process.btagging) #djs - for full reconstruction #this part is from nate #djs commented out as a test - the above works and should be used when done testing
import FWCore.ParameterSet.Config as cms from RecoJets.Configuration.GenJetParticles_cff import genParticlesForJetsNoMuNoNu from RecoJets.Configuration.GenJetParticles_cff import genParticlesForJetsNoNu from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets cms2antikt5GenJets = ak5GenJets.clone( src=cms.InputTag("genParticlesForJetsNoMuNoNu")) cms2antikt5PFGenJets = ak5GenJets.clone( src=cms.InputTag("genParticlesForJetsNoNu")) genJetSequence = cms.Sequence(genParticlesForJetsNoMuNoNu * cms2antikt5GenJets * genParticlesForJetsNoNu * cms2antikt5PFGenJets)
), src=cms.InputTag(pvSrc) ) ############################### ########## Gen Setup ########## ############################### process.load("RecoJets.Configuration.GenJetParticles_cff") from RecoJets.JetProducers.ca4GenJets_cfi import ca4GenJets from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.ca8GenJetsNoNu = ca4GenJets.clone( rParam = cms.double(0.8), src = cms.InputTag("genParticlesForJetsNoNu")) process.ak8GenJetsNoNu = ak5GenJets.clone( rParam = cms.double(0.8), src = cms.InputTag("genParticlesForJetsNoNu")) process.load("TopQuarkAnalysis.TopEventProducers.sequences.ttGenEvent_cff") # add the flavor history process.load("PhysicsTools.HepMCCandAlgos.flavorHistoryPaths_cfi") # prune gen particles process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi") process.prunedGenParticles = cms.EDProducer("GenParticlePruner", src = cms.InputTag("genParticles"), select = cms.vstring( "drop *" ,"keep status = 3" #keeps particles from the hard matrix element
process.ZMMProducer = cms.EDProducer('ZMMFilter', PFCandidateInputTag = cms.InputTag("particleFlow"), VertexTag = cms.InputTag("goodOfflinePrimaryVertices"), TrackTag = cms.InputTag("generalTracks"), MuonTag = cms.InputTag("muons"), ) from RecoParticleFlow.PFProducer.particleFlowTmpPtrs_cfi import particleFlowTmpPtrs process.particleFlowNoMuonTmpPtrs = particleFlowTmpPtrs.clone( src = cms.InputTag('ZMMProducer', 'PFCandidatesNoMuon') ) ## GenJet from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.load('RecoJets.Configuration.GenJetParticles_cff') process.ak4GenJets = ak5GenJets.clone(rParam = 0.4, jetPtMin= 0.0, doAreaFastjet=True, src = 'genParticlesForJetsNoMuNoNu') ## AK4 PFJet from RecoJets.JetProducers.ak5PFJets_cfi import ak5PFJets process.ak4PFJets = ak5PFJets.clone(rParam = 0.4, jetPtMin= 0.0, doAreaFastjet=True, src = cms.InputTag("ZMMProducer", "PFCandidatesNoMuon") ) ## AK4 CHSJet from CommonTools.ParticleFlow.pfNoPileUpJME_cff import * process.load('CommonTools.ParticleFlow.pfNoPileUpJME_cff') process.PFPileUpZMM = pfPileUp.clone( PFCandidates='particleFlowNoMuonTmpPtrs', Vertices = 'goodOfflinePrimaryVertices') process.PFNoPileUpZMM = pfNoPileUp.clone(topCollection = 'PFPileUpZMM', bottomCollection = 'particleFlowNoMuonTmpPtrs' ) process.ak4PFJetsCHS = ak5PFJets.clone(rParam = 0.4, doAreaFastjet=True, jetPtMin = 0.0, src = cms.InputTag("PFNoPileUpZMM")) ## AK4 PuppiJet from Dummy.Puppi.Puppi_PhaseII_SLHC23p1_v2 import *
# process.load("Configuration.StandardSequences.Reconstruction_cff") ############################################################################### ############################################################################### if not isData: ############################ # Add ReReco of ak5GenJets # ############################ process.load('RecoJets.Configuration.GenJetParticles_cff') process.load('RecoJets.Configuration.RecoGenJets_cff') from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets process.lowPtak5GenJets = ak5GenJets.clone(jetPtMin = 0.5 ) process.LowPtGenJetsReCluster = cms.Path(process.lowPtak5GenJets) ############################################################################### ############################################################################### # ############################################################################### # ############################################################################### # ########################################################################### # # Add PAT calo Jets # !!! Still with old JEC (Jet Energy Corrections) !!! # # ########################################################################### # ## switch to uncheduled mode # process.options.allowUnscheduled = cms.untracked.bool(True) # #process.Tracer = cms.Service("Tracer")