def miniAOD_customizeAllMCFastSim(process):
miniAOD_customizeCommon(process)
miniAOD_customizeMC(process)
from PhysicsTools.PatAlgos.slimming.metFilterPaths_cff import miniAOD_customizeMETFiltersFastSim
process = miniAOD_customizeMETFiltersFastSim(process)
from PhysicsTools.PatAlgos.slimming.isolatedTracks_cfi import miniAOD_customizeIsolatedTracksFastSim
process = miniAOD_customizeIsolatedTracksFastSim(process)
process.patMuons.addTriggerMatching = False
# Disable pixelClusterTagInfos in FastSim (no siPixelCluster available)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(process.patJets, addTagInfos = cms.bool(False) )
fastSim.toModify(process.slimmedJetsNoDeepFlavour, dropTagInfos = cms.string('1') )
fastSim.toModify(process.updatedPatJetsSlimmedDeepFlavour, addTagInfos = cms.bool(False) )
fastSim.toModify(process.updatedPatJetsTransientCorrectedSlimmedDeepFlavour, addTagInfos = cms.bool(False) )
return process
KFSmootherForRefitOutsideIn = KFTrajectorySmoother.clone(
ComponentName='KFSmootherForRefitOutsideIn',
Propagator='SmartPropagatorAnyRKOpposite',
Updator='KFUpdator',
Estimator='Chi2EstimatorForRefit',
errorRescaling=100.0,
minHits=3)
#
KFFitterForRefitInsideOut = KFTrajectoryFitter.clone(
ComponentName='KFFitterForRefitInsideOut',
Propagator='SmartPropagatorAnyRK',
Updator='KFUpdator',
Estimator='Chi2EstimatorForRefit',
minHits=3)
KFSmootherForRefitInsideOut = KFTrajectorySmoother.clone(
ComponentName='KFSmootherForRefitInsideOut',
Propagator='SmartPropagatorAnyRK',
Updator='KFUpdator',
Estimator='Chi2EstimatorForRefit',
errorRescaling=100.0,
minHits=3)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
# FastSim doesn't use Runge Kute for propagation
# the following propagators are not used in FastSim, but just to be sure...
fastSim.toModify(KFFitterForRefitOutsideIn, Propagator='SmartPropagatorAny')
fastSim.toModify(KFSmootherForRefitOutsideIn, Propagator='SmartPropagator')
fastSim.toModify(KFFitterForRefitInsideOut, Propagator="SmartPropagatorAny")
fastSim.toModify(KFSmootherForRefitInsideOut, Propagator="SmartPropagatorAny")
from RecoTracker.FinalTrackSelectors.trackTfClassifier_cfi import *
from RecoTracker.FinalTrackSelectors.trackSelectionTf_cfi import *
from RecoTracker.FinalTrackSelectors.trackSelectionTf_CKF_cfi import *
trackdnn.toReplaceWith(
jetCoreRegionalStep,
trackTfClassifier.clone(
src='jetCoreRegionalStepTracks',
qualityCuts=qualityCutDictionary.JetCoreRegionalStep.value()))
trackdnn.toReplaceWith(
jetCoreRegionalStepBarrel,
trackTfClassifier.clone(
src='jetCoreRegionalStepBarrelTracks',
qualityCuts=qualityCutDictionary.JetCoreRegionalStep.value()))
fastSim.toModify(jetCoreRegionalStep,
vertices='firstStepPrimaryVerticesBeforeMixing')
jetCoreRegionalStepEndcap = jetCoreRegionalStep.clone(
src='jetCoreRegionalStepEndcapTracks', )
# Final sequence
JetCoreRegionalStepTask = cms.Task(
jetsForCoreTracking,
firstStepGoodPrimaryVertices,
#jetCoreRegionalStepClusters,
jetCoreRegionalStepSeedLayers,
jetCoreRegionalStepTrackingRegions,
jetCoreRegionalStepHitDoublets,
jetCoreRegionalStepSeeds,
jetCoreRegionalStepTrackCandidates,
jetCoreRegionalStepTracks,
import FWCore.ParameterSet.Config as cms # Modifier for FastSim from Configuration.Eras.Modifier_fastSim_cff import fastSim # PFRecTracks from generalTracks from RecoParticleFlow.PFTracking.pfTrack_cfi import * lowPtGsfElePfTracks = pfTrack.clone() lowPtGsfElePfTracks.TkColList = ['generalTracks'] lowPtGsfElePfTracks.GsfTracksInEvents = False lowPtGsfElePfTracks.GsfTrackModuleLabel = '' fastSim.toModify(lowPtGsfElePfTracks, TkColList=['generalTracksBeforeMixing']) # Low pT ElectronSeeds from RecoEgamma.EgammaElectronProducers.lowPtGsfElectronSeeds_cfi import * # Electron track candidates from TrackingTools.GsfTracking.CkfElectronCandidateMaker_cff import * lowPtGsfEleTrajectoryFilter = TrajectoryFilterForElectrons.clone() lowPtGsfEleTrajectoryFilter.minPt = 0. lowPtGsfEleTrajectoryFilter.minimumNumberOfHits = 3 lowPtGsfEleTrajectoryBuilder = TrajectoryBuilderForElectrons.clone() lowPtGsfEleTrajectoryBuilder.trajectoryFilter.refToPSet_ = 'lowPtGsfEleTrajectoryFilter' lowPtGsfEleCkfTrackCandidates = electronCkfTrackCandidates.clone() lowPtGsfEleCkfTrackCandidates.TrajectoryBuilderPSet.refToPSet_ = 'lowPtGsfEleTrajectoryBuilder' lowPtGsfEleCkfTrackCandidates.src = 'lowPtGsfElectronSeeds' import FastSimulation.Tracking.electronCkfTrackCandidates_cff
)
)
simSiStripDigis = cms.EDAlias(
mix = cms.VPSet(
cms.PSet(type = cms.string('SiStripDigiedmDetSetVector')),
cms.PSet(type = cms.string('SiStripRawDigiedmDetSetVector')),
cms.PSet(type = cms.string('StripDigiSimLinkedmDetSetVector'))
)
)
#mergedtruth = cms.EDAlias(
# mix = cms.VPSet(
# cms.PSet(type = cms.string('TrackingParticles')),
# cms.PSet(type = cms.string('TrackingVertexs'))
# )
#)
genPUProtons = cms.EDAlias(
mixData = cms.VPSet(
cms.PSet( type = cms.string('recoGenParticles') )
)
)
# no castor,pixel,strip digis in fastsim
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(simCastorDigis, mix = None)
fastSim.toModify(simSiPixelDigis, mix = None)
fastSim.toModify(simSiStripDigis, mix = None)
from Configuration.Eras.Modifier_run3_common_cff import run3_common
run3_common.toModify(simCastorDigis, mix = None)
ptMin = 0,
)
# Select B-hadron TPs
trackingParticlesBHadron = _trackingParticleBHadronRefSelector.clone()
## MTV instances
trackValidator = Validation.RecoTrack.MultiTrackValidator_cfi.multiTrackValidator.clone(
useLogPt = cms.untracked.bool(True),
dodEdxPlots = True,
doPVAssociationPlots = True
#,minpT = cms.double(-1)
#,maxpT = cms.double(3)
#,nintpT = cms.int32(40)
)
fastSim.toModify(trackValidator,
dodEdxPlots = False)
for _eraName, _postfix, _era in _relevantEras:
_setForEra(trackValidator, _eraName, _era,
label = ["generalTracks", locals()["_generalTracksHp"+_postfix]] +
locals()["_selectorsByAlgo"+_postfix] + locals()["_selectorsByAlgoHp"+_postfix] +
locals()["_selectorsByOriginalAlgo"+_postfix] +
["generalTracksPt09"] + locals()["_selectorsPt09"+_postfix] +
[
"cutsRecoTracksBtvLike",
"cutsRecoTracksAK4PFJets"
],
doResolutionPlotsForLabels = [
"generalTracks",
locals()["_generalTracksHp"+_postfix],
"generalTracksPt09",
cms.InputTag('g4SimHits', 'MuonDTHits'),
cms.InputTag('g4SimHits', 'MuonCSCHits'),
cms.InputTag('g4SimHits', 'MuonRPCHits'),
cms.InputTag('g4SimHits', 'TrackerHitsTIBLowTof'),
cms.InputTag('g4SimHits', 'TrackerHitsTIBHighTof'),
cms.InputTag('g4SimHits', 'TrackerHitsTIDLowTof'),
cms.InputTag('g4SimHits', 'TrackerHitsTIDHighTof'),
cms.InputTag('g4SimHits', 'TrackerHitsTOBLowTof'),
cms.InputTag('g4SimHits', 'TrackerHitsTOBHighTof'),
cms.InputTag('g4SimHits', 'TrackerHitsTECLowTof'),
cms.InputTag('g4SimHits', 'TrackerHitsTECHighTof'),
cms.InputTag('g4SimHits', 'TrackerHitsPixelBarrelLowTof'),
cms.InputTag('g4SimHits', 'TrackerHitsPixelBarrelHighTof'),
cms.InputTag('g4SimHits', 'TrackerHitsPixelEndcapLowTof'),
cms.InputTag('g4SimHits', 'TrackerHitsPixelEndcapHighTof')),
trackingParticleSrc=cms.InputTag('mix', 'MergedTrackTruth'))
from Configuration.Eras.Modifier_phase2_tracker_cff import phase2_tracker
phase2_tracker.toModify(simHitTPAssocProducer,
simHitSrc=[
"g4SimHits:TrackerHitsPixelBarrelLowTof",
"g4SimHits:TrackerHitsPixelEndcapLowTof"
])
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(simHitTPAssocProducer,
simHitSrc=[
"famosSimHits:TrackerHits", "MuonSimHits:MuonCSCHits",
"MuonSimHits:MuonDTHits", "MuonSimHits:MuonRPCHits"
])
import FWCore.ParameterSet.Config as cms
from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
HcalSimHitsAnalyser = DQMEDAnalyzer('HcalSimHitsValidation',
outputFile = cms.untracked.string('')
)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify( HcalSimHitsAnalyser, ModuleLabel = cms.untracked.string("fastSimProducer") )
from Configuration.Eras.Modifier_run2_HCAL_2017_cff import run2_HCAL_2017
run2_HCAL_2017.toModify( HcalSimHitsAnalyser, TestNumber = cms.untracked.bool(True), EEHitCollection = cms.untracked.string("") )
DQMOffline_DoubleMu4_Mass8_PFHT350 = DQMOffline_DoubleMu4_Mass8_DZ_PFHT350.clone()
DQMOffline_DoubleMu4_Mass8_PFHT350.folderName = cms.string('HLT_DoubleMu4_Mass8_PFHT350')
DQMOffline_DoubleMu4_Mass8_PFHT350.numGenericTriggerEventPSet.hltPaths = cms.vstring("HLT_DoubleMu4_Mass8_PFHT350_v*")
DQMOffline_DoubleEle8_CaloIdM_TrackIdM_Mass8_PFHT350 = DQMOffline_DoubleEle8_CaloIdM_TrackIdM_Mass8_DZ_PFHT350.clone()
DQMOffline_DoubleEle8_CaloIdM_TrackIdM_Mass8_PFHT350.folderName = cms.string('HLT_DoubleEle8_CaloIdM_TrackIdM_Mass8_PFHT350')
DQMOffline_DoubleEle8_CaloIdM_TrackIdM_Mass8_PFHT350.numGenericTriggerEventPSet.hltPaths = cms.vstring("HLT_DoubleEle8_CaloIdM_TrackIdM_Mass8_PFHT350_v*")
DQMOffline_Mu8_Ele8_CaloIdM_TrackIdM_Mass8_PFHT350 = DQMOffline_Mu8_Ele8_CaloIdM_TrackIdM_Mass8_PFHT350_DZ.clone()
DQMOffline_Mu8_Ele8_CaloIdM_TrackIdM_Mass8_PFHT350.folderName = cms.string('HLT_Mu8_Ele8_CaloIdM_TrackIdM_Mass8_PFHT350')
DQMOffline_Mu8_Ele8_CaloIdM_TrackIdM_Mass8_PFHT350.numGenericTriggerEventPSet.hltPaths = cms.vstring("HLT_Mu8_Ele8_CaloIdM_TrackIdM_Mass8_PFHT350_v*")
# fastsim has no conversion collection
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(DQMOffline_Ele15_HT600,conversionCollection=cms.InputTag(''))
fastSim.toModify(DQMOffline_Ele15_HT450,conversionCollection=cms.InputTag(''))
fastSim.toModify(DQMOffline_Ele50_HT450,conversionCollection=cms.InputTag(''))
fastSim.toModify(DQMOffline_DoubleEle8_CaloIdM_TrackIdM_Mass8_DZ_PFHT350,conversionCollection=cms.InputTag(''))
fastSim.toModify(DQMOffline_Mu8_Ele8_CaloIdM_TrackIdM_Mass8_PFHT350_DZ,conversionCollection=cms.InputTag(''))
DQMOffline_LepHT_POSTPROCESSING = DQMEDHarvester("DQMGenericClient",
subDirs = cms.untracked.vstring('HLT/SUSY/LepHT/*'),
efficiency = cms.vstring(
"lepPtTurnOn_eff ';Offline lepton p_{T} [GeV];#epsilon' lepPtTurnOn_num lepPtTurnOn_den",
"lepEtaTurnOn_eff ';Offline lepton #eta;#epsilon' lepEtaTurnOn_num lepEtaTurnOn_den",
"lepPhiTurnOn_eff ';Offline lepton #phi;#epsilon' lepPhiTurnOn_num lepPhiTurnOn_den",
"lepEtaPhiTurnOn_eff ';Offline lepton #eta;Offline lepton #phi;#epsilon' lepEtaPhiTurnOn_num lepEtaPhiTurnOn_den",
"pfHTTurnOn_eff ';Offline PF H_{T} [GeV];#epsilon' pfHTTurnOn_num pfHTTurnOn_den",
import FWCore.ParameterSet.Config as cms
muondtdigianalyzer = cms.EDAnalyzer("MuonDTDigis",
# Label to retrieve Digis from the event
DigiLabel = cms.InputTag('simMuonDTDigis'),
# Label to retrieve SimHits from the event
SimHitLabel = cms.InputTag('g4SimHits', "MuonDTHits"),
# Name of the root file which will contain the histos
outputFile = cms.untracked.string(''),
# Switch on/off the verbosity
verbose = cms.untracked.bool(False)
)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(muondtdigianalyzer, SimHitLabel = "MuonSimHits:MuonDTHits")
import FWCore.ParameterSet.Config as cms
cscDigiValidation = cms.EDAnalyzer("CSCDigiValidation",
simHitsTag = cms.InputTag("mix", "g4SimHitsMuonCSCHits"),
wireDigiTag = cms.InputTag("simMuonCSCDigis","MuonCSCWireDigi"),
outputFile = cms.string(''),
stripDigiTag = cms.InputTag("simMuonCSCDigis","MuonCSCStripDigi"),
comparatorDigiTag = cms.InputTag("simMuonCSCDigis","MuonCSCComparatorDigi"),
alctDigiTag = cms.InputTag("simCscTriggerPrimitiveDigis"),
clctDigiTag = cms.InputTag("simCscTriggerPrimitiveDigis"),
doSim = cms.bool(False)
)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(cscDigiValidation, simHitsTag = "mix:MuonSimHitsMuonCSCHits")
import FWCore.ParameterSet.Config as cms
from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
hcalSimHitStudy = DQMEDAnalyzer('HcalSimHitStudy',
ModuleLabel = cms.untracked.string('g4SimHits'),
outputFile = cms.untracked.string(''),
Verbose = cms.untracked.bool(False),
HitCollection = cms.untracked.string('HcalHits'),
TestNumber = cms.bool(False),
hep17 = cms.bool(False)
)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify( hcalSimHitStudy, ModuleLabel = cms.untracked.string('fastSimProducer') )
from Configuration.Eras.Modifier_run2_HCAL_2017_cff import run2_HCAL_2017
(run2_HCAL_2017 & ~fastSim).toModify( hcalSimHitStudy, TestNumber = cms.bool(True) )
from Configuration.Eras.Modifier_run2_HEPlan1_2017_cff import run2_HEPlan1_2017
run2_HEPlan1_2017.toModify( hcalSimHitStudy, hep17 = cms.bool(True) )
),
mcPartons = cms.InputTag("hltBtagJetsbyValAlgo")
)
#put all in a path
hltbtagValidationSequence = cms.Sequence(
# remove noisy warnings
# hltBtagTriggerSelection +
hltBtagJetMCTools +
HltVertexValidationVertices +
hltbTagValidation
)
# fastsim customs
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(HltVertexValidationVertices, SimVertexCollection = "famosSimHits")
# are these customs actually needed?
#HltVertexValidationVertices.HLTPathNames =cms.vstring(
#'HLT_PFMET120_NoiseCleaned_BTagCSV07_v',
#'HLT_PFMET120_NoiseCleaned_BTagCSV07_v',
# 'HLT_CaloMHTNoPU90_PFMET90_PFMHT90_IDLoose_',
# 'HLT_CaloMHTNoPU90_PFMET90_PFMHT90_IDLoose_',
# 'HLT_QuadPFJet_VBF',
# 'HLT_QuadPFJet_VBF',
# 'HLT_Ele32_eta2p1_',
# 'HLT_IsoMu24_eta2p1_')
#HltVertexValidationVertices.Vertex = cms.VInputTag(
# cms.InputTag("hltVerticesL3"),
# cms.InputTag("hltFastPVPixelVertices"),
# cms.InputTag("hltVerticesPF"),
#)
import FastSimulation.Tracking.TrackCandidateProducer_cfi
_fastSim_detachedQuadStepTrackCandidates = FastSimulation.Tracking.TrackCandidateProducer_cfi.trackCandidateProducer.clone(
src = cms.InputTag("detachedQuadStepSeeds"),
MinNumberOfCrossedLayers = 4,
hitMasks = cms.InputTag("detachedQuadStepMasks")
)
fastSim.toReplaceWith(detachedQuadStepTrackCandidates,_fastSim_detachedQuadStepTrackCandidates)
# TRACK FITTING
import RecoTracker.TrackProducer.TrackProducer_cfi
detachedQuadStepTracks = RecoTracker.TrackProducer.TrackProducer_cfi.TrackProducer.clone(
AlgorithmName = 'detachedQuadStep',
src = 'detachedQuadStepTrackCandidates',
Fitter = 'FlexibleKFFittingSmoother',
)
fastSim.toModify(detachedQuadStepTracks,TTRHBuilder = 'WithoutRefit')
# TRACK SELECTION AND QUALITY FLAG SETTING.
from RecoTracker.FinalTrackSelectors.TrackMVAClassifierDetached_cfi import *
detachedQuadStep = TrackMVAClassifierDetached.clone(
src = 'detachedQuadStepTracks',
mva = dict(GBRForestLabel = 'MVASelectorDetachedQuadStep_Phase1'),
qualityCuts = [-0.5,0.0,0.5],
)
fastSim.toModify(detachedQuadStep,vertices = "firstStepPrimaryVerticesBeforeMixing")
highBetaStar_2018.toModify(detachedQuadStep,qualityCuts = [-0.7,0.0,0.5])
# For Phase2PU140
import RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi
detachedQuadStepSelector = RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.multiTrackSelector.clone(
src = 'detachedQuadStepTracks',
from SimGeneral.TrackingAnalysis.trackingParticleNumberOfLayersProducer_cfi import * from Configuration.Eras.Modifier_fastSim_cff import fastSim fastSim.toModify(trackingParticleNumberOfLayersProducer, simHits=['fastSimProducer:TrackerHits']) from Configuration.Eras.Modifier_phase2_tracker_cff import phase2_tracker phase2_tracker.toModify(trackingParticleNumberOfLayersProducer, simHits = ["g4SimHits:TrackerHitsPixelBarrelLowTof", "g4SimHits:TrackerHitsPixelEndcapLowTof"]) from Configuration.ProcessModifiers.premix_stage2_cff import premix_stage2 premix_stage2.toModify(trackingParticleNumberOfLayersProducer, trackingParticles = "mixData:MergedTrackTruth")
muoncosmicreco = cms.Sequence(muoncosmicrecoTask)
muoncosmichighlevelrecoTask = cms.Task(muoncosmicreco2legsHighLevelTask,
muoncosmicreco1legHighLevelTask,
cosmicsMuonIdTask)
muoncosmichighlevelreco = cms.Sequence(muoncosmichighlevelrecoTask)
#### High level sequence (i.e., post PF reconstruction) ###
from RecoMuon.MuonIdentification.muons_cfi import *
from RecoMuon.MuonIdentification.displacedMuons_cfi import *
from RecoMuon.MuonIsolation.muonPFIsolation_cff import *
from RecoMuon.MuonIdentification.muonReducedTrackExtras_cfi import *
from RecoMuon.MuonIdentification.displacedMuonReducedTrackExtras_cfi import *
# clusters are not present in fastsim
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(muonReducedTrackExtras, outputClusters=False)
# cluster collections are different in phase 2
from Configuration.Eras.Modifier_phase2_tracker_cff import phase2_tracker
phase2_tracker.toModify(muonReducedTrackExtras, outputClusters=False)
phase2_tracker.toModify(displacedMuonReducedTrackExtras, outputClusters=False)
muonshighlevelrecoTask = cms.Task(muonPFIsolationTask,
displacedMuonPFIsolationTask, muons,
displacedMuons, muonReducedTrackExtras,
displacedMuonReducedTrackExtras)
muonshighlevelreco = cms.Sequence(muonshighlevelrecoTask)
# displaced sequences do not run in fastsim
fastSim.toReplaceWith(
muonshighlevelrecoTask,
# CMS NOTE 2003-17: (0.244) m/ns
vPropWire=cms.double(24.4),
# The Synchronization Algo Name. Algos Type= [DTDigiSyncTOFCorr,DTDigiSyncFromTable]
SyncName=cms.string('DTDigiSyncTOFCorr'),
# configure the creation of Digi-Sim links:
# one-to-one correspondence between digis and SimHits (MultipleLinks=false)
# or association of SimHits within a given time window (of the order of the resolution, in ns)
MultipleLinks=cms.bool(True),
# constant drift velocity for the IdealModel (cm/ns)
IdealModelConstantDriftVelocity=cms.double(0.0055),
interpolate=cms.bool(True),
# further configurable smearing
Smearing=cms.double(2.4),
# Dead time for signals on the same wire (number from M. Pegoraro)
deadTime=cms.double(150.0),
#Name of Collection used for create the XF
mixLabel=cms.string('mix'),
InputCollection=cms.string('g4SimHitsMuonDTHits'),
debug=cms.untracked.bool(False),
# Its parameters
pset=cms.PSet(TOFCorrection=cms.int32(3), offset=cms.double(500.0)),
# ideal model. Used for debug.
# It uses constant drift velocity and it hasn't any external delay
IdealModel=cms.bool(False),
LinksTimeWindow=cms.double(10.0),
onlyMuHits=cms.bool(False),
GeometryType=cms.string('idealForDigi'))
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(simMuonDTDigis, InputCollection='MuonSimHitsMuonDTHits')
photonPostprocessing.batch = cms.bool(False)
photonPostprocessing.standalone = cms.bool(False)
photonPostprocessing.isRunCentrally = cms.bool(True)
#
pfPhotonPostprocessing = Validation.RecoEgamma.photonPostprocessing_cfi.photonPostprocessing.clone()
pfPhotonPostprocessing.ComponentName = cms.string('pfPhotonPostprocessing')
pfPhotonPostprocessing.analyzerName = cms.string('pfPhotonValidator')
pfPhotonPostprocessing.batch = cms.bool(False)
pfPhotonPostprocessing.standalone = cms.bool(False)
pfPhotonPostprocessing.isRunCentrally = cms.bool(True)
#
oldpfPhotonPostprocessing = Validation.RecoEgamma.photonPostprocessing_cfi.photonPostprocessing.clone()
oldpfPhotonPostprocessing.ComponentName = cms.string('oldpfPhotonPostprocessing')
oldpfPhotonPostprocessing.analyzerName = cms.string('oldpfPhotonValidator')
oldpfPhotonPostprocessing.batch = cms.bool(False)
oldpfPhotonPostprocessing.standalone = cms.bool(False)
oldpfPhotonPostprocessing.isRunCentrally = cms.bool(True)
#
from Validation.RecoEgamma.conversionPostprocessing_cfi import *
conversionPostprocessing.batch = cms.bool(False)
conversionPostprocessing.standalone = cms.bool(False)
photonPostProcessor = cms.Sequence(photonPostprocessing*pfPhotonPostprocessing*conversionPostprocessing)
#photonPostProcessor = cms.Sequence(photonPostprocessing*conversionPostprocessing)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(photonPostprocessing, fastSim = True)
fastSim.toModify(oldpfPhotonPostprocessing, fastSim = True)
import FWCore.ParameterSet.Config as cms
from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
HcalSimHitsAnalyser = DQMEDAnalyzer('HcalSimHitsValidation',
outputFile = cms.untracked.string('')
)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify( HcalSimHitsAnalyser, ModuleLabel = cms.untracked.string("famosSimHits") )
from Configuration.Eras.Modifier_run2_HCAL_2017_cff import run2_HCAL_2017
run2_HCAL_2017.toModify( HcalSimHitsAnalyser, TestNumber = cms.untracked.bool(True), EEHitCollection = cms.untracked.string("") )
pdgId = cms.vint32(),
signalOnly = cms.bool(True),
intimeOnly = cms.bool(False),
stableOnly = cms.bool(False),
minRapidity = cms.double(-2.5),
minHit = cms.int32(0),
ptMin = cms.double(0.9),
ptMax = cms.double(1e100),
maxRapidity = cms.double(2.5),
tip = cms.double(3.5),
minPhi = cms.double(-3.2),
maxPhi = cms.double(3.2),
)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(generalTpSelectorBlock, stableOnly = True)
TpSelectorForEfficiencyVsEtaBlock = generalTpSelectorBlock.clone()
TpSelectorForEfficiencyVsPhiBlock = generalTpSelectorBlock.clone()
TpSelectorForEfficiencyVsPtBlock = generalTpSelectorBlock.clone(ptMin = 0.050 )
TpSelectorForEfficiencyVsVTXRBlock = generalTpSelectorBlock.clone(tip = 60.0)
TpSelectorForEfficiencyVsVTXZBlock = generalTpSelectorBlock.clone()
def _modifyForPhase1(pset):
pset.minRapidity = -3
pset.maxRapidity = 3
from Configuration.Eras.Modifier_phase1Pixel_cff import phase1Pixel
phase1Pixel.toModify(generalTpSelectorBlock, _modifyForPhase1)
phase1Pixel.toModify(TpSelectorForEfficiencyVsEtaBlock, _modifyForPhase1)
from Configuration.Eras.Modifier_phase2_tracker_cff import phase2_tracker
phase2_tracker.toModify(generalTpSelectorBlock, minRapidity=-4.5, maxRapidity=4.5)
tobTecFlexibleKFFittingSmoother = TrackingTools.TrackFitters.FlexibleKFFittingSmoother_cfi.FlexibleKFFittingSmoother.clone(
ComponentName = cms.string('tobTecFlexibleKFFittingSmoother'),
standardFitter = cms.string('tobTecStepFitterSmoother'),
looperFitter = cms.string('tobTecStepFitterSmootherForLoopers'),
)
# TRACK FITTING
import RecoTracker.TrackProducer.TrackProducer_cfi
tobTecStepTracks = RecoTracker.TrackProducer.TrackProducer_cfi.TrackProducer.clone(
src = 'tobTecStepTrackCandidates',
AlgorithmName = cms.string('tobTecStep'),
#Fitter = 'tobTecStepFitterSmoother',
Fitter = 'tobTecFlexibleKFFittingSmoother',
)
fastSim.toModify(tobTecStepTracks, TTRHBuilder = 'WithoutRefit')
# TRACK SELECTION AND QUALITY FLAG SETTING.
from RecoTracker.FinalTrackSelectors.TrackMVAClassifierPrompt_cfi import *
from RecoTracker.FinalTrackSelectors.TrackMVAClassifierDetached_cfi import *
tobTecStepClassifier1 = TrackMVAClassifierDetached.clone()
tobTecStepClassifier1.src = 'tobTecStepTracks'
tobTecStepClassifier1.mva.GBRForestLabel = 'MVASelectorIter6_13TeV'
tobTecStepClassifier1.qualityCuts = [-0.6,-0.45,-0.3]
fastSim.toModify(tobTecStepClassifier1, vertices = "firstStepPrimaryVerticesBeforeMixing")
tobTecStepClassifier2 = TrackMVAClassifierPrompt.clone()
tobTecStepClassifier2.src = 'tobTecStepTracks'
tobTecStepClassifier2.mva.GBRForestLabel = 'MVASelectorIter0_13TeV'
tobTecStepClassifier2.qualityCuts = [0.0,0.0,0.0]
collectionDM = cms.string('MuonCSCWireDigiSimLinks'),
),
),
)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
from FastSimulation.Tracking.recoTrackAccumulator_cfi import recoTrackAccumulator as _recoTrackAccumulator
fastSim.toModify(mixData,
# from signal: mix tracks not strip or pixel digis
workers = dict(
pixel = None,
strip = None,
pixelSimLink = None,
stripSimLink = None,
tracks = cms.PSet(
workerType = cms.string("PreMixingDigiAccumulatorWorker"),
accumulator = _recoTrackAccumulator.clone(
pileUpTracks = "mix:generalTracks"
)
)
),
)
from Configuration.Eras.Modifier_phase2_common_cff import phase2_common
from Configuration.Eras.Modifier_phase2_tracker_cff import phase2_tracker
from Configuration.Eras.Modifier_phase2_ecal_cff import phase2_ecal
from Configuration.Eras.Modifier_phase2_hcal_cff import phase2_hcal
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
from Configuration.Eras.Modifier_phase2_muon_cff import phase2_muon
phase2_common.toModify(mixData, input = dict(producers = [])) # we use digis directly, no need for raw2digi producers
from Configuration.ProcessModifiers.premix_stage1_cff import premix_stage1
(premix_stage1 & phase2_hgcal).toModify(simHGCalUnsuppressedDigis,
mix={
0: dict(type="PHGCSimAccumulator"),
1: dict(type="PHGCSimAccumulator"),
2: dict(type="PHGCSimAccumulator"),
})
from Configuration.Eras.Modifier_phase2_hfnose_cff import phase2_hfnose
(~phase2_hfnose).toModify(simHFNoseUnsuppressedDigis, mix=None)
from Configuration.Eras.Modifier_phase1Pixel_cff import phase1Pixel
phase1Pixel.toModify(
simSiPixelDigis,
mix=_pixelCommon +
[cms.PSet(type=cms.string('PixelFEDChanneledmNewDetSetVector'))])
from Configuration.Eras.Modifier_phase2_tracker_cff import phase2_tracker
phase2_tracker.toModify(simSiStripDigis, mix=None)
phase2_tracker.toModify(simAPVsaturation, mix=None)
# no castor,pixel,strip digis in fastsim
from Configuration.Eras.Modifier_fastSim_cff import fastSim
(fastSim & ~run3_common).toModify(simCastorDigis, mix=None)
fastSim.toModify(simSiPixelDigis, mix=None)
fastSim.toModify(simSiStripDigis, mix=None)
fastSim.toModify(simAPVsaturation, mix=None)
castor = cms.PSet(
castorDigitizer
),
puVtx = cms.PSet(
pileupVtxDigitizer
),
mergedtruth = cms.PSet(
trackingParticles
)
)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(theDigitizers,
# fastsim does not model castor
castor = None,
# fastsim does not digitize pixel and strip hits
pixel = None,
strip = None,
tracks = recoTrackAccumulator
)
from Configuration.ProcessModifiers.premix_stage2_cff import premix_stage2
(fastSim & premix_stage2).toModify(theDigitizers,
tracks = None
)
from SimCalorimetry.HGCalSimProducers.hgcalDigitizer_cfi import hgceeDigitizer, hgchebackDigitizer, hgchefrontDigitizer, HGCAL_noise_fC, HGCAL_noise_MIP, HGCAL_chargeCollectionEfficiencies, HGCAL_noises
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toModify( theDigitizers,
hgceeDigitizer = cms.PSet(hgceeDigitizer),
hgchebackDigitizer = cms.PSet(hgchebackDigitizer),
import FastSimulation.Tracking.TrackCandidateProducer_cfi
_fastSim_detachedQuadStepTrackCandidates = FastSimulation.Tracking.TrackCandidateProducer_cfi.trackCandidateProducer.clone(
src=cms.InputTag("detachedQuadStepSeeds"),
MinNumberOfCrossedLayers=4,
hitMasks=cms.InputTag("detachedQuadStepMasks"))
fastSim.toReplaceWith(detachedQuadStepTrackCandidates,
_fastSim_detachedQuadStepTrackCandidates)
# TRACK FITTING
import RecoTracker.TrackProducer.TrackProducer_cfi
detachedQuadStepTracks = RecoTracker.TrackProducer.TrackProducer_cfi.TrackProducer.clone(
AlgorithmName='detachedQuadStep',
src='detachedQuadStepTrackCandidates',
Fitter='FlexibleKFFittingSmoother',
)
fastSim.toModify(detachedQuadStepTracks, TTRHBuilder='WithoutRefit')
# TRACK SELECTION AND QUALITY FLAG SETTING.
from RecoTracker.FinalTrackSelectors.TrackMVAClassifierDetached_cfi import *
detachedQuadStep = TrackMVAClassifierDetached.clone(
src='detachedQuadStepTracks',
mva=dict(GBRForestLabel='MVASelectorDetachedQuadStep_Phase1'),
qualityCuts=[-0.5, 0.0, 0.5],
)
fastSim.toModify(detachedQuadStep,
vertices="firstStepPrimaryVerticesBeforeMixing")
highBetaStar_2018.toModify(detachedQuadStep, qualityCuts=[-0.7, 0.0, 0.5])
# For Phase2PU140
import RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi
detachedQuadStepSelector = RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.multiTrackSelector.clone(
MinApproachHigh = cms.double(1.0), #Track pair min distance at approaching point on X-Y
rCut = cms.double(2.0),#analytical track cross point
dz = cms.double(5.0),#track pair inner position difference
# kinematic vertex fit parameters
maxDelta = cms.double(0.01),#delta of parameters
maxReducedChiSq = cms.double(225.),#maximum chi^2 per degree of freedom before fit is terminated
minChiSqImprovement = cms.double(50.),#threshold for "significant improvement" in the fit termination logic
maxNbrOfIterations = cms.int32(40),#maximum number of convergence iterations
UsePvtx = cms.bool(True),
AllowD0 = cms.bool(True), #Allow d0*charge cut
AllowDeltaPhi = cms.bool(False),
AllowTrackBC = cms.bool(False), #Allow to match track-basic cluster
AllowDeltaCot = cms.bool(True), #Allow pairing using delta cot theta cut
AllowMinApproach = cms.bool(True), #Allow pairing using min approach cut
AllowOppCharge = cms.bool(True), #use opposite charge tracks to pair
AllowVertex = cms.bool(True),
bypassPreselGsf = cms.bool(True), #bypass preselection for gsf + X pairs
bypassPreselEcal = cms.bool(False), #bypass preselection for ecal-seeded + X pairs
bypassPreselEcalEcal = cms.bool(True), #bypass preselection for ecal-seeded + ecal-seeded pairs
AllowSingleLeg = cms.bool(False), #Allow single track conversion
AllowRightBC = cms.bool(False) #Require second leg matching basic cluster
)
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toModify( allConversions, bypassPreselGsf = cms.bool(False) )
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(allConversions, src = 'gsfGeneralConversionTrackMerger')
from RecoMuon.MuonIdentification.displacedMuonReducedTrackExtras_cfi import displacedMuonReducedTrackExtras
slimmedDisplacedMuonTrackExtras = displacedMuonReducedTrackExtras.clone(
muonTag="selectedPatDisplacedMuons",
cut="pt > 4.5",
trackExtraTags=[
"displacedMuonReducedTrackExtras", "displacedStandAloneMuons"
],
trackExtraAssocs=["displacedMuonReducedTrackExtras"],
pixelClusterTag="displacedMuonReducedTrackExtras",
stripClusterTag="displacedMuonReducedTrackExtras")
# no clusters in fastsim
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(slimmedMuonTrackExtras, outputClusters=False)
# cluster collections are different in phase 2, so skip this for now
from Configuration.Eras.Modifier_phase2_tracker_cff import phase2_tracker
phase2_tracker.toModify(slimmedMuonTrackExtras, outputClusters=False)
phase2_tracker.toModify(slimmedDisplacedMuonTrackExtras, outputClusters=False)
# lower minimum pt for B-parking
from Configuration.Eras.Modifier_bParking_cff import bParking
bParking.toModify(slimmedMuonTrackExtras, cut="pt > 3.0")
# full set of track extras not available in existing AOD
from Configuration.Eras.Modifier_run2_miniAOD_80XLegacy_cff import run2_miniAOD_80XLegacy
from Configuration.Eras.Modifier_run2_miniAOD_94XFall17_cff import run2_miniAOD_94XFall17
jetCoreRegionalStep.mva.maxChi2 = [9999.,9999.,9999.]
jetCoreRegionalStep.mva.maxChi2n = [1.6,1.0,0.7]
jetCoreRegionalStep.mva.minLayers = [3,5,5]
jetCoreRegionalStep.mva.min3DLayers = [1,2,3]
jetCoreRegionalStep.mva.maxLostLayers = [4,3,2]
jetCoreRegionalStep.mva.maxDz = [0.5,0.35,0.2];
jetCoreRegionalStep.mva.maxDr = [0.3,0.2,0.1];
jetCoreRegionalStep.vertices = 'firstStepGoodPrimaryVertices'
from RecoTracker.FinalTrackSelectors.TrackMVAClassifierPrompt_cfi import *
trackingPhase1.toReplaceWith(jetCoreRegionalStep, TrackMVAClassifierPrompt.clone(
src = 'jetCoreRegionalStepTracks',
mva = dict(GBRForestLabel = 'MVASelectorJetCoreRegionalStep_Phase1'),
qualityCuts = [-0.2,0.0,0.4],
))
fastSim.toModify(jetCoreRegionalStep,vertices = "firstStepPrimaryVerticesBeforeMixing")
# Final sequence
JetCoreRegionalStepTask = cms.Task(jetsForCoreTracking,
firstStepGoodPrimaryVertices,
#jetCoreRegionalStepClusters,
jetCoreRegionalStepSeedLayers,
jetCoreRegionalStepTrackingRegions,
jetCoreRegionalStepHitDoublets,
jetCoreRegionalStepSeeds,
jetCoreRegionalStepTrackCandidates,
jetCoreRegionalStepTracks,
# jetCoreRegionalStepClassifier1,jetCoreRegionalStepClassifier2,
jetCoreRegionalStep)
JetCoreRegionalStep = cms.Sequence(JetCoreRegionalStepTask)
fastSim.toReplaceWith(JetCoreRegionalStepTask,
import FWCore.ParameterSet.Config as cms
# SimHits Validation Analyzer after Simulation
from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
validSimHit = DQMEDAnalyzer('MuonSimHitsValidAnalyzer',
# Name of the root file which will contain the histos
DT_outputFile = cms.untracked.string(''),
Name = cms.untracked.string('MuonSimHitsValidAnalyzer'),
RPCHitsSrc = cms.InputTag("g4SimHits","MuonRPCHits"),
Verbosity = cms.untracked.int32(0),
ProvenanceLookup = cms.PSet(
PrintProvenanceInfo = cms.untracked.bool(False),
GetAllProvenances = cms.untracked.bool(False)
),
DTHitsSrc = cms.InputTag("g4SimHits","MuonDTHits"),
CSCHitsSrc = cms.InputTag("g4SimHits","MuonCSCHits"),
Label = cms.string('Hits')
)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(validSimHit,
DTHitsSrc = "MuonSimHits:MuonDTHits",
CSCHitsSrc = "MuonSimHits:MuonCSCHits",
RPCHitsSrc = "MuonSimHits:MuonRPCHits"
)
import FastSimulation.Tracking.TrackCandidateProducer_cfi
fastSim.toReplaceWith(initialStepTrackCandidates,
FastSimulation.Tracking.TrackCandidateProducer_cfi.trackCandidateProducer.clone(
src = cms.InputTag("initialStepSeeds"),
MinNumberOfCrossedLayers = 3
))
# fitting
import RecoTracker.TrackProducer.TrackProducer_cfi
initialStepTracks = RecoTracker.TrackProducer.TrackProducer_cfi.TrackProducer.clone(
src = 'initialStepTrackCandidates',
AlgorithmName = cms.string('initialStep'),
Fitter = cms.string('FlexibleKFFittingSmoother')
)
fastSim.toModify(initialStepTracks, TTRHBuilder = 'WithoutRefit')
#vertices
from RecoVertex.PrimaryVertexProducer.OfflinePrimaryVertices_cfi import offlinePrimaryVertices as _offlinePrimaryVertices
firstStepPrimaryVerticesUnsorted = _offlinePrimaryVertices.clone()
firstStepPrimaryVerticesUnsorted.TrackLabel = cms.InputTag("initialStepTracks")
firstStepPrimaryVerticesUnsorted.vertexCollections = [_offlinePrimaryVertices.vertexCollections[0].clone()]
# 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
light=cms.vuint32(1, 2, 3, 21), # udsg
c=cms.vuint32(4),
b=cms.vuint32(5),
g=cms.vuint32(21),
uds=cms.vuint32(1, 2, 3)),
mcPartons=cms.InputTag("hltBtagJetsbyValAlgo"))
#put all in a path
hltbtagValidationSequence = cms.Sequence(
# remove noisy warnings
# hltBtagTriggerSelection +
hltBtagJetMCTools + HltVertexValidationVertices + hltbTagValidation)
# fastsim customs
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(HltVertexValidationVertices,
SimVertexCollection="fastSimProducer")
# are these customs actually needed?
#HltVertexValidationVertices.HLTPathNames =cms.vstring(
#'HLT_PFMET120_NoiseCleaned_BTagCSV07_v',
#'HLT_PFMET120_NoiseCleaned_BTagCSV07_v',
# 'HLT_CaloMHTNoPU90_PFMET90_PFMHT90_IDLoose_',
# 'HLT_CaloMHTNoPU90_PFMET90_PFMHT90_IDLoose_',
# 'HLT_QuadPFJet_VBF',
# 'HLT_QuadPFJet_VBF',
# 'HLT_Ele32_eta2p1_',
# 'HLT_IsoMu24_eta2p1_')
#HltVertexValidationVertices.Vertex = cms.VInputTag(
# cms.InputTag("hltVerticesL3"),
# cms.InputTag("hltFastPVPixelVertices"),
# cms.InputTag("hltVerticesPF"),
#)
#
# Modify for running with the Stage 1 or Stage 2 trigger
#
from Configuration.Eras.Modifier_stage1L1Trigger_cff import stage1L1Trigger
from Configuration.Eras.Modifier_stage2L1Trigger_cff import stage2L1Trigger
_caloStage1LegacyFormatDigis = "caloStage1LegacyFormatDigis"
_params = dict(
etTotalSource = cms.InputTag(_caloStage1LegacyFormatDigis),
nonIsolatedEmSource = cms.InputTag(_caloStage1LegacyFormatDigis,"nonIsoEm"),
etMissSource = cms.InputTag(_caloStage1LegacyFormatDigis),
htMissSource = cms.InputTag(_caloStage1LegacyFormatDigis),
forwardJetSource = cms.InputTag(_caloStage1LegacyFormatDigis,"forJets"),
centralJetSource = cms.InputTag(_caloStage1LegacyFormatDigis,"cenJets"),
tauJetSource = cms.InputTag(_caloStage1LegacyFormatDigis,"tauJets"),
isoTauJetSource = cms.InputTag(_caloStage1LegacyFormatDigis,"isoTauJets"),
isolatedEmSource = cms.InputTag(_caloStage1LegacyFormatDigis,"isoEm"),
etHadSource = cms.InputTag(_caloStage1LegacyFormatDigis),
hfRingEtSumsSource = cms.InputTag(_caloStage1LegacyFormatDigis),
hfRingBitCountsSource = cms.InputTag(_caloStage1LegacyFormatDigis),
muonSource = cms.InputTag("gtDigis"),
centralBxOnly = True)
stage1L1Trigger.toModify( l1extraParticles, **_params)
stage2L1Trigger.toModify( l1extraParticles, **_params)
# fastsim runs L1Reco and HLT in one step
# this requires to set :
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(l1extraParticles, centralBxOnly = True)
import FWCore.ParameterSet.Config as cms
from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
ecalDigisValidation = DQMEDAnalyzer(
'EcalDigisValidation',
outputFile=cms.untracked.string(''),
verbose=cms.untracked.bool(False),
EBdigiCollection=cms.InputTag("simEcalDigis", "ebDigis"),
ESdigiCollection=cms.InputTag("simEcalPreshowerDigis"),
moduleLabelMC=cms.string('generatorSmeared'),
EEdigiCollection=cms.InputTag("simEcalDigis", "eeDigis"),
moduleLabelG4=cms.string('g4SimHits'))
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(ecalDigisValidation, moduleLabelG4='fastSimProducer')
HFDarkening = cms.bool(False),
minFCToDelay=cms.double(5.), # old TC model! set to 5 for the new one
debugCaloSamples=cms.bool(False),
ignoreGeantTime=cms.bool(False),
# settings for SimHit test injection
injectTestHits = cms.bool(False),
# if no time is specified for injected hits, t = 0 will be used
# (recommendation: enable "ignoreGeantTime" in that case to set t = tof)
# otherwise, need 1 time value per energy value
injectTestHitsEnergy = cms.vdouble(),
injectTestHitsTime = cms.vdouble(),
# format for cells: subdet, ieta, iphi, depth
# multiple quadruplets can be specified
# if instead only 1 value is given,
# it will be interpreted as an entire subdetector
injectTestHitsCells = cms.vint32()
)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify( hcalSimBlock, hitsProducer=cms.string('famosSimHits') )
from Configuration.Eras.Modifier_run2_HCAL_2017_cff import run2_HCAL_2017
run2_HCAL_2017.toModify( hcalSimBlock, TestNumbering = cms.bool(True) )
# remove HE processing for phase 2, completely put in HGCal land
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toModify(hcalSimBlock, killHE = cms.bool(True) )
from Configuration.Eras.Modifier_phase2_common_cff import phase2_common
phase2_common.toModify(hcalSimBlock, doNeutralDensityFilter = cms.bool(False))
from RecoEgamma.EgammaElectronProducers.gsfElectrons_cfi import ecalDrivenGsfElectrons lowPtGsfElectronsPreRegression = ecalDrivenGsfElectrons.clone(gsfElectronCoresTag = "lowPtGsfElectronCores") from Configuration.Eras.Modifier_fastSim_cff import fastSim fastSim.toModify(lowPtGsfElectronsPreRegression,ctfTracksTag = "generalTracksBeforeMixing")
import FWCore.ParameterSet.Config as cms
from RecoMuon.TrackingTools.MuonServiceProxy_cff import *
from RecoMuon.TrackingTools.MuonTrackLoader_cff import *
from RecoMuon.GlobalTrackingTools.GlobalTrajectoryBuilderCommon_cff import *
from RecoMuon.GlobalTrackingTools.GlobalMuonRefitter_cff import *
tevMuons = cms.EDProducer("TevMuonProducer",
MuonTrackLoaderForGLB,
# InputTag MuonCollectionLabel = standAloneMuons:UpdatedAtVtx
MuonServiceProxy,
RefitIndex = cms.vint32(1, 2, 3, 4),
Refits = cms.vstring('default',
'firstHit',
'picky',
'dyt'),
MuonCollectionLabel = cms.InputTag("globalMuons"),
RefitterParameters = cms.PSet(
GlobalMuonRefitter
)
)
# FastSim has no template fit on tracker hits
# FastSim doesn't use Runge Kute for propagation
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(tevMuons,
RefitterParameters = dict(TrackerRecHitBuilder = 'WithoutRefit',
Propagator = "SmartPropagatorAny")
)
import FWCore.ParameterSet.Config as cms
hcalSimHitStudy = cms.EDAnalyzer(
"HcalSimHitStudy",
ModuleLabel=cms.untracked.string('g4SimHits'),
outputFile=cms.untracked.string(''),
Verbose=cms.untracked.bool(False),
HitCollection=cms.untracked.string('HcalHits'),
TestNumber=cms.bool(False))
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(hcalSimHitStudy,
ModuleLabel=cms.untracked.string('famosSimHits'))
from RecoTracker.TkTrackingRegions.globalTrackingRegionWithVertices_cff import globalTrackingRegionWithVertices as _globalTrackingRegionWithVertices
pixelPairStepTrackingRegions = _globalTrackingRegionWithVertices.clone(RegionPSet = dict(
ptMin = 0.6,
originRadius = 0.015,
fixedError = 0.03,
useMultipleScattering = True,
))
from Configuration.Eras.Modifier_trackingLowPU_cff import trackingLowPU
trackingLowPU.toModify(pixelPairStepTrackingRegions, RegionPSet=dict(useMultipleScattering=False))
_region_Phase1 = dict(
useMultipleScattering = False,
maxNVertices = 5,
)
trackingPhase1.toModify(pixelPairStepTrackingRegions, RegionPSet=_region_Phase1)
trackingPhase2PU140.toModify(pixelPairStepTrackingRegions, RegionPSet=_region_Phase1)
fastSim.toModify(pixelPairStepTrackingRegions, RegionPSet=dict(VertexCollection = "firstStepPrimaryVerticesBeforeMixing"))
# SEEDS
from RecoTracker.TkHitPairs.hitPairEDProducer_cfi import hitPairEDProducer as _hitPairEDProducer
pixelPairStepHitDoublets = _hitPairEDProducer.clone(
seedingLayers = "pixelPairStepSeedLayers",
trackingRegions = "pixelPairStepTrackingRegions",
produceSeedingHitSets = True,
)
from RecoTracker.TkSeedGenerator.seedCreatorFromRegionConsecutiveHitsEDProducer_cff import seedCreatorFromRegionConsecutiveHitsEDProducer as _seedCreatorFromRegionConsecutiveHitsEDProducer
pixelPairStepSeeds = _seedCreatorFromRegionConsecutiveHitsEDProducer.clone(
seedingHitSets = "pixelPairStepHitDoublets",
SeedComparitorPSet = dict(# FIXME: is this defined in any cfi that could be imported instead of copy-paste?
ComponentName = 'PixelClusterShapeSeedComparitor',
FilterAtHelixStage = cms.bool(True),
FilterPixelHits = cms.bool(True),
leptonPtThreshold = cms.untracked.double(55.0),
htThreshold = cms.untracked.double(500.0),
metThreshold = cms.untracked.double(-1.0),
csvThreshold = cms.untracked.double(-1.0)
)
SUSY_HLT_Ele50_HT400_SingleLepton_POSTPROCESSING = cms.EDAnalyzer('DQMGenericClient',
subDirs = cms.untracked.vstring('HLT/SUSYBSM/HLT_Ele50_IsoVVVL_PFHT400'),
efficiency = cms.vstring(
"leptonTurnOn_eff ';Offline Electron p_{T} [GeV];#epsilon' leptonTurnOn_num leptonTurnOn_den",
"pfHTTurnOn_eff ';Offline PF H_{T} [GeV];#epsilon' pfHTTurnOn_num pfHTTurnOn_den"
),
resolution = cms.vstring('')
)
# fastsim has no conversion collection (yet)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(SUSY_HLT_Ele15_HT600_SingleLepton,conversionCollection=cms.InputTag(''))
fastSim.toModify(SUSY_HLT_Ele15_HT400_SingleLepton,conversionCollection=cms.InputTag(''))
fastSim.toModify(SUSY_HLT_Ele50_HT400_SingleLepton,conversionCollection=cms.InputTag(''))
SUSY_HLT_Ele_HT_SingleLepton = cms.Sequence( SUSY_HLT_Ele15_HT600_SingleLepton
+ SUSY_HLT_Ele15_HT400_SingleLepton
+ SUSY_HLT_Ele50_HT400_SingleLepton
)
SUSY_HLT_Ele_HT_SingleLepton_POSTPROCESSING = cms.Sequence( SUSY_HLT_Ele15_HT600_SingleLepton_POSTPROCESSING
+ SUSY_HLT_Ele15_HT400_SingleLepton_POSTPROCESSING
+ SUSY_HLT_Ele50_HT400_SingleLepton_POSTPROCESSING
)
#
# NANOAOD
#
#
from PhysicsTools.NanoAOD.NanoAODEDMEventContent_cff import *
#
#
# FastSim
#
#
from FastSimulation.Configuration.EventContent_cff import FASTPUEventContent
import FastSimulation.Configuration.EventContent_cff as fastSimEC
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(RecoLocalTrackerRECO, outputCommands = fastSimEC.RecoLocalTracker.outputCommands)
fastSim.toModify(RecoLocalTrackerFEVT, outputCommands = fastSimEC.RecoLocalTracker.outputCommands)
fastSim.toReplaceWith(SimG4CoreRAW, fastSimEC.SimRAW)
fastSim.toReplaceWith(SimG4CoreRECO, fastSimEC.SimRECO)
#
#
# Top level additional keep statements
#
#
CommonEventContent = cms.PSet(
outputCommands = cms.untracked.vstring('keep *_logErrorHarvester_*_*')
)
#
#
from FastSimulation.Tracking.recoTrackAccumulator_cfi import recoTrackAccumulator as _recoTrackAccumulator
fastSim.toModify(mixData,
# from signal: mix tracks not strip or pixel digis
workers = dict(
pixel = None,
strip = None,
pixelSimLink = None,
stripSimLink = None,
tracks = cms.PSet(
workerType = cms.string("PreMixingDigiAccumulatorWorker"),
accumulator = _recoTrackAccumulator.clone(
pileUpTracks = "mix:generalTracks"
)
),
dt = dict(
labelSig = "mix:MuonSimHitsMuonDTHits",
pileInputTag = "mix:MuonSimHitsMuonDTHits",
collectionDM = "MuonSimHitsMuonDTHits",
),
rpc = dict(
labelSig = "mix:MuonSimHitsMuonRPCHits",
pileInputTag = "mix:MuonSimHitsMuonRPCHits",
collectionDM = "MuonSimHitsMuonRPCHits",
),
csc = dict(
labelSig = "mix:MuonSimHitsMuonCSCHits",
pileInputTag = "mix:MuonSimHitsMuonCSCHits",
collectionDM = "MuonSimHitsMuonCSCHits",
),
),
)
'TrackerHitsPixelEndcapLowTof',
'TrackerHitsPixelEndcapHighTof'),
#
# to associate to reco::Muon segments
inputDTRecSegment4DCollection = cms.InputTag("dt4DSegments"),
inputCSCSegmentCollection = cms.InputTag("cscSegments"),
)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(muonAssociatorByHitsCommonParameters,
simtracksTag = "fastSimProducer",
DTsimhitsTag = "MuonSimHits:MuonDTHits",
CSCsimHitsTag = "MuonSimHits:MuonCSCHits",
RPCsimhitsTag = "MuonSimHits:MuonRPCHits",
simtracksXFTag = "mix:fastSimProducer",
DTsimhitsXFTag = "mix:MuonSimHitsMuonDTHits",
CSCsimHitsXFTag = "mix:MuonSimHitsMuonCSCHits",
RPCsimhitsXFTag = "mix:MuonSimHitsMuonRPCHits",
ROUList = ['fastSimProducerTrackerHits']
)
muonAssociatorByHits = cms.EDProducer("MuonAssociatorEDProducer",
# COMMON CONFIGURATION
muonAssociatorByHitsCommonParameters,
# for Muon Track association
#
# input collections
#
# ... reco::Track collection
# tracksTag = cms.InputTag("standAloneMuons"),
SiStripAPVListDM = cms.string('SiStripAPVList'),
#
# Calorimeter Digis
#
EBDigiCollectionDM = cms.string(''),
EEDigiCollectionDM = cms.string(''),
ESDigiCollectionDM = cms.string(''),
HBHEDigiCollectionDM = cms.string(''),
HODigiCollectionDM = cms.string(''),
HFDigiCollectionDM = cms.string(''),
QIE10DigiCollectionDM = cms.string(''),
QIE11DigiCollectionDM = cms.string(''),
ZDCDigiCollectionDM = cms.string('')
)
mixData.doEB = cms.bool(True)
mixData.doEE = cms.bool(True)
mixData.doES = cms.bool(True)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
from FastSimulation.Tracking.recoTrackAccumulator_cfi import recoTrackAccumulator as _recoTrackAccumulator
fastSim.toModify(mixData,
# from signal: mix tracks not strip or pixel digis
TrackerMergeType = "tracks",
tracker = _recoTrackAccumulator.clone(
pileUpTracks = "mix:generalTracks"
),
hitsProducer = "famosSimHits"
)
pdgId=cms.vint32(),
signalOnly=cms.bool(True),
intimeOnly=cms.bool(False),
stableOnly=cms.bool(False),
minRapidity=cms.double(-2.5),
minHit=cms.int32(0),
ptMin=cms.double(0.9),
ptMax=cms.double(1e100),
maxRapidity=cms.double(2.5),
tip=cms.double(3.5),
minPhi=cms.double(-3.2),
maxPhi=cms.double(3.2),
)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(generalTpSelectorBlock, stableOnly=True)
TpSelectorForEfficiencyVsEtaBlock = generalTpSelectorBlock.clone()
TpSelectorForEfficiencyVsPhiBlock = generalTpSelectorBlock.clone()
TpSelectorForEfficiencyVsPtBlock = generalTpSelectorBlock.clone(ptMin=0.050)
TpSelectorForEfficiencyVsVTXRBlock = generalTpSelectorBlock.clone(tip=60.0)
TpSelectorForEfficiencyVsVTXZBlock = generalTpSelectorBlock.clone()
def _modifyForPhase1(pset):
pset.minRapidity = -3
pset.maxRapidity = 3
from Configuration.Eras.Modifier_phase1Pixel_cff import phase1Pixel
phase1Pixel.toModify(generalTpSelectorBlock, _modifyForPhase1)
NEWtpToTevFirstMuonAssociation+NEWtpToTevPickyMuonAssociation+NEWtpToTevDytMuonAssociation
)
NewMuonAssociationDisplaced_seq = cms.Sequence(
NEWseedsOfDisplacedSTAmuons_seq+NEWtpToDisplacedStaSeedAssociation+NEWtpToDisplacedStaMuonAssociation
+NEWtpToDisplacedTrkMuonAssociation+NEWtpToDisplacedGlbMuonAssociation
)
NewMuonAssociationRefit_seq = cms.Sequence(
NEWtpToStaRefitMuonAssociation+NEWtpToStaRefitUpdMuonAssociation
)
NewMuonAssociationCosmic_seq = cms.Sequence(
NEWtpToTkCosmicSelMuonAssociation+ NEWtpToTkCosmic1LegSelMuonAssociation
+NEWtpToStaCosmicSelMuonAssociation+NEWtpToStaCosmic1LegSelMuonAssociation
+NEWtpToGlbCosmicSelMuonAssociation+NEWtpToGlbCosmic1LegSelMuonAssociation
)
NewMuonAssociationHLT_seq = cms.Sequence(
NEWtpToL2MuonAssociation+NEWtpToL2UpdMuonAssociation+NEWtpToL3TkMuonAssociation+NEWtpToL3MuonAssociation
)
# fastsim has no hlt specific dt hit collection
from Configuration.Eras.Modifier_fastSim_cff import fastSim
_DTrechitTag = SimMuon.MCTruth.NewMuonAssociatorByHits_cfi.NewMuonAssociatorByHits.DTrechitTag
fastSim.toModify(NEWtpToL3TkMuonAssociation, DTrechitTag = _DTrechitTag)
fastSim.toModify(NEWtpToL2MuonAssociation, DTrechitTag = _DTrechitTag)
fastSim.toModify(NEWtpToL2UpdMuonAssociation, DTrechitTag = _DTrechitTag)
fastSim.toModify(NEWtpToL3MuonAssociation, DTrechitTag = _DTrechitTag)
import FWCore.ParameterSet.Config as cms
rpcRecHitV = cms.EDAnalyzer("RPCRecHitValid",
subDir = cms.string("RPC/RPCRecHitV/SimVsReco"),
simHit = cms.InputTag("g4SimHits", "MuonRPCHits"),
recHit = cms.InputTag("rpcRecHits"),
simTrack = cms.InputTag("mix", "MergedTrackTruth"),
simHitAssoc = cms.InputTag("simHitTPAssocProducer"),
muon = cms.InputTag("muons"),
)
rpcRecHitValidation_step = cms.Sequence(rpcRecHitV)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(rpcRecHitV, simHit = "MuonSimHits:MuonRPCHits")
import FWCore.ParameterSet.Config as cms from RecoEgamma.EgammaElectronProducers.lowPtGsfElectronCores_cfi import lowPtGsfElectronCores from Configuration.Eras.Modifier_fastSim_cff import fastSim fastSim.toModify(lowPtGsfElectronCores, ctfTracks="generalTracksBeforeMixing")
RecHitSeverityToBeExcluded,
recHitSeverityToBeExcludedEndcaps=cleanedHybridSuperClusters.
RecHitSeverityToBeExcluded,
# Isolation algos configuration
trkIsol03Cfg=trkIsol03CfgV1,
trkIsol04Cfg=trkIsol04CfgV1,
trkIsolHEEP03Cfg=trkIsol03CfgV2,
trkIsolHEEP04Cfg=trkIsol04CfgV2,
# regression. The labels are needed in all cases.
ecalRefinedRegressionWeightLabels=cms.vstring(),
#"gedelectron_EBCorrection_offline_v1"
#"gedelectron_EECorrection_offline_v1"
#"gedelectron_EBUncertainty_offline_v1"
#"gedelectron_EEUncertainty_offline_v1"
combinationRegressionWeightLabels=cms.vstring(),
#"gedelectron_p4combination_offline"
# Iso values
useIsolationValues=cms.bool(False),
)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(lowPtGsfElectrons,
ctfTracksTag=cms.InputTag("generalTracksBeforeMixing"))
from Configuration.Eras.Modifier_pp_on_AA_2018_cff import pp_on_AA_2018
pp_on_AA_2018.toModify(lowPtGsfElectrons.preselection, minSCEtBarrel=15.0)
pp_on_AA_2018.toModify(lowPtGsfElectrons.preselection, minSCEtEndcaps=15.0)
from Configuration.Eras.Modifier_run3_GEM_cff import run3_GEM
from Configuration.Eras.Modifier_phase2_tracker_cff import phase2_tracker
run3_GEM.toModify(muonAssociatorByHitsCommonParameters, useGEMs=True)
phase2_tracker.toModify(
muonAssociatorByHitsCommonParameters,
usePhase2Tracker=True,
pixelSimLinkSrc="simSiPixelDigis:Pixel",
)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(muonAssociatorByHitsCommonParameters,
simtracksTag="fastSimProducer",
DTsimhitsTag="MuonSimHits:MuonDTHits",
CSCsimHitsTag="MuonSimHits:MuonCSCHits",
RPCsimhitsTag="MuonSimHits:MuonRPCHits",
simtracksXFTag="mix:fastSimProducer",
DTsimhitsXFTag="mix:MuonSimHitsMuonDTHits",
CSCsimHitsXFTag="mix:MuonSimHitsMuonCSCHits",
RPCsimhitsXFTag="mix:MuonSimHitsMuonRPCHits",
ROUList=['fastSimProducerTrackerHits'])
muonAssociatorByHits = cms.EDProducer(
"MuonAssociatorEDProducer",
# COMMON CONFIGURATION
muonAssociatorByHitsCommonParameters,
# for Muon Track association
#
# input collections
#
# ... reco::Track collection
# tracksTag = cms.InputTag("standAloneMuons"),
theDigitizersMixPreMix = cms.PSet(pixel=cms.PSet(pixelDigitizer),
strip=cms.PSet(stripDigitizer),
castor=cms.PSet(castorDigitizer),
puVtx=cms.PSet(pileupVtxDigitizer),
mergedtruth=cms.PSet(trackingParticles))
theDigitizersMixPreMix.strip.Noise = cms.bool(
False) # will be added in DataMixer
theDigitizersMixPreMix.strip.PreMixingMode = cms.bool(
True) #Special mode to save all hit strips
theDigitizersMixPreMix.strip.FedAlgorithm = cms.int32(
5) # special ZS mode: accept adc>0
theDigitizersMixPreMix.pixel.AddPixelInefficiency = cms.bool(
False) # will be added in DataMixer
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(
theDigitizersMixPreMix,
# fastsim does not model castor
castor=None,
# fastsim does not digitize pixel and strip hits
pixel=None,
strip=None)
from Configuration.Eras.Modifier_run3_common_cff import run3_common
run3_common.toModify(theDigitizersMixPreMix, castor=None)
theDigitizersMixPreMixValid = cms.PSet(theDigitizersMixPreMix)
theDigitizersMixPreMix.mergedtruth.select.signalOnlyTP = cms.bool(True)
pileInputTag=cms.InputTag("simMuonCSCDigis",
"MuonCSCWireDigiSimLinks"),
collectionDM=cms.string('MuonCSCWireDigiSimLinks'),
),
),
)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
from FastSimulation.Tracking.recoTrackAccumulator_cfi import recoTrackAccumulator as _recoTrackAccumulator
fastSim.toModify(
mixData,
# from signal: mix tracks not strip or pixel digis
workers=dict(pixel=None,
strip=None,
pixelSimLink=None,
stripSimLink=None,
tracks=cms.PSet(
workerType=cms.string("PreMixingDigiAccumulatorWorker"),
accumulator=_recoTrackAccumulator.clone(
pileUpTracks="mix:generalTracks"))),
)
from Configuration.Eras.Modifier_phase2_common_cff import phase2_common
from Configuration.Eras.Modifier_phase2_tracker_cff import phase2_tracker
from Configuration.Eras.Modifier_phase2_ecal_cff import phase2_ecal
from Configuration.Eras.Modifier_phase2_hcal_cff import phase2_hcal
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
from Configuration.Eras.Modifier_phase2_hfnose_cff import phase2_hfnose
from Configuration.Eras.Modifier_phase2_muon_cff import phase2_muon
phase2_common.toModify(mixData, input=dict(
import FWCore.ParameterSet.Config as cms
from Validation.CSCRecHits.cscRecHitPSet import *
from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
cscRecHitValidation = DQMEDAnalyzer('CSCRecHitValidation',
cscRecHitPSet,
doSim=cms.bool(True),
useGEMs=cms.bool(False),
simHitsTag=cms.InputTag(
"mix", "g4SimHitsMuonCSCHits"))
from Configuration.Eras.Modifier_run3_GEM_cff import run3_GEM
run3_GEM.toModify(cscRecHitValidation, useGEMs=True)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(cscRecHitValidation, simHitsTag="mix:MuonSimHitsMuonCSCHits")
import FastSimulation.Tracking.TrackCandidateProducer_cfi
_fastSim_highPtTripletStepTrackCandidates = FastSimulation.Tracking.TrackCandidateProducer_cfi.trackCandidateProducer.clone(
src='highPtTripletStepSeeds',
MinNumberOfCrossedLayers=3,
hitMasks=cms.InputTag('highPtTripletStepMasks'))
fastSim.toReplaceWith(highPtTripletStepTrackCandidates,
_fastSim_highPtTripletStepTrackCandidates)
# TRACK FITTING
import RecoTracker.TrackProducer.TrackProducer_cfi
highPtTripletStepTracks = RecoTracker.TrackProducer.TrackProducer_cfi.TrackProducer.clone(
src='highPtTripletStepTrackCandidates',
AlgorithmName='highPtTripletStep',
Fitter='FlexibleKFFittingSmoother',
)
fastSim.toModify(highPtTripletStepTracks, TTRHBuilder='WithoutRefit')
# Final selection
from RecoTracker.FinalTrackSelectors.TrackMVAClassifierPrompt_cfi import *
highPtTripletStep = TrackMVAClassifierPrompt.clone(
mva=dict(GBRForestLabel='MVASelectorHighPtTripletStep_Phase1'),
src='highPtTripletStepTracks',
qualityCuts=[0.2, 0.3, 0.4])
from RecoTracker.FinalTrackSelectors.TrackLwtnnClassifier_cfi import *
from RecoTracker.FinalTrackSelectors.trackSelectionLwtnn_cfi import *
trackdnn.toReplaceWith(
highPtTripletStep,
TrackLwtnnClassifier.clone(
src='highPtTripletStepTracks',
qualityCuts=[0.75, 0.775, 0.8],
from TrackingTools.TrajectoryCleaning.TrajectoryCleanerBySharedHits_cfi import trajectoryCleanerBySharedHits
pixelLessStepTrajectoryCleanerBySharedHits = trajectoryCleanerBySharedHits.clone(
ComponentName='pixelLessStepTrajectoryCleanerBySharedHits',
fractionShared=0.11,
allowSharedFirstHit=True)
trackingLowPU.toModify(pixelLessStepTrajectoryCleanerBySharedHits,
fractionShared=0.19)
# TRACK FITTING
import RecoTracker.TrackProducer.TrackProducer_cfi
pixelLessStepTracks = RecoTracker.TrackProducer.TrackProducer_cfi.TrackProducer.clone(
src='pixelLessStepTrackCandidates',
AlgorithmName='pixelLessStep',
Fitter='FlexibleKFFittingSmoother')
fastSim.toModify(pixelLessStepTracks, TTRHBuilder='WithoutRefit')
# TRACK SELECTION AND QUALITY FLAG SETTING.
from RecoTracker.FinalTrackSelectors.TrackMVAClassifierPrompt_cfi import *
from RecoTracker.FinalTrackSelectors.TrackMVAClassifierDetached_cfi import *
pixelLessStepClassifier1 = TrackMVAClassifierPrompt.clone(
src='pixelLessStepTracks',
mva=dict(GBRForestLabel='MVASelectorIter5_13TeV'),
qualityCuts=[-0.4, 0.0, 0.4])
fastSim.toModify(pixelLessStepClassifier1,
vertices='firstStepPrimaryVerticesBeforeMixing')
pixelLessStepClassifier2 = TrackMVAClassifierPrompt.clone(
src='pixelLessStepTracks',
mva=dict(GBRForestLabel='MVASelectorIter0_13TeV'),
qualityCuts=[-0.0, 0.0, 0.0])
import TrackingTools.TrackFitters.FlexibleKFFittingSmoother_cfi
tobTecFlexibleKFFittingSmoother = TrackingTools.TrackFitters.FlexibleKFFittingSmoother_cfi.FlexibleKFFittingSmoother.clone(
ComponentName=cms.string('tobTecFlexibleKFFittingSmoother'),
standardFitter=cms.string('tobTecStepFitterSmoother'),
looperFitter=cms.string('tobTecStepFitterSmootherForLoopers'),
)
# TRACK FITTING
import RecoTracker.TrackProducer.TrackProducer_cfi
tobTecStepTracks = RecoTracker.TrackProducer.TrackProducer_cfi.TrackProducer.clone(
src='tobTecStepTrackCandidates',
AlgorithmName=cms.string('tobTecStep'),
#Fitter = 'tobTecStepFitterSmoother',
Fitter='tobTecFlexibleKFFittingSmoother',
)
fastSim.toModify(tobTecStepTracks, TTRHBuilder='WithoutRefit')
# TRACK SELECTION AND QUALITY FLAG SETTING.
from RecoTracker.FinalTrackSelectors.TrackMVAClassifierPrompt_cfi import *
from RecoTracker.FinalTrackSelectors.TrackMVAClassifierDetached_cfi import *
tobTecStepClassifier1 = TrackMVAClassifierDetached.clone()
tobTecStepClassifier1.src = 'tobTecStepTracks'
tobTecStepClassifier1.mva.GBRForestLabel = 'MVASelectorIter6_13TeV'
tobTecStepClassifier1.qualityCuts = [-0.6, -0.45, -0.3]
fastSim.toModify(tobTecStepClassifier1,
vertices="firstStepPrimaryVerticesBeforeMixing")
tobTecStepClassifier2 = TrackMVAClassifierPrompt.clone()
tobTecStepClassifier2.src = 'tobTecStepTracks'
tobTecStepClassifier2.mva.GBRForestLabel = 'MVASelectorIter0_13TeV'
tobTecStepClassifier2.qualityCuts = [0.0, 0.0, 0.0]
### Allow switching off particular histograms
doSummaryPlots=cms.untracked.bool(True),
doSimPlots=cms.untracked.bool(True),
doSimTrackPlots=cms.untracked.bool(True),
doRecoTrackPlots=cms.untracked.bool(True),
dodEdxPlots=cms.untracked.bool(False),
doPVAssociationPlots=cms.untracked.bool(
False
), # do plots that require true PV, if True, label_vertex and vertexAssociator are read
doSeedPlots=cms.untracked.bool(
False), # input comes from TrackFromSeedProducer
doMVAPlots=cms.untracked.bool(
False), # needs input from track MVA selectors
### do resolution plots only for these labels (or all if empty)
doResolutionPlotsForLabels=cms.VInputTag(),
cores=cms.InputTag("highPtJetsForTrk"), #ak4CaloJets with pt>1 TeV
)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(multiTrackValidator, sim=['fastSimProducer:TrackerHits'])
from Configuration.ProcessModifiers.premix_stage2_cff import premix_stage2
premix_stage2.toModify(
multiTrackValidator,
label_tp_effic="mixData:MergedTrackTruth",
label_tp_fake="mixData:MergedTrackTruth",
label_tv="mixData:MergedTrackTruth",
)
doSeedingRegionRebuilding=True,
useHitsSplitting=True)
import FastSimulation.Tracking.TrackCandidateProducer_cfi
fastSim.toReplaceWith(
initialStepTrackCandidates,
FastSimulation.Tracking.TrackCandidateProducer_cfi.trackCandidateProducer.
clone(src=cms.InputTag("initialStepSeeds"), MinNumberOfCrossedLayers=3))
# fitting
import RecoTracker.TrackProducer.TrackProducer_cfi
initialStepTracks = RecoTracker.TrackProducer.TrackProducer_cfi.TrackProducer.clone(
src='initialStepTrackCandidates',
AlgorithmName=cms.string('initialStep'),
Fitter=cms.string('FlexibleKFFittingSmoother'))
fastSim.toModify(initialStepTracks, TTRHBuilder='WithoutRefit')
#vertices
from RecoVertex.PrimaryVertexProducer.OfflinePrimaryVertices_cfi import offlinePrimaryVertices as _offlinePrimaryVertices
firstStepPrimaryVerticesUnsorted = _offlinePrimaryVertices.clone()
firstStepPrimaryVerticesUnsorted.TrackLabel = cms.InputTag("initialStepTracks")
firstStepPrimaryVerticesUnsorted.vertexCollections = [
_offlinePrimaryVertices.vertexCollections[0].clone()
]
from Configuration.Eras.Modifier_pp_on_XeXe_2017_cff import pp_on_XeXe_2017
from Configuration.Eras.Modifier_pp_on_AA_2018_cff import pp_on_AA_2018
(pp_on_XeXe_2017 | pp_on_AA_2018).toModify(
firstStepPrimaryVerticesUnsorted,
TkFilterParameters=dict(trackQuality="any"))
import FWCore.ParameterSet.Config as cms
ecalSimHitsValidation = cms.EDAnalyzer("EcalSimHitsValidation",
ESHitsCollection = cms.string('EcalHitsES'),
outputFile = cms.untracked.string(''),
verbose = cms.untracked.bool(False),
moduleLabelMC = cms.string('generatorSmeared'),
EBHitsCollection = cms.string('EcalHitsEB'),
EEHitsCollection = cms.string('EcalHitsEE'),
moduleLabelG4 = cms.string('g4SimHits')
)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(ecalSimHitsValidation, moduleLabelG4 = "famosSimHits")
