AsgSelectionAlg.cxx

Go to the documentation of this file.

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

 
 
 
//
// includes
//
 
#include <AsgAnalysisAlgorithms/AsgSelectionAlg.h>
 
#include <PATInterfaces/ISystematicsTool.h>
 
//
// method implementations
//
 
namespace CP
{
 
  StatusCode AsgSelectionAlg ::
  initialize ()
  {
    if (!m_selectionTool.empty())
    {
      ANA_CHECK (m_selectionTool.retrieve());
      m_systematicsTool = dynamic_cast<ISystematicsTool*>(&*m_selectionTool);
      if (m_systematicsTool)
        ANA_CHECK (m_systematicsList.addSystematics (*m_systematicsTool));
    }
 
    ANA_CHECK (m_particlesHandle.initialize (m_systematicsList));
    ANA_CHECK (m_preselection.initialize (m_systematicsList, m_particlesHandle, SG::AllowEmpty));
    ANA_CHECK (m_selectionHandle.initialize (m_systematicsList, m_particlesHandle));
    ANA_CHECK (m_systematicsList.initialize());
 
    if (!m_selectionTool.empty())
    {
      asg::AcceptData blankAccept (&m_selectionTool->getAcceptInfo());
      // Just in case this isn't initially set up as a failure clear it this one
      // time. This only calls reset on the bitset
      blankAccept.clear();
      m_setOnFail = selectionFromAccept(blankAccept);
 
      if (!m_nameSvc.empty())
      {
        ANA_CHECK (m_nameSvc.retrieve());
        ANA_CHECK (m_nameSvc->addAcceptInfo (m_particlesHandle.getNamePattern(), m_selectionHandle.getLabel(),
            m_selectionTool->getAcceptInfo()));
      }
    }
 
    return StatusCode::SUCCESS;
  }
 
 
 
  StatusCode AsgSelectionAlg ::
  execute ()
  {
    for (const auto& sys : m_systematicsList.systematicsVector())
    {
      if (m_systematicsTool)
        ANA_CHECK (m_systematicsTool->applySystematicVariation (sys));
 
      const xAOD::IParticleContainer *particles = nullptr;
      ANA_CHECK (m_particlesHandle.retrieve (particles, sys));
      for (const xAOD::IParticle *particle : *particles)
      {
        if (m_preselection.getBool (*particle, sys))
        {
          if (!m_selectionTool.empty())
          {
            m_selectionHandle.setBits
              (*particle, selectionFromAccept (m_selectionTool->accept (particle)), sys);
          }
          else
          {
            m_selectionHandle.setBool (*particle, true, sys);
          }
        }
        else
        {
          if (!m_selectionTool.empty())
          {
            m_selectionHandle.setBits (*particle, m_setOnFail, sys);
          }
          else
          {
            m_selectionHandle.setBool (*particle, false, sys);
          }
        }
      }
    }
 
    return StatusCode::SUCCESS;
  }
}