MCDataHolder.h

Go to the documentation of this file.

#ifndef MCDATAHOLDER_H
#define MCDATAHOLDER_H

#include "lardataobj/MCBase/MCLimits.h"
#include <algorithm>
#include <array>

namespace sim {

  class MCEnDep {

  public:
    MCEnDep() { Reset(); }
    ~MCEnDep() {}

  protected:
    std::array<float, 3> fVertex;
    float fEnergy;
    unsigned int fTrackId;

  public:
    void Reset()
    {
      fVertex.fill(::sim::kINVALID_FLOAT);
      fEnergy = ::sim::kINVALID_FLOAT;
      fTrackId = ::sim::kINVALID_UINT;
    }

    void SetVertex(float x, float y, float z)
    {
      fVertex.at(0) = x;
      fVertex.at(1) = y;
      fVertex.at(2) = z;
    }

    void SetEnergy(float e) { fEnergy = e; }

    void SetTrackId(unsigned int id) { fTrackId = id; }

    const std::array<float, 3>& Vertex() const { return fVertex; }

    float Energy() const { return fEnergy; }

    unsigned int TrackId() const { return fTrackId; }

    inline bool operator<(const MCEnDep& e) const
    {
      if (fTrackId < e.fTrackId) return true;
      if (fTrackId > e.fTrackId) return false;

      if (fVertex.at(0) < e.fVertex.at(0)) return true;
      if (fVertex.at(0) > e.fVertex.at(0)) return false;

      if (fVertex.at(1) < e.fVertex.at(1)) return true;
      if (fVertex.at(1) > e.fVertex.at(1)) return false;

      if (fVertex.at(2) < e.fVertex.at(2)) return true;
      if (fVertex.at(2) > e.fVertex.at(2)) return false;

      if (fEnergy < e.fEnergy) return true;
      if (fEnergy > e.fEnergy) return false;

      return false;
    }
  };
}

namespace std {
  template <>
  class less<sim::MCEnDep*> {
  public:
    bool operator()(const sim::MCEnDep* lhs, const sim::MCEnDep* rhs) { return (*lhs) < (*rhs); }
  };
}

#endif