eFexTauRoI_v1.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
 
// System include(s):
#include <stdexcept>
#include <numbers>
 
// xAOD include(s):
#include "xAODCore/AuxStoreAccessorMacros.h"
 
// Local include(s):
#include "xAODTrigger/versions/eFexTauRoI_v1.h"
 
namespace xAOD {

  const float eFexTauRoI_v1::s_tobEtScale = 100.;
  const float eFexTauRoI_v1::s_xTobEtScale = 25.;
  const float eFexTauRoI_v1::s_towerEtaWidth = 0.1;
 
   eFexTauRoI_v1::eFexTauRoI_v1()
      : SG::AuxElement() {
 
   }
 
   void eFexTauRoI_v1::initialize( unsigned int eFexNumber, unsigned int shelf, uint32_t word0 ) {
      using std::numbers::pi;
 
      // xTOBs will have eFEX and Shelf numbers in word 1 
      // To save space, use the second word of this object, which is not part of a TOB, to store these values
      uint32_t word1 = 0;
      word1 |= (eFexNumber&s_eFexMask)<<s_eFexBit;
      word1 |= (shelf&s_shelfMask)<<s_shelfBit;
      setWord0( word0 );
      setWord1( word1 );

      setEt( etTOB()*s_tobEtScale );
      float etaVal = computeEta();
      setEta( etaVal );
      float phiVal = iPhi() * pi/32. + pi/64.;
      if (phiVal > pi) phiVal = phiVal - 2.*pi;
      setPhi( phiVal );

      if (type() == TOB) setIsTOB(1);
      else               setIsTOB(0);
 
      return;
   }
 
 
   void eFexTauRoI_v1::initialize( uint32_t word0, uint32_t word1 ) {
 
      using std::numbers::pi;
 
      // xTOBs will have eFEX and Shelf numbers in word 1 
      // So all we need to do is set the TOB words
      setWord0( word0 );
      setWord1( word1 );

      setEt( etTOB()*s_tobEtScale );
      float etaVal = computeEta();
      setEta( etaVal );
      float phiVal = iPhi() * pi/32. + pi/64.;
      if (phiVal > pi) phiVal = phiVal - 2.*pi;
      setPhi( phiVal );

      if (type() == TOB) setIsTOB(1);
      else               setIsTOB(0);
 
      return;
   }

   AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( eFexTauRoI_v1, uint32_t, word0,
                                         setWord0 )
   AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( eFexTauRoI_v1, uint32_t, word1,
                                         setWord1 )
 
   
   AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( eFexTauRoI_v1, uint16_t, rCoreNumerator,
                                         setRCoreNumerator )
   AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( eFexTauRoI_v1, uint16_t, rCoreDenominator,
                                         setRCoreDenominator )
   AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( eFexTauRoI_v1, uint16_t, rHadNumerator,
                                         setRHadNumerator )
   AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( eFexTauRoI_v1, uint16_t, rHadDenominator,
                                         setRHadDenominator )

   AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( eFexTauRoI_v1, char, isTOB,
                                         setIsTOB )

   AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( eFexTauRoI_v1, float, et,
                                         setEt )
   AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( eFexTauRoI_v1, float, eta,
                                         setEta )
   AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( eFexTauRoI_v1, float, phi,
                                         setPhi )
 


   unsigned int eFexTauRoI_v1::tobVersion() const {
     return (word0() >> s_versionBit) & s_versionMask;
   }

   unsigned int eFexTauRoI_v1::eFexNumber() const {
     return (word1() >> s_eFexBit) & s_eFexMask;
   }
 
   // Shelf number.  Usually 0 or 1, but can be 12 for Surface Test Facility.
   unsigned int eFexTauRoI_v1::shelfNumber() const {
     return (word1() >> s_shelfBit) & s_shelfMask;
   }

   eFexTauRoI_v1::ObjectType eFexTauRoI_v1::type() const {
     if (etXTOB() == 0) return TOB;
     else              return xTOB;
   }
    
   unsigned int eFexTauRoI_v1::fpga() const {
     return (word0() >> s_fpgaBit) & s_fpgaMask;
   }
    
   unsigned int eFexTauRoI_v1::fpgaEta() const {
     return (word0() >> s_etaBit) & s_etaMask;
   }
    
   unsigned int eFexTauRoI_v1::fpgaPhi() const {
     return (word0() >> s_phiBit) & s_phiMask;
   }
    
   bool eFexTauRoI_v1::upNotDown() const {
     return (word0() >> s_updownBit) & s_updownMask;
   }
    
   unsigned int eFexTauRoI_v1::seed() const {
     return (word0() >> s_seedBit) & s_seedMask;
   }
    
   bool eFexTauRoI_v1::seedMax() const {
     return (word0() >> s_maxBit) & s_maxMask;
   }

   unsigned int eFexTauRoI_v1::etTOB() const {
     // Data content = TOB
     if (etXTOB() == 0) {
       return (word0() >> s_etBit) & s_etMask;
     }
     // Data Content = xTOB. Need to remove lower bits and cap range
     else {
       unsigned int etWord = (etXTOB() >> s_etBitXTOB);
       if (etWord > s_etMask) etWord = s_etMask; 
       return etWord;
     }
   }
   
   unsigned int eFexTauRoI_v1::etXTOB() const {
     return (word1() >> s_etBit) & s_etFullMask;
   }
    
   
   unsigned int eFexTauRoI_v1::bdtScore() const {
     if (tobVersion() == Heuristic || type() != xTOB) return 0;
     return (word0() >> s_bdtScoreBit) & s_bdtScoreMask;
   }
    

   unsigned int eFexTauRoI_v1::rCoreThresholds() const {
     // Return rCore results if used, otherwise return 0
    return (tobVersion() == Heuristic ? tauOneThresholds() : 0);
   }
   
   unsigned int eFexTauRoI_v1::bdtThresholds() const {
     // Return BDT results if used, otherwise return 0
    return (tobVersion() == Heuristic ? 0 : tauOneThresholds());
   }
    
   unsigned int eFexTauRoI_v1::rHadThresholds() const {
     return tauTwoThresholds();
   }
    
   unsigned int eFexTauRoI_v1::tauOneThresholds() const {
     return (word0() >> s_veto1Bit) & s_veto1Mask;
   }
    
   unsigned int eFexTauRoI_v1::tauTwoThresholds() const {
     return (word0() >> s_veto2Bit) & s_veto2Mask;
   }
 
   unsigned int eFexTauRoI_v1::tauThreeThresholds() const {
     return (word0() >> s_veto3Bit) & s_veto3Mask;
   }
    
   unsigned int eFexTauRoI_v1::bcn4() const {
     return (word1() >> s_bcn4Bit) & s_bcn4Mask;
   }

   uint32_t eFexTauRoI_v1::tobWord() const {
     // Do something sensible if called for a TOB
     if (etXTOB() == 0) return word0();
     // When called for xTOB
     else {
       uint32_t word = word0() + etTOB();
       return word;
     }
   }

   float eFexTauRoI_v1::rCore() const {
     if (rCoreDenominator() > 0) return (static_cast<float>(rCoreNumerator())/static_cast<float>(rCoreDenominator()));
     else         return 0.;
   }
 
   float eFexTauRoI_v1::rHad() const {
     if (rHadDenominator() > 0) return (static_cast<float>(rHadNumerator())/static_cast<float>(rHadDenominator()));
     else         return 0.;
   }


  int eFexTauRoI_v1::iPhi() const {

    // only use LSB of shelfNumber to determine the octant
    // (shelfNumber is 0 or 1 in P1, but can be 12 in the STF)
    unsigned int octant = int(eFexNumber()/3) + (shelfNumber()%2)*s_shelfPhiWidth;

    unsigned int index = s_eFexPhiWidth*octant + fpgaPhi() + s_eFexPhiOffset;
    if (index >= s_numPhi) index -= s_numPhi;
 
    return index;
  }

   int eFexTauRoI_v1::iEta() const {

     int index = s_minEta + (eFexNumber()%3)*s_eFexEtaWidth + fpga()*s_fpgaEtaWidth + fpgaEta();

     return index;
 
   }

  int eFexTauRoI_v1::iPhiTopo() const {

     return iPhi()*2;
 
   }

   int eFexTauRoI_v1::iEtaTopo() const {

     return iEta()*4 + seed();
 
   }
 
   float eFexTauRoI_v1::computeEta() const {
     if (tobVersion() == Heuristic) {
       return iEta()*s_towerEtaWidth + (seed()+0.5)*s_towerEtaWidth/4;
     } 
     return iEta()*s_towerEtaWidth + s_towerEtaWidth/2;
   }
 
 
} // namespace xAOD