Ignore:
Timestamp:
Dec 22, 2010, 3:52:27 PM (13 years ago)
Author:
garnier
Message:

geant4 tag 9.4

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/source/processes/hadronic/models/util/src/G4NuclearFermiDensity.cc

    r819 r1347  
    2929#include "G4NuclearFermiDensity.hh"
    3030
    31 G4NuclearFermiDensity::G4NuclearFermiDensity(G4double anA, G4double aZ)
    32   :  a(0.545 * fermi)
     31G4NuclearFermiDensity::G4NuclearFermiDensity(G4int anA, G4int aZ)
     32  :  theA(anA), theZ(aZ), a(0.545 * fermi)
    3333{
    3434//        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;
    3836        theR= r0 * std::pow(anA, 1./3. );
    3937        Setrho0(3./ (4. * pi * std::pow(r0,3.) * theA * ( 1. + sqr(a/theR)*pi2 )));
Note: See TracChangeset for help on using the changeset viewer.