eFexEMRoI_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/eFexEMRoI_v1.h"
 
namespace xAOD {

  const float eFexEMRoI_v1::s_tobEtScale = 100.;
  const float eFexEMRoI_v1::s_xTobEtScale = 25.;
  const float eFexEMRoI_v1::s_towerEtaWidth = 0.1;
 
   eFexEMRoI_v1::eFexEMRoI_v1()
      : SG::AuxElement() {
 
   }
 
   void eFexEMRoI_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 = iEta()*s_towerEtaWidth + (seed()+0.5)*s_towerEtaWidth/4;
      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 eFexEMRoI_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 = iEta()*s_towerEtaWidth + (seed()+0.5)*s_towerEtaWidth/4;
      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( eFexEMRoI_v1, uint32_t, word0,
                                         setWord0 )
   AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( eFexEMRoI_v1, uint32_t, word1,
                                         setWord1 )

   AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( eFexEMRoI_v1, uint16_t, RetaCore,
                                         setRetaCore )
   AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( eFexEMRoI_v1, uint16_t, RetaEnv,
                                         setRetaEnv )
   AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( eFexEMRoI_v1, uint16_t, RhadEM,
                                         setRhadEM )
   AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( eFexEMRoI_v1, uint16_t, RhadHad,
                                         setRhadHad )
   AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( eFexEMRoI_v1, uint16_t, WstotNumerator,
                                         setWstotNumerator )
   AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( eFexEMRoI_v1, uint16_t, WstotDenominator,
                                         setWstotDenominator )
 

   AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( eFexEMRoI_v1, char, isTOB,
                                         setIsTOB )

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


   unsigned int eFexEMRoI_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 eFexEMRoI_v1::shelfNumber() const {
     return (word1() >> s_shelfBit) & s_shelfMask;
   }

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

   unsigned int eFexEMRoI_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 eFexEMRoI_v1::etXTOB() const {
     return (word1() >> s_etBit) & s_etFullMask;
   }
    
   unsigned int eFexEMRoI_v1::RetaThresholds() const {
     return (word0() >> s_veto1Bit) & s_veto1Mask;
   }
    
   unsigned int eFexEMRoI_v1::RhadThresholds() const {
     return (word0() >> s_veto3Bit) & s_veto3Mask;
   }
    
   unsigned int eFexEMRoI_v1::WstotThresholds() const {
     return (word0() >> s_veto2Bit) & s_veto2Mask;
   }
    
   unsigned int eFexEMRoI_v1::bcn4() const {
     return (word1() >> s_bcn4Bit) & s_bcn4Mask;
   }

   uint32_t eFexEMRoI_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 eFexEMRoI_v1::Reta() const {
     float sum = RetaCore() + RetaEnv();
     if (sum > 0) return (1. - RetaCore()/sum);
     else         return -1.;
   }
 
   float eFexEMRoI_v1::Rhad() const {
     float sum = RhadEM() + RhadHad();
     if (sum > 0) return (RhadHad()/sum);
     else         return -1.;
   }
 
   float eFexEMRoI_v1::Wstot() const {
     if (WstotDenominator() > 0) return (static_cast<float>(WstotNumerator())/WstotDenominator());
     else         return -1.;
   }


  int eFexEMRoI_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;

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

   int eFexEMRoI_v1::iEta() const {

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

     return index;
 
   }

  int eFexEMRoI_v1::iPhiTopo() const {

     return iPhi()*2;
 
   }

   int eFexEMRoI_v1::iEtaTopo() const {

     return iEta()*4 + seed();
 
   }
 
 
} // namespace xAOD