AmBeSource.hh

Go to the documentation of this file.

 
#ifndef __RAT_AmBeSource__
#define __RAT_AmBeSource__
 
#include <CLHEP/Vector/LorentzVector.h>
 
#include <cmath>
#include <iostream>
#include <vector>
 
namespace RAT {
 
class AmBeSource {
 public:
  AmBeSource();
 
  ~AmBeSource();
 
  AmBeSource(const AmBeSource &AmBeSource);
 
  AmBeSource &operator=(const AmBeSource &rhs);
 
  int GetNumNeutron() const { return Nneutron; }
 
  int GetNumGamma() const { return Ngamma; }
 
  CLHEP::HepLorentzVector GetAmBeNeutronMomentum(int n) const { return neutronE[n]; }
  double GetAmBeNeutronTime(int n) const { return Tneutron[n]; }
 
  CLHEP::HepLorentzVector GetAmBeGammaMomentum(int n) const { return gammaE[n]; }
  double GetAmBeGammaTime(int n) const { return Tgamma[n]; }
 
 private:
  static const int maxNeutron = 1;
  static const int maxGamma = 1;
  int Nneutron, Ngamma;
 
  std::vector<CLHEP::HepLorentzVector> neutronE;
  std::vector<double> Tneutron;
 
  std::vector<CLHEP::HepLorentzVector> gammaE;
  std::vector<double> Tgamma;
 
  // G4 particle definitions.
  static double massNeutron;
 
  static double AmBeNeutronSpectrum(const double &x);
};
 
}  // namespace RAT
 
#endif