GLG4HitPhoton.hh

Go to the documentation of this file.

#ifndef __GLG4HitPhoton_hh__
#define __GLG4HitPhoton_hh__
#include <iostream>
 
class GLG4HitPhoton {
 public:
  GLG4HitPhoton() { fPrepulse = false; }
 
  void SetPMTID(int id) { fPMTID = id; }
  void SetTime(double t) { fTime = t; }
  void SetCreationTime(double t) { fCreationTime = t; }
  void SetKineticEnergy(double KE);
  void SetWavelength(double wl);
  void SetPosition(double x, double y, double z);
  void SetMomentum(double x, double y, double z);
  void SetPolarization(double x, double y, double z);
  void SetCount(int count) { fCount = count; }
  void AddCount(int dcount) { fCount += dcount; }
  void SetTrackID(int trackID) { fTrackID = trackID; }
  void SetCreatorProcess(std::string process) { fCreatorProcess = process; }
  void SetPrepulse(bool prepulse) { fPrepulse = prepulse; }
 
  int GetPMTID() const { return fPMTID; }
  double GetTime() const { return fTime; }
  double GetCreationTime() const { return fCreationTime; }
  double GetKineticEnergy() const;
  double GetWavelength() const;
  template <class T>
  inline void GetPosition(T &x, T &y, T &z) const;
  template <class T>
  inline void GetMomentum(T &x, T &y, T &z) const;
  template <class T>
  inline void GetPolarization(T &x, T &y, T &z) const;
  int GetCount() const { return fCount; }
  int GetTrackID() const { return fTrackID; }
  std::string GetCreatorProcess() const { return fCreatorProcess; }
  bool GetPrepulse() const { return fPrepulse; }
 
  void Print(std::ostream &) const;
 
 private:
  double fTime;                 // time of hit
  double fCreationTime;         
  int fPMTID;                   // ID number of PMT the HitPhoton hit
  double fKE;                   // kinetic energy
  double fPosition[3];          // x,y,z components of position
  double fMomentum[3];          // x,y,z components of momentum (normalized?)
  double fPolarization[3];      // x,y,z components of polarization
  int fCount;                   // count of photons, often 1
  int fTrackID;                 // ID number of track which generated this photoelectron
  std::string fCreatorProcess;  // Process that created the photon
  bool fPrepulse;               // if this photon was transmitted through the
                                // the photocathode and a photoelectron was created
                                // at the first dynode, we will choose from the
                                // time and charge distributions for prepulses
};
 
template <class T>
inline void GLG4HitPhoton::GetPosition(T &x, T &y, T &z) const {
  x = fPosition[0];
  y = fPosition[1];
  z = fPosition[2];
}
 
template <class T>
inline void GLG4HitPhoton::GetMomentum(T &x, T &y, T &z) const {
  x = fMomentum[0];
  y = fMomentum[1];
  z = fMomentum[2];
}
 
template <class T>
inline void GLG4HitPhoton::GetPolarization(T &x, T &y, T &z) const {
  x = fPolarization[0];
  y = fPolarization[1];
  z = fPolarization[2];
}
 
inline bool Compare_HitPhotonPtr_TimeAscending(const GLG4HitPhoton *a, const GLG4HitPhoton *b) {
  return a->GetTime() < b->GetTime();
}
 
#endif  // __GLG4HitPhoton_hh__