SpacePointDumpers.h

Go to the documentation of this file.

#ifndef LARDATA_RECOBASE_DUMPERS_SPACEPOINTDUMPERS_H
#define LARDATA_RECOBASE_DUMPERS_SPACEPOINTDUMPERS_H 1


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

// C/C++ standard libraries
#include <string>
#include <utility> // std::forward<>()
#include <type_traits> // std::decay<>


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


namespace recob {
  namespace dumper {
    
    struct SpacePointPrintOptions_t {
      IndentOptions_t indent; 
      bool hexFloats = false; 
      
      SpacePointPrintOptions_t() = default;
      
      SpacePointPrintOptions_t
        (IndentOptions_t indentOptions, bool bHexFloats)
        : indent(indentOptions), hexFloats(bHexFloats)
        {}
      
    }; // SpacePointPrintOptions_t
    
    
    template
      <typename Stream, typename NewLineRef = recob::dumper::NewLine<Stream>>
    auto DumpSpacePoint(
      Stream&& out,
      recob::SpacePoint const& sp,
      SpacePointPrintOptions_t const& options = {}
      ) -> std::enable_if_t
      <std::is_same<NewLine<std::decay_t<Stream>>, std::decay_t<NewLineRef>>::value>;
    
  } // namespace dumper
} // namespace recob


//==============================================================================
//=== template implementation
//===
//------------------------------------------------------------------------------
//--- recob::dumper::DumpSpacePoint
//---
template <typename Stream, typename NewLineRef>
auto recob::dumper::DumpSpacePoint(
  Stream&& out,
  recob::SpacePoint const& sp,
  SpacePointPrintOptions_t const& options /* = {} */
) -> std::enable_if_t<
  std::is_same<
    NewLine<std::decay_t<Stream>>,
    std::decay_t<NewLineRef>
  >::value>
{
  double const* pos = sp.XYZ();
  double const* err = sp.ErrXYZ();
  
  NewLineRef nl(out, options.indent);
  lar::OptionalHexFloat hexfloat(options.hexFloats);
  
  nl()
    << "ID=" << sp.ID() << " at (" << hexfloat(pos[0])
    << ", " << hexfloat(pos[1]) << ", " << hexfloat(pos[2])
    << ") cm, chi^2/NDF=" << hexfloat(sp.Chisq());
  
  nl()
    << "variances { x^2=" << hexfloat(err[0]) << " y^2=" << hexfloat(err[2])
    << " z^2=" << hexfloat(err[5])
    << " xy=" << hexfloat(err[1]) << " xz=" << hexfloat(err[3])
    << " yz=" << hexfloat(err[4]) << " }";
  
} // recob::dumper::DumpSpacePoint()

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

#endif // LARDATA_RECOBASE_DUMPERS_SPACEPOINTDUMPERS_H