AsgEventScaleFactorAlg.cxx

Go to the documentation of this file.

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

 
 
//
// includes
//
 
#include <AsgAnalysisAlgorithms/AsgEventScaleFactorAlg.h>
 
//
// method implementations
//
 
namespace CP
{
  StatusCode AsgEventScaleFactorAlg ::
  initialize ()
  {
    if (m_scaleFactorInputDecoration.empty() || m_scaleFactorOutputDecoration.empty())
    {
      ANA_MSG_ERROR ("no scale factor decoration name set");
      return StatusCode::FAILURE;
    }
 
    ANA_CHECK (m_eventInfoHandle.initialize (m_systematicsList));
    ANA_CHECK (m_particleHandle.initialize (m_systematicsList));
    ANA_CHECK (m_preselection.initialize (m_systematicsList, m_particleHandle, SG::AllowEmpty));
    ANA_CHECK (m_scaleFactorInputDecoration.initialize (m_systematicsList, m_particleHandle));
    ANA_CHECK (m_scaleFactorOutputDecoration.initialize (m_systematicsList, m_eventInfoHandle));
    ANA_CHECK (m_systematicsList.initialize());
 
    return StatusCode::SUCCESS;
  }
 
 
 
  StatusCode AsgEventScaleFactorAlg ::
  execute (const EventContext& ctx)
  {
    for (const auto& sys : m_systematicsList.systematicsVector())
    {
      const xAOD::EventInfo *eventInfo = nullptr;
      ANA_CHECK (m_eventInfoHandle.retrieve (eventInfo, sys, ctx));
 
      const xAOD::IParticleContainer *particles = nullptr;
      ANA_CHECK (m_particleHandle.retrieve (particles, sys, ctx));
 
      float scaleFactor = 1;
      for (const xAOD::IParticle *particle : *particles)
      {
        if (m_preselection.getBool (*particle, sys))
        {
          scaleFactor *= m_scaleFactorInputDecoration.get (*particle, sys);
        }
      }
 
      m_scaleFactorOutputDecoration.set (*eventInfo, scaleFactor, sys);
    }
 
    return StatusCode::SUCCESS;
  }
}