Shower.cxx

Go to the documentation of this file.

//
// \brief Definition of shower object for LArSoft
//
// \author brebel@fnal.gov
//

#include "lardataobj/RecoBase/Shower.h"

#include <utility> // std::move()
#include <iomanip>
#include <iostream>

namespace recob{
  
  // the value is assigned in the class definition:
  constexpr double Shower::InvalidOpeningAngle;

  //----------------------------------------------------------------------
  Shower::Shower()
  {
  }

  //----------------------------------------------------------------------
  Shower::Shower(TVector3 const& dcosVtx,
           TVector3 const& dcosVtxErr,
           TVector3 const& xyz,
           TVector3 const& xyzErr,
           std::vector< double >  TotalEnergy,
           std::vector< double >  TotalEnergyErr,
           std::vector< double >  dEdx,
           std::vector< double >  dEdxErr,
           int bestplane,
           int     id,
           double length,
           double openAngle)
    : fID(id)
    , fDCosStart(dcosVtx)
    , fSigmaDCosStart(dcosVtxErr)
    , fXYZstart(xyz)
    , fSigmaXYZstart(xyzErr)
    , fTotalEnergy(std::move(TotalEnergy))
    , fSigmaTotalEnergy(std::move(TotalEnergyErr))
    , fdEdx(std::move(dEdx))
    , fSigmadEdx(std::move(dEdxErr))
    , fBestPlane(bestplane)
    , fLength(length)
    , fOpenAngle(openAngle)
  {
   
  }

  //----------------------------------------------------------------------
  std::ostream& operator<< (std::ostream& o, Shower const& a)
  {
    o << std::setiosflags(std::ios::fixed) << std::setprecision(3);
    o << " Shower ID "        << std::setw(4)  << std::right << a.ID();
    o << " Energy    "        << std::setw(4)  << std::right << a.Energy()[a.best_plane()];
    o << " dEdx    "         << std::setw(4)  << std::right << a.dEdx()[a.best_plane()];
    return o;
  }

  //----------------------------------------------------------------------
  // < operator.
  //
  bool operator < (const Shower & a, const Shower & b)
  {
    if(a.ID() != b. ID())
      return a.ID() < b.ID();

    return false; //They are equal
  }


}