egamma add lowMu mode e.g for UPC

egamma add lowMu mode e.g for UPC

Reconstruction/egamma/egammaAlgs/python/egammaSuperClusterBuilderConfig.py

@@ -16,14 +16,14 @@ from egammaMVACalib.egammaMVACalibConfig import egammaMVASvcCfg
 
 def electronSuperClusterBuilderCfg(flags,
                                    name='electronSuperClusterBuilder',
-                                   sequenceName = None,
+                                   sequenceName=None,
                                    **kwargs):
 
     mlog = logging.getLogger(name)
     mlog.debug('Start configuration')
 
     seqkw = {'sequence': sequenceName} if sequenceName else {}
-    acc = ComponentAccumulator (**seqkw)
+    acc = ComponentAccumulator(**seqkw)
 
     if "TrackMatchBuilderTool" not in kwargs:
         emtrkmatch = EMTrackMatchBuilderCfg(flags)
@@ -47,7 +47,7 @@ def electronSuperClusterBuilderCfg(flags,
         "egammaCheckEnergyDepositTool",
         CompFactory.egammaCheckEnergyDepositTool())
     kwargs.setdefault("EtThresholdCut", 1000)
-    kwargs.setdefault("UseExtendedTG3",flags.GeoModel.Run is LHCPeriod.Run3)
+    kwargs.setdefault("UseExtendedTG3", flags.GeoModel.Run is LHCPeriod.Run3)
     elscAlg = CompFactory.electronSuperClusterBuilder(name, **kwargs)
 
     acc.addEventAlgo(elscAlg)
@@ -57,11 +57,11 @@ def electronSuperClusterBuilderCfg(flags,
 def photonSuperClusterBuilderCfg(
         flags,
         name='photonSuperClusterBuilder',
-        sequenceName = None,
+        sequenceName=None,
         **kwargs):
 
     seqkw = {'sequence': sequenceName} if sequenceName else {}
-    acc = ComponentAccumulator (**seqkw)
+    acc = ComponentAccumulator(**seqkw)
 
     photonSuperClusterBuilder = CompFactory.photonSuperClusterBuilder
     egammaCheckEnergyDepositTool = CompFactory.egammaCheckEnergyDepositTool
@@ -87,7 +87,11 @@ def photonSuperClusterBuilderCfg(
     kwargs.setdefault(
         "egammaCheckEnergyDepositTool",
         egammaCheckEnergyDepositTool())
-    kwargs.setdefault("UseExtendedTG3",flags.GeoModel.Run is LHCPeriod.Run3)
+    kwargs.setdefault("UseExtendedTG3", flags.GeoModel.Run is LHCPeriod.Run3)
+
+    kwargs.setdefault(
+        "EtThresholdCut",
+        1500 if not flags.Egamma.doLowMu else 300)
     phscAlg = photonSuperClusterBuilder(name, **kwargs)
 
     acc.addEventAlgo(phscAlg)

Reconstruction/egamma/egammaAlgs/python/egammaTopoClusterCopierConfig.py

@@ -25,6 +25,10 @@ def egammaTopoClusterCopierCfg(
         "OutputTopoCollectionShallow",
         "tmp_"+egtopocluster)
 
+    kwargs.setdefault(
+        "ECut",
+        700 if not flags.Egamma.doLowMu else 300)
+
     egcopierAlg = CompFactory.egammaTopoClusterCopier(name, **kwargs)
 
     acc.addEventAlgo(egcopierAlg)

Reconstruction/egamma/egammaConfig/python/egammaConfigFlags.py

@@ -47,6 +47,9 @@ def createEgammaConfigFlags():
     egcf.addFlag("Egamma.slimGSFTrkTracks",
                  lambda prevFlags: prevFlags.Egamma.doTracking)
 
+    # Egamma runs in low <mu> mode (e.g UPC )
+    egcf.addFlag("Egamma.doLowMu", False)
+
     # The cluster corrections/calib
     egcf.addFlag("Egamma.Calib.ClusterCorrectionVersion",
                  'v12phiflip_noecorrnogap')
@@ -91,8 +94,11 @@ def createEgammaConfigFlags():
                  '-topoetconeCorrBitset')
 
     egcf.addFlag("Egamma.Keys.Output.Electrons", 'Electrons')
-    #Remove GlobalFELinks - these are links between FlowElement (FE) containers created in jet finding and electrons/photons. Since these transient FE containers are not in the ESD/AOD, we should not write out these links.
-    gpf_vars='-chargedGlobalFELinks.-neutralGlobalFELinks.'
+    # Remove GlobalFELinks - these are links between FlowElement (FE)
+    # containers created in jet finding and electrons/photons.
+    # Since these transient FE containers are not in the ESD/AOD,
+    # we should not write out these links.
+    gpf_vars = '-chargedGlobalFELinks.-neutralGlobalFELinks.'
     egcf.addFlag("Egamma.Keys.Output.ElectronsSuppESD", gpf_vars)
     egcf.addFlag("Egamma.Keys.Output.ElectronsSuppAOD",
                  lambda prevFlags: (
@@ -137,7 +143,7 @@ def createEgammaConfigFlags():
                  lambda prevFlags: (
                      prevFlags.Egamma.Keys.Output.TruthParticlesSuppESD))
 
-    # not xAOD
+    # Not written in AOD
     egcf.addFlag("Egamma.Keys.Output.GSFTracks", 'GSFTracks')
 
     # These are the clusters that are used to determine
@@ -154,8 +160,6 @@ def createEgammaConfigFlags():
 
 
 # self test
-
-
 class TestEgammaConfigFlags(unittest.TestCase):
 
     def runTest(self):