- Timestamp:
- Nov 5, 2010, 3:45:55 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/models/cascade/cascade/src/nucleiLevelDensity.cc
r819 r1340 23 23 // * acceptance of all terms of the Geant4 Software license. * 24 24 // ******************************************************************** 25 // $Id: nucleiLevelDensity.cc,v 1.8 2010/09/14 17:51:36 mkelsey Exp $ 25 26 // 27 // 20100914 M. Kelsey -- Migrate to integer A and Z; simplify code 28 26 29 #include "G4InuclSpecialFunctions.hh" 27 30 28 G4double G4InuclSpecialFunctions::nucleiLevelDensity(G4double a) {29 G4int verboseLevel = 2;30 31 31 if (verboseLevel > 3) { 32 G4cout << " >>> G4InuclSpecialFunctions::nucleiLevelDensity" << G4endl; 33 } 34 32 G4double G4InuclSpecialFunctions::nucleiLevelDensity(G4int a) { 35 33 const G4double NLD[226] = { 36 34 // 20 - 29 … … 81 79 29.87, 30.25, 30.5, 29.8, 29.17, 28.67}; 82 80 83 G4int ia = G4int(a + 0.1) - 20; 84 85 if (ia > 0) { 86 return NLD[ia]; 87 88 } else { 89 return 0.1 * a; 90 }; 81 return (a>=20) ? NLD[a-20] : 0.1*a; 91 82 }
Note: See TracChangeset
for help on using the changeset viewer.