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.caloJetsForTrk_cff import *
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', 'BPix1+BPix3', 'BPix2+BPix3',
'BPix1+FPix1_pos', 'BPix1+FPix1_neg',
'BPix2+FPix1_pos', 'BPix2+FPix1_neg',
'FPix1_pos+FPix2_pos', 'FPix1_neg+FPix2_neg',
#'BPix2+TIB1','BPix2+TIB2',
'BPix3+TIB1','BPix3+TIB2'),
TIB = cms.PSet(
matchedRecHits = cms.InputTag("siStripMatchedRecHits","matchedRecHit"),
TTRHBuilder = cms.string('WithTrackAngle'), clusterChargeCut = cms.PSet(refToPSet_ = cms.string('SiStripClusterChargeCutNone'))
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
hiInitialStepTrackRefsForJets = trackRefsForJets.clone(
src=cms.InputTag('hiGlobalPrimTracks'))
#change this to import Bkg substracted Heavy Ion jets:
from RecoHI.HiJetAlgos.hiCaloJetsForTrk_cff import *
hiJetsForCoreTracking = cms.EDFilter(
"CandPtrSelector",
src=cms.InputTag("akPu4CaloJetsSelected"),
cut=cms.string("pt > 30 && abs(eta) < 2.5"))
# care only at tracks from main PV
hiFirstStepGoodPrimaryVertices = cms.EDFilter(
"PrimaryVertexObjectFilter",
filterParams=cms.PSet(minNdof=cms.double(25.0),
maxZ=cms.double(15.0),
maxRho=cms.double(2.0)),
src=cms.InputTag('hiSelectedPixelVertex'))
# SEEDING LAYERS
hiJetCoreRegionalStepSeedLayers = cms.EDProducer(
"SeedingLayersEDProducer",
layerList=cms.vstring(
'BPix1+BPix2+BPix3',
'BPix1+BPix2+FPix1_pos',
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
hiInitialStepTrackRefsForJets = trackRefsForJets.clone(
src='hiGlobalPrimTracks')
#change this to import Bkg substracted Heavy Ion jets:
from RecoHI.HiJetAlgos.hiCaloJetsForTrk_cff import *
hiJetsForCoreTracking = cms.EDFilter(
"CandPtrSelector",
src=cms.InputTag("akPu4CaloJetsSelected"),
cut=cms.string("pt > 30 && abs(eta) < 2.5"))
# care only at tracks from main PV
hiFirstStepGoodPrimaryVertices = cms.EDFilter(
"PrimaryVertexObjectFilter",
filterParams=cms.PSet(minNdof=cms.double(25.0),
maxZ=cms.double(15.0),
maxRho=cms.double(2.0)),
src=cms.InputTag('hiSelectedPixelVertex'))
# SEEDING LAYERS
hiJetCoreRegionalStepSeedLayers = cms.EDProducer(
"SeedingLayersEDProducer",
layerList=cms.vstring(
'BPix1+BPix2+BPix3',
'BPix1+BPix2+FPix1_pos',
TrackLabel = 'initialStepTracks',
vertexCollections = [_offlinePrimaryVertices.vertexCollections[0].clone()]
)
(pp_on_XeXe_2017 | pp_on_AA).toModify(firstStepPrimaryVerticesUnsorted, TkFilterParameters = dict(trackQuality = 'any'))
# we need a replacment for the firstStepPrimaryVerticesUnsorted
# that includes tracker information of signal and pile up
# after mixing there is no such thing as initialStepTracks,
# so we replace the input collection for firstStepPrimaryVerticesUnsorted with generalTracks
firstStepPrimaryVerticesBeforeMixing = firstStepPrimaryVerticesUnsorted.clone()
fastSim.toModify(firstStepPrimaryVerticesUnsorted, TrackLabel = 'generalTracks')
from RecoJets.JetProducers.TracksForJets_cff import trackRefsForJets
initialStepTrackRefsForJets = trackRefsForJets.clone(
src = 'initialStepTracks'
)
fastSim.toModify(initialStepTrackRefsForJets, src = 'generalTracks')
from RecoJets.JetProducers.caloJetsForTrk_cff import *
from CommonTools.RecoAlgos.sortedPrimaryVertices_cfi import sortedPrimaryVertices as _sortedPrimaryVertices
firstStepPrimaryVertices = _sortedPrimaryVertices.clone(
vertices = 'firstStepPrimaryVerticesUnsorted',
particles = 'initialStepTrackRefsForJets',
)
# Final selection
from RecoTracker.FinalTrackSelectors.TrackMVAClassifierPrompt_cfi import *
from RecoTracker.FinalTrackSelectors.TrackMVAClassifierDetached_cfi import *
initialStepClassifier1 = TrackMVAClassifierPrompt.clone(
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
hiInitialStepTrackRefsForJets = trackRefsForJets.clone(src = cms.InputTag('hiGlobalPrimTracks'))
#change this to import Bkg substracted Heavy Ion jets:
from RecoHI.HiJetAlgos.hiCaloJetsForTrk_cff import *
hiJetsForCoreTracking = cms.EDFilter("CandPtrSelector", src = cms.InputTag("akPu4CaloJetsSelected"), cut = cms.string("pt > 30 && abs(eta) < 2.5"))
# care only at tracks from main PV
hiFirstStepGoodPrimaryVertices = cms.EDFilter("PrimaryVertexObjectFilter",
filterParams = cms.PSet(
minNdof = cms.double(25.0),
maxZ = cms.double(15.0),
maxRho = cms.double(2.0)
),
src=cms.InputTag('hiSelectedVertex')
)
# SEEDING LAYERS
hiJetCoreRegionalStepSeedLayers = cms.EDProducer("SeedingLayersEDProducer",
layerList = cms.vstring('BPix1+BPix2+BPix3',
'BPix1+BPix2+FPix1_pos',
'BPix1+BPix2+FPix1_neg',
'BPix1+FPix1_pos+FPix2_pos',
'BPix1+FPix1_neg+FPix2_neg',
'BPix1+BPix2+TIB1',
