| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // 26 //
>> 27 // $Id: G4StatMFParameters.cc,v 1.5 2008/07/25 11:20:47 vnivanch Exp $
>> 28 // GEANT4 tag $Name: geant4-09-02-patch-03 $
27 // 29 //
28 // Hadronic Process: Nuclear De-excitations 30 // Hadronic Process: Nuclear De-excitations
29 // by V. Lara 31 // by V. Lara
30 32
31 33
32 #include "G4StatMFParameters.hh" 34 #include "G4StatMFParameters.hh"
33 #include "G4SystemOfUnits.hh" <<
34 #include "G4PhysicalConstants.hh" <<
35 35
36 const G4double G4StatMFParameters::fKappa = 1. <<
37 36
38 const G4double G4StatMFParameters::fKappaCoulo << 37 const G4double G4StatMFParameters::_Kappa = 1.0; // dimensionless
39 38
40 const G4double G4StatMFParameters::fEpsilon0 = << 39 const G4double G4StatMFParameters::_KappaCoulomb = 2.0; // dimensionless
>> 40
>> 41 const G4double G4StatMFParameters::_Epsilon0 = 16.0*MeV;
41 42
42 // Bethe-Weizsacker coefficients 43 // Bethe-Weizsacker coefficients
43 const G4double G4StatMFParameters::fE0 = 16.0* << 44 const G4double G4StatMFParameters::_E0 = 16.0*MeV;
44 45
45 const G4double G4StatMFParameters::fBeta0 = 18 << 46 const G4double G4StatMFParameters::_Beta0 = 18.0*MeV;
46 47
47 const G4double G4StatMFParameters::fGamma0 = 2 << 48 const G4double G4StatMFParameters::_Gamma0 = 25.0*MeV;
48 49
49 // Critical temperature (for liquid-gas phase 50 // Critical temperature (for liquid-gas phase transitions)
50 const G4double G4StatMFParameters::fCriticalTe << 51 const G4double G4StatMFParameters::_CriticalTemp = 18.0*MeV;
51 52
52 // Nuclear radius 53 // Nuclear radius
53 const G4double G4StatMFParameters::fr0 = 1.17* << 54 const G4double G4StatMFParameters::_r0 = 1.17*fermi;
54 <<
55 const G4double G4StatMFParameters::fCoulomb = <<
56 (1.0 - 1.0/std::pow(1.0+fKappaCoulomb,1./3.) <<
57 <<
58 G4StatMFParameters::G4StatMFParameters() <<
59 {} <<
60 <<
61 G4StatMFParameters::~G4StatMFParameters() <<
62 {} <<
63 <<
64 G4double G4StatMFParameters::GetKappa() <<
65 { <<
66 return fKappa; <<
67 } <<
68 <<
69 G4double G4StatMFParameters::GetKappaCoulomb() <<
70 { <<
71 return fKappaCoulomb; <<
72 } <<
73 <<
74 G4double G4StatMFParameters::GetEpsilon0() <<
75 { <<
76 return fEpsilon0; <<
77 } <<
78 <<
79 G4double G4StatMFParameters::GetE0() <<
80 { <<
81 return fE0; <<
82 } <<
83 <<
84 G4double G4StatMFParameters::GetBeta0() <<
85 { <<
86 return fBeta0; <<
87 } <<
88 <<
89 G4double G4StatMFParameters::GetGamma0() <<
90 { <<
91 return fGamma0; <<
92 } <<
93 <<
94 G4double G4StatMFParameters::GetCriticalTemp() <<
95 { <<
96 return fCriticalTemp; <<
97 } <<
98 <<
99 G4double G4StatMFParameters::Getr0() <<
100 { <<
101 return fr0; <<
102 } <<
103 <<
104 G4double G4StatMFParameters::GetCoulomb() <<
105 { <<
106 return fCoulomb; <<
107 } <<
108 55
109 G4double G4StatMFParameters::Beta(G4double T) << 56 G4double G4StatMFParameters::Beta(const G4double T)
110 { 57 {
111 G4double res = 0.0; << 58 if (T > _CriticalTemp) return 0.0;
112 if (T < fCriticalTemp) { << 59 else {
113 G4double CriticalTempSqr = fCriticalTemp*f << 60 G4double CriticalTempSqr = _CriticalTemp*_CriticalTemp;
114 G4double TempSqr = T*T; 61 G4double TempSqr = T*T;
115 G4double tmp = (CriticalTempSqr-TempSqr)/( 62 G4double tmp = (CriticalTempSqr-TempSqr)/(CriticalTempSqr+TempSqr);
116 63
117 res = fBeta0*tmp*std::pow(tmp,0.25); << 64 return _Beta0*tmp*std::pow(tmp,1.0/4.0);
118 } 65 }
119 return res; <<
120 } 66 }
121 67
122 G4double G4StatMFParameters::DBetaDT(G4double << 68 G4double G4StatMFParameters::DBetaDT(const G4double T)
123 { 69 {
124 G4double res = 0.0; << 70 if (T > _CriticalTemp) return 0.0;
125 if (T < fCriticalTemp) { << 71 else {
126 G4double CriticalTempSqr = fCriticalTemp*f << 72 G4double CriticalTempSqr = _CriticalTemp*_CriticalTemp;
127 G4double TempSqr = T*T; 73 G4double TempSqr = T*T;
128 G4double tmp = (CriticalTempSqr-TempSqr)/( 74 G4double tmp = (CriticalTempSqr-TempSqr)/(CriticalTempSqr+TempSqr);
129 75
130 res = -5.0*fBeta0*std::pow(tmp,0.25)*(Crit << 76 return -5.0*_Beta0*std::pow(tmp,1.0/4.0)*(CriticalTempSqr*T)/
131 ((CriticalTempSqr+TempSqr)*(CriticalTemp 77 ((CriticalTempSqr+TempSqr)*(CriticalTempSqr+TempSqr));
132 } 78 }
133 return res; <<
134 } 79 }
135 80
136 G4double << 81 G4double G4StatMFParameters::GetMaxAverageMultiplicity(const G4int A)
137 G4StatMFParameters::GetMaxAverageMultiplicity( <<
138 { 82 {
139 // Maximun average multiplicity: M_0 = 2.6 f 83 // Maximun average multiplicity: M_0 = 2.6 for A ~ 200
140 // and M_0 = 3.3 for A <= 110 84 // and M_0 = 3.3 for A <= 110
141 G4double MaxAverageMultiplicity = 2.6; 85 G4double MaxAverageMultiplicity = 2.6;
142 if (A <= 110) { MaxAverageMultiplicity = 3.3 << 86 if (A <= 110) MaxAverageMultiplicity = 3.3;
143 return MaxAverageMultiplicity; 87 return MaxAverageMultiplicity;
144 } 88 }
>> 89
>> 90 G4StatMFParameters G4StatMFParameters::theStatMFParameters;
>> 91
>> 92
>> 93 G4StatMFParameters * G4StatMFParameters::GetAddress()
>> 94 { return &theStatMFParameters; }
145 95
146 96