示例#1
0
import FWCore.ParameterSet.Config as cms

from RecoJets.JetProducers.ak4CaloJets_cfi import ak4CaloJets

from RecoJets.JetProducers.fixedGridRhoProducerFastjet_cfi import fixedGridRhoFastjetAllCalo

doPileup = cms.bool(True)

ak4CaloJetsPUCorr = ak4CaloJets.clone(doPUOffsetCorr=doPileup)
ak4CaloJets.doAreaFastjet = True

fixedGridRhoFastjetCentralCalo = fixedGridRhoFastjetAllCalo.clone(
    maxRapidity=cms.double(2.5))

recoJets = cms.Sequence(fixedGridRhoFastjetAllCalo +
                        fixedGridRhoFastjetCentralCalo + ak4CaloJets)

recoAllJets = cms.Sequence(fixedGridRhoFastjetAllCalo +
                           fixedGridRhoFastjetCentralCalo +
                           fixedGridRhoFastjetAllCalo + ak4CaloJets)

recoAllJetsPUOffsetCorr = cms.Sequence(fixedGridRhoFastjetAllCalo +
                                       fixedGridRhoFastjetCentralCalo +
                                       ak4CaloJetsPUCorr)
示例#2
0
import FWCore.ParameterSet.Config as cms

from RecoJets.JetProducers.ak4CaloJets_cfi import ak4CaloJets as _ak4CaloJets
from RecoLocalCalo.CaloTowersCreator.calotowermaker_cfi import calotowermaker
caloTowerForTrk = calotowermaker.clone(hbheInput=cms.InputTag('hbheprereco'))
ak4CaloJetsForTrk = _ak4CaloJets.clone(
    srcPVs=cms.InputTag('firstStepPrimaryVerticesUnsorted'),
    src=cms.InputTag('caloTowerForTrk'))
from Configuration.Eras.Modifier_trackingLowPU_cff import trackingLowPU
trackingLowPU.toModify(ak4CaloJetsForTrk, srcPVs="pixelVertices")
from Configuration.Eras.Modifier_trackingPhase2PU140_cff import trackingPhase2PU140
trackingPhase2PU140.toModify(ak4CaloJetsForTrk, srcPVs="pixelVertices")

caloJetsForTrk = cms.Sequence(caloTowerForTrk * ak4CaloJetsForTrk)

from Configuration.Eras.Modifier_phase2_hcal_cff import phase2_hcal
phase2_hcal.toModify(caloTowerForTrk, hbheInput=cms.InputTag("hbhereco"))
示例#3
0
import FWCore.ParameterSet.Config as cms

# This step runs over all clusters

# run only if there are high pT jets
from RecoJets.JetProducers.TracksForJets_cff import trackRefsForJets

initialStepTrackRefsForJets = trackRefsForJets.clone(src=cms.InputTag("initialStepTracks"))
from RecoJets.JetProducers.ak4CaloJets_cfi import ak4CaloJets
from RecoLocalCalo.CaloTowersCreator.calotowermaker_cfi import calotowermaker

caloTowerForTrk = calotowermaker.clone(hbheInput=cms.InputTag("hbheprereco"))
ak4CaloJetsForTrk = ak4CaloJets.clone(
    srcPVs=cms.InputTag("firstStepPrimaryVertices"), src=cms.InputTag("caloTowerForTrk")
)
jetsForCoreTracking = cms.EDFilter(
    "CandPtrSelector", src=cms.InputTag("ak4CaloJetsForTrk"), cut=cms.string("pt > 100 && abs(eta) < 2.5")
)

# care only at tracks from main PV
firstStepGoodPrimaryVertices = cms.EDFilter(
    "PrimaryVertexObjectFilter",
    filterParams=cms.PSet(minNdof=cms.double(25.0), maxZ=cms.double(15.0), maxRho=cms.double(2.0)),
    src=cms.InputTag("firstStepPrimaryVertices"),
)

# SEEDING LAYERS
jetCoreRegionalStepSeedLayers = cms.EDProducer(
    "SeedingLayersEDProducer",
    layerList=cms.vstring(
        "BPix1+BPix2",
示例#4
0
import FWCore.ParameterSet.Config as cms

from RecoJets.JetProducers.ak4CaloJets_cfi import ak4CaloJets
from RecoLocalCalo.CaloTowersCreator.calotowermaker_cfi import calotowermaker
caloTowerForTrk = calotowermaker.clone(hbheInput=cms.InputTag('hbheprereco'))
ak4CaloJetsForTrk = ak4CaloJets.clone(srcPVs = cms.InputTag('firstStepPrimaryVertices'), src= cms.InputTag('caloTowerForTrk'))
from Configuration.Eras.Modifier_trackingLowPU_cff import trackingLowPU
trackingLowPU.toModify(ak4CaloJetsForTrk,
    srcPVs = "pixelVertices"
)
from Configuration.Eras.Modifier_trackingPhase2PU140_cff import trackingPhase2PU140
trackingPhase2PU140.toModify(ak4CaloJetsForTrk,
    srcPVs = "pixelVertices"
)

caloJetsForTrk = cms.Sequence(caloTowerForTrk*ak4CaloJetsForTrk)

from Configuration.Eras.Modifier_phase2_hcal_cff import phase2_hcal
phase2_hcal.toModify( caloTowerForTrk, hbheInput = cms.InputTag("hbhereco") )

示例#5
0
# Restrict SISCone algorithm to 1000 towers input
sisCone5CaloJets.restrictInputs = cms.bool(True)
sisCone5CaloJets.maxInputs = cms.uint32(1000)

sisCone7CaloJets.restrictInputs = cms.bool(True)
sisCone7CaloJets.maxInputs = cms.uint32(1000)

doPileup = cms.bool(True)

sisCone5CaloJetsPUCorr      =sisCone5CaloJets.clone      (doPUOffsetCorr = doPileup)
sisCone7CaloJetsPUCorr      =sisCone7CaloJets.clone      (doPUOffsetCorr = doPileup)
kt4CaloJetsPUCorr           =kt4CaloJets.clone           (doPUOffsetCorr = doPileup)
kt6CaloJetsPUCorr           =kt6CaloJets.clone           (doPUOffsetCorr = doPileup)
iterativeCone5CaloJetsPUCorr=iterativeCone5CaloJets.clone(doPUOffsetCorr = doPileup)
ak5CaloJetsPUCorr           =ak5CaloJets.clone           (doPUOffsetCorr = doPileup)
ak4CaloJetsPUCorr           =ak4CaloJets.clone           (doPUOffsetCorr = doPileup)
ak7CaloJetsPUCorr           =ak7CaloJets.clone           (doPUOffsetCorr = doPileup)
gk5CaloJetsPUCorr           =gk5CaloJets.clone           (doPUOffsetCorr = doPileup)
gk7CaloJetsPUCorr           =gk7CaloJets.clone           (doPUOffsetCorr = doPileup)
ca4CaloJetsPUCorr           =ca4CaloJets.clone           (doPUOffsetCorr = doPileup)
ca6CaloJetsPUCorr           =ca6CaloJets.clone           (doPUOffsetCorr = doPileup)

#compute areas for Fastjet PU subtraction  
kt6CaloJets.doRhoFastjet = True
kt6CaloJets.doAreaFastjet = True
#use active areas and not Voronoi tessellation for the moment
#kt6CaloJets.voronoiRfact = 0.9
ak5CaloJets.doAreaFastjet = True
ak4CaloJets.doAreaFastjet = True
ak7CaloJets.doAreaFastjet = True
示例#6
0
import FWCore.ParameterSet.Config as cms

from RecoJets.JetProducers.ak4CaloJets_cfi import ak4CaloJets as _ak4CaloJets
from RecoHI.HiJetAlgos.HiRecoJets_cff import akPu4CaloJets as _akPu4CaloJets
from RecoLocalCalo.CaloTowersCreator.calotowermaker_cfi import calotowermaker
caloTowerForTrk = calotowermaker.clone(hbheInput='hbheprereco',
                                       missingHcalRescaleFactorForEcal=1.0)

ak4CaloJetsForTrk = _ak4CaloJets.clone(
    srcPVs='firstStepPrimaryVerticesUnsorted', src='caloTowerForTrk')

from Configuration.ProcessModifiers.pp_on_AA_cff import pp_on_AA
pp_on_AA.toReplaceWith(
    ak4CaloJetsForTrk,
    _akPu4CaloJets.clone(srcPVs='firstStepPrimaryVerticesUnsorted',
                         src='caloTowerForTrk'))

from Configuration.Eras.Modifier_trackingLowPU_cff import trackingLowPU
trackingLowPU.toModify(ak4CaloJetsForTrk, srcPVs="pixelVertices")

caloJetsForTrkTask = cms.Task(caloTowerForTrk, ak4CaloJetsForTrk)
caloJetsForTrk = cms.Sequence(caloJetsForTrkTask)

from Configuration.Eras.Modifier_run3_HB_cff import run3_HB
run3_HB.toModify(caloTowerForTrk, hbheInput="hbhereco")

from Configuration.Eras.Modifier_pf_badHcalMitigationOff_cff import pf_badHcalMitigationOff
pf_badHcalMitigationOff.toModify(caloTowerForTrk,
                                 missingHcalRescaleFactorForEcal=0.0)
示例#7
0
import FWCore.ParameterSet.Config as cms

# This step runs over all clusters

# run only if there are high pT jets
from RecoJets.JetProducers.TracksForJets_cff import trackRefsForJets

iter0TrackRefsForJets = trackRefsForJets.clone(
    src=cms.InputTag('initialStepTracks'))
from RecoJets.JetProducers.ak4CaloJets_cfi import ak4CaloJets
from RecoLocalCalo.CaloTowersCreator.calotowermaker_cfi import calotowermaker

caloTowerForTrk = calotowermaker.clone(hbheInput=cms.InputTag('hbheprereco'))
ak4CaloJetsForTrk = ak4CaloJets.clone(srcPVs=cms.InputTag('pixelVertices'),
                                      src=cms.InputTag('caloTowerForTrk'))
jetsForCoreTracking = cms.EDFilter(
    "CandPtrSelector",
    src=cms.InputTag("ak4CaloJetsForTrk"),
    cut=cms.string("pt > 100 && abs(eta) < 2.5"))

# care only at tracks from main PV
import RecoVertex.PrimaryVertexProducer.OfflinePrimaryVertices_cfi

firstStepPrimaryVertices = RecoVertex.PrimaryVertexProducer.OfflinePrimaryVertices_cfi.offlinePrimaryVertices.clone(
)
firstStepPrimaryVertices.TrackLabel = cms.InputTag("initialStepTracks")
firstStepGoodPrimaryVertices = cms.EDFilter(
    "PrimaryVertexObjectFilter",
    filterParams=cms.PSet(minNdof=cms.double(25.0),
                          maxZ=cms.double(15.0),
                          maxRho=cms.double(2.0)),