import FWCore.ParameterSet.Config as cms from RecoJets.JetProducers.sc5TrackJets_cfi import sisCone5TrackJets from RecoJets.JetProducers.ic5TrackJets_cfi import iterativeCone5TrackJets from RecoJets.JetProducers.ak5TrackJets_cfi import ak5TrackJets from RecoJets.JetProducers.ak4TrackJets_cfi import ak4TrackJets from RecoJets.JetProducers.gk5TrackJets_cfi import gk5TrackJets from RecoJets.JetProducers.kt4TrackJets_cfi import kt4TrackJets from RecoJets.JetProducers.ca4TrackJets_cfi import ca4TrackJets from CommonTools.RecoAlgos.TrackWithVertexRefSelector_cfi import * from RecoJets.JetProducers.TracksForJets_cff import * sisCone7TrackJets = sisCone5TrackJets.clone(rParam=0.7) ak7TrackJets = ak5TrackJets.clone(rParam=0.7) gk7TrackJets = gk5TrackJets.clone(rParam=0.7) kt6TrackJets = kt4TrackJets.clone(rParam=0.6) ca6TrackJets = ca4TrackJets.clone(rParam=0.6) recoTrackJets = cms.Sequence(trackWithVertexRefSelector + trackRefsForJets + ak4TrackJets) recoAllTrackJets = cms.Sequence(trackWithVertexRefSelector + trackRefsForJets + sisCone5TrackJets + sisCone7TrackJets + kt4TrackJets + kt6TrackJets + iterativeCone5TrackJets + ak4TrackJets + ak5TrackJets + ak7TrackJets + gk5TrackJets + gk7TrackJets + ca4TrackJets + ca6TrackJets)
ueAntiKt5ChgGenJet700 = ak5GenJets.clone( src=cms.InputTag("chargeParticles"), jetPtMin=cms.double(1.0), inputEtMin=cms.double(0.7) ) ueAntiKt5GenJet1100 = ak5GenJets.clone( src=cms.InputTag("goodParticles"), jetPtMin=cms.double(1.0), inputEtMin=cms.double(1.1) ) ueAntiKt5ChgGenJet1100 = ak5GenJets.clone( src=cms.InputTag("chargeParticles"), jetPtMin=cms.double(1.0), inputEtMin=cms.double(1.1) ) # RECO jet Tracks ueAntiKt5TracksJet500 = ak5TrackJets.clone( src=cms.InputTag("goodTracks"), jetPtMin=cms.double(1.0), inputEtMin=cms.double(0.5) ) ueAntiKt5TracksJet = ak5TrackJets.clone( src=cms.InputTag("goodTracks"), jetPtMin=cms.double(1.0), inputEtMin=cms.double(0.9) ) ueAntiKt5TracksJet1500 = ak5TrackJets.clone( src=cms.InputTag("goodTracks"), jetPtMin=cms.double(1.0), inputEtMin=cms.double(1.5) ) ueAntiKt5TracksJet700 = ak5TrackJets.clone( src=cms.InputTag("goodTracks"), jetPtMin=cms.double(1.0), inputEtMin=cms.double(0.7) ) ueAntiKt5TracksJet1100 = ak5TrackJets.clone(
import FWCore.ParameterSet.Config as cms from RecoJets.JetProducers.sc5TrackJets_cfi import sisCone5TrackJets from RecoJets.JetProducers.ic5TrackJets_cfi import iterativeCone5TrackJets from RecoJets.JetProducers.ak5TrackJets_cfi import ak5TrackJets from RecoJets.JetProducers.ak4TrackJets_cfi import ak4TrackJets from RecoJets.JetProducers.gk5TrackJets_cfi import gk5TrackJets from RecoJets.JetProducers.kt4TrackJets_cfi import kt4TrackJets from RecoJets.JetProducers.ca4TrackJets_cfi import ca4TrackJets from CommonTools.RecoAlgos.TrackWithVertexRefSelector_cfi import * from RecoJets.JetProducers.TracksForJets_cff import * sisCone7TrackJets = sisCone5TrackJets.clone( rParam = 0.7 ) ak7TrackJets = ak5TrackJets.clone( rParam = 0.7 ) gk7TrackJets = gk5TrackJets.clone( rParam = 0.7 ) kt6TrackJets = kt4TrackJets.clone( rParam = 0.6 ) ca6TrackJets = ca4TrackJets.clone( rParam = 0.6 ) recoTrackJets =cms.Sequence(trackWithVertexRefSelector+ trackRefsForJets+ ak4TrackJets+ak5TrackJets+kt4TrackJets) recoAllTrackJets=cms.Sequence(trackWithVertexRefSelector+ trackRefsForJets+ sisCone5TrackJets+sisCone7TrackJets+ kt4TrackJets+kt6TrackJets+ iterativeCone5TrackJets+ ak4TrackJets+ak5TrackJets+ak7TrackJets+ gk5TrackJets+gk7TrackJets+
process.TFileService = cms.Service("TFileService", fileName = cms.string("Phase1_140PU_aged_AK4_SLHC23p1.root"), closeFileFast = cms.untracked.bool(True) ) ## 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 = 'genParticlesForJetsNoNu') process.load('CommonTools.RecoAlgos.TrackWithVertexRefSelector_cfi') process.load('RecoJets.JetProducers.TracksForJets_cff') from RecoJets.JetProducers.ak5TrackJets_cfi import ak5TrackJets process.ak4TrackJets = ak5TrackJets.clone(rParam = 0.4, jetPtMin= 0.0) ## PFJet from RecoJets.Configuration.RecoPFJets_cff import * process.ak4PFJets = ak5PFJets.clone(rParam = 0.4, jetPtMin= 0.0, doAreaFastjet=True) process.neuPar = cms.EDFilter("CandPtrSelector", src = cms.InputTag("particleFlow"), cut = cms.string("charge == 0")) process.ak4NeuJets = ak5PFJets.clone(rParam = 0.4, jetPtMin= 0.0, doAreaFastjet=True, src = cms.InputTag('neuPar')) process.ChaPar = cms.EDFilter("CandPtrSelector", src = cms.InputTag("particleFlow"), cut = cms.string("charge != 0")) process.ak4ChaJets = ak5PFJets.clone(rParam = 0.4, jetPtMin= 0.0, doAreaFastjet=True, src = cms.InputTag('ChaPar')) ## CHSJet process.load('CommonTools.ParticleFlow.pfNoPileUpJME_cff') process.ak4PFJetsCHS = ak5PFJetsCHS.clone(rParam = 0.4, jetPtMin= 0.0, doAreaFastjet=True) ## CHS no good vertex
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.sc5TrackJets_cfi import sisCone5TrackJets from RecoJets.JetProducers.ak5TrackJets_cfi import ak5TrackJets process.sisCone5TrackJets = sisCone5TrackJets.clone( rParam = 0.5, jetPtMin=1.0, UseOnlyVertexTracks=True, UseOnlyOnePV=True, src = cms.InputTag("trackRefsForJets"), inputEtMin = cms.double(0.5)) process.sisCone7TrackJets = sisCone5TrackJets.clone(rParam = 0.7, jetPtMin=1.0, UseOnlyVertexTracks=True, UseOnlyOnePV=True, src = cms.InputTag("trackRefsForJets"), inputEtMin = cms.double(0.5)) process.ak5TrackJets = ak5TrackJets.clone( rParam = 0.5, jetPtMin=1.0, UseOnlyVertexTracks=True, UseOnlyOnePV=True, src = cms.InputTag("trackRefsForJets"), inputEtMin = cms.double(0.5)) process.ak4TrackJets = ak5TrackJets.clone( rParam = 0.4, jetPtMin=1.0, UseOnlyVertexTracks=True, UseOnlyOnePV=True, src = cms.InputTag("trackRefsForJets"), inputEtMin = cms.double(0.5)) process.ak7TrackJets = ak5TrackJets.clone( rParam = 0.7, jetPtMin=1.0, UseOnlyVertexTracks=True, UseOnlyOnePV=True, src = cms.InputTag("trackRefsForJets"), inputEtMin = cms.double(0.5)) process.ak10TrackJets = ak5TrackJets.clone( rParam = 1.0, jetPtMin=1.0, UseOnlyVertexTracks=True, UseOnlyOnePV=True, src = cms.InputTag("trackRefsForJets"), inputEtMin = cms.double(0.5)) process.chargedjets = cms.Path(process.sisCone5TrackJets* process.sisCone7TrackJets*process.ak4TrackJets *process.ak5TrackJets*process.ak7TrackJets*process.ak10TrackJets) ##### stop from Diego # Here starts the CFF specific part
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)
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 charged gen jet ak5ChgGenJets500 = ak5GenJets.clone( src = cms.InputTag("chargeParticles"), jetPtMin = cms.double(0.3), inputEtMin = cms.double(0.3) ) # DATA and MC reco track jet ak5TrackJets500 = ak5TrackJets.clone( src = cms.InputTag("goodTracks"), jetPtMin = cms.double(0.3), inputEtMin = cms.double(0.3) ) UEAnalysisNeededJets = cms.Sequence(ak5ChgGenJets500*ak5TrackJets500)
# Make the jets import FWCore.ParameterSet.Config as cms 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 charged gen jet ak5ChgGenJets500 = ak5GenJets.clone(src=cms.InputTag("chargeParticles"), jetPtMin=cms.double(0.3), inputEtMin=cms.double(0.3)) # DATA and MC reco track jet ak5TrackJets500 = ak5TrackJets.clone(src=cms.InputTag("goodTracks"), jetPtMin=cms.double(0.3), inputEtMin=cms.double(0.3)) UEAnalysisNeededJets = cms.Sequence(ak5ChgGenJets500 * ak5TrackJets500)