- Timestamp:
- Dec 22, 2010, 3:52:27 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/models/util/src/G4NuclearFermiDensity.cc
r819 r1347 29 29 #include "G4NuclearFermiDensity.hh" 30 30 31 G4NuclearFermiDensity::G4NuclearFermiDensity(G4 double anA, G4doubleaZ)32 : a(0.545 * fermi)31 G4NuclearFermiDensity::G4NuclearFermiDensity(G4int anA, G4int aZ) 32 : theA(anA), theZ(aZ), a(0.545 * fermi) 33 33 { 34 34 // const G4double r0=1.14*fermi; 35 const G4double r0=1.16 * ( 1. - 1.16 * std::pow(anA, -2./3.)) * fermi; 36 theA = G4int(anA); 37 theZ = G4int(aZ); 35 const G4double r0=1.16 * ( 1. - 1.16 * std::pow(G4double(anA), -2./3.)) * fermi; 38 36 theR= r0 * std::pow(anA, 1./3. ); 39 37 Setrho0(3./ (4. * pi * std::pow(r0,3.) * theA * ( 1. + sqr(a/theR)*pi2 )));
Note: See TracChangeset
for help on using the changeset viewer.