00001
00006 #ifndef MYHESOURCE_HH
00007 #define MYHESOURCE_HH
00008 #include <math.h>
00009 #include <iostream.h>
00010 #include "ReactorConstants.hh"
00011
00018 class MyHeSource{
00019
00020 public:
00021
00026 MyHeSource(int newIsotope=ISOTOPE,
00027 double newRmaxgen=RMAXGEN);
00028
00030 ~MyHeSource();
00031
00033 MyHeSource(const MyHeSource& HeSource);
00034
00036 MyHeSource& operator=(const MyHeSource& rhs);
00037
00039 double GetHeNeutronEnergy() const {return neutronE;}
00040
00042 double GetHeBetaEnergy() const {return betaE;}
00043
00045 double GetHeGammaEnergy() const {return gammaE;}
00046
00047
00048 static float BetaHighQ(const float& x);
00049 static float BetaMiddleQ(const float& x);
00050 static float BetaLowQ(const float& x);
00051 static double Fermi(const int& Z, const int& A, const float& x);
00052
00054 double GetHeWeight() const {return Weight;}
00055
00057 double GetHeDecayTime() const {return DecayTime;}
00058
00059 private:
00060
00061 int Isotope;
00062 double Rmaxgen;
00063
00064 int Nneutron, Nbeta, Ngamma;
00065
00066 float neutronE, betaE, gammaE;
00067
00068 double Weight;
00069 double DecayTime;
00070 };
00071 #endif