12 process.patMuons.isoDeposits = cms.PSet()
13 process.patElectrons.isoDeposits = cms.PSet()
14 process.patTaus.isoDeposits = cms.PSet()
15 process.patPhotons.isoDeposits = cms.PSet()
17 process.patMuons.embedTrack =
True 18 process.patMuons.embedCombinedMuon =
True 19 process.patMuons.embedMuonBestTrack =
True 20 process.patMuons.embedStandAloneMuon =
True 21 process.patMuons.embedPickyMuon =
False 22 process.patMuons.embedTpfmsMuon =
False 23 process.patMuons.embedDytMuon =
False 24 process.patMuons.addPuppiIsolation = cms.bool(
True)
25 process.patMuons.puppiIsolationChargedHadrons = cms.InputTag(
"muonPUPPIIsolation",
"h+-DR040-ThresholdVeto000-ConeVeto000")
26 process.patMuons.puppiIsolationNeutralHadrons = cms.InputTag(
"muonPUPPIIsolation",
"h0-DR040-ThresholdVeto000-ConeVeto001")
27 process.patMuons.puppiIsolationPhotons = cms.InputTag(
"muonPUPPIIsolation",
"gamma-DR040-ThresholdVeto000-ConeVeto001")
28 process.patMuons.puppiNoLeptonsIsolationChargedHadrons = cms.InputTag(
"muonPUPPINoLeptonsIsolation",
"h+-DR040-ThresholdVeto000-ConeVeto000")
29 process.patMuons.puppiNoLeptonsIsolationNeutralHadrons = cms.InputTag(
"muonPUPPINoLeptonsIsolation",
"h0-DR040-ThresholdVeto000-ConeVeto001")
30 process.patMuons.puppiNoLeptonsIsolationPhotons = cms.InputTag(
"muonPUPPINoLeptonsIsolation",
"gamma-DR040-ThresholdVeto000-ConeVeto001")
32 process.patMuons.computeMiniIso =
True 33 process.patMuons.computeMuonIDMVA =
True 34 process.patMuons.computeSoftMuonMVA =
True 36 process.patMuons.addTriggerMatching =
True 37 from Configuration.Eras.Modifier_run2_muon_2016_cff
import run2_muon_2016
38 from Configuration.Eras.Modifier_run2_muon_2017_cff
import run2_muon_2017
39 from Configuration.Eras.Modifier_run2_muon_2018_cff
import run2_muon_2018
40 from Configuration.Eras.Modifier_run3_muon_cff
import run3_muon
41 run2_muon_2016.toModify( process.patMuons, effectiveAreaVec = [0.0735,0.0619,0.0465,0.0433,0.0577])
42 run2_muon_2017.toModify( process.patMuons, effectiveAreaVec = [0.0566, 0.0562, 0.0363, 0.0119, 0.0064])
43 run2_muon_2018.toModify( process.patMuons, effectiveAreaVec = [0.0566, 0.0562, 0.0363, 0.0119, 0.0064])
44 run3_muon.toModify( process.patMuons, effectiveAreaVec = [0.0566, 0.0562, 0.0363, 0.0119, 0.0064])
46 process.patMuons.computePuppiCombinedIso =
True 49 process.patElectrons.embedGsfElectronCore =
False 50 process.patElectrons.embedSuperCluster =
False 51 process.patElectrons.embedPflowSuperCluster =
False 52 process.patElectrons.embedSeedCluster =
False 53 process.patElectrons.embedBasicClusters =
False 54 process.patElectrons.embedPreshowerClusters =
False 55 process.patElectrons.embedPflowBasicClusters =
False 56 process.patElectrons.embedPflowPreshowerClusters =
False 57 process.patElectrons.embedRecHits =
False 58 process.patElectrons.electronSource = cms.InputTag(
"reducedEgamma",
"reducedGedGsfElectrons")
59 process.patElectrons.usePfCandidateMultiMap =
True 60 process.patElectrons.pfCandidateMultiMap = cms.InputTag(
"reducedEgamma",
"reducedGsfElectronPfCandMap")
61 process.patElectrons.electronIDSources = cms.PSet()
64 process.patElectrons.addPuppiIsolation = cms.bool(
True)
65 process.patElectrons.puppiIsolationChargedHadrons = cms.InputTag(
"egmElectronPUPPIIsolation",
"h+-DR030-BarVeto000-EndVeto001")
66 process.patElectrons.puppiIsolationNeutralHadrons = cms.InputTag(
"egmElectronPUPPIIsolation",
"h0-DR030-BarVeto000-EndVeto000")
67 process.patElectrons.puppiIsolationPhotons = cms.InputTag(
"egmElectronPUPPIIsolation",
"gamma-DR030-BarVeto000-EndVeto008")
68 process.patElectrons.puppiNoLeptonsIsolationChargedHadrons = cms.InputTag(
"egmElectronPUPPINoLeptonsIsolation",
"h+-DR030-BarVeto000-EndVeto001")
69 process.patElectrons.puppiNoLeptonsIsolationNeutralHadrons = cms.InputTag(
"egmElectronPUPPINoLeptonsIsolation",
"h0-DR030-BarVeto000-EndVeto000")
70 process.patElectrons.puppiNoLeptonsIsolationPhotons = cms.InputTag(
"egmElectronPUPPINoLeptonsIsolation",
"gamma-DR030-BarVeto000-EndVeto008")
72 process.patElectrons.computeMiniIso = cms.bool(
True)
74 process.elPFIsoDepositChargedPAT.src = cms.InputTag(
"reducedEgamma",
"reducedGedGsfElectrons")
75 process.elPFIsoDepositChargedAllPAT.src = cms.InputTag(
"reducedEgamma",
"reducedGedGsfElectrons")
76 process.elPFIsoDepositNeutralPAT.src = cms.InputTag(
"reducedEgamma",
"reducedGedGsfElectrons")
77 process.elPFIsoDepositGammaPAT.src = cms.InputTag(
"reducedEgamma",
"reducedGedGsfElectrons")
78 process.elPFIsoDepositPUPAT.src = cms.InputTag(
"reducedEgamma",
"reducedGedGsfElectrons")
80 process.patPhotons.embedSuperCluster =
False 81 process.patPhotons.embedSeedCluster =
False 82 process.patPhotons.embedBasicClusters =
False 83 process.patPhotons.embedPreshowerClusters =
False 84 process.patPhotons.embedRecHits =
False 87 process.patPhotons.addPuppiIsolation = cms.bool(
True)
88 process.patPhotons.puppiIsolationChargedHadrons = cms.InputTag(
"egmPhotonPUPPIIsolation",
"h+-DR030-")
89 process.patPhotons.puppiIsolationNeutralHadrons = cms.InputTag(
"egmPhotonPUPPIIsolation",
"h0-DR030-")
90 process.patPhotons.puppiIsolationPhotons = cms.InputTag(
"egmPhotonPUPPIIsolation",
"gamma-DR030-")
92 process.patPhotons.photonSource = cms.InputTag(
"reducedEgamma",
"reducedGedPhotons")
93 process.patPhotons.electronSource = cms.InputTag(
"reducedEgamma",
"reducedGedGsfElectrons")
95 process.phPFIsoDepositChargedPAT.src = cms.InputTag(
"reducedEgamma",
"reducedGedPhotons")
96 process.phPFIsoDepositChargedAllPAT.src = cms.InputTag(
"reducedEgamma",
"reducedGedPhotons")
97 process.phPFIsoDepositNeutralPAT.src = cms.InputTag(
"reducedEgamma",
"reducedGedPhotons")
98 process.phPFIsoDepositGammaPAT.src = cms.InputTag(
"reducedEgamma",
"reducedGedPhotons")
99 process.phPFIsoDepositPUPAT.src = cms.InputTag(
"reducedEgamma",
"reducedGedPhotons")
101 process.patOOTPhotons.photonSource = cms.InputTag(
"reducedEgamma",
"reducedOOTPhotons")
102 process.patOOTPhotons.electronSource = cms.InputTag(
"reducedEgamma",
"reducedGedGsfElectrons")
104 process.selectedPatJets.cut = cms.string(
"pt > 15")
105 process.selectedPatMuons.cut = cms.string(
"pt > 5 || isPFMuon || (pt > 3 && (isGlobalMuon || isStandAloneMuon || numberOfMatches > 0 || muonID('RPCMuLoose')))")
107 from Configuration.Eras.Modifier_phase2_muon_cff
import phase2_muon
108 phase2_muon.toModify(process.selectedPatMuons, cut =
"pt > 5 || isPFMuon || (pt > 3 && (isGlobalMuon || isStandAloneMuon || numberOfMatches > 0 || muonID('RPCMuLoose') || muonID('ME0MuonArbitrated') || muonID('GEMMuonArbitrated')) )")
109 from Configuration.ProcessModifiers.pp_on_AA_cff
import pp_on_AA
110 pp_on_AA.toModify(process.selectedPatMuons, cut =
"pt > 5 || isPFMuon || (pt > 1.2 && (isGlobalMuon || isStandAloneMuon) )")
111 from Configuration.Eras.Modifier_run3_upc_cff
import run3_upc
112 run3_upc.toModify(process.selectedPatMuons, cut =
"")
113 from Configuration.Eras.Modifier_ppRef_2024_cff
import ppRef_2024
114 ppRef_2024.toModify(process.selectedPatMuons, cut =
"pt > 5 || isPFMuon || isTrackerMuon || (pt > 1.2 && (isGlobalMuon || isStandAloneMuon))")
116 process.selectedPatElectrons.cut = cms.string(
"")
117 process.selectedPatTaus.cut = cms.string(
"pt > 18. && tauID('decayModeFindingNewDMs')> 0.5")
118 process.selectedPatPhotons.cut = cms.string(
"")
120 _dummyPatJets = process.selectedPatJets.clone(cut =
"pt < 0")
123 def _applySubstructure(process):
126 from PhysicsTools.PatAlgos.slimming.applySubstructure_cff
import applySubstructure
128 (~pp_on_AA).toModify(process, _applySubstructure)
130 pp_on_AA.toModify(process, func =
lambda p:
addToProcessAndTask(
'slimmedJets', p.selectedPatJets.clone(), p, task))
131 pp_on_AA.toModify(process, func =
lambda p:
addToProcessAndTask(
'slimmedJetsAK8', _dummyPatJets.clone(), p, task))
134 from PhysicsTools.PatAlgos.tools.trigTools
import switchOnTriggerStandAlone
136 process.patTrigger.packTriggerPathNames = cms.bool(
True)
141 from PhysicsTools.PatUtils.tools.runMETCorrectionsAndUncertainties
import runMetCorAndUncForMiniAODProduction
143 jetCollUnskimmed=
"patJets")
146 from PhysicsTools.PatAlgos.tools.metTools
import addMETCollection
147 addMETCollection(process,
148 labelName =
"patCaloMet",
149 metSource =
"caloMetM" 154 process.noHFCands = cms.EDFilter(
"GenericPFCandidateSelector",
155 src=cms.InputTag(
"particleFlow"),
156 cut=cms.string(
"abs(pdgId)!=1 && abs(pdgId)!=2 && abs(eta)<3.0")
158 task.add(process.noHFCands)
161 pfCandColl=cms.InputTag(
"noHFCands"),
162 recoMetFromPFCs=
True,
163 jetSelection=
"pt>15 && abs(eta)<3.",
167 process.load(
'PhysicsTools.PatAlgos.slimming.slimmedMETs_cfi')
168 task.add(process.slimmedMETs)
169 (~pp_on_AA).toModify(process.slimmedMETs, addDeepMETs =
True)
174 process.load(
"CommonTools.ParticleFlow.pfCHS_cff")
175 task.add(process.pfCHS)
178 process.pfMetCHS = pfMet.clone(src =
'pfCHS')
179 task.add(process.pfMetCHS)
181 addMETCollection(process,
182 labelName =
"patCHSMet",
183 metSource =
"pfMetCHS" 186 process.patCHSMet.computeMETSignificance = cms.bool(
False)
191 process.TrkCands = chargedPackedCandsForTkMet.clone()
192 task.add(process.TrkCands)
194 process.pfMetTrk = pfMet.clone(src =
'TrkCands')
195 task.add(process.pfMetTrk)
197 addMETCollection(process,
198 labelName =
"patTrkMet",
199 metSource =
"pfMetTrk" 202 process.patTrkMet.computeMETSignificance = cms.bool(
False)
207 process.load(
'RecoJets.JetProducers.hfJetShowerShape_cfi')
208 task.add(process.hfJetShowerShape)
210 process.patJets.userData.userFloats.src += [
'hfJetShowerShape:sigmaEtaEta',
'hfJetShowerShape:sigmaPhiPhi']
211 process.patJets.userData.userInts.src += [
'hfJetShowerShape:centralEtaStripSize',
'hfJetShowerShape:adjacentEtaStripsSize']
214 def _add_deepFlavour(process):
215 process.load(
'RecoBTag.Combined.deepFlavour_cff')
216 task.add(process.pfDeepCSVDiscriminatorsJetTags)
217 process.patJets.discriminatorSources.extend([
218 'pfDeepCSVDiscriminatorsJetTags:BvsAll',
219 'pfDeepCSVDiscriminatorsJetTags:CvsB',
220 'pfDeepCSVDiscriminatorsJetTags:CvsL',
224 process.caloJetMap = cms.EDProducer(
"RecoJetDeltaRValueMapProducer",
225 src = process.patJets.jetSource,
226 matched = cms.InputTag(
"ak4CaloJets"),
227 distMax = cms.double(0.4),
228 values = cms.vstring(
'pt',
'emEnergyFraction'),
229 valueLabels = cms.vstring(
'pt',
'emEnergyFraction'),
230 lazyParser = cms.bool(
True) )
231 task.add(process.caloJetMap)
232 process.patJets.userData.userFloats.src += [
'caloJetMap:pt',
'caloJetMap:emEnergyFraction' ]
234 pp_on_AA.toModify(process.patJets.userData.userInts, src = [] )
235 pp_on_AA.toModify(process.patJets.userData.userFloats, src = [] )
238 from PhysicsTools.PatAlgos.slimming.muonIsolationsPUPPI_cfi
import makeInputForPUPPIIsolationMuon
242 from PhysicsTools.PatAlgos.slimming.egmIsolationsPUPPI_cfi
import makeInputForPUPPIIsolationEgm
244 from RecoEgamma.EgammaTools.egammaObjectModificationsInMiniAOD_cff
import egamma_modifications
245 process.slimmedElectrons.modifierConfig.modifications = egamma_modifications
246 process.slimmedPhotons.modifierConfig.modifications = egamma_modifications
249 process.patElectrons.addElectronID = cms.bool(
True)
250 electron_ids = [
'RecoEgamma.ElectronIdentification.Identification.heepElectronID_HEEPV70_cff',
251 'RecoEgamma.ElectronIdentification.Identification.heepElectronID_HEEPV71_cff',
252 'RecoEgamma.ElectronIdentification.Identification.cutBasedElectronID_Fall17_94X_V1_cff',
253 'RecoEgamma.ElectronIdentification.Identification.cutBasedElectronID_Fall17_94X_V2_cff',
254 'RecoEgamma.ElectronIdentification.Identification.cutBasedElectronID_Winter22_122X_V1_cff',
255 'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_Fall17_noIso_V1_cff',
256 'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_Fall17_iso_V1_cff',
257 'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_Fall17_noIso_V2_cff',
258 'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_Fall17_iso_V2_cff',
259 'RecoEgamma.ElectronIdentification.Identification.cutBasedElectronID_Summer16_80X_V1_cff',
260 'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_Spring16_GeneralPurpose_V1_cff',
261 'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_Spring16_HZZ_V1_cff',
262 'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_Summer16UL_ID_ISO_cff',
263 'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_Summer17UL_ID_ISO_cff',
264 'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_Summer18UL_ID_ISO_cff',
265 'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_RunIIIWinter22_noIso_V1_cff',
266 'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_RunIIIWinter22_iso_V1_cff',
267 'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_Winter22_HZZ_V1_cff' 270 process.egmGsfElectronIDs.physicsObjectSrc = cms.InputTag(
"reducedEgamma",
"reducedGedGsfElectrons")
271 process.electronMVAValueMapProducer.src = cms.InputTag(
'reducedEgamma',
'reducedGedGsfElectrons')
276 process.load(
"RecoEgamma.EgammaTools.gedGsfElectronsTo106X_cff")
277 from Configuration.Eras.Modifier_pp_on_AA_2018_cff
import pp_on_AA_2018
278 pp_on_AA_2018.toModify(task, func=
lambda t: t.add(process.gedGsfElectronsFrom94XTo106XTask))
279 pp_on_AA_2018.toModify(process.electronMVAValueMapProducer,
280 keysForValueMaps = cms.InputTag(
'reducedEgamma',
'reducedGedGsfElectrons'),
281 src =
"gedGsfElectronsFrom94XTo106X")
283 for idmod
in electron_ids:
287 process.patPhotons.addPhotonID = cms.bool(
True)
288 photon_ids = [
'RecoEgamma.PhotonIdentification.Identification.cutBasedPhotonID_Fall17_94X_V1_TrueVtx_cff',
289 'RecoEgamma.PhotonIdentification.Identification.cutBasedPhotonID_Fall17_94X_V2_cff',
290 'RecoEgamma.PhotonIdentification.Identification.mvaPhotonID_Fall17_94X_V1p1_cff',
291 'RecoEgamma.PhotonIdentification.Identification.mvaPhotonID_Fall17_94X_V2_cff',
292 'RecoEgamma.PhotonIdentification.Identification.cutBasedPhotonID_Spring16_V2p2_cff',
293 'RecoEgamma.PhotonIdentification.Identification.mvaPhotonID_Spring16_nonTrig_V1_cff',
294 'RecoEgamma.PhotonIdentification.Identification.mvaPhotonID_Winter22_122X_V1_cff',
295 'RecoEgamma.PhotonIdentification.Identification.cutBasedPhotonID_RunIIIWinter22_122X_V1_cff']
298 process.egmPhotonIDs.physicsObjectSrc = cms.InputTag(
"reducedEgamma",
"reducedGedPhotons")
299 process.photonMVAValueMapProducer.src = cms.InputTag(
'reducedEgamma',
'reducedGedPhotons')
300 for idmod
in photon_ids:
304 from RecoEgamma.EgammaTools.egammaObjectModifications_tools
import makeVIDBitsModifier
305 egamma_modifications.append(
makeVIDBitsModifier(process,
"egmGsfElectronIDs",
"egmPhotonIDs"))
308 from RecoTauTag.Configuration.boostedHPSPFTaus_cfi
import addBoostedTaus
311 import RecoTauTag.RecoTau.tools.runTauIdMVA
as tauIdConfig
312 _updatedBoostedTauName =
'slimmedTausBoostedDeepIDs' 313 _noUpdatedBoostedTauName =
'slimmedTausBoostedNoDeepIDs' 314 boostedTauIdEmbedder = tauIdConfig.TauIDEmbedder(
315 process, debug =
False,
316 originalTauName = _noUpdatedBoostedTauName,
317 updatedTauName = _updatedBoostedTauName,
318 postfix =
'BoostedForMini',
319 toKeep = [
'boostedDeepTauRunIIv2p0']
321 boostedTauIdEmbedder.runTauID()
322 addToProcessAndTask(_noUpdatedBoostedTauName, process.slimmedTausBoosted.clone(),process,task)
323 delattr(process,
'slimmedTausBoosted')
324 process.slimmedTausBoosted = getattr(process, _updatedBoostedTauName).
clone()
325 process.rerunMvaIsolationTaskBoostedForMini.add(process.slimmedTausBoosted)
326 task.add(process.rerunMvaIsolationTaskBoostedForMini)
327 process.load(
"RecoTauTag.Configuration.RecoPFTauTag_cff")
328 process.load(
"RecoTauTag.Configuration.HPSPFTaus_cff")
330 _makePatTausTaskWithRetrainedMVATauID = process.makePatTausTask.copy()
331 _makePatTausTaskWithRetrainedMVATauID.add(process.hpsPFTauBasicDiscriminatorsTask,
332 process.hpsPFTauDiscriminationByIsolationMVArun2v1DBoldDMwLTTask,
333 process.hpsPFTauDiscriminationByIsolationMVArun2v1DBnewDMwLTTask,
334 process.hpsPFTauBasicDiscriminatorsdR03Task,
335 process.hpsPFTauDiscriminationByIsolationMVArun2v1DBdR03oldDMwLTTask,
336 process.hpsPFTauDiscriminationByMVA6rawElectronRejection,
337 process.hpsPFTauDiscriminationByMVA6ElectronRejection,
338 process.hpsPFTauDiscriminationByMuonRejection3)
339 from Configuration.ProcessModifiers.run2_miniAOD_UL_cff
import run2_miniAOD_UL
340 (run2_miniAOD_UL).toReplaceWith(
341 process.makePatTausTask, _makePatTausTaskWithRetrainedMVATauID
344 _updatedTauName =
'slimmedTausDeepIDs' 345 _noUpdatedTauName =
'slimmedTausNoDeepIDs' 346 tauIdEmbedder = tauIdConfig.TauIDEmbedder(
347 process, debug =
False,
348 originalTauName = _noUpdatedTauName,
349 updatedTauName = _updatedTauName,
351 toKeep = [
'deepTau2018v2p5']
353 from Configuration.Eras.Modifier_phase2_common_cff
import phase2_common
354 _tauIds_phase2 = [
'deepTau2026v2p5']
355 phase2_common.toModify(tauIdEmbedder.toKeep, func=
lambda t:t.extend(_tauIds_phase2))
356 tauIdEmbedder.runTauID()
358 delattr(process,
'slimmedTaus')
359 process.slimmedTaus = getattr(process, _updatedTauName).
clone()
360 process.rerunMvaIsolationTaskForMini.add(process.slimmedTaus)
361 task.add(process.rerunMvaIsolationTaskForMini)
365 _makePatTausTaskWithDeadECalVeto = process.makePatTausTask.copy()
366 _makePatTausTaskWithDeadECalVeto.add(
367 process.hpsPFTauDiscriminationByDeadECALElectronRejection
369 _run2_miniAOD_ANY = (run2_miniAOD_UL)
370 _run2_miniAOD_ANY.toReplaceWith(
371 process.makePatTausTask, _makePatTausTaskWithDeadECalVeto
375 _makePatTausTaskWithTauReReco = process.makePatTausTask.copy()
376 _makePatTausTaskWithTauReReco.add(process.PFTauTask)
377 (pp_on_AA).toReplaceWith(
378 process.makePatTausTask, _makePatTausTaskWithTauReReco
382 process.load(
'CommonTools.PileupAlgos.Puppi_cff')
383 process.load(
'RecoJets.JetProducers.ak4PFJets_cfi')
384 from Configuration.Eras.Modifier_pA_2016_cff
import pA_2016
385 _rerun_puppijets_task = task.copy()
386 _rerun_puppijets_task.add(process.puppi, process.ak4PFJetsPuppi)
387 (_run2_miniAOD_ANY | pA_2016 | pp_on_AA).toReplaceWith(task, _rerun_puppijets_task)
390 process.ak4PFJetsPuppiTracksAssociatorAtVertex = cms.EDProducer(
"JetTracksAssociatorAtVertex",
392 jets = cms.InputTag(
"ak4PFJetsPuppi")
394 task.add(process.ak4PFJetsPuppiTracksAssociatorAtVertex)
395 process.patJetPuppiCharge = cms.EDProducer(
"JetChargeProducer",
396 src = cms.InputTag(
"ak4PFJetsPuppiTracksAssociatorAtVertex"),
397 var = cms.string(
'Pt'),
398 exp = cms.double(1.0)
400 task.add(process.patJetPuppiCharge)
403 process.load(
"RecoJets.JetProducers.PileupJetID_cfi")
404 task.add(process.pileUpJetIDPuppiTask)
406 def _add_jetsPuppi(process):
408 noDeepFlavourDiscriminators = [x.value()
if isinstance(x, cms.InputTag)
else x
for x
in process.patJets.discriminatorSources
409 if (
not "DeepFlavour" in str(x)
and not "Transformer" in str(x))]
410 addJetCollection(process, postfix =
"", labelName =
'Puppi', jetSource = cms.InputTag(
'ak4PFJetsPuppi'),
411 jetCorrections = (
'AK4PFPuppi', [
'L2Relative',
'L3Absolute'],
''),
412 pfCandidates = cms.InputTag(
"particleFlow"),
413 algo=
'AK', rParam = 0.4, btagDiscriminators = noDeepFlavourDiscriminators
416 process.patJetGenJetMatchPuppi.matched =
'slimmedGenJets' 418 process.patJetsPuppi.jetChargeSource = cms.InputTag(
"patJetPuppiCharge")
421 process.patJetsPuppi.userData.userFloats.src += [cms.InputTag(
"pileupJetIdPuppi:fullDiscriminant")]
422 process.patJetsPuppi.userData.userInts.src += [cms.InputTag(
"pileupJetIdPuppi:fullId")]
424 process.selectedPatJetsPuppi.cut = cms.string(
"pt > 10")
426 from PhysicsTools.PatAlgos.slimming.applyDeepBtagging_cff
import applyDeepBtagging
429 process.slimmedJetsNoDeepFlavour.dropTagInfos =
'0' 430 process.updatedPatJetsTransientCorrectedSlimmedDeepFlavour.addTagInfos =
True 431 process.updatedPatJetsTransientCorrectedSlimmedDeepFlavour.tagInfoSources = [
"pixelClusterTagInfos"]
432 _run2_miniAOD_ANY.toModify(process.updatedPatJetsTransientCorrectedSlimmedDeepFlavour, addTagInfos =
False )
433 (~pp_on_AA).toModify(process, _add_jetsPuppi)
435 pp_on_AA.toModify(process, func =
lambda p:
addToProcessAndTask(
'slimmedJetsPuppi', _dummyPatJets.clone(), p, task))
438 process.patJets.addTagInfos =
True 439 process.patJets.tagInfoSources = [
"pixelClusterTagInfos"]
441 _run2_miniAOD_ANY.toModify(process.patJets, addTagInfos =
False )
443 from Configuration.Eras.Modifier_pp_on_AA_2018_cff
import pp_on_AA_2018
444 pp_on_AA_2018.toModify(process.patJets, tagInfoSources = cms.VInputTag([
"impactParameterTagInfos",
"secondaryVertexTagInfos"]) )
447 def _add_metPuppi(process):
448 process.load(
'RecoMET.METProducers.pfMetPuppi_cfi')
449 _rerun_puppimet_task = task.copy()
450 _rerun_puppimet_task.add(process.puppiNoLep, process.pfMetPuppi)
451 (_run2_miniAOD_ANY | pA_2016 | pp_on_AA ).toReplaceWith(task, _rerun_puppimet_task)
454 jetCollUnskimmed=
"slimmedJetsPuppi",
455 recoMetFromPFCs=
True,
456 jetFlavor=
"AK4PFPuppi",
459 (~pp_on_AA).toModify(process, _add_metPuppi)
461 process.load(
'PhysicsTools.PatAlgos.slimming.slimmedMETs_cfi')
462 task.add(process.slimmedMETs)
464 def _add_slimmedMETsPuppi(process):
466 process.slimmedMETsPuppi.src = cms.InputTag(
"patMETsPuppi")
467 process.slimmedMETsPuppi.rawVariation = cms.InputTag(
"patPFMetPuppi")
468 process.slimmedMETsPuppi.t1Uncertainties = cms.InputTag(
"patPFMetT1%sPuppi")
469 process.slimmedMETsPuppi.t01Variation = cms.InputTag(
"patPFMetT0pcT1Puppi")
470 process.slimmedMETsPuppi.t1SmearedVarsAndUncs = cms.InputTag(
"patPFMetT1Smear%sPuppi")
471 process.slimmedMETsPuppi.tXYUncForRaw = cms.InputTag(
"patPFMetTxyPuppi")
472 process.slimmedMETsPuppi.tXYUncForT1 = cms.InputTag(
"patPFMetT1TxyPuppi")
473 process.slimmedMETsPuppi.tXYUncForT01 = cms.InputTag(
"patPFMetT0pcT1TxyPuppi")
474 process.slimmedMETsPuppi.tXYUncForT1Smear = cms.InputTag(
"patPFMetT1SmearTxyPuppi")
475 process.slimmedMETsPuppi.tXYUncForT01Smear = cms.InputTag(
"patPFMetT0pcT1SmearTxyPuppi")
476 del process.slimmedMETsPuppi.caloMET
477 (~pp_on_AA).toModify(process, _add_slimmedMETsPuppi)
479 def _add_deepMET(process):
480 from RecoMET.METPUSubtraction.deepMETProducer_cff
import deepMETsResolutionTune, deepMETsResponseTune
484 (~pp_on_AA).toModify(process, _add_deepMET)
487 process.load(
"TrackingTools.TrackAssociator.DetIdAssociatorESProducer_cff")
491 from Configuration.Eras.Modifier_phase2_hgcal_cff
import phase2_hgcal
492 process.load(
"RecoEgamma.EgammaTools.slimmedEgammaHGC_cff")
493 phase2_hgcal.toModify(task, func=
lambda t: t.add(process.slimmedEgammaHGCTask))
497 def _addUTagToTaus(process, task,
498 storePNetCHSjets = False,
499 storeUParTPUPPIjets = False,
501 if not (storePNetCHSjets
or storeUParTPUPPIjets):
return process
502 noUpdatedTauName =
'slimmedTausNoUTag' 505 process.pfParticleNetFromMiniAODAK4CHSCentralTagInfosSlimmedDeepFlavour.taus = noUpdatedTauName
506 process.pfParticleNetFromMiniAODAK4PuppiCentralTagInfosSlimmedPuppiWithDeepTags.taus = noUpdatedTauName
507 process.pfParticleNetFromMiniAODAK4PuppiForwardTagInfosSlimmedPuppiWithDeepTags.taus = noUpdatedTauName
508 from PhysicsTools.PatAlgos.patTauHybridProducer_cfi
import patTauHybridProducer
510 jetCollection =
'slimmedJets' 511 TagName =
'pfParticleNetFromMiniAODAK4CHSCentralJetTags' 512 tag_prefix =
'byUTagCHS' 513 updatedTauName =
'slimmedTausWithUTagCHS' 515 from RecoBTag.ONNXRuntime.pfParticleNetFromMiniAODAK4_cff
import pfParticleNetFromMiniAODAK4CHSCentralJetTags
516 Discriminators = [TagName+
":"+tag
for tag
in pfParticleNetFromMiniAODAK4CHSCentralJetTags.flav_names.value()]
518 setattr(process, updatedTauName,
519 patTauHybridProducer.clone(
520 src = noUpdatedTauName,
521 jetSource = jetCollection,
526 UTagScoreNames = Discriminators,
527 tagPrefix = tag_prefix,
530 checkTauScoreIsBest =
False,
531 chargeAssignmentProbMin = 0.2,
532 addGenJetMatch = addGenJet,
537 process.tauGenJetMatch.clone(src = jetCollection),
539 getattr(process,updatedTauName).genJetMatch =
'tauGenJetMatchCHSJet' 540 if storeUParTPUPPIjets: task.add(getattr(process,updatedTauName))
541 if storeUParTPUPPIjets:
542 jetCollection =
'slimmedJetsPuppi' 543 TagName =
'pfUnifiedParticleTransformerAK4JetTags' 544 tag_prefix =
'byUTagPUPPI' 545 updatedTauName =
'slimmedTausWithUTagPUPPI' if not storePNetCHSjets
else 'slimmedTausWithUTagCHSAndUTagPUPPI' 547 from RecoBTag.ONNXRuntime.pfUnifiedParticleTransformerAK4JetTags_cfi
import pfUnifiedParticleTransformerAK4JetTags
548 Discriminators = [TagName+
":"+tag
for tag
in pfUnifiedParticleTransformerAK4JetTags.flav_names.value()]
550 setattr(process, updatedTauName,
551 patTauHybridProducer.clone(
552 src = noUpdatedTauName
if not storePNetCHSjets
else 'slimmedTausWithUTagCHS',
553 jetSource = jetCollection,
558 UTagScoreNames = Discriminators,
559 tagPrefix = tag_prefix,
562 checkTauScoreIsBest =
False,
563 chargeAssignmentProbMin = 0.2,
564 addGenJetMatch = addGenJet,
569 process.tauGenJetMatch.clone(src = jetCollection),
571 getattr(process,updatedTauName).genJetMatch =
'tauGenJetMatchPUPPIJet' 573 delattr(process,
'slimmedTaus')
574 process.slimmedTaus = getattr(process, updatedTauName).
clone()
575 task.add(process.slimmedTaus)
578 _uTagToTaus_switches = cms.PSet(
579 storePNetCHSjets = cms.bool(
True),
580 storeUParTPUPPIjets = cms.bool(
True),
581 addGenJet = cms.bool(
True)
583 pp_on_AA.toModify(_uTagToTaus_switches, storePNetCHSjets=
False, storeUParTPUPPIjets=
False)
584 _addUTagToTaus(process, task,
585 storePNetCHSjets = _uTagToTaus_switches.storePNetCHSjets.value(),
586 storeUParTPUPPIjets = _uTagToTaus_switches.storePNetCHSjets.value(),
587 addGenJet = _uTagToTaus_switches.addGenJet.value()
591 from Configuration.Eras.Modifier_run2_L1prefiring_cff
import run2_L1prefiring
592 from Configuration.Eras.Modifier_stage2L1Trigger_cff
import stage2L1Trigger
593 from Configuration.Eras.Modifier_stage2L1Trigger_2017_cff
import stage2L1Trigger_2017
594 from Configuration.Eras.Modifier_stage2L1Trigger_2018_cff
import stage2L1Trigger_2018
595 from Configuration.Eras.Modifier_tracker_apv_vfp30_2016_cff
import tracker_apv_vfp30_2016
596 if not hasattr(process,
'prefiringweight'):
597 process.load(
"PhysicsTools.PatUtils.L1PrefiringWeightProducer_cff")
598 (stage2L1Trigger & tracker_apv_vfp30_2016).toModify(process.prefiringweight, DataEraECAL =
"UL2016preVFP", DataEraMuon =
"2016preVFP" )
599 (stage2L1Trigger & ~tracker_apv_vfp30_2016).toModify(process.prefiringweight, DataEraECAL =
"UL2016postVFP", DataEraMuon =
"2016postVFP" )
600 stage2L1Trigger_2017.toModify(process.prefiringweight, DataEraECAL =
"UL2017BtoF", DataEraMuon =
"20172018")
601 stage2L1Trigger_2018.toModify(process.prefiringweight, DataEraECAL =
"None", DataEraMuon =
"20172018")
602 run2_L1prefiring.toModify(task, func=
lambda t: t.add(process.prefiringweight))
604 from PhysicsTools.PatAlgos.producersHeavyIons.heavyIonJetSetup
import removeL1FastJetJECs
605 pp_on_AA.toModify(process, removeL1FastJetJECs)
def makeInputForPUPPIIsolationMuon(process)
def applySubstructure(process, postfix="")
def addToProcessAndTask(label, module, process, task)
def makeInputForPUPPIIsolationEgm(process)
def runMetCorAndUncForMiniAODProduction(process, metType="PF", jetCollUnskimmed="patJets", photonColl="selectedPatPhotons", electronColl="selectedPatElectrons", muonColl="selectedPatMuons", tauColl="selectedPatTaus", pfCandColl="particleFlow", jetCleaning="LepClean", jetSelection="pt>15 && abs(eta)<9.9", jecUnFile="", jetFlavor="AK4PFchs", recoMetFromPFCs=False, postfix="")
def applyDeepBtagging(process, postfix="")
TEveGeoShape * clone(const TEveElement *element, TEveElement *parent)
def getPatAlgosToolsTask(process)
def addBoostedTaus(process)