Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer << 3 // * DISCLAIMER * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers * 6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which * 7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: * 8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license * 9 // * include a list of copyright holders. << 10 // * 9 // * * 11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file << 14 // * use. * 16 // * for the full disclaimer and the limitatio << 17 // * 15 // * * 18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the * 19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. * 20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work * 21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this * 22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. * 23 // * acceptance of all terms of the Geant4 Sof << 24 // ******************************************* 21 // ******************************************************************** 25 // 22 // 26 // 23 // >> 24 // $Id: G4StatMFMacroTetraNucleon.cc,v 1.4 2005/06/04 13:27:49 jwellisc Exp $ >> 25 // GEANT4 tag $Name: geant4-08-00-patch-01 $ 27 // 26 // 28 // Hadronic Process: Nuclear De-excitations 27 // Hadronic Process: Nuclear De-excitations 29 // by V. Lara 28 // by V. Lara 30 29 31 #include "G4StatMFMacroTetraNucleon.hh" 30 #include "G4StatMFMacroTetraNucleon.hh" 32 #include "G4PhysicalConstants.hh" << 33 #include "G4SystemOfUnits.hh" << 34 #include "G4Log.hh" << 35 #include "G4Exp.hh" << 36 #include "G4Pow.hh" << 37 << 38 G4StatMFMacroTetraNucleon::G4StatMFMacroTetraN << 39 : G4VStatMFMacroCluster(4) << 40 {} << 41 << 42 G4StatMFMacroTetraNucleon::~G4StatMFMacroTetra << 43 {} << 44 << 45 G4double << 46 G4StatMFMacroTetraNucleon::CalcMeanMultiplicit << 47 const G4double mu, << 48 const G4double nu, << 49 const G4double T) << 50 { << 51 G4double ThermalWaveLenght = 16.15*fermi/std << 52 G4double lambda3 = ThermalWaveLenght*Thermal << 53 static const G4double degeneracy = 1; // He << 54 31 55 //old value was 30.11*MeV << 32 // Copy constructor 56 G4double BindingE = G4NucleiProperties::GetB << 33 G4StatMFMacroTetraNucleon:: >> 34 G4StatMFMacroTetraNucleon(const G4StatMFMacroTetraNucleon & ) : >> 35 G4VStatMFMacroCluster(0) // Beacuse the def. constr. of base class is private >> 36 { >> 37 throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMacroTetraNucleon::copy_constructor meant to not be accessable"); >> 38 } >> 39 >> 40 // Operators >> 41 >> 42 G4StatMFMacroTetraNucleon & G4StatMFMacroTetraNucleon:: >> 43 operator=(const G4StatMFMacroTetraNucleon & ) >> 44 { >> 45 throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMacroTetraNucleon::operator= meant to not be accessable"); >> 46 return *this; >> 47 } >> 48 >> 49 >> 50 G4bool G4StatMFMacroTetraNucleon::operator==(const G4StatMFMacroTetraNucleon & ) const >> 51 { >> 52 throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMacroTetraNucleon::operator== meant to not be accessable"); >> 53 return false; >> 54 } >> 55 >> 56 >> 57 G4bool G4StatMFMacroTetraNucleon::operator!=(const G4StatMFMacroTetraNucleon & ) const >> 58 { >> 59 throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMacroTetraNucleon::operator!= meant to not be accessable"); >> 60 return true; >> 61 } >> 62 >> 63 >> 64 >> 65 G4double G4StatMFMacroTetraNucleon::CalcMeanMultiplicity(const G4double FreeVol, const G4double mu, >> 66 const G4double nu, const G4double T) >> 67 { >> 68 const G4double ThermalWaveLenght = 16.15*fermi/std::sqrt(T); >> 69 >> 70 const G4double lambda3 = ThermalWaveLenght*ThermalWaveLenght*ThermalWaveLenght; >> 71 >> 72 const G4double degeneracy = 1; // He4 >> 73 >> 74 const G4double Coulomb = (3./5.)*(elm_coupling/G4StatMFParameters::Getr0())* >> 75 (1.0 - 1.0/std::pow(1.0+G4StatMFParameters::GetKappaCoulomb(),1./3.)); >> 76 >> 77 const G4double BindingE = G4NucleiPropertiesTable::GetBindingEnergy(2,theA); //old value was 30.11*MeV 57 78 58 G4double exponent = (BindingE + theA*(mu+nu* << 79 G4double exponent = (BindingE + theA*(mu+nu*theZARatio+T*T/_InvLevelDensity) - 59 - G4StatMFParameters::GetCoulomb()* << 80 Coulomb*theZARatio*theZARatio*std::pow(static_cast<G4double>(theA),5./3.))/T; 60 theA*G4Pow::GetInstance()->Z23(theA << 81 if (exponent > 700.0) exponent = 700.0; 61 if (exponent > 300.0) exponent = 300.0; << 62 82 63 _MeanMultiplicity = ( degeneracy*FreeVol*the << 83 _MeanMultiplicity = ( degeneracy*FreeVol* static_cast<G4double>(theA)* 64 G4Exp(exponent); << 84 std::sqrt(static_cast<G4double>(theA))/lambda3)* >> 85 std::exp(exponent); 65 86 66 return _MeanMultiplicity; << 87 return _MeanMultiplicity; 67 } 88 } 68 89 >> 90 69 G4double G4StatMFMacroTetraNucleon::CalcEnergy 91 G4double G4StatMFMacroTetraNucleon::CalcEnergy(const G4double T) 70 { 92 { 71 return _Energy = -G4NucleiProperties::GetBin << 93 const G4double Coulomb = (3./5.)*(elm_coupling/G4StatMFParameters::Getr0())* 72 G4StatMFParameters::GetCoulomb() * theZARa << 94 (1.0 - 1.0/std::pow(1.0+G4StatMFParameters::GetKappaCoulomb(),1./3.)); 73 theA*G4Pow::GetInstance()->Z23(theA) + << 95 74 1.5 * T + theA * T*T/_InvLevelDensity; << 96 return _Energy = -G4NucleiPropertiesTable::GetBindingEnergy(2,theA) + 75 } << 97 Coulomb * theZARatio * theZARatio * std::pow(static_cast<G4double>(theA),5./3.) + 76 << 98 (3./2.) * T + 77 G4double << 99 theA * T*T/_InvLevelDensity; 78 G4StatMFMacroTetraNucleon::CalcEntropy(const G << 100 79 const G4double FreeVol) << 101 } 80 { << 102 81 G4double Entropy = 0.0; << 103 82 if (_MeanMultiplicity > 0.0) { << 104 83 G4double ThermalWaveLenght = 16.15*fermi/s << 105 G4double G4StatMFMacroTetraNucleon::CalcEntropy(const G4double T, const G4double FreeVol) 84 G4double lambda3 = ThermalWaveLenght*Therm << 106 { 85 Entropy = _MeanMultiplicity*(2.5 + G4Log(8 << 107 const G4double ThermalWaveLenght = 16.15*fermi/std::sqrt(T); 86 8.0*T/_InvLevelDensity; << 108 const G4double lambda3 = ThermalWaveLenght*ThermalWaveLenght*ThermalWaveLenght; 87 } << 109 88 return Entropy; << 110 G4double Entropy = 0.0; >> 111 if (_MeanMultiplicity > 0.0) >> 112 Entropy = _MeanMultiplicity*(5./2.+ >> 113 std::log(8.0*FreeVol/(lambda3*_MeanMultiplicity)))+ // 8 = theA*std::sqrt(theA) >> 114 8.0*T/_InvLevelDensity; >> 115 >> 116 return Entropy; 89 } 117 } 90 118