) sisCone5PFJets = cms.EDProducer( "FastjetJetProducer", PFJetParameters, AnomalousCellParameters, jetAlgorithm = cms.string("SISCone"), rParam = cms.double(0.4) ) #MC jet ueSisCone5GenJet = sisCone5GenJets.clone( src = cms.InputTag("goodParticles"), jetPtMin = cms.double(1.0), inputEtMin = cms.double(0.9) ) ueSisCone5ChgGenJet = sisCone5GenJets.clone( src = cms.InputTag("chargeParticles"), jetPtMin = cms.double(1.0), inputEtMin = cms.double(0.9) ) ueSisCone5GenJet500 = sisCone5GenJets.clone( src = cms.InputTag("goodParticles"), jetPtMin = cms.double(1.0), inputEtMin = cms.double(0.5) )
# 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)
##### from Diego # 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"),
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.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"))
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 jet ueSisCone5ChgGenJet500 = sisCone5GenJets.clone( src=cms.InputTag("chargeParticles"), jetPtMin=cms.double(1.0), inputEtMin=cms.double(0.5) ) # RECO jet Tracks ueSisCone5TracksJet500 = sisCone5TrackJets.clone( src=cms.InputTag("goodTracks"), jetPtMin=cms.double(1.0), inputEtMin=cms.double(0.5) ) UEAnalysisNeededJets = cms.Sequence(ueSisCone5ChgGenJet500 * ueSisCone5TracksJet500)
process.load("Configuration.StandardSequences.MagneticField_cff") ##### from Diego # 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()