MCHit.h

Go to the documentation of this file.

#ifndef MCHIT_H
#define MCHIT_H

// C++ includes
#include "lardataobj/MCBase/MCLimits.h"
#include <functional>

namespace sim {

  class MCHit {

  public:
    MCHit() { Reset(); }

    void Reset()
    {
      fSignalTime = fSignalWidth = ::sim::kINVALID_FLOAT;
      fPeakAmp = fCharge = ::sim::kINVALID_FLOAT;
      for (int i = 0; i < 3; ++i)
        fPartVertex[i] = ::sim::kINVALID_FLOAT;
      fPartEnergy = ::sim::kINVALID_FLOAT;
      fPartTrackId = ::sim::kINVALID_INT;
    }

  private:
    //
    // MCHit core information
    //
    float fSignalTime;  
    float fSignalWidth; 

    float fPeakAmp; 
    float fCharge;  

    //
    // Particle information that caused this MCHit
    //

    float fPartVertex[3]; 
    float fPartEnergy;    
    int fPartTrackId;     

  public:
    void SetCharge(float qsum, float amp)
    {
      fCharge = qsum;
      fPeakAmp = amp;
    }

    void SetTime(const float peak, const float width)
    {
      fSignalTime = peak;
      fSignalWidth = width;
    }

    void SetParticleInfo(const float vtx[], const float energy, const int trackId)
    {
      for (size_t i = 0; i < 3; ++i)
        fPartVertex[i] = vtx[i];
      fPartEnergy = energy;
      fPartTrackId = trackId;
    }

    float PeakTime() const { return fSignalTime; }

    float PeakWidth() const { return fSignalWidth; }

    float Charge(bool max = false) const { return (max ? fPeakAmp : fCharge); }

    const float* PartVertex() const { return fPartVertex; }

    float PartEnergy() const { return fPartEnergy; }

    int PartTrackId() const { return fPartTrackId; }

    inline bool operator<(const MCHit& rhs) const { return fSignalTime < rhs.fSignalTime; }

    inline bool operator<(const float& rhs) const { return fSignalTime < rhs; }
  };

}

// Define a pointer comparison
namespace std {
  template <>
  class less<sim::MCHit*> {
  public:
    bool operator()(const sim::MCHit* lhs, const sim::MCHit* rhs) { return (*lhs) < (*rhs); }
  };
}

#endif