CfSource.hh

Go to the documentation of this file.

 
#ifndef __RAT_CfSource__
#define __RAT_CfSource__
 
#include <CLHEP/Vector/LorentzVector.h>
 
#include <cmath>
#include <iostream>
#include <vector>
 
namespace RAT {
 
class CfSource {
 public:
  CfSource(int newIsotope = 252);
 
  ~CfSource();
 
  CfSource(const CfSource &CfSource);
 
  CfSource &operator=(const CfSource &rhs);
 
  int GetNumNeutron() const { return Nneutron; }
 
  int GetNumGamma() const { return Ngamma; }
 
  CLHEP::HepLorentzVector GetCfNeutronMomentum(int n) const { return neutronE[n]; }
  double GetCfNeutronTime(int n) const { return Tneutron[n]; }
 
  CLHEP::HepLorentzVector GetCfGammaMomentum(int n) const { return gammaE[n]; }
  double GetCfGammaTime(int n) const { return Tgamma[n]; }
 
 private:
  int Isotope;
 
  static const int maxNeutron = 8;
  static const int maxGamma = 25;
  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 Cf252NeutronSpectrum(const double &x);
 
  static double Cf252GammaMultiplicity(const int &x);
 
  static double Cf252GammaMultiplicityFit(const double &x);
 
  static double Cf252GammaSpectrum(const double &x);
};
 
}  // namespace RAT
 
#endif