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