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=k.Isotope,
00027 int newRmaxgen=k.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 static ReactorConstants k;
00063
00064 int Isotope;
00065 int Rmaxgen;
00066
00067 int Nneutron, Nbeta, Ngamma;
00068
00069 float neutronE, betaE, gammaE;
00070
00071 double Weight;
00072 double DecayTime;
00073 };
00074 #endif