EGammaIsoConfig.py

Go to the documentation of this file.

# Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
 
# ComponentAccumulator version of EGammaIso.py
 
# listOfScheduledAlgs = [ 'JetAlgorithm/jetalg_ConstitModCorrectPFOCHS_EMPFlow',
#                        'PseudoJetAlgorithm/PseudoJetAlgForIsoNFlow',
#                        'EventDensityAthAlg/CentralDensityForNFlowIsoAlg',
#                        'EventDensityAthAlg/ForwardDensityForNFlowIsoAlg',
#                        'IsolationBuilder/PFlowIsolationBuilder',
#                        'JetAlgorithm/jetalg_ConstitModCorrectPFOCSSKCHS_EMPFlowCSSK',
#                        'PseudoJetAlgorithm/PseudoJetAlgForIsoCSSKNFlow',
#                        'EventDensityAthAlg/CentralDensityForCSSKNFlowIsoAlg',
#                        'EventDensityAthAlg/ForwardDensityForCSSKNFlowIsoAlg',
#                        'IsolationBuilder/CSSKPFlowIsolationBuilder' ]
 
from AthenaConfiguration.ComponentAccumulator import ComponentAccumulator
from IsolationAlgs.IsolationSteeringDerivConfig import IsolationSteeringDerivCfg
 
 
def makeEGammaCommonIsoCfg(ConfigFlags):
    acc = ComponentAccumulator()
 
    addIsoVar = ""
    densityList = []
    densityDict = {}
    for inT in ["CSSK"]:
        acc.merge(IsolationSteeringDerivCfg(ConfigFlags, inType="EMPFlow" + inT))
 
        for eta in ["Central", "Forward"]:
            # I do not understand why this is needed
            densityDict.update(
                {
                    "NeutralParticle"
                    + inT
                    + "FlowIso"
                    + eta
                    + "EventShape": "xAOD::EventShape",
                    "NeutralParticle"
                    + inT
                    + "FlowIso"
                    + eta
                    + "EventShapeAux": "xAOD::EventShapeAuxInfo",
                }
            )
            densityList += [
                "NeutralParticle" + inT + "FlowIso" + eta + "EventShape.Density"
            ]
 
        suff = "" if len(inT) == 0 else "_" + inT
        addIsoVar += (
            f".neflowisol20{suff}.neflowisol30{suff}.neflowisol40{suff}"
            + f".neflowisolcoreConeEnergyCorrection{suff}"
            + f".neflowisolcoreConeSCEnergyCorrection{suff}"
        )
 
    return addIsoVar, densityList, densityDict, acc