EgammaCalibrationAndSmearingAlg.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/

 
 
 
//
// includes
//
 
#include <EgammaAnalysisAlgorithms/EgammaCalibrationAndSmearingAlg.h>
 
//
// method implementations
//
 
namespace CP
{
 
  StatusCode EgammaCalibrationAndSmearingAlg ::
  initialize ()
  {
    ANA_CHECK (m_calibrationAndSmearingTool.retrieve());
    ANA_CHECK (m_egammaHandle.initialize (m_systematicsList));
    ANA_CHECK (m_preselection.initialize (m_systematicsList, m_egammaHandle, SG::AllowEmpty));
    ANA_CHECK (m_caloClusterEnergyResoHandle.initialize (m_systematicsList, m_egammaHandle));
    if (!m_noToolSystematics.value())
      ANA_CHECK (m_systematicsList.addSystematics (*m_calibrationAndSmearingTool));
    ANA_CHECK (m_systematicsList.initialize());
    ANA_CHECK (m_outOfValidity.initialize());
    return StatusCode::SUCCESS;
  }
 
 
 
  StatusCode EgammaCalibrationAndSmearingAlg ::
  execute (const EventContext& ctx)
  {
    for (const auto& sys : m_systematicsList.systematicsVector())
    {
      // always need to call `getCopy` first to ensure that the shallow copies
      // are all there if requested
      xAOD::EgammaContainer *egammas = nullptr;
      ANA_CHECK (m_egammaHandle.getCopy (egammas, sys, ctx));
 
      if (sys.empty() && m_skipNominal.value())
        continue;
 
      if (!m_noToolSystematics.value())
        ANA_CHECK (m_calibrationAndSmearingTool->applySystematicVariation (sys));
      for (xAOD::Egamma *egamma : *egammas)
      {
        if (m_preselection.getBool (*egamma, sys))
        {
          ANA_CHECK_CORRECTION (m_outOfValidity, *egamma, m_calibrationAndSmearingTool->applyCorrection (*egamma));
          float caloClusterEnergyReso = m_calibrationAndSmearingTool->getResolution (*egamma);
          m_caloClusterEnergyResoHandle.set(*egamma, caloClusterEnergyReso, sys);
        }
      }
    }
    return StatusCode::SUCCESS;
  }
}