---

MyCfSource.hh

Go to the documentation of this file.

#ifndef MYCFSOURCE_HH
#define MYCFSOURCE_HH
#include <math.h>
#include <iostream.h>
#include "ReactorConstants.hh"

class MyCfSource{

public:

  MyCfSource(int newIsotope=ISOTOPE,
             double newRmaxgen=RMAXGEN);

  ~MyCfSource();

  MyCfSource(const MyCfSource& CfSource);

  MyCfSource& operator=(const MyCfSource& rhs);

  static float Cf252NeutronSpectrum(const float& x);

  static float Cf252GammaMultiplicity(const int& x);

  static float Cf252GammaMultiplicityFit(const float& x);

  static float Cf252GammaSpectrum(const float& x);

  int GetNumCfSources() const {return Nsources;}

  int GetNumNeutron(int& n) const {return Nneutron[n];}

  int GetNumGamma(int& n) const {return Ngamma[n];}

  double GetCfNeutronEnergy(int& n,int& nn) const {return neutronE[n][nn];}
  double GetCfNeutronTime(int& n,int& nn) const {return Tneutron[n][nn];}

  double GetCfGammaEnergy(int& n,int& nn) const {return gammaE[n][nn];}
  double GetCfGammaTime(int& n,int& nn) const {return Tgamma[n][nn];}

  double GetCfVertexR(int& n) const {return R[n];}

  double GetCfVertexCosTheta(int& n) const {return zR[n];}

  double GetCfVertexPhi(int& n) const {return phiR[n];}

private:

  int Isotope;
  double Rmaxgen;

  static const int Sdim = 1121; // entries in cf252_neutron_spectrum.txt
  float neutronEnergy[Sdim];
  float neutronProb[Sdim];

  int Nsources;

  double R[Ndim],zR[Ndim],phiR[Ndim];

  static const int maxNeutron = 8;
  static const int maxGamma = 25;
  int Nneutron[Ndim], Ngamma[Ndim];

  float neutronE[Ndim][maxNeutron];
  float Tneutron[Ndim][maxNeutron];

  float gammaE[Ndim][maxGamma];
  float Tgamma[Ndim][maxGamma];
};
#endif

---