PCAxisDumpers.h

Go to the documentation of this file.

#ifndef LARDATA_RECOBASE_DUMPERS_PCAXISDUMPERS_H
#define LARDATA_RECOBASE_DUMPERS_PCAXISDUMPERS_H 1

// LArSoft libraries
#include "lardataobj/RecoBase/PCAxis.h"

// C/C++ standard libraries
#include <string>
#include <iomanip>
#include <type_traits> // std::decay<>


// --- for the implementation ---
// LArSoft libraries
#include "lardata/ArtDataHelper/Dumpers/NewLine.h"


namespace recob {
  namespace dumper {
    
    template <typename Stream, typename NewLineRef>
    std::enable_if_t
      <std::is_same<recob::dumper::NewLine<std::decay_t<Stream>>, std::decay_t<NewLineRef>>::value>
    DumpPCAxis
      (Stream&& out, recob::PCAxis const& pca, NewLineRef&& nl);
    
    template <typename Stream>
    void DumpPCAxis(Stream&& out, recob::PCAxis const& pca,
      std::string indent = "",
      bool indentFirst = true
      )
      {
        DumpPCAxis(
          std::forward<Stream>(out), pca, makeNewLine(out, indent, !indentFirst)
          );
      }
    
    
  } // namespace dumper
} // namespace lar


//==============================================================================
//=== template implementation
//===
//------------------------------------------------------------------------------
//--- recob::dumper::DumpPCAxis
//---
template <typename Stream, typename NewLineRef>
std::enable_if_t
  <std::is_same<recob::dumper::NewLine<std::decay_t<Stream>>, std::decay_t<NewLineRef>>::value>
recob::dumper::DumpPCAxis
  (Stream&& out, recob::PCAxis const& pca, NewLineRef&& nl)
{
  
  if (!pca.getSvdOK()) {
    nl() << "<not valid>";
    return;
  }

  nl() << std::setiosflags(std::ios::fixed) << std::setprecision(2)
    << " ID " << pca.getID()
    << " run on " << pca.getNumHitsUsed() << " space points";
  nl()
    << "  - center position: " << std::setw(6) << pca.getAvePosition()[0]
    << ", " << pca.getAvePosition()[1]
    << ", " << pca.getAvePosition()[2];
  nl()
    << "  - eigen values: " << std::setw(8) << std::right
    << pca.getEigenValues()[0] << ", "
    << pca.getEigenValues()[1] << ", " << pca.getEigenValues()[2];
  nl()
    << "  - average doca: " << pca.getAveHitDoca();
  nl()
    << "  - principle axis: "
    << std::setw(7) << std::setprecision(4) << pca.getEigenVectors()[0][0]
    << ", " << pca.getEigenVectors()[0][1]
    << ", " << pca.getEigenVectors()[0][2];
  nl()
    << "  - second axis: "
    << std::setw(7) << std::setprecision(4) << pca.getEigenVectors()[1][0]
    << ", " << pca.getEigenVectors()[1][1]
    << ", " << pca.getEigenVectors()[1][2];
  nl()
    << "  - third axis: "
    << std::setw(7) << std::setprecision(4) << pca.getEigenVectors()[2][0]
    << ", " << pca.getEigenVectors()[2][1]
    << ", " << pca.getEigenVectors()[2][2];
  
} // recob::dumper::DumpPCAxis()

//------------------------------------------------------------------------------

#endif // LARDATA_RECOBASE_DUMPERS_PCAXISDUMPERS_H