00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033 #include "G4StatMFParameters.hh"
00034 #include "G4SystemOfUnits.hh"
00035
00036 const G4double G4StatMFParameters::_Kappa = 1.0;
00037
00038 const G4double G4StatMFParameters::_KappaCoulomb = 2.0;
00039
00040 const G4double G4StatMFParameters::_Epsilon0 = 16.0*MeV;
00041
00042
00043 const G4double G4StatMFParameters::_E0 = 16.0*MeV;
00044
00045 const G4double G4StatMFParameters::_Beta0 = 18.0*MeV;
00046
00047 const G4double G4StatMFParameters::_Gamma0 = 25.0*MeV;
00048
00049
00050 const G4double G4StatMFParameters::_CriticalTemp = 18.0*MeV;
00051
00052
00053 const G4double G4StatMFParameters::_r0 = 1.17*fermi;
00054
00055 G4double G4StatMFParameters::Beta(const G4double T)
00056 {
00057 if (T > _CriticalTemp) return 0.0;
00058 else {
00059 G4double CriticalTempSqr = _CriticalTemp*_CriticalTemp;
00060 G4double TempSqr = T*T;
00061 G4double tmp = (CriticalTempSqr-TempSqr)/(CriticalTempSqr+TempSqr);
00062
00063 return _Beta0*tmp*std::pow(tmp,1.0/4.0);
00064 }
00065 }
00066
00067 G4double G4StatMFParameters::DBetaDT(const G4double T)
00068 {
00069 if (T > _CriticalTemp) return 0.0;
00070 else {
00071 G4double CriticalTempSqr = _CriticalTemp*_CriticalTemp;
00072 G4double TempSqr = T*T;
00073 G4double tmp = (CriticalTempSqr-TempSqr)/(CriticalTempSqr+TempSqr);
00074
00075 return -5.0*_Beta0*std::pow(tmp,1.0/4.0)*(CriticalTempSqr*T)/
00076 ((CriticalTempSqr+TempSqr)*(CriticalTempSqr+TempSqr));
00077 }
00078 }
00079
00080 G4double G4StatMFParameters::GetMaxAverageMultiplicity(const G4int A)
00081 {
00082
00083
00084 G4double MaxAverageMultiplicity = 2.6;
00085 if (A <= 110) MaxAverageMultiplicity = 3.3;
00086 return MaxAverageMultiplicity;
00087 }
00088
00089 G4StatMFParameters G4StatMFParameters::theStatMFParameters;
00090
00091
00092 G4StatMFParameters * G4StatMFParameters::GetAddress()
00093 { return &theStatMFParameters; }
00094