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)
# Make the jets
import FWCore.ParameterSet.Config as cms
from RecoJets.JetProducers.FastjetParameters_cfi import *
from RecoJets.JetProducers.sc5GenJets_cfi import sisCone5GenJets
from RecoJets.JetProducers.sc5TrackJets_cfi import sisCone5TrackJets
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
SisCone5ChgGenJet500 = sisCone5GenJets.clone(
src=cms.InputTag("chargeParticles"),
jetPtMin=cms.double(1.0),
inputEtMin=cms.double(0.5))
# DATA and MC reco track jet
SisCone5TracksJet500 = sisCone5TrackJets.clone(src=cms.InputTag("goodTracks"),
jetPtMin=cms.double(1.0),
inputEtMin=cms.double(0.5))
UEAnalysisNeededJets = cms.Sequence(SisCone5ChgGenJet500 *
SisCone5TracksJet500)
src = cms.InputTag("goodParticles"),
jetPtMin = cms.double(1.0),
inputEtMin = cms.double(1.1)
)
ueSisCone5ChgGenJet1100 = sisCone5GenJets.clone(
src = cms.InputTag("chargeParticles"),
jetPtMin = cms.double(1.0),
inputEtMin = cms.double(1.1)
)
#RECO jet Tracks
ueSisCone5TracksJet500 = sisCone5TrackJets.clone(
src = cms.InputTag("goodTracks"),
jetPtMin = cms.double(1.0),
inputEtMin = cms.double(0.5)
)
ueSisCone5TracksJet = sisCone5TrackJets.clone(
src = cms.InputTag("goodTracks"),
jetPtMin = cms.double(1.0),
inputEtMin = cms.double(0.9)
)
ueSisCone5TracksJet1500 = sisCone5TrackJets.clone(
src = cms.InputTag("goodTracks"),
jetPtMin = cms.double(1.0),
inputEtMin = cms.double(1.5)
)
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.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.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
