# input collections
gsfElectronCoresTag="gedGsfElectronCores",
# steering
resetMvaValuesUsingPFCandidates=True,
applyPreselection=True,
ecalDrivenEcalEnergyFromClassBasedParameterization=False,
ecalDrivenEcalErrorFromClassBasedParameterization=False,
useEcalRegression=True,
useCombinationRegression=True,
# regression. The labels are needed in all cases.
ecalRefinedRegressionWeightLabels=[
"gedelectron_EBCorrection_offline_v1",
"gedelectron_EECorrection_offline_v1",
"gedelectron_EBUncertainty_offline_v1",
"gedelectron_EEUncertainty_offline_v1"
],
combinationRegressionWeightLabels=["gedelectron_p4combination_offline"],
)
from Configuration.Eras.Modifier_pp_on_AA_2018_cff import pp_on_AA_2018
pp_on_AA_2018.toModify(gedGsfElectronsTmp.preselection, minSCEtBarrel=15.0)
pp_on_AA_2018.toModify(gedGsfElectronsTmp.preselection, minSCEtEndcaps=15.0)
from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
egamma_lowPt_exclusive.toModify(gedGsfElectronsTmp.preselection,
minSCEtBarrel=1.0,
minSCEtEndcaps=1.0)
egamma_lowPt_exclusive.toModify(gedGsfElectronsTmp, applyPreselection=False)
nTrackSolidConeBarrel=cms.double(999999999),
nTrackSolidConeEndcap=cms.double(999999999),
nTrackHollowConeBarrel=cms.double(999999999),
nTrackHollowConeEndcap=cms.double(999999999),
trackPtSumSolidConeBarrel=cms.double(999999999),
trackPtSumSolidConeEndcap=cms.double(999999999),
trackPtSumHollowConeBarrel=cms.double(999999999),
trackPtSumHollowConeEndcap=cms.double(999999999),
sigmaIetaIetaCutBarrel=cms.double(999999999),
sigmaIetaIetaCutEndcap=cms.double(999999999),
posCalcParameters=cms.PSet(T0_barl=cms.double(7.4),
T0_endc=cms.double(6.3),
T0_endcPresh=cms.double(3.6),
LogWeighted=cms.bool(True),
W0=cms.double(4.2),
X0=cms.double(0.89)),
RecHitFlagToBeExcludedEB=cleanedHybridSuperClusters.RecHitFlagToBeExcluded,
RecHitSeverityToBeExcludedEB=cleanedHybridSuperClusters.
RecHitSeverityToBeExcluded,
RecHitFlagToBeExcludedEE=multi5x5BasicClustersCleaned.
RecHitFlagToBeExcluded,
RecHitSeverityToBeExcludedEE=cleanedHybridSuperClusters.
RecHitSeverityToBeExcluded,
)
from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
egamma_lowPt_exclusive.toModify(
photons,
minSCEtBarrel=1.0, #default 10
minSCEtEndcap=1.0) #default 10
_pp_on_AA_reducedEcalRecHitsTask = reducedEcalRecHitsTask.copy()
_pp_on_AA_reducedEcalRecHitsTask.remove(interestingEcalDetIdOOTPFEB)
_pp_on_AA_reducedEcalRecHitsTask.remove(interestingEcalDetIdOOTPFEE)
_pp_on_AA_reducedEcalRecHitsTask.remove(interestingEcalDetIdOOTPFES)
from Configuration.ProcessModifiers.pp_on_AA_cff import pp_on_AA
pp_on_AA.toReplaceWith(reducedEcalRecHitsTask,
_pp_on_AA_reducedEcalRecHitsTask)
pp_on_AA.toModify(
reducedEcalRecHitsEB.interestingDetIdCollections,
func=lambda list: list.remove(cms.InputTag("interestingEcalDetIdOOTPFEB")))
pp_on_AA.toModify(
reducedEcalRecHitsEB.interestingDetIdCollections,
func=lambda list: list.remove(cms.InputTag("interestingOotGamIsoDetIdEB")))
pp_on_AA.toModify(
reducedEcalRecHitsEE.interestingDetIdCollections,
func=lambda list: list.remove(cms.InputTag("interestingEcalDetIdOOTPFEE")))
pp_on_AA.toModify(
reducedEcalRecHitsEE.interestingDetIdCollections,
func=lambda list: list.remove(cms.InputTag("interestingOotGamIsoDetIdEE")))
pp_on_AA.toModify(
reducedEcalRecHitsES.interestingDetIds,
func=lambda list: list.remove(cms.InputTag("interestingEcalDetIdOOTPFES")))
pp_on_AA.toModify(reducedEcalRecHitsES.interestingDetIdsNotToClean,
func=lambda list: list.remove(
cms.InputTag("interestingOotEgammaIsoESDetId")))
from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
egamma_lowPt_exclusive.toModify(reducedEcalRecHitsES, scEtThreshold=1.0)
PreIdLabel=cms.string('preid'),
ProducePreId=cms.untracked.bool(True),
PtThresholdSavePreId=cms.untracked.double(1.0),
Min_dr=cms.double(0.2))
from Configuration.Eras.Modifier_pp_on_XeXe_2017_cff import pp_on_XeXe_2017
from Configuration.ProcessModifiers.pp_on_AA_cff import pp_on_AA
for e in [pp_on_XeXe_2017, pp_on_AA]:
e.toModify(trackerDrivenElectronSeeds, MinPt=5.0)
# tracker driven electron seeds depend on the generalTracks trajectory collection
# However, in FastSim jobs, trajectories are only available for the 'before mixing' track collections
# Therefore we let the seeds depend on the 'before mixing' generalTracks collection
# TODO: investigate whether the dependence on trajectories can be avoided
from Configuration.Eras.Modifier_fastSim_cff import fastSim
trackerDrivenElectronSeedsTmp = trackerDrivenElectronSeeds.clone(
TkColList=["generalTracksBeforeMixing"])
import FastSimulation.Tracking.ElectronSeedTrackRefFix_cfi
_fastSim_trackerDrivenElectronSeeds = FastSimulation.Tracking.ElectronSeedTrackRefFix_cfi.fixedTrackerDrivenElectronSeeds.clone(
)
_fastSim_trackerDrivenElectronSeeds.seedCollection.setModuleLabel(
"trackerDrivenElectronSeedsTmp")
_fastSim_trackerDrivenElectronSeeds.idCollection = [
"trackerDrivenElectronSeedsTmp:preid",
]
fastSim.toReplaceWith(trackerDrivenElectronSeeds,
_fastSim_trackerDrivenElectronSeeds)
from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
egamma_lowPt_exclusive.toModify(trackerDrivenElectronSeeds, MinPt=1.0)
# calibHF_b_EMHAD = cms.vdouble(10.00,10.00,10.00,10.00,10.00,10.00,10.00,10.00,10.00,10.00,10.00)
calibHF_a_EMonly=cms.vdouble(0.96945, 0.96701, 0.76309, 0.82268, 0.87583,
0.89718, 0.98674, 1.4681, 1.4580, 1.4580),
calibHF_b_HADonly=cms.vdouble(1.27541, 0.85361, 0.86333, 0.89091, 0.94348,
0.94348, 0.94370, 1.0034, 1.0444, 1.0444),
calibHF_a_EMHAD=cms.vdouble(1.42215, 1.00496, 0.68961, 0.81656, 0.98504,
0.98504, 1.00802, 1.0593, 1.4576, 1.4576),
calibHF_b_EMHAD=cms.vdouble(1.27541, 0.85361, 0.86333, 0.89091, 0.94348,
0.94348, 0.94370, 1.0034, 1.0444, 1.0444)
# toRead = cms.untracked.vstring("PFfa_BARREL",
# "PFfa_ENDCAP",
# "PFfb_BARREL",
# "PFfb_ENDCAP",
# "PFfc_BARREL",
# "PFfc_ENDCAP",
# "PFfaEta_BARREL",
# "PFfaEta_ENDCAP",
# "PFfbEta_BARREL",
# "PFfbEta_ENDCAP") # same strings as fType
)
from Configuration.Eras.Modifier_pf_badHcalMitigation_cff import pf_badHcalMitigation
pf_badHcalMitigation.toModify(
particleFlowTmp,
electron_protectionsForBadHcal=dict(enableProtections=True),
photon_protectionsForBadHcal=dict(enableProtections=True))
from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
egamma_lowPt_exclusive.toModify(particleFlowTmp, photon_MinEt=1.)
gedPhotonsTmp.outputPhotonCollection = cms.string("")
gedPhotonsTmp.reconstructionStep = cms.string("tmp")
gedPhotonTaskTmp = cms.Task(gedPhotonCore, gedPhotonsTmp)
gedPhotonSequenceTmp = cms.Sequence(gedPhotonTaskTmp)
gedPhotons = RecoEgamma.EgammaPhotonProducers.gedPhotons_cfi.gedPhotons.clone()
gedPhotons.photonProducer = cms.InputTag("gedPhotonsTmp")
gedPhotons.outputPhotonCollection = cms.string("")
gedPhotons.reconstructionStep = cms.string("final")
gedPhotons.pfECALClusIsolation = cms.InputTag("photonEcalPFClusterIsolationProducer")
gedPhotons.pfHCALClusIsolation = cms.InputTag("photonHcalPFClusterIsolationProducer")
gedPhotons.pfIsolCfg = cms.PSet(
chargedHadronIso = cms.InputTag("photonIDValueMaps","phoChargedIsolation"),
neutralHadronIso = cms.InputTag("photonIDValueMaps","phoNeutralHadronIsolation"),
photonIso = cms.InputTag("photonIDValueMaps","phoPhotonIsolation"),
chargedHadronWorstVtxIso = cms.InputTag("photonIDValueMaps","phoWorstChargedIsolation"),
chargedHadronWorstVtxGeomVetoIso = cms.InputTag("photonIDValueMaps","phoWorstChargedIsolationConeVeto"),
chargedHadronPFPVIso = cms.InputTag("egmPhotonIsolationCITK:h+-DR030-"),
)
gedPhotonTask = cms.Task(gedPhotons)
gedPhotonSequence = cms.Sequence(gedPhotonTask)
from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
egamma_lowPt_exclusive.toModify(gedPhotons,
minSCEtBarrel = 1.0,
minSCEtEndcap = 1.0)
egamma_lowPt_exclusive.toModify(gedPhotonsTmp,
minSCEtBarrel = 1.0,
minSCEtEndcap = 1.0)
MinDPhiGsfSC = cms.double(0.15),
MinSCEnergy = cms.double(4.0),
TTRHBuilder = cms.string('WithTrackAngle'),
GsfTrackModuleLabel = cms.InputTag("electronGsfTracks"),
Propagator = cms.string('fwdElectronPropagator'),
PFRecTrackLabel = cms.InputTag("pfTrack"),
PFEcalClusters = cms.InputTag("particleFlowClusterECAL"),
PrimaryVertexLabel = cms.InputTag("offlinePrimaryVertices"),
useConvBremFinder = cms.bool(True),
PFNuclear = cms.InputTag("pfDisplacedTrackerVertex"),
PFConversions = cms.InputTag("pfConversions"),
PFV0 = cms.InputTag("pfV0"),
useNuclear = cms.bool(False),
useV0 = cms.bool(False),
useConversions = cms.bool(False),
debugGsfCleaning = cms.bool(False),
AbsEtaBarrelEndcapsSeparation = cms.double(1.479),
PtLowHighSeparation = cms.double(20),
pf_convBremFinderID_mvaCutBarrelLowPt = cms.double(0.6),
pf_convBremFinderID_mvaCutBarrelHighPt = cms.double(0.97),
pf_convBremFinderID_mvaCutEndcapsLowPt = cms.double(0.9),
pf_convBremFinderID_mvaCutEndcapsHighPt = cms.double(0.995),
pf_convBremFinderID_mvaWeightFileBarrelLowPt = cms.FileInPath('RecoParticleFlow/PFTracking/data/TMVAClassification_ConvBremFinder_Testetlt20absetalt1_479_BDT.weights.xml'),
pf_convBremFinderID_mvaWeightFileBarrelHighPt = cms.FileInPath('RecoParticleFlow/PFTracking/data/TMVAClassification_ConvBremFinder_Testetgt20absetalt1_479_BDT.weights.xml'),
pf_convBremFinderID_mvaWeightFileEndcapsLowPt = cms.FileInPath('RecoParticleFlow/PFTracking/data/TMVAClassification_ConvBremFinder_Testetlt20absetagt1_479_BDT.weights.xml'),
pf_convBremFinderID_mvaWeightFileEndcapsHighPt = cms.FileInPath('RecoParticleFlow/PFTracking/data/TMVAClassification_ConvBremFinder_Testetgt20absetagt1_479_BDT.weights.xml')
)
from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
egamma_lowPt_exclusive.toModify(pfTrackElec,MinSCEnergy = 1.0)
import FWCore.ParameterSet.Config as cms
selectDigi = cms.EDProducer("EcalDigiSelector",
EcalEBDigiTag = cms.InputTag("ecalDigis","ebDigis"),
nclus_sel = cms.int32(2),
EcalEEDigiTag = cms.InputTag("ecalDigis","eeDigis"),
barrelSuperClusterProducer = cms.InputTag('uncleanedHybridSuperClusters'),
endcapSuperClusterProducer = cms.InputTag('correctedMulti5x5SuperClustersWithPreshower'),
EcalEBRecHitTag = cms.InputTag("ecalRecHit","EcalRecHitsEB"),
EcalEERecHitTag = cms.InputTag("ecalRecHit","EcalRecHitsEE"),
cluster_pt_thresh = cms.double(10.0),
single_cluster_thresh = cms.double(15.0),
selectedEcalEBDigiCollection = cms.string('selectedEcalEBDigiCollection'),
selectedEcalEEDigiCollection = cms.string('selectedEcalEEDigiCollection')
)
from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
egamma_lowPt_exclusive.toModify(selectDigi,
cluster_pt_thresh = 1.0,
single_cluster_thresh = 1.0)
'keep recoCaloClusters_hgcalLayerClustersHFNose_*_*',
'keep *_hgcalLayerClustersHFNose_timeLayerCluster_*',
'keep *_hgcalLayerClustersHFNose_InitialLayerClustersMask_*'
])
(pA_2016 | pp_on_AA_2018).toModify(
RecoLocalCaloAOD.outputCommands,
func=lambda outputCommands: outputCommands.extend([
'keep *_zdcreco_*_*', 'keep ZDCDataFramesSorted_hcalDigis_*_*',
'keep ZDCDataFramesSorted_castorDigis_*_*',
'keep QIE10DataFrameHcalDataFrameContainer_hcalDigis_ZDC_*'
]))
from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
egamma_lowPt_exclusive.toModify(
RecoLocalCaloAOD,
outputCommands=RecoLocalCaloAOD.outputCommands + [
'keep *_towerMaker_*_*', 'keep *_zdcreco_*_*',
'keep ZDCDataFramesSorted_hcalDigis_*_*',
'keep ZDCDataFramesSorted_castorDigis_*_*',
'keep QIE10DataFrameHcalDataFrameContainer_hcalDigis_ZDC_*'
])
#RECO content
RecoLocalCaloRECO = cms.PSet(outputCommands=cms.untracked.vstring(
'keep *_hbhereco_*_*', 'keep *_hbheprereco_*_*', 'keep *_hfprereco_*_*',
'keep *_hfreco_*_*', 'keep *_horeco_*_*',
'keep HBHERecHitsSorted_hbherecoMB_*_*',
'keep HORecHitsSorted_horecoMB_*_*', 'keep HFRecHitsSorted_hfrecoMB_*_*',
'keep ZDCDataFramesSorted_hcalDigis_*_*',
'keep ZDCDataFramesSorted_castorDigis_*_*',
'keep QIE10DataFrameHcalDataFrameContainer_hcalDigis_ZDC_*',
'keep ZDCRecHitsSorted_zdcreco_*_*'))
RecoLocalCaloRECO.outputCommands.extend(RecoLocalCaloAOD.outputCommands)
import FWCore.ParameterSet.Config as cms
#AOD content
RecoLocalTrackerAOD = cms.PSet(outputCommands=cms.untracked.vstring(
'keep ClusterSummary_clusterSummaryProducer_*_*'))
from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
egamma_lowPt_exclusive.toModify(
RecoLocalTrackerAOD,
outputCommands=RecoLocalTrackerAOD.outputCommands +
['keep *_siPixelRecHits_*_*', 'keep *_siPixelClusters_*_*'])
#RECO content
RecoLocalTrackerRECO = cms.PSet(outputCommands=cms.untracked.vstring(
'keep DetIdedmEDCollection_siStripDigis_*_*',
'keep DetIdedmEDCollection_siPixelDigis_*_*',
'keep PixelFEDChanneledmNewDetSetVector_siPixelDigis_*_*',
'keep *_siPixelClusters_*_*', 'keep *_siStripClusters_*_*'))
RecoLocalTrackerRECO.outputCommands.extend(RecoLocalTrackerAOD.outputCommands)
from Configuration.Eras.Modifier_phase2_tracker_cff import phase2_tracker
phase2_tracker.toModify(RecoLocalTrackerRECO,
outputCommands=RecoLocalTrackerRECO.outputCommands +
['keep *_siPhase2Clusters_*_*'])
#Full Event content
RecoLocalTrackerFEVT = cms.PSet(
outputCommands=cms.untracked.vstring('keep *_clusterSummaryProducer_*_*'))
RecoLocalTrackerFEVT.outputCommands.extend(RecoLocalTrackerRECO.outputCommands)
phase2_tracker.toModify(RecoLocalTrackerFEVT,
outputCommands=RecoLocalTrackerFEVT.outputCommands +
['keep *_siPhase2Clusters_*_*'])
ePhiMin1 = cms.double(-0.125),
ePhiMax1 = cms.double(0.075),
pPhiMin1 = cms.double(-0.075),
pPhiMax1 = cms.double(0.125),
PhiMin2B = cms.double(-0.002), ## barrel
PhiMax2B = cms.double(0.002), ## barrel
PhiMin2F = cms.double(-0.003), ## forward
PhiMax2F = cms.double(0.003), ## forward
)
)
from Configuration.Eras.Modifier_pp_on_AA_2018_cff import pp_on_AA_2018
pp_on_AA_2018.toModify(ecalDrivenElectronSeeds.SeedConfiguration, SCEtCut = 15.0)
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toModify(
ecalDrivenElectronSeeds,
endcapSuperClusters = 'particleFlowSuperClusterHGCal',
SeedConfiguration = dict( allowHGCal = True )
)
from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
egamma_lowPt_exclusive.toModify(ecalDrivenElectronSeeds.SeedConfiguration,
LowPtThreshold =1.0,
applyHOverECut = False)
# create ecal driven seeds for electron using HGCal Multiclusters
ecalDrivenElectronSeedsFromMultiCl = ecalDrivenElectronSeeds.clone(
endcapSuperClusters = 'particleFlowSuperClusterHGCalFromMultiCl')
from RecoEgamma.PhotonIdentification.isolationCalculator_cfi import *
#
# producer for photonCore collection
#
photonCore = cms.EDProducer(
"PhotonCoreProducer",
conversionProducer=cms.InputTag(""),
# conversionCollection = cms.string(''),
scHybridBarrelProducer=cms.InputTag("correctedHybridSuperClusters"),
scIslandEndcapProducer=cms.InputTag(
"correctedMulti5x5SuperClustersWithPreshower"),
photonCoreCollection=cms.string(''),
pixelSeedProducer=cms.InputTag('electronMergedSeeds'),
minSCEt=cms.double(10.0),
risolveConversionAmbiguity=cms.bool(True),
# MVA_weights_location = cms.string('RecoEgamma/EgammaTools/data/TMVAnalysis_Likelihood.weights.txt')
)
photonCoreFromMultiCl = photonCore.clone(
scIslandEndcapProducer='particleFlowSuperClusterHGCalFromMultiCl',
pixelSeedProducer='electronMergedSeedsFromMultiCl')
islandPhotonCore = photonCore.clone(
scHybridBarrelProducer="correctedIslandBarrelSuperClusters",
scIslandEndcapProducer="correctedIslandEndcapSuperClusters",
minSCEt=8.0)
from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
egamma_lowPt_exclusive.toModify(photonCore, minSCEt=0) #
egamma_lowPt_exclusive.toModify(islandPhotonCore, minSCEt=1.0) #default 8
interestingEgammaIsoDetIds = cms.Sequence(interestingEgammaIsoDetIdsTask)
_pp_on_AA_interestingEgammaIsoDetIdsTask = interestingEgammaIsoDetIdsTask.copy(
)
_pp_on_AA_interestingEgammaIsoDetIdsTask.remove(interestingOotGamIsoDetIdEB)
_pp_on_AA_interestingEgammaIsoDetIdsTask.remove(interestingOotGamIsoDetIdEE)
_pp_on_AA_interestingEgammaIsoDetIdsTask.remove(
interestingOotEgammaIsoHCALDetId)
_pp_on_AA_interestingEgammaIsoDetIdsTask.remove(interestingOotEgammaIsoESDetId)
from Configuration.ProcessModifiers.pp_on_AA_cff import pp_on_AA
pp_on_AA.toReplaceWith(interestingEgammaIsoDetIdsTask,
_pp_on_AA_interestingEgammaIsoDetIdsTask)
from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
egamma_lowPt_exclusive.toModify(
interestingGedEgammaIsoESDetId,
minSCEt=1.0, #default 500
minEleEt=1.0, #default 20
minPhoEt=1.0 #default 20
)
egamma_lowPt_exclusive.toModify(
interestingGedEgammaIsoHCALDetId,
minSCEt=1.0, #default 20
minEleEt=1.0, #default 20
minPhoEt=1.0 #default 20
)
from Configuration.Eras.Modifier_run3_HB_cff import run3_HB
run3_HB.toModify(interestingEgammaIsoHCALSel, minEnergyHB=0.1)
import FWCore.ParameterSet.Config as cms
from RecoEcal.EgammaClusterProducers.particleFlowSuperClusterECALMustache_cfi import particleFlowSuperClusterECALMustache as _particleFlowSuperClusterECALMustache
# define the default ECAL clustering (Mustache or Box)
particleFlowSuperClusterECAL = _particleFlowSuperClusterECALMustache.clone()
from Configuration.Eras.Modifier_pp_on_AA_2018_cff import pp_on_AA_2018
pp_on_AA_2018.toModify(particleFlowSuperClusterECAL,
useDynamicDPhiWindow=False)
pp_on_AA_2018.toModify(particleFlowSuperClusterECAL,
phiwidth_SuperClusterBarrel=0.20)
pp_on_AA_2018.toModify(particleFlowSuperClusterECAL,
phiwidth_SuperClusterEndcap=0.20)
from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
egamma_lowPt_exclusive.toModify(particleFlowSuperClusterECAL,
thresh_SCEt=1.0,
thresh_PFClusterSeedBarrel=0.5,
thresh_PFClusterSeedEndcap=0.5)
linkType = cms.string("HCAL:BREM"),
useKDTree = cms.bool(False) ),
cms.PSet( linkerName = cms.string("SCAndECALLinker"),
linkType = cms.string("SC:ECAL"),
useKDTree = cms.bool(False),
SuperClusterMatchByRef = cms.bool(True) )
)
)
for imp in particleFlowBlock.elementImporters:
if imp.importerName.value() == "SuperClusterImporter":
_scImporter = imp
from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
egamma_lowPt_exclusive.toModify(_scImporter,
minSuperClusterPt = 1.0,
minPTforBypass = 0.0)
def _findIndicesByModule(name):
ret = []
for i, pset in enumerate(particleFlowBlock.elementImporters):
if pset.importerName.value() == name:
ret.append(i)
return ret
#
# kill tracks in the HGCal
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
_insertTrackImportersWithVeto = {}
_trackImporters = ['GeneralTracksImporter','ConvBremTrackImporter',
'ConversionTrackImporter','NuclearInteractionTrackImporter']
electronDnnHighPtEndcapThr=cms.double(0.056),
electronDnnLowPtThr=cms.double(0.075),
electronDnnExtEta1Thr=cms.double(0.113604),
electronDnnExtEta2Thr=cms.double(0.118044))
# Thresholds for electron: Bkg_nonIsolated
particleFlowTmp.PFEGammaFiltersParameters.electronDnnBkgThresholds = cms.PSet(
electronDnnBkgHighPtBarrelThr=cms.double(0.8),
electronDnnBkgHighPtEndcapThr=cms.double(0.75),
electronDnnBkgLowPtThr=cms.double(0.75),
electronDnnBkgExtEta1Thr=cms.double(0.7),
electronDnnBkgExtEta2Thr=cms.double(0.85))
# Thresholds for photons
particleFlowTmp.PFEGammaFiltersParameters.photonDnnThresholds = cms.PSet(
photonDnnBarrelThr=cms.double(0.22), photonDnnEndcapThr=cms.double(0.35))
from Configuration.Eras.Modifier_pf_badHcalMitigationOff_cff import pf_badHcalMitigationOff
pf_badHcalMitigationOff.toModify(
particleFlowTmp.PFEGammaFiltersParameters,
electron_protectionsForBadHcal=dict(enableProtections=False),
photon_protectionsForBadHcal=dict(enableProtections=False))
from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
egamma_lowPt_exclusive.toModify(particleFlowTmp.PFEGammaFiltersParameters,
photon_MinEt=1.)
# Activate Egamma PFID with DNN for Run3
from Configuration.Eras.Modifier_run3_common_cff import run3_common
run3_common.toModify(particleFlowTmp.PFEGammaFiltersParameters,
useElePFidDnn=True,
usePhotonPFidDnn=True)
from Configuration.Eras.Modifier_peripheralPbPb_cff import peripheralPbPb
peripheralPbPb.toModify(
trackerClusterCheck,
doClusterCheck=True, #FIXMETOO
cut=
"strip < 400000 && pixel < 40000 && (strip < 60000 + 7.0*pixel) && (pixel < 8000 + 0.14*strip)"
)
from Configuration.Eras.Modifier_pp_on_XeXe_2017_cff import pp_on_XeXe_2017
pp_on_XeXe_2017.toModify(
trackerClusterCheck,
doClusterCheck=True, #FIXMETOO
cut=
"strip < 1000000 && pixel < 100000 && (strip < 50000 + 10*pixel) && (pixel < 5000 + strip/2.)",
MaxNumberOfPixelClusters=100000)
from Configuration.Eras.Modifier_pp_on_AA_2018_cff import pp_on_AA_2018
pp_on_AA_2018.toModify(
trackerClusterCheck,
doClusterCheck=True,
cut=
"strip < 1000000 && pixel < 150000 && (strip < 50000 + 10*pixel) && (pixel < 5000 + strip/2.)",
MaxNumberOfPixelClusters=150000,
MaxNumberOfCosmicClusters=500000)
from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
egamma_lowPt_exclusive.toModify(trackerClusterCheck,
doClusterCheck=True,
cut="strip < 1000 && pixel < 300 ",
MaxNumberOfPixelClusters=300,
MaxNumberOfCosmicClusters=1000)