gFexJetRoI_v1.cxx

Go to the documentation of this file.

 
/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
 
 
// System include(s):
#include <stdexcept>
 
// xAOD include(s):
#include "xAODCore/AuxStoreAccessorMacros.h"
 
// Local include(s):
#include "xAODTrigger/versions/gFexJetRoI_v1.h"
 
namespace xAOD {
 
  const float gFexJetRoI_v1::s_PhiWidth = (2*M_PI)/32; //In central region, gFex has 32 bins in phi
  const std::vector<float> gFexJetRoI_v1::s_EtaEdge     = { -4.9, -4.45, -4.0, -3.5, -3.3, -3.1, 
                                                            -2.9, -2.7, -2.5, -2.2, -2.0, -1.8, -1.6, -1.4, -1.2, -1.0,  
                                                            -0.8, -0.6, -0.4, -0.2, 0.0, 0.2, 0.4, 0.6, 0.8, 1.0,                                                 
                                                            1.2, 1.4, 1.6, 1.8, 2.0, 2.2, 2.5, 2.7, 2.9,
                                                            3.1, 3.3, 3.5, 4.0, 4.45, 4.9 }; //Indices 0-3 and 37-40 are niot currently used in hardware (don't receive energy info)
  
  const std::vector<float> gFexJetRoI_v1::s_EtaCenter   = { -4.7, -4.2, -3.7, -3.4, -3.2, -3, 
                                                            -2.8, -2.6, -2.35, -2.1, -1.9, -1.7, -1.5, -1.3, -1.1, -0.9,  
                                                            -0.7, -0.5, -0.3, -0.1, 0.1, 0.3, 0.5, 0.7, 0.9, 1.1,                                                 
                                                            1.3, 1.5, 1.7, 1.9, 2.1, 2.35, 2.6, 2.8, 3.0,
                                                            3.2, 3.4, 3.7, 4.2, 4.7};
 
 
  //vector<float> gFexJetRoI_v1::s_maxEta = {}; 
 
  gFexJetRoI_v1::gFexJetRoI_v1()
    : SG::AuxElement() {
 
  }
 
  void gFexJetRoI_v1::initialize( uint32_t word, int tobEtScale ) {
 
    setWord( word );
    setScale( tobEtScale );
    setgFexType(unpackType());
    setTobEt( unpackEt() );
    setEta( unpackEtaIndex() );
    setPhi( unpackPhiIndex() );
    setStatus( unpackStatus());
    setSaturated(unpackSaturated());
 
  }
 
    int gFexJetRoI_v1::menuEta() const {
      return static_cast<int>(eta()*40)/4;
    }
 
 
  AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( gFexJetRoI_v1, uint32_t, word,
                                         setWord )
  AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( gFexJetRoI_v1, int, tobEtScale,
                                         setScale )
  AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( gFexJetRoI_v1, uint8_t, saturated,
                                         setSaturated )
 
   
  AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( gFexJetRoI_v1, int, gFexType,
                                         setgFexType )
  AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( gFexJetRoI_v1, int16_t, gFexTobEt,
                                         setTobEt )
  AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( gFexJetRoI_v1, uint8_t, iEta,
                                         setEta )
  AUXSTORE_PRIMITIVE_SETTER_AND_GETTER( gFexJetRoI_v1, uint8_t, iPhi,
                                         setPhi )
 
  uint8_t gFexJetRoI_v1::status() const {
  static const Accessor< uint8_t > acc( "gFexJetStatus" );
  return acc( *this );                                   
  }
  void gFexJetRoI_v1::setStatus ( uint8_t value ) {
     static const Accessor< uint8_t > acc( "gFexJetStatus" );
     acc( *this ) = value;
    return;
  }
 
 
  // Methods to decode data from the TOB/RoI and return to the user
 
 
  // Object disambiguation ()
  int gFexJetRoI_v1::unpackType() const {
    auto tobID = (word() >> s_tobIDBit) & s_tobIDMask;
    if (tobID == 0 ){
      return gRho;
    }  
    else if (tobID == 1 || tobID == 2) {
      return gBlockLead;
    } 
    else if (tobID == 3 || tobID == 4 ) {
      return gBlockSub;
    }      
    else if (tobID == 5 || tobID == 6) {
      return gJet;
    }  
    else return -999;
  }
 
  bool gFexJetRoI_v1::isgBlockLead() const {
    return gFexType() == gBlockLead;
  }
 
  bool gFexJetRoI_v1::isgBlockSub() const {
    return gFexType() == gBlockSub;
  }
 
  bool gFexJetRoI_v1::isgJet() const {
    return gFexType() == gJet;
  }
 
  bool gFexJetRoI_v1::isgRho() const {
    return gFexType() == gRho;
  }
 
 
  unsigned int gFexJetRoI_v1::unpackStatus() const{
    return (word() >> s_statusBit) & s_statusMask;
  }
 
  unsigned int gFexJetRoI_v1::unpackSaturated() const{
    return (word() >> s_saturBit) & s_saturMask;
  }
    
 
  // Raw ET on TOB scale 
  int16_t gFexJetRoI_v1::unpackEt() const {
  // Data content = TOB
    int16_t energy = (word() >> s_etBit) & s_etMask; 
    return energy; 
  }
 
  // Return an eta index in the range 0-49
  unsigned int gFexJetRoI_v1::unpackEtaIndex() const {
    return (word() >> s_etaBit) & s_etaMask;
  }
   
  // Return an eta index in the range 0-32
  unsigned int gFexJetRoI_v1::unpackPhiIndex() const {
    return (word() >> s_phiBit) & s_phiMask;
 
  }
 
  // Methods that require combining results or applying scales
  
  float gFexJetRoI_v1::et() const {
      return gFexTobEt()*tobEtScale();
  }
 
 
  // Floating point coordinates. Return the center of Eta. 
  float gFexJetRoI_v1::eta() const {
    if (gFexType() != gRho){
      return s_EtaCenter[iEta()];
    } 
    return 0;
   }
 
  float gFexJetRoI_v1::etaMin() const {
    if (gFexType() != gRho){
      return s_EtaEdge[iEta()];
    } 
    return 0;
  }
   
  float gFexJetRoI_v1::etaMax() const {
    if (gFexType() != gRho){
      return s_EtaEdge[iEta()+1];
    } 
    return 0;
  }
 
 
 
  // Floating point coordinates using gFex convention [0, 2pi].
  // Returns the centre phi coordinate.
  float gFexJetRoI_v1::phi_gFex() const {
    float phi_out = 0;
    if (gFexType() != gRho){
        phi_out = (iPhi() * s_PhiWidth) + (s_PhiWidth/2);
    } 
    return phi_out; 
  }
 
 
  // Floating point coordinates using gFex convention [0, 2pi]. 
  // Returns the low edge of the phi bin.
  float gFexJetRoI_v1::phiMin_gFex() const {
    float phi_out = 0;
    if (gFexType() != gRho){
        phi_out = iPhi() * s_PhiWidth;
    } 
    return phi_out; 
  }
 
  // Floating point coordinates using gFex convention [0, 2pi].
  // Returns the high edge of the phi bin.
  float gFexJetRoI_v1::phiMax_gFex() const {
    float phi_out = 0;
    if (gFexType() != gRho){
        phi_out = iPhi() * s_PhiWidth + s_PhiWidth;
    } 
    return phi_out; 
  }
 
  // Floating point coordinates using ATLAS convention [-pi, pi]. 
  // Returns the centre phi coordinate.
  float gFexJetRoI_v1::phi() const {
    float phi_out = 0;
    if (gFexType() != gRho){
      if (phi_gFex() < M_PI) phi_out = phi_gFex();
      else phi_out = (phi_gFex() - 2*M_PI);
    }
    return phi_out; 
  }
 
  // Floating point coordinates using ATLAS convention [-pi, pi].
  // Returns the high edge of the phi bin. 
  float gFexJetRoI_v1::phiMax() const {
    float phi_out = 0;
    if (gFexType() != gRho){
      if (phiMax_gFex() < M_PI) phi_out = phiMax_gFex();
      else phi_out = (phiMax_gFex() - 2*M_PI);
    }
    return phi_out; 
  }
 
  // Floating point coordinates using ATLAS convention [-pi, pi]. 
  // Returns the low edge of the phi bin.
  float gFexJetRoI_v1::phiMin() const {
    float phi_out = 0;
    if (gFexType() != gRho){
      if (phiMin_gFex() < M_PI) phi_out = phiMin_gFex();
      else phi_out = (phiMin_gFex() - 2*M_PI);
    }
    return phi_out; 
  }
 
  // Returns phi index in the range used by L1Topo (0->127).
  int gFexJetRoI_v1::iPhiTopo() const {
    int phi_out = 0;
    if (gFexType() != gRho){
      if (( iEta() <= 3 ) || ( (iEta() >= 36) )){
        phi_out = ( 8 * iPhi() ) + 4;
      }
      else if ( iEta() >3  && iEta() < 36 ){
        phi_out = ( 4 * iPhi() ) + 2;
      } 
      else return -999; 
    } 
    return phi_out; 
  }
 
 
} // namespace xAOD