IParticleRetrievalTool.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
// Package includes
#include "TriggerMatchingTool/IParticleRetrievalTool.h"
#include <ostream>
 
// Anonymous namespace contains helper functions
namespace {
  using CLIDTuple_t =
    std::tuple<HLT::class_id_type, xAOD::Type::ObjectType, std::string>;
 
  template <std::size_t FROM, std::size_t TO>
    bool CLIDLookUp(
      typename std::tuple_element<TO, CLIDTuple_t>::type& to,
      const typename std::tuple_element<FROM, CLIDTuple_t>::type& from)
    {
      static const std::vector<CLIDTuple_t> CLIDVector {
        {1178459224, xAOD::Type::Muon, "xAOD::MuonContainer"},
        {1087532415, xAOD::Type::Electron, "xAOD::ElectronContainer"},
        {1219821989, xAOD::Type::CaloCluster, "xAOD::CaloClusterContainer"},
        {1105575213, xAOD::Type::Photon, "xAOD::PhotonContainer"},
        {1177172564, xAOD::Type::Tau, "xAOD::TauJetContainer"}
      };
      auto itr = std::find_if(CLIDVector.begin(), CLIDVector.end(),
          [from] (const CLIDTuple_t& tup) {
            return std::get<FROM>(tup) == from;
          });
      if (itr == CLIDVector.end() )
        return false;
      else {
        to = std::get<TO>(*itr);
        return true;
      }
    }
  bool CLIDToObjectType(xAOD::Type::ObjectType& type, HLT::class_id_type clid)
  {
    return CLIDLookUp<0, 1>(type, clid);
  }
  bool ObjectTypeToCLID(HLT::class_id_type& clid, xAOD::Type::ObjectType type)
  {
    return CLIDLookUp<1, 0>(clid, type);
  }
  bool CLIDToContainerType(std::string& name, HLT::class_id_type clid)
  {
    return CLIDLookUp<0, 2>(name, clid);
  }
 
} //> end anonymous namespace
 
namespace Trig {
  IParticleRetrievalTool::IParticleRetrievalTool(const std::string& name) :
    asg::AsgTool(name)
  {
    declareProperty("TrigDecisionTool", m_tdt,
        "The trigger decision tool to use.");
    declareProperty("WarnOnNavigationError", m_warnOnNavigationFailure = true,
        "Only spit a warning on a navigation error, don't cause a job failure.");
 
  }
 
  IParticleRetrievalTool::~IParticleRetrievalTool() {}
 
  StatusCode IParticleRetrievalTool::initialize()
  {
    ATH_MSG_INFO("Initializing " << name() );
    ATH_CHECK( m_tdt.retrieve() );
    return StatusCode::SUCCESS;
  }
 
  StatusCode IParticleRetrievalTool::retrieveParticles(
      std::vector<std::vector<const xAOD::IParticle*>>& combinations,
      const std::string& chain,
      bool rerun) const
  {
    // Make sure what we're getting is cleared.
    combinations.clear();
    // Start by getting the chain group
    const ChainGroup* cg = m_tdt->getChainGroup(chain);
    // Make sure that this group actually contains triggers. If it doesn't then
    // this is very wrong and the job shouldn't be allowed to continue
    if (cg->getListOfTriggers().size() == 0) {
      ATH_MSG_ERROR("Chain group " << chain 
          << " is empty! This means that no matching chains were found!");
      return StatusCode::FAILURE;
    }
    unsigned int condition = TrigDefs::Physics;
    if (rerun)
      condition |= TrigDefs::allowResurrectedDecision;
    if (!cg->isPassed(condition) ) {
      ATH_MSG_DEBUG("Chain: " << chain << " was not passed!");
      return StatusCode::SUCCESS;
    }
    
    FeatureContainer features = cg->features(condition);
    for (const Combination& combo : features.getCombinations() ) {
      // The assumption here is that each combination represents a *single* way
      // in which the trigger could have been passed. This should be true for
      // the types of trigger that this tool expects to examine. It won't be
      // true for (e.g.) jet triggers. 
      std::vector<const xAOD::IParticle*> currentCombination;
      bool navFailure = false;
      for (const HLT::TriggerElement* te : combo.tes() ) {
        if (!retrieveParticles(currentCombination, te, navFailure).isSuccess() ) {
          // Interpret a failure this way so that we can report the chain name
          ATH_MSG_ERROR("Failed to retrieve particles for chain " << chain );
          return StatusCode::FAILURE;
        }
        // If the navigation failed for this combination ignore it completely
        if (navFailure)
          break;
      }
      if (navFailure)
        ATH_MSG_WARNING("Skipping combination for chain " << chain 
            << " due to navigation failure");
      else
        combinations.push_back(currentCombination);
    }
    return StatusCode::SUCCESS;
  }
 
  StatusCode IParticleRetrievalTool::retrieveParticles(
      std::vector<const xAOD::IParticle*>& combination,
      const HLT::TriggerElement* te,
      bool& navFailure) const
  {
    ATH_MSG_DEBUG( "Processing TE " << Trig::getTEName(*te) );
    // Keep track of whether or not we found a particle here.
    const xAOD::IParticle* part = nullptr;
    for (const auto& feature : te->getFeatureAccessHelpers() ) {
      xAOD::Type::ObjectType type;
      if (!CLIDToObjectType(type, feature.getCLID() ) )
        // Skip any features that don't hold 'final' state particles
        continue;
      if (type == xAOD::Type::Electron || type == xAOD::Type::Photon) {
        // If this is an egamma type we have to check if we are meant to be
        // looking for a CaloCluster instead!
        xAOD::Type::ObjectType egType = getEGammaTEType(te);
        if (egType == xAOD::Type::Other) {
          std::string message = 
            "Unable to determine the correct type for TE " + Trig::getTEName(*te);
          navFailure = true;
          if (m_warnOnNavigationFailure) {
            ATH_MSG_WARNING( message );
            return StatusCode::SUCCESS;
          }
          else {
            ATH_MSG_ERROR( message );
            return StatusCode::FAILURE;
          }
        }
        else if (egType == xAOD::Type::CaloCluster) {
          HLT::class_id_type clid = 0;
          ObjectTypeToCLID(clid, egType);
          // This will be where we store the TE containing the calo cluster
          // feature
          const HLT::TriggerElement* sourceTE = nullptr;
          const HLT::TrigNavStructure* navigation = 
            m_tdt->ExperimentalAndExpertMethods().getNavigation();
          HLT::TriggerElement::FeatureAccessHelper egFeature = 
            navigation->getFeatureRecursively(te, clid, "", sourceTE);
          if (!sourceTE) {
            std::ostringstream os;
            os << "Unable to retrieve feature of type " << egType
                << " from TE " << Trig::getTEName(*te);
            navFailure = true;
            if (m_warnOnNavigationFailure) {
              ATH_MSG_WARNING( os.str() );
              return StatusCode::SUCCESS;
            }
            else {
              ATH_MSG_ERROR( os.str() );
              return StatusCode::FAILURE;
            }
          }
          ATH_CHECK( retrieveFeatureParticle(part, egFeature, sourceTE, navFailure) );
          if (navFailure)
            return StatusCode::SUCCESS;
          // If it's a calo-cluster like TE then stop here, otherwise we'll
          // encounter the other EG type again and fail...
          break;
        }
        else if (egType != type) {
          // This is the wrong feature to be looking at
          continue;
        }
        // Otherwise we continue as before
      }
      // If we found a particle from another feature access helper then this is
      // a problem. Our assumption is that there is only one particle per leg!
      if (part) {
        std::string message = "TE" + Trig::getTEName(*te) + "has multiple " +
          "'final' particles attached to it! This breaks this tool's asumptions!";
        navFailure = true;
        if (m_warnOnNavigationFailure) {
          ATH_MSG_WARNING(message);
          return StatusCode::SUCCESS;
        }
        else {
          ATH_MSG_ERROR(message);
          return StatusCode::FAILURE;
        }
      }
      ATH_CHECK( retrieveFeatureParticle(part, feature, te, navFailure) );
      if (navFailure)
        return StatusCode::SUCCESS;
    } //> end loop over features
 
    // If we found a particle then we can stop going through this branch of the
    // tree
    if (part) {
      combination.push_back(part);
      return StatusCode::SUCCESS;
    }
 
    // Otherwise look at each of the next TEs separately
    for (const HLT::TriggerElement* nextTE : 
        HLT::TrigNavStructure::getDirectPredecessors(te) )
      ATH_CHECK( retrieveParticles(combination, nextTE, navFailure) );
 
    return StatusCode::SUCCESS;
  }
 
  xAOD::Type::ObjectType IParticleRetrievalTool::getEGammaTEType(
      const HLT::TriggerElement* te) const
  {
    // We need to check the TE name for this
    std::string teName = Trig::getTEName(*te);
    if (teName.find("etcut") != std::string::npos &&
        teName.find("trkcut") == std::string::npos)
      return xAOD::Type::CaloCluster;
    else if (teName.starts_with( "EF_e"))
      return xAOD::Type::Electron;
    else if (teName.starts_with( "EF_g"))
      return xAOD::Type::Photon;
    else 
      return xAOD::Type::Other;
  }
 
  StatusCode IParticleRetrievalTool::retrieveFeatureParticle(
      const xAOD::IParticle*& particle,
      const HLT::TriggerElement::FeatureAccessHelper& feature,
      const HLT::TriggerElement* te,
      bool& navFailure) const
  {
    // Get the right container type name
    std::string containerType;
    if (!CLIDToContainerType(containerType, feature.getCLID() ) ) {
      ATH_MSG_ERROR("Unrecognised CLID " << feature.getCLID() << " received!");
      // This isn't a navigation error - this is a failure in the tool's
      // internal logic!
      return StatusCode::FAILURE;
    }
    // Build a vector of typeless features so that we can use the central
    // functions.
    const HLT::TrigNavStructure* navigation = 
      m_tdt->ExperimentalAndExpertMethods().getNavigation();
 

    auto typelessHolder = navigation->getHolder(feature);
 
    if (!typelessHolder) {
      std::string message = "Typeless holder for feature not present!";
      navFailure = true;
      if (m_warnOnNavigationFailure) {
        ATH_MSG_WARNING(message);
        return StatusCode::SUCCESS;
      }
      else {
        ATH_MSG_ERROR(message);
        return StatusCode::FAILURE;
      }
    }
    const xAOD::IParticleContainer* cont(nullptr);
    // Get the name used in the event store
    std::string key = HLTNavDetails::formatSGkey(
        "HLT", containerType, typelessHolder->label() );
    // Now things are *much* more familiar
    if (!evtStore()->contains<xAOD::IParticleContainer>(key) ) {
      std::string message = "Store does not contain " + key + "!";
      navFailure = true;
      if (m_warnOnNavigationFailure) {
        ATH_MSG_WARNING(message);
        return StatusCode::SUCCESS;
      }
      else {
        ATH_MSG_ERROR(message);
        return StatusCode::FAILURE;
      }
    }
    ATH_CHECK( evtStore()->retrieve(cont, key) );
    HLT::TriggerElement::ObjectIndex idx = feature.getIndex();
    if (cont->size() < idx.objectsEnd() ) {
      std::ostringstream os; 
      os << "Featured object end " << idx.objectsEnd()
        << " is *after* the end of container " << key;
      navFailure = true;
      if (m_warnOnNavigationFailure) {
        ATH_MSG_WARNING( os.str() );
        return StatusCode::SUCCESS;
      }
      else {
        ATH_MSG_ERROR( os.str() );
        return StatusCode::FAILURE;
      }
    }
    std::vector<const xAOD::IParticle*> particleFeatures;
    particleFeatures.reserve(idx.objectsEnd() - idx.objectsBegin() );
    auto begin = cont->begin();
    auto end = cont->begin();
    std::advance(begin, idx.objectsBegin() );
    std::advance(end, idx.objectsEnd() );
    particleFeatures.insert(particleFeatures.end(), begin, end);
 
    // Make sure the answer is what we expect
    std::ostringstream os;
    switch (particleFeatures.size() ) {
      case 0:
        os << "No particles retrieved from feature "
            << navigation->label(feature.getCLID(), feature.getIndex().subTypeIndex() )
            << " from TE " << Trig::getTEName(*te);
        navFailure = true;
        if (m_warnOnNavigationFailure) {
          ATH_MSG_WARNING(os.str() );
          return StatusCode::SUCCESS;
        }
        else {
          ATH_MSG_ERROR(os.str() );
          return StatusCode::FAILURE;
        }
      case 1:
        // Set the output.
        particle = particleFeatures.at(0);
        break;
      default:
        // Some TEs can end up reporting multiple outputs within the same RoI.
        // AFAIK this only happens within EGamma TEs but I don't know that for
        // sure. In any case this shouldn't matter too much for the matching
        // given that they will be nearby each other. Just return the highest pT
        // object.
        particle = *(std::max_element(
              particleFeatures.begin(), particleFeatures.end(),
              [] (const xAOD::IParticle* lhs, const xAOD::IParticle* rhs)
              { return lhs->pt() < rhs->pt(); }) );
    }
    return StatusCode::SUCCESS;
  }
 
} //> end namespace Trig