SimEnergyDepositLite.h

Go to the documentation of this file.

#ifndef LARDATAOBJ_SIMULATION_SIMENERGYDEPOSITLITE_H
#define LARDATAOBJ_SIMULATION_SIMENERGYDEPOSITLITE_H

// LArSoft includes
// Define the LArSoft standard geometry types and methods.
#include "larcoreobj/SimpleTypesAndConstants/geo_vectors.h"
#include "lardataobj/Simulation/SimEnergyDeposit.h"

// C++ includes
#include <vector>

namespace sim {
  class SimEnergyDepositLite {
  public:
    using Length_t = float;
    using Point_t = geo::Point_t;

    SimEnergyDepositLite(double e = 0.,
                         geo::Point_t middle = {0., 0., 0.},
                         double t = 0.,
                         int id = 0)
      : edep(e), middlePos(middle), middleTime(t), trackID(id)
    {}

    double Energy() const { return edep; }
    geo::Point_t const& Position() const noexcept { return middlePos; }
    geo::Point_t const& MidPoint() const noexcept
    {
      return Position();
    } 
    double Time() const { return middleTime; }
    int TrackID() const { return trackID; }

    geo::Length_t X() const { return middlePos.X(); }
    geo::Length_t Y() const { return middlePos.Y(); }
    geo::Length_t Z() const { return middlePos.Z(); }
    double T() const { return middleTime; }
    double E() const { return edep; }

    bool operator<(const SimEnergyDepositLite& rhs) const
    {
      if (trackID < rhs.trackID) return true;
      if (trackID > rhs.trackID) return false;
      if (middleTime < rhs.middleTime) return true;
      if (middleTime > rhs.middleTime) return false;
      if (middlePos.Z() < rhs.Z()) return true;
      if (middlePos.Z() > rhs.Z()) return false;
      if (middlePos.Y() < rhs.Y()) return true;
      if (middlePos.Y() > rhs.Y()) return false;
      if (middlePos.X() < rhs.X()) return true;
      if (middlePos.X() > rhs.X()) return false;
      return (edep > rhs.edep); // sort by _decreasing_ energy
    }

    // Need to be aware that information is obviously not complete
    // (e.g. no exact start/end point or time, only middle point or time)
    operator sim::SimEnergyDeposit() const
    {
      return sim::SimEnergyDeposit(
        0, 0, 0, edep, middlePos, middlePos, middleTime, middleTime, trackID, 0);
    }

  private:
    float edep;             
    geo::Point_t middlePos; 
    double middleTime;      
    int trackID;            
  };

  typedef std::vector<SimEnergyDepositLite> SimEnergyDepositLiteCollection;
} // namespace sim
#endif // LARDATAOBJ_SIMULATION_SIMENERGYDEPOSITLITE_H