Ignore:
Timestamp:
Nov 5, 2010, 3:45:55 PM (14 years ago)
Author:
garnier
Message:

update ti head

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/source/processes/hadronic/models/cascade/cascade/src/nucleiLevelDensity.cc

    r819 r1340  
    2323// * acceptance of all terms of the Geant4 Software license.          *
    2424// ********************************************************************
     25// $Id: nucleiLevelDensity.cc,v 1.8 2010/09/14 17:51:36 mkelsey Exp $
    2526//
     27// 20100914  M. Kelsey -- Migrate to integer A and Z; simplify code
     28
    2629#include "G4InuclSpecialFunctions.hh"
    2730
    28 G4double G4InuclSpecialFunctions::nucleiLevelDensity(G4double a) {
    29   G4int verboseLevel = 2;
    3031
    31   if (verboseLevel > 3) {
    32     G4cout << " >>> G4InuclSpecialFunctions::nucleiLevelDensity" << G4endl;
    33   }
    34 
     32G4double G4InuclSpecialFunctions::nucleiLevelDensity(G4int a) {
    3533  const G4double NLD[226] = {
    3634    // 20 - 29
     
    8179    29.87, 30.25, 30.5, 29.8, 29.17, 28.67};
    8280
    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;
    9182}
Note: See TracChangeset for help on using the changeset viewer.