MissingETPlusTransverseMassSelectorAlg.cxx

Go to the documentation of this file.

/*
   Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
 
#include "EventSelectionAlgorithms/MissingETPlusTransverseMassSelectorAlg.h"
 
namespace CP {
 
  MissingETPlusTransverseMassSelectorAlg::MissingETPlusTransverseMassSelectorAlg(const std::string &name, ISvcLocator *pSvcLocator)
  : EL::AnaAlgorithm(name, pSvcLocator)
  {}
 
  StatusCode MissingETPlusTransverseMassSelectorAlg::initialize() {
    ANA_CHECK(m_metHandle.initialize(m_systematicsList));
    ANA_CHECK(m_electronsHandle.initialize(m_systematicsList, SG::AllowEmpty));
    ANA_CHECK(m_electronSelection.initialize(m_systematicsList, m_electronsHandle, SG::AllowEmpty));
    ANA_CHECK(m_muonsHandle.initialize(m_systematicsList, SG::AllowEmpty));
    ANA_CHECK(m_muonSelection.initialize(m_systematicsList, m_muonsHandle, SG::AllowEmpty));
    ANA_CHECK(m_eventInfoHandle.initialize(m_systematicsList));
 
    ANA_CHECK(m_preselection.initialize(m_systematicsList, m_eventInfoHandle, SG::AllowEmpty));
    ANA_CHECK(m_decoration.initialize(m_systematicsList, m_eventInfoHandle));
    ANA_CHECK(m_systematicsList.initialize());
 
    m_signEnum = SignEnum::stringToOperator.at( m_sign );
 
    return StatusCode::SUCCESS;
  }
 
  StatusCode MissingETPlusTransverseMassSelectorAlg::execute() {
    // accessors
    static const SG::AuxElement::ConstAccessor<float> acc_pt_dressed("pt_dressed");
    static const SG::AuxElement::ConstAccessor<float> acc_phi_dressed("phi_dressed");
 
    for (const auto &sys : m_systematicsList.systematicsVector()) {
      // retrieve the EventInfo
      const xAOD::EventInfo *evtInfo = nullptr;
      ANA_CHECK(m_eventInfoHandle.retrieve(evtInfo, sys));
 
      // default-decorate EventInfo
      m_decoration.setBool(*evtInfo, 0, sys);
 
      // check the preselection
      if (m_preselection && !m_preselection.getBool(*evtInfo, sys))
        continue;
 
      // retrieve the MET container
      const xAOD::MissingETContainer *met = nullptr;
      ANA_CHECK(m_metHandle.retrieve(met, sys));
      // retrieve the electron container
      const xAOD::IParticleContainer *electrons = nullptr;
      if (m_electronsHandle)
    ANA_CHECK(m_electronsHandle.retrieve(electrons, sys));
      // retrieve the electron container
      const xAOD::IParticleContainer *muons = nullptr;
      if (m_muonsHandle)
    ANA_CHECK(m_muonsHandle.retrieve(muons, sys));
 
      // compute the W boson transverse mass
      if ((*met)[m_metTerm.value()] == nullptr) {
        ANA_MSG_ERROR("The MET term " << m_metTerm.value() << " doesn't exist! Aborting.");
        return StatusCode::FAILURE;
      }
      float etmiss_pt = (*met)[m_metTerm.value()]->met();
      float etmiss_phi = (*met)[m_metTerm.value()]->phi();
      float lep_pt, lep_phi;
      int lep_count = 0;
 
      if (m_electronsHandle) {
    for (const xAOD::IParticle *el : *electrons) {
      if (!m_electronSelection || m_electronSelection.getBool(*el, sys)){
        if (m_useDressedProperties) {
          lep_pt = acc_pt_dressed(*el);
          lep_phi = acc_phi_dressed(*el);
        } else {
          lep_pt = el->pt();
          lep_phi = el->phi();
        }
        lep_count++;
        break;
      }
    }
      }
      if (m_muonsHandle) {
    for (const xAOD::IParticle *mu : *muons) {
      if (!m_muonSelection || m_muonSelection.getBool(*mu, sys)) {
        if (m_useDressedProperties) {
          lep_pt = acc_pt_dressed(*mu);
          lep_phi = acc_phi_dressed(*mu);
        } else {
          lep_pt = mu->pt();
          lep_phi = mu->phi();
        }
        lep_count++;
        break;
      }
    }
      }
 
      if (lep_count == 0){
        ATH_MSG_ERROR("No charged lepton in the event, cannot compute MWT!");
        return StatusCode::FAILURE;
      }
 
      if (lep_count > 1){
        ATH_MSG_ERROR("More than one lepton is associated with this event - cannot determine which is the W decay lepton.");
        return StatusCode::FAILURE;
      }
 
      float mwt = sqrt(2. * lep_pt * etmiss_pt * (1. - std::cos(lep_phi - etmiss_phi)));
 
      // compute the sum
      float sum = etmiss_pt + mwt;
 
      // calculate decision
      bool decision = SignEnum::checkValue(m_sumref.value(), m_signEnum, sum);
      m_decoration.setBool(*evtInfo, decision, sys);
    }
    return StatusCode::SUCCESS;
  }
} // namespace CP