PCAxis.cxx

Go to the documentation of this file.

//
// \file  PCAxis.cxx
// \brief Definition of Principal Components Axis object for LArSoft
//
// \author usher@slac.stanford.edu
//
// \see   PCAxis.h
//

#include "lardataobj/RecoBase/PCAxis.h"

#include <iomanip>
#include <ostream>

namespace recob {

  //----------------------------------------------------------------------
  PCAxis::PCAxis() {}

  //----------------------------------------------------------------------
  PCAxis::PCAxis(bool ok,
                 int nHits,
                 const double* eigenValues,
                 const recob::PCAxis::EigenVectors& eigenVecs,
                 const double* avePos,
                 const double aveHitDoca,
                 size_t id)
    : fSvdOK(ok), fNumHitsUsed(nHits), fEigenVectors(eigenVecs), fAveHitDoca(aveHitDoca), fID(id)
  {
    fEigenValues[0] = eigenValues[0];
    fEigenValues[1] = eigenValues[1];
    fEigenValues[2] = eigenValues[2];
    fAvePosition[0] = avePos[0];
    fAvePosition[1] = avePos[1];
    fAvePosition[2] = avePos[2];
  }

  //----------------------------------------------------------------------
  // ostream operator.
  //
  std::ostream& operator<<(std::ostream& o, const PCAxis& a)
  {
    if (a.fSvdOK) {
      o << std::setiosflags(std::ios::fixed) << std::setprecision(2);
      o << " PCAxis ID " << a.fID << " run with " << a.fNumHitsUsed << " space points" << std::endl;
      o << "   - center position: " << std::setw(6) << a.fAvePosition[0] << ", "
        << a.fAvePosition[1] << ", " << a.fAvePosition[2] << std::endl;
      o << "   - eigen values: " << std::setw(8) << std::right << a.fEigenValues[0] << ", "
        << a.fEigenValues[1] << ", " << a.fEigenValues[2] << std::endl;
      o << "   - average doca: " << a.fAveHitDoca << std::endl;
      o << "   - Principle axis: " << std::setw(7) << std::setprecision(4) << a.fEigenVectors[0][0]
        << ", " << a.fEigenVectors[0][1] << ", " << a.fEigenVectors[0][2] << std::endl;
      o << "   - second axis:    " << std::setw(7) << std::setprecision(4) << a.fEigenVectors[1][0]
        << ", " << a.fEigenVectors[1][1] << ", " << a.fEigenVectors[1][2] << std::endl;
      o << "   - third axis:     " << std::setw(7) << std::setprecision(4) << a.fEigenVectors[2][0]
        << ", " << a.fEigenVectors[2][1] << ", " << a.fEigenVectors[2][2] << std::endl;
    }
    else
      o << " Principal Components Axis is not valid" << std::endl;

    return o;
  }

  //----------------------------------------------------------------------------
  // comparison operator
  //
  bool operator<(const PCAxis& a, const PCAxis& b)
  {
    if (a.getID() != b.getID()) return a.getID() < b.getID();

    return false; //They are equal
  }

} // End of namespace