Geant4  10.03
G4StatMFParameters.cc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 // $Id: G4StatMFParameters.cc 91834 2015-08-07 07:24:22Z gcosmo $
28 //
29 // Hadronic Process: Nuclear De-excitations
30 // by V. Lara
31 
32 
33 #include "G4StatMFParameters.hh"
34 #include "G4SystemOfUnits.hh"
35 #include "G4PhysicalConstants.hh"
36 
37 const G4double G4StatMFParameters::fKappa = 1.0; // dimensionless
38 
39 const G4double G4StatMFParameters::fKappaCoulomb = 2.0; // dimensionless
40 
42 
43 // Bethe-Weizsacker coefficients
45 
47 
49 
50 // Critical temperature (for liquid-gas phase transitions)
52 
53 // Nuclear radius
55 
56 const G4double G4StatMFParameters::fCoulomb = 0.6*(CLHEP::elm_coupling/fr0)*
57  (1.0 - 1.0/std::pow(1.0+fKappaCoulomb,1./3.));
58 
60 {}
61 
63 {}
64 
66 {
67  return fKappa;
68 }
69 
71 {
72  return fKappaCoulomb;
73 }
74 
76 {
77  return fEpsilon0;
78 }
79 
81 {
82  return fE0;
83 }
84 
86 {
87  return fBeta0;
88 }
89 
91 {
92  return fGamma0;
93 }
94 
96 {
97  return fCriticalTemp;
98 }
99 
101 {
102  return fr0;
103 }
104 
106 {
107  return fCoulomb;
108 }
109 
111 {
112  G4double res = 0.0;
113  if (T < fCriticalTemp) {
114  G4double CriticalTempSqr = fCriticalTemp*fCriticalTemp;
115  G4double TempSqr = T*T;
116  G4double tmp = (CriticalTempSqr-TempSqr)/(CriticalTempSqr+TempSqr);
117 
118  res = fBeta0*tmp*std::pow(tmp,0.25);
119  }
120  return res;
121 }
122 
124 {
125  G4double res = 0.0;
126  if (T < fCriticalTemp) {
127  G4double CriticalTempSqr = fCriticalTemp*fCriticalTemp;
128  G4double TempSqr = T*T;
129  G4double tmp = (CriticalTempSqr-TempSqr)/(CriticalTempSqr+TempSqr);
130 
131  res = -5.0*fBeta0*std::pow(tmp,0.25)*(CriticalTempSqr*T)/
132  ((CriticalTempSqr+TempSqr)*(CriticalTempSqr+TempSqr));
133  }
134  return res;
135 }
136 
137 G4double
139 {
140  // Maximun average multiplicity: M_0 = 2.6 for A ~ 200
141  // and M_0 = 3.3 for A <= 110
142  G4double MaxAverageMultiplicity = 2.6;
143  if (A <= 110) { MaxAverageMultiplicity = 3.3; }
144  return MaxAverageMultiplicity;
145 }
146 
static G4double GetGamma0()
static const G4double fr0
static const G4double fGamma0
static const G4double fCriticalTemp
static G4double GetKappaCoulomb()
static const G4double fEpsilon0
int G4int
Definition: G4Types.hh:78
static const G4double fE0
static G4double GetMaxAverageMultiplicity(G4int A)
static G4double Getr0()
double A(double temperature)
static G4double GetE0()
static G4double DBetaDT(G4double T)
static G4double GetCoulomb()
static const G4double fCoulomb
static const G4double fKappa
static constexpr double MeV
Definition: G4SIunits.hh:214
static G4double GetEpsilon0()
static G4double GetBeta0()
double G4double
Definition: G4Types.hh:76
static const G4double fKappaCoulomb
static constexpr double fermi
Definition: G4SIunits.hh:103
static G4double Beta(G4double T)
static G4double GetKappa()
static G4double GetCriticalTemp()
static const G4double fBeta0