eTauTob.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
#include "eTauTob.h"
 
namespace GlobalSim {
 
  eTauTob::eTauTob(const GepAlgoHypothesisPortsIn& ports_in) {
 
    const auto& w_tob = ports_in.m_I_eTauTobs;
 
    std::size_t i;
    std::size_t j;
 
    for (i = 0; i != AlgoConstants::eFexEtBitWidth; ++i) {
      m_Et[i] = (*w_tob)[i];
    }
 
    for (i = 16, j=0;
     i != 16+ AlgoConstants::eFexPhiBitWidth;
     ++i, ++j) {
      m_Phi[j] = (*w_tob)[i];
    }
 
    for (i = 32, j=0;
     i != 32+ AlgoConstants::eFexEtaBitWidth;
     ++i, ++j) {
      m_Eta[j] = (*w_tob)[i];
    }
 
 
    for (i = 48, j=0;
     i != 48+ AlgoConstants::eFexDiscriminantBitWidth;
     ++i, ++j)
      {
    m_RCore[j] = (*w_tob)[i];
      }
    
 
    for (i = 48+ AlgoConstants::eFexDiscriminantBitWidth, j=0;
     i != 48 + 2*AlgoConstants::eFexDiscriminantBitWidth;
     ++i, ++j)
      {
    m_RHad[j] = (*w_tob)[i];
      }
 
    
    m_Overflow[0] = (*w_tob)[63];
 
  }
 
 
  std::bitset<32>  eTauTob::as_bits() const  {
    auto result = std::bitset<32>();
    std::size_t r_ptr{0};
 
    for (std::size_t s_ptr=0; s_ptr < m_Et.size(); ++s_ptr, ++r_ptr) {
      if(m_Et.test(s_ptr)) {result.set(r_ptr);}
    }
    
    for (std::size_t s_ptr=0; s_ptr < m_RCore.size(); ++s_ptr, ++r_ptr) {
      if(m_RCore.test(s_ptr)) {result.set(r_ptr);}
    }
 
    for (std::size_t s_ptr=0; s_ptr < m_RHad.size(); ++s_ptr, ++r_ptr) {
      if(m_RHad.test(s_ptr)) {result.set(r_ptr);}
    }
 
    for (std::size_t s_ptr=0; s_ptr < m_Eta.size(); ++s_ptr, ++r_ptr) {
      if(m_Eta.test(s_ptr)) {result.set(r_ptr);}
    }
 
    for (std::size_t s_ptr=0; s_ptr < m_Phi.size(); ++s_ptr, ++r_ptr) {
      if(m_Phi.test(s_ptr)) {result.set(r_ptr);}
    }
 
    if(m_Overflow.test(0)) {result.set(r_ptr++);}
    
 
 
    return result;
  }
 
  using eTauTobPtr = std::shared_ptr<eTauTob>;
 
 
  std::ostream& operator << (std::ostream& os, const GlobalSim::eTauTob& tob) {
 
    os << "GlobalSim::eTauTob\n"
       << "Et: " << std::hex << tob.Et_bits().to_ulong() << '\n'
       << "RCore: " << std::hex <<  tob.RCore_bits().to_ulong() << '\n'
       << "RHad: " << std::hex<< tob.RHad_bits().to_ulong() << '\n'
       << "Eta: " << std::hex<< tob.Eta_bits().to_ulong() << '\n'
       << "Phi: " << std::hex<< tob.Phi_bits().to_ulong() << '\n'
       << tob.as_bits() << ' ' << std::hex << tob.as_bits().to_ulong()
       << '\n';
    return os;
  }
}