MCHit.h

Go to the documentation of this file.

#ifndef MCHIT_H
#define MCHIT_H

// C++ includes
#include <iostream>
#include <sstream>
#include <algorithm>
#include "MCBaseException.h"
#include "MCLimits.h"

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