Event.cxx

Go to the documentation of this file.

//
// \brief Definition of event object for LArSoft
//
// \author brebel@fnal.gov
//
#include "lardataobj/RecoBase/Event.h"

#include "larcoreobj/SimpleTypesAndConstants/PhysicalConstants.h"
#include "messagefacility/MessageLogger/MessageLogger.h"

#include <algorithm> // for std::sort
#include <iomanip>
#include <ostream>

namespace recob {

  //----------------------------------------------------------------------
  Event::Event() : Event(util::kBogusI) {}

  //----------------------------------------------------------------------
  Event::Event(int id) : fID(id) {}

  //----------------------------------------------------------------------
  double Event::Energy() const
  {
    // loop over all vertex objects and get the
    mf::LogWarning("Event") << "Event::Energy() is not yet defined.  Need to decide "
                            << " how to calculate energy of Vertex"
                            << " Return util:kBogusD for now.";

    return util::kBogusD;
  }

  //----------------------------------------------------------------------
  double Event::SigmaEnergy() const
  {
    // loop over all vertex objects and get the
    mf::LogWarning("Event") << "Event::SigmaEnergy() is not yet defined.  Need to decide "
                            << " how to calculate uncertainty in energy of Prong/Vertex"
                            << " Return util:kBogusD for now.";

    return util::kBogusD;
  }

  //----------------------------------------------------------------------
  const recob::Vertex* Event::PrimaryVertex(std::vector<const recob::Vertex*>& vtxs) const
  {
    // sort the vertices, set the first one in z to be the primary
    std::sort(vtxs.begin(), vtxs.end());
    return vtxs.front();
  }

  //----------------------------------------------------------------------
  // ostream operator.
  //
  std::ostream& operator<<(std::ostream& o, const Event& a)
  {

    o << std::setprecision(5);
    o << "Event " << a.fID << std::setw(5) << " Energy = " << a.Energy() << " +/- "
      << a.SigmaEnergy() << std::endl;

    return o;
  }

  //----------------------------------------------------------------------
  // < operator.
  //
  bool operator<(const Event& a, const Event& b)
  {
    return a.Energy() < b.Energy();
  }

}