EgammaDetails.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef XAODEGAMMA_EGAMMADETAILS_H
#define XAODEGAMMA_EGAMMADETAILS_H
 
#include "xAODEgamma/PhotonxAODHelpers.h"
#include "xAODEgamma/ElectronxAODHelpers.h"
 
namespace xAOD
{
  namespace EgammaDetails
  {
 
    using EgammaParameters::ConversionType;

    inline ConversionType conversionType(const bool hasTrk1, const bool hasTrk2, const std::uint8_t nSiHits1, const std::uint8_t nSiHits2) {
 
      if (!hasTrk1) {return xAOD::EgammaParameters::unconverted;}
 
      if (!hasTrk2)
        {return nSiHits1 ? xAOD::EgammaParameters::singleSi : xAOD::EgammaParameters::singleTRT;}
      
      if (nSiHits1 && nSiHits2){
        return xAOD::EgammaParameters::doubleSi;
      }
      if (nSiHits1 || nSiHits2){
        return xAOD::EgammaParameters::doubleSiTRT;
      }  
      else{
        return xAOD::EgammaParameters::doubleTRT;
      }
    }

    inline bool isConvertedPhoton (const bool excludeTRT, const float eta, const std::size_t nVertices, const ConversionType conversionType) {
      const bool hasVertices = nVertices > 0;
      if (excludeTRT) {
        // special case for Run3: consider unconv if TRT Conv in the barrel
        using enum EgammaParameters::ConversionType;
        const bool isTRTConv = (conversionType == singleTRT) || (conversionType == doubleTRT);
        return hasVertices && (std::abs(eta) > 0.8 || !isTRTConv);
      }
      return hasVertices;
    }
  }
}
 
#endif