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 <type_traits> // std::decay<>
#include <utility>     // std::forward<>()

// --- for the implementation ---
// LArSoft libraries
#include "lardata/ArtDataHelper/Dumpers/NewLine.h"
#include "lardata/ArtDataHelper/Dumpers/hexfloat.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