)


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)
示例#3
0
##### 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"),
示例#4
0
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") )
示例#5
0
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()