TrigBmumuxComboHypoTool.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
#include "TrigBmumuxComboHypoTool.h"
 
 
using TrigCompositeUtils::Decision;
using TrigCompositeUtils::DecisionID;
using TrigCompositeUtils::DecisionIDContainer;
 
TrigBmumuxComboHypoTool::TrigBmumuxComboHypoTool(const std::string& type, const std::string& name, const IInterface* parent)
    : ComboHypoToolBase(type, name, parent) {}
 
 
StatusCode TrigBmumuxComboHypoTool::initialize()
{
  if (!m_monTool.empty()) {
    ATH_CHECK( m_monTool.retrieve() );
    ATH_MSG_DEBUG( "GenericMonitoringTool name:" << m_monTool );
  }
  else {
    ATH_MSG_DEBUG( "No GenericMonitoringTool configured: no monitoring histograms will be available" );
  }
 
  ATH_MSG_DEBUG( "Initialization completed successfully" );
 
  return StatusCode::SUCCESS;
}
 
 
bool TrigBmumuxComboHypoTool::passed(const xAOD::TrigBphys* trigBphys) const {
 
  using namespace Monitored;
 
  auto mon_chi2 = Monitored::Scalar<float>("Chi2", -1.);
  auto mon_fitmass = Monitored::Scalar<float>("Fitmass", -1.);
  auto mon_mass = Monitored::Scalar<float>("Mass", -1.);
  auto mon_pt = Monitored::Scalar<float>("Pt", -1.);
  auto mon_eta = Monitored::Scalar<float>("Eta", -100.);
  auto group = Monitored::Group(m_monTool, mon_chi2, mon_fitmass, mon_mass, mon_pt, mon_eta);
 
  ATH_MSG_DEBUG( "in TrigBmumuxComboHypoTool::decideOnSingleObject(), looking at TrigBphys object");
 
  bool result = false;
 
  if (static_cast<int>(trigBphys->particleType()) == m_decay) {
    result = true;
    ATH_MSG_DEBUG("accepting event");
  }
 
  if (result) {
    mon_chi2 = trigBphys->fitchi2();
    mon_fitmass = trigBphys->fitmass();
    mon_mass = trigBphys->mass();
    mon_pt = trigBphys->pt();
    mon_eta = trigBphys->eta();
  }
 
  return result;
}
 
 
StatusCode TrigBmumuxComboHypoTool::decideOnSingleObject(Decision* decision, const std::vector<const DecisionIDContainer*>& previousDecisionIDs) const {
 
  ATH_CHECK( decision->hasObjectLink(TrigCompositeUtils::featureString()) );
 
  auto trigBphysEL = decision->objectLink<xAOD::TrigBphysContainer>(TrigCompositeUtils::featureString());
  ATH_CHECK( trigBphysEL.isValid() );
 
  const std::vector<HLT::Identifier>& legDecisionIDs = legDecisionIds();
  if (m_isBmux) {  // bBmux chain: only one previousDecision is available, in case of assymmetric chain it should match to one arbitrary leg
    ATH_CHECK( previousDecisionIDs.size() == 1 );
    bool passed = false;
    for (size_t i = 0; i < legDecisionIDs.size(); ++i) {
      const DecisionID id = legDecisionIDs[i].numeric();
      if (TrigCompositeUtils::passed(id, *previousDecisionIDs[0])) {
        passed = true;
        break;
      }
    }
    if (!passed) {
      ATH_MSG_DEBUG( "bBmux chain did not match to the previous decisions" );
      return StatusCode::SUCCESS;
    }
  }
  else if (legDecisionIDs.size() == 1) {  // bBmumux chain with symmetric legs
    auto n = static_cast<size_t>(legMultiplicity().at(0));
    ATH_CHECK( previousDecisionIDs.size() == n );
    const DecisionID id = legDecisionIDs[0].numeric();
    for (size_t i = 0; i < n; ++i) {
      if (!TrigCompositeUtils::passed(id, *previousDecisionIDs[i])) {
        ATH_MSG_DEBUG( "bBmumux chain with symmetric legs did not match to the previous decisions" );
        return StatusCode::SUCCESS;
      }
    }
  }
  else if (legDecisionIDs.size() == 2) {  // bBmumux chain with asymmetric legs
    ATH_CHECK( previousDecisionIDs.size() == 2 );
    bool direct = true;
    bool inverse = true;
    for (size_t i = 0; i < 2; ++i) {
      if (direct && !TrigCompositeUtils::passed(legDecisionIDs.at(i).numeric(), *previousDecisionIDs.at(i))) direct = false;
      if (inverse && !TrigCompositeUtils::passed(legDecisionIDs.at(i).numeric(), *previousDecisionIDs.at(1-i))) inverse = false;
    }
    if (!direct && !inverse) {
      ATH_MSG_DEBUG( "bBmumux chain with asymmetric legs matched to the previous decisions neither direct nor inverse way" );
      return StatusCode::SUCCESS;
    }
  }
  else {
    ATH_MSG_ERROR( "TrigBmumuxComboHypoTool can not check previous decisions for " << name() );
    return StatusCode::FAILURE;
  }
 
  if (passed(*trigBphysEL)) {
    TrigCompositeUtils::addDecisionID(decisionId(), decision);
  }
 
  return StatusCode::SUCCESS;
}
 
 
bool TrigBmumuxComboHypoTool::executeAlg(const std::vector<Combo::LegDecision>&) const {
  ATH_MSG_ERROR("executeAlg not supported for TrigBmumuxComboHypoTool.");
  return true;
}