Update the egamma regression for UPC PbPb

cms-sw · Oct 31, 2024 · 0aa2b12 · 0aa2b12
1 parent ee08414
commit 0aa2b12
Show file tree

Hide file tree

Showing 3 changed files with 47 additions and 8 deletions.
diff --git a/RecoEcal/EgammaClusterProducers/python/particleFlowSuperClusterECAL_cfi.py b/RecoEcal/EgammaClusterProducers/python/particleFlowSuperClusterECAL_cfi.py
@@ -24,6 +24,12 @@
 
 from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
 egamma_lowPt_exclusive.toModify(particleFlowSuperClusterECAL,
+                                regressionConfig = dict(
+                                    regressionKeyEB  = 'pfscecal_ebCorrection_offline_v2',
+                                    uncertaintyKeyEB = 'pfscecal_ebUncertainty_offline_v2',
+                                    regressionKeyEE  = 'pfscecal_eeCorrection_offline_v2',
+                                    uncertaintyKeyEE = 'pfscecal_eeUncertainty_offline_v2'
+                                ),
                                 thresh_SCEt = 1.0,
                                 thresh_PFClusterSeedBarrel = 0.5,
                                 thresh_PFClusterSeedEndcap = 0.5)

diff --git a/RecoEgamma/EgammaElectronProducers/python/lowPtGsfElectrons_cfi.py b/RecoEgamma/EgammaElectronProducers/python/lowPtGsfElectrons_cfi.py
@@ -39,6 +39,41 @@
     ),
 )
 
+from RecoEgamma.EgammaTools.regressionModifier_cfi import regressionModifier
+_lowPtRegressionModifierUPC = regressionModifier.clone(
+    eleRegs = dict(
+        ecalOnlyMean = dict(
+            ebLowEtForestName = ":lowPtElectron_eb_ecalOnly_1To20_0p2To2_mean",
+            ebHighEtForestName = ":lowPtElectron_eb_ecalOnly_1To20_0p2To2_mean",
+            eeLowEtForestName = ":lowPtElectron_ee_ecalOnly_1To20_0p2To2_mean",
+            eeHighEtForestName = ":lowPtElectron_ee_ecalOnly_1To20_0p2To2_mean",
+            ),
+        ecalOnlySigma = dict(
+            ebLowEtForestName = ":lowPtElectron_eb_ecalOnly_1To20_0p0002To0p5_sigma",
+            ebHighEtForestName = ":lowPtElectron_eb_ecalOnly_1To20_0p0002To0p5_sigma",
+            eeLowEtForestName = ":lowPtElectron_ee_ecalOnly_1To20_0p0002To0p5_sigma",
+            eeHighEtForestName = ":lowPtElectron_ee_ecalOnly_1To20_0p0002To0p5_sigma",
+            ),
+        epComb = dict(
+            ecalTrkRegressionConfig = dict(
+                ebLowEtForestName = ":lowPtElectron_eb_ecalTrk_1To20_0p2To2_mean",
+                ebHighEtForestName = ":lowPtElectron_eb_ecalTrk_1To20_0p2To2_mean",
+                eeLowEtForestName = ":lowPtElectron_ee_ecalTrk_1To20_0p2To2_mean",
+                eeHighEtForestName = ":lowPtElectron_ee_ecalTrk_1To20_0p2To2_mean",
+                ),
+            ecalTrkRegressionUncertConfig = dict(
+                ebLowEtForestName = ":lowPtElectron_eb_ecalTrk_1To20_0p0002To0p5_sigma",
+                ebHighEtForestName = ":lowPtElectron_eb_ecalTrk_1To20_0p0002To0p5_sigma",
+                eeLowEtForestName = ":lowPtElectron_ee_ecalTrk_1To20_0p0002To0p5_sigma",
+                eeHighEtForestName = ":lowPtElectron_ee_ecalTrk_1To20_0p0002To0p5_sigma",
+                ),
+        )
+    ),
+)
+from Configuration.Eras.Modifier_run3_egamma_cff import run3_egamma
+from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
+(egamma_lowPt_exclusive & run3_egamma).toReplaceWith(lowPtRegressionModifier,_lowPtRegressionModifierUPC)
+
 from RecoEgamma.EgammaElectronProducers.lowPtGsfElectronFinalizer_cfi import lowPtGsfElectronFinalizer
 lowPtGsfElectrons = lowPtGsfElectronFinalizer.clone(
     previousGsfElectronsTag = "lowPtGsfElectronsPreRegression",

diff --git a/RecoEgamma/EgammaTools/python/regressionModifier_cfi.py b/RecoEgamma/EgammaTools/python/regressionModifier_cfi.py
@@ -245,10 +245,14 @@
 from Configuration.Eras.Modifier_run2_egamma_2016_cff import run2_egamma_2016
 from Configuration.Eras.Modifier_run2_egamma_2017_cff import run2_egamma_2017
 from Configuration.Eras.Modifier_run2_egamma_2018_cff import run2_egamma_2018
+from Configuration.Eras.Modifier_run3_egamma_cff import run3_egamma
 
-(run2_egamma_2016 | run2_egamma_2017 | run2_egamma_2018).toReplaceWith(regressionModifier,regressionModifierRun2)
+(run2_egamma_2016 | run2_egamma_2017 | run2_egamma_2018 | run3_egamma).toReplaceWith(regressionModifier,regressionModifierRun2)
 
-regressionModifierRun3 = regressionModifierRun2.clone(
+from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
+egamma_lowPt_exclusive.toReplaceWith(regressionModifier,regressionModifier103XLowPtPho)
+
+run3_egamma.toModify(regressionModifier,
     eleRegs = dict(
         ecalOnlyMean = dict(
             rangeMinHighEt = 0.2,
@@ -262,9 +266,3 @@
         )
     )
 )
-
-from Configuration.Eras.Modifier_run3_egamma_cff import run3_egamma
-run3_egamma.toReplaceWith(regressionModifier,regressionModifierRun3)
-
-from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
-egamma_lowPt_exclusive.toReplaceWith(regressionModifier,regressionModifier103XLowPtPho)