---

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=k.Isotope,
             int newRmaxgen=k.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:

  static ReactorConstants k;

  int Isotope;
  int Rmaxgen;

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

  int Nsources;

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

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

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

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

---