TrkrClusterv3.cc

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file TrkrClusterv3.cc

#include "TrkrClusterv3.h"

#include <cmath>
#include <utility>          // for swap

namespace
{
  // square convenience function
  template<class T> inline constexpr T square( const T& x ) 
    { return x*x; }
}

TrkrClusterv3::TrkrClusterv3()
  : m_cluskey(TrkrDefs::CLUSKEYMAX)
  , m_subsurfkey(TrkrDefs::SUBSURFKEYMAX)
  , m_adc(0xFFFFFFFF)
{
  for (int i = 0; i < 2; i++)
    {
      m_local[i] = NAN;
      for(int j = 0; j < 2; j ++)
        {
          m_actsLocalErr[i][j] = NAN;
        }
    }
}

void TrkrClusterv3::identify(std::ostream& os) const
{
  os << "---TrkrClusterv3--------------------" << std::endl;

  os << " (rphi,z) =  (" << getLocalX();
  os << ", " << getLocalY() << ") cm ";

  os << " adc = " << getAdc() << std::endl;

  os << "Error " << std::endl;
  os << "         ( ";
  os << getActsLocalError(0, 0) << " , ";
  os << getActsLocalError(0, 1) << " ) " << std::endl;
  os << "         ( ";
  os << getActsLocalError(1, 0) << ", ";
  os << getActsLocalError(1, 1) << " ) " << std::endl;

  os << std::endl;
  os << "-----------------------------------------------" << std::endl;

  return;
}

int TrkrClusterv3::isValid() const
{
  if (m_cluskey == TrkrDefs::CLUSKEYMAX) { return 0; }
  for (int i = 0; i < 2; ++i)
  {
    if (std::isnan(getPosition(i))) { return 0; }
    for (int j = 0; j < 2; ++j)
      { if (std::isnan(getActsLocalError(i,j))) { return 0; }}
  }
  if (m_adc == 0xFFFFFFFF) { return 0; }

  return 1;
}

void TrkrClusterv3::CopyFrom( const TrkrCluster& source )
{
  // do nothing if copying onto oneself
  if( this == &source ) return;
 
  // parent class method
  TrkrCluster::CopyFrom( source );
 
  setLocalX( source.getLocalX() );
  setLocalY( source.getLocalY() );
  
  for (int j = 0; j < 2; ++j)
    for (int i = 0; i < 2; ++i)
  { setActsLocalError(i, j, source.getActsLocalError(i, j)); }
  
  setSubSurfKey( source.getSubSurfKey() );
  setAdc( source.getAdc() );
}

float TrkrClusterv3::getRPhiError() const
{ return std::sqrt(m_actsLocalErr[0][0]); }

float TrkrClusterv3::getZError() const
{ return std::sqrt(m_actsLocalErr[1][1]); }

void TrkrClusterv3::setActsLocalError(unsigned int i, unsigned int j,
                                      float value)
{
  m_actsLocalErr[i][j] = value;
}