MetSignificanceAlg.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
 
 
 
//
// includes
//
 
#include <MetAnalysisAlgorithms/MetSignificanceAlg.h>
 
#include <xAODMissingET/MissingETAuxContainer.h>
#include "xAODEventInfo/EventInfo.h"
 
//
// method implementations
//
namespace CP
{
 
  StatusCode MetSignificanceAlg ::
  initialize ()
  {
    if (m_significanceDecoration.empty())
    {
      ANA_MSG_ERROR ("no significance decoration name set");
      return StatusCode::FAILURE;
    }
    m_significanceAccessor = std::make_unique<SG::AuxElement::Accessor<float> > (m_significanceDecoration);
 
    ANA_CHECK (m_significanceTool.retrieve());
    ANA_CHECK (m_metHandle.initialize (m_systematicsList));
    ANA_CHECK (m_systematicsList.initialize());
    return StatusCode::SUCCESS;
  }
 
 
 
  StatusCode MetSignificanceAlg ::
  execute ()
  {
    for (const auto& sys : m_systematicsList.systematicsVector())
    {
      // I'm not sure why this can't be const, but the interface
      // requires a non-const object
      xAOD::MissingETContainer *met {};
      ANA_CHECK (m_metHandle.getCopy (met, sys));
    
      const xAOD::EventInfo* evtInfo = 0;
      ANA_CHECK( evtStore()->retrieve( evtInfo, "EventInfo" ) );
 
      ANA_CHECK (m_significanceTool->varianceMET (met, evtInfo->averageInteractionsPerCrossing(), m_jetTermName, m_softTermName, m_totalMETName));
      const float significance = m_significanceTool->GetSignificance();
      (*m_significanceAccessor) (*(*met)[m_totalMETName]) = significance;
    }
 
    return StatusCode::SUCCESS;
  }
}