Esempio n. 1
0
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)
Esempio n. 2
0
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(
Esempio n. 3
0
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
Esempio n. 5
0
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
Esempio n. 6
0
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)
Esempio n. 7
0
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)