G4Nucleus Class Reference

#include <G4Nucleus.hh>

Public Member Functions
	G4Nucleus ()
	G4Nucleus (const G4double A, const G4double Z)
	G4Nucleus (const G4int A, const G4int Z)
	G4Nucleus (const G4Material *aMaterial)
	~G4Nucleus ()
	G4Nucleus (const G4Nucleus &right)
G4Nucleus &	operator= (const G4Nucleus &right)
G4bool	operator== (const G4Nucleus &right) const
G4bool	operator!= (const G4Nucleus &right) const
void	ChooseParameters (const G4Material *aMaterial)
void	SetParameters (const G4double A, const G4double Z)
void	SetParameters (const G4int A, const G4int Z)
G4int	GetA_asInt () const
G4int	GetN_asInt () const
G4int	GetZ_asInt () const
const G4Isotope *	GetIsotope ()
void	SetIsotope (const G4Isotope *iso)
G4DynamicParticle *	ReturnTargetParticle () const
G4double	AtomicMass (const G4double A, const G4double Z) const
G4double	AtomicMass (const G4int A, const G4int Z) const
G4double	GetThermalPz (const G4double mass, const G4double temp) const
G4ReactionProduct	GetThermalNucleus (G4double aMass, G4double temp=-1) const
G4ReactionProduct	GetBiasedThermalNucleus (G4double aMass, G4ThreeVector aVelocity, G4double temp=-1) const
G4double	Cinema (G4double kineticEnergy)
G4double	EvaporationEffects (G4double kineticEnergy)
G4double	AnnihilationEvaporationEffects (G4double kineticEnergy, G4double ekOrg)
G4double	GetPNBlackTrackEnergy () const
G4double	GetDTABlackTrackEnergy () const
G4double	GetAnnihilationPNBlackTrackEnergy () const
G4double	GetAnnihilationDTABlackTrackEnergy () const
G4ThreeVector	GetFermiMomentum ()
G4ReactionProductVector *	Fragmentate ()
void	AddExcitationEnergy (G4double anEnergy)
void	AddMomentum (const G4ThreeVector aMomentum)
G4double	GetEnergyDeposit ()

Detailed Description

Definition at line 50 of file G4Nucleus.hh.

Constructor & Destructor Documentation

G4Nucleus::G4Nucleus

(

)

Definition at line 54 of file G4Nucleus.cc.

  : theA(0), theZ(0), aEff(0.0), zEff(0)
{
  pnBlackTrackEnergy = 0.0;
  dtaBlackTrackEnergy = 0.0;
  pnBlackTrackEnergyfromAnnihilation = 0.0;
  dtaBlackTrackEnergyfromAnnihilation = 0.0;
  excitationEnergy = 0.0;
  momentum = G4ThreeVector(0.,0.,0.);
  fermiMomentum = 1.52*hbarc/fermi;
  theTemp = 293.16*kelvin;
  fIsotope = 0;
}

G4Nucleus::G4Nucleus	(	const G4double	A,
		const G4double	Z
	)

Definition at line 68 of file G4Nucleus.cc.

{
  SetParameters( A, Z );
  pnBlackTrackEnergy = 0.0;
  dtaBlackTrackEnergy = 0.0;
  pnBlackTrackEnergyfromAnnihilation = 0.0;
  dtaBlackTrackEnergyfromAnnihilation = 0.0;
  excitationEnergy = 0.0;
  momentum = G4ThreeVector(0.,0.,0.);
  fermiMomentum = 1.52*hbarc/fermi;
  theTemp = 293.16*kelvin;
  fIsotope = 0;
}

Here is the call graph for this function:

G4Nucleus::G4Nucleus	(	const G4int	A,
		const G4int	Z
	)

Definition at line 82 of file G4Nucleus.cc.

{
  SetParameters( A, Z );
  pnBlackTrackEnergy = 0.0;
  dtaBlackTrackEnergy = 0.0;
  pnBlackTrackEnergyfromAnnihilation = 0.0;
  dtaBlackTrackEnergyfromAnnihilation = 0.0;
  excitationEnergy = 0.0;
  momentum = G4ThreeVector(0.,0.,0.);
  fermiMomentum = 1.52*hbarc/fermi;
  theTemp = 293.16*kelvin;
  fIsotope = 0;
}

Here is the call graph for this function:

G4Nucleus::G4Nucleus

(

const G4Material *

aMaterial

)

Definition at line 96 of file G4Nucleus.cc.

{
  ChooseParameters( aMaterial );
  pnBlackTrackEnergy = 0.0;
  dtaBlackTrackEnergy = 0.0;
  pnBlackTrackEnergyfromAnnihilation = 0.0;
  dtaBlackTrackEnergyfromAnnihilation = 0.0;
  excitationEnergy = 0.0;
  momentum = G4ThreeVector(0.,0.,0.);
  fermiMomentum = 1.52*hbarc/fermi;
  theTemp = aMaterial->GetTemperature();
  fIsotope = 0;
}

Here is the call graph for this function:

G4Nucleus::~G4Nucleus

(

)

Definition at line 110 of file G4Nucleus.cc.

{}

G4Nucleus::G4Nucleus ( const G4Nucleus &  right )

inline

Definition at line 61 of file G4Nucleus.hh.

    { *this = right; }

Member Function Documentation

void G4Nucleus::AddExcitationEnergy

(

G4double

anEnergy

)

Definition at line 449 of file G4Nucleus.cc.

  {
    excitationEnergy+=anEnergy;
  }

void G4Nucleus::AddMomentum

(

const G4ThreeVector

aMomentum

)

Definition at line 444 of file G4Nucleus.cc.

  {
    momentum+=(aMomentum);
  }

G4double G4Nucleus::AnnihilationEvaporationEffects	(	G4double	kineticEnergy,
		G4double	ekOrg
	)

Definition at line 337 of file G4Nucleus.cc.

  {
    // Nuclear evaporation as a function of atomic number and kinetic 
    // energy (MeV) of primary particle.  Modified for annihilation effects. 
    //
    if( aEff < 1.5 || ekOrg < 0.)
    {
      pnBlackTrackEnergyfromAnnihilation = 0.0;
      dtaBlackTrackEnergyfromAnnihilation = 0.0;
      return 0.0;
    }
    G4double ek = kineticEnergy/GeV;
    G4float ekin = std::min( 4.0, std::max( 0.1, ek ) );
    const G4float atno = std::min( 120., aEff ); 
    const G4float gfa = 2.0*((aEff-1.0)/70.)*G4Exp(-(aEff-1.0)/70.);

    G4float cfa = std::max( 0.15, 0.35 + ((0.35-0.05)/2.3)*G4Log(ekin) );
    G4float exnu = 7.716 * cfa * G4Exp(-cfa)
      * ((atno-1.0)/120.)*G4Exp(-(atno-1.0)/120.);
    G4float fpdiv = std::max( 0.5, 1.0-0.25*ekin*ekin );

    pnBlackTrackEnergyfromAnnihilation = exnu*fpdiv;
    dtaBlackTrackEnergyfromAnnihilation = exnu*(1.0-fpdiv);
    
    G4double ran1 = -6.0;
    G4double ran2 = -6.0;
    for( G4int i=0; i<12; ++i ) {
      ran1 += G4UniformRand();
      ran2 += G4UniformRand();
    }
    pnBlackTrackEnergyfromAnnihilation *= 1.0 + ran1*gfa;
    dtaBlackTrackEnergyfromAnnihilation *= 1.0 + ran2*gfa;

    pnBlackTrackEnergyfromAnnihilation = std::max( 0.0, pnBlackTrackEnergyfromAnnihilation);
    dtaBlackTrackEnergyfromAnnihilation = std::max( 0.0, dtaBlackTrackEnergyfromAnnihilation);
    G4double blackSum = pnBlackTrackEnergyfromAnnihilation+dtaBlackTrackEnergyfromAnnihilation;
    if (blackSum >= ekOrg/GeV) {
      pnBlackTrackEnergyfromAnnihilation *= ekOrg/GeV/blackSum;
      dtaBlackTrackEnergyfromAnnihilation *= ekOrg/GeV/blackSum;
    }

    return (pnBlackTrackEnergyfromAnnihilation+dtaBlackTrackEnergyfromAnnihilation)*GeV;
  }

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4Nucleus::AtomicMass	(	const G4double	A,
		const G4double	Z
	)		const

Definition at line 254 of file G4Nucleus.cc.

  {
    // Now returns (atomic mass - electron masses) 
    return G4NucleiProperties::GetNuclearMass(A, Z);
  }

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4Nucleus::AtomicMass	(	const G4int	A,
		const G4int	Z
	)		const

Definition at line 261 of file G4Nucleus.cc.

  {
    // Now returns (atomic mass - electron masses) 
    return G4NucleiProperties::GetNuclearMass(A, Z);
  }

Here is the call graph for this function:

void G4Nucleus::ChooseParameters

(

const G4Material *

aMaterial

)

Definition at line 172 of file G4Nucleus.cc.

{
  G4double random = G4UniformRand();
  G4double sum = aMaterial->GetTotNbOfAtomsPerVolume();
  const G4ElementVector* theElementVector = aMaterial->GetElementVector();
  G4double running(0);
  //  G4Element* element(0);
  G4Element* element = (*theElementVector)[aMaterial->GetNumberOfElements()-1];

  for (unsigned int i = 0; i < aMaterial->GetNumberOfElements(); ++i) {
    running += aMaterial->GetVecNbOfAtomsPerVolume()[i];
    if (running > random*sum) {
      element = (*theElementVector)[i];
      break;
    }
  }

  if (element->GetNumberOfIsotopes() > 0) {
    G4double randomAbundance = G4UniformRand();
    G4double sumAbundance = element->GetRelativeAbundanceVector()[0];
    unsigned int iso=0;
    while (iso < element->GetNumberOfIsotopes() &&  /* Loop checking, 02.11.2015, A.Ribon */
           sumAbundance < randomAbundance) {
      ++iso;
      sumAbundance += element->GetRelativeAbundanceVector()[iso];
    }
    theA=element->GetIsotope(iso)->GetN();
    theZ=element->GetIsotope(iso)->GetZ();
    aEff=theA;
    zEff=theZ;
  } else {   
    aEff = element->GetN();
    zEff = element->GetZ();
    theZ = G4int(zEff + 0.5);
    theA = G4int(aEff + 0.5);   
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4Nucleus::Cinema

(

G4double

kineticEnergy

)

Definition at line 382 of file G4Nucleus.cc.

  {
    // derived from original FORTRAN code CINEMA by H. Fesefeldt (14-Oct-1987)
    //
    // input: kineticEnergy (MeV)
    // returns modified kinetic energy (MeV)
    //
    static const G4double expxu =  82.;           // upper bound for arg. of exp
    static const G4double expxl = -expxu;         // lower bound for arg. of exp
    
    G4double ek = kineticEnergy/GeV;
    G4double ekLog = G4Log( ek );
    G4double aLog = G4Log( aEff );
    G4double em = std::min( 1.0, 0.2390 + 0.0408*aLog*aLog );
    G4double temp1 = -ek * std::min( 0.15, 0.0019*aLog*aLog*aLog );
    G4double temp2 = G4Exp( std::max( expxl, std::min( expxu, -(ekLog-em)*(ekLog-em)*2.0 ) ) );
    G4double result = 0.0;
    if( std::abs( temp1 ) < 1.0 )
    {
      if( temp2 > 1.0e-10 )result = temp1*temp2;
    }
    else result = temp1*temp2;
    if( result < -ek )result = -ek;
    return result*GeV;
  }

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4Nucleus::EvaporationEffects

(

G4double

kineticEnergy

)

Definition at line 278 of file G4Nucleus.cc.

  {
    // derived from original FORTRAN code EXNU by H. Fesefeldt (10-Dec-1986)
    //
    // Nuclear evaporation as function of atomic number
    // and kinetic energy (MeV) of primary particle
    //
    // returns kinetic energy (MeV)
    //
    if( aEff < 1.5 )
    {
      pnBlackTrackEnergy = dtaBlackTrackEnergy = 0.0;
      return 0.0;
    }
    G4double ek = kineticEnergy/GeV;
    G4float ekin = std::min( 4.0, std::max( 0.1, ek ) );
    const G4float atno = std::min( 120., aEff ); 
    const G4float gfa = 2.0*((aEff-1.0)/70.)*G4Exp(-(aEff-1.0)/70.);
    //
    // 0.35 value at 1 GeV
    // 0.05 value at 0.1 GeV
    //
    G4float cfa = std::max( 0.15, 0.35 + ((0.35-0.05)/2.3)*G4Log(ekin) );
    G4float exnu = 7.716 * cfa * G4Exp(-cfa)
      * ((atno-1.0)/120.)*G4Exp(-(atno-1.0)/120.);
    G4float fpdiv = std::max( 0.5, 1.0-0.25*ekin*ekin );
    //
    // pnBlackTrackEnergy  is the kinetic energy (in GeV) available for
    //                     proton/neutron black track particles
    // dtaBlackTrackEnergy is the kinetic energy (in GeV) available for
    //                     deuteron/triton/alpha black track particles
    //
    pnBlackTrackEnergy = exnu*fpdiv;
    dtaBlackTrackEnergy = exnu*(1.0-fpdiv);
    
    if( G4int(zEff+0.1) != 82 )
    { 
      G4double ran1 = -6.0;
      G4double ran2 = -6.0;
      for( G4int i=0; i<12; ++i )
      {
        ran1 += G4UniformRand();
        ran2 += G4UniformRand();
      }
      pnBlackTrackEnergy *= 1.0 + ran1*gfa;
      dtaBlackTrackEnergy *= 1.0 + ran2*gfa;
    }
    pnBlackTrackEnergy = std::max( 0.0, pnBlackTrackEnergy );
    dtaBlackTrackEnergy = std::max( 0.0, dtaBlackTrackEnergy );
    while( pnBlackTrackEnergy+dtaBlackTrackEnergy >= ek )  /* Loop checking, 02.11.2015, A.Ribon */
    {
      pnBlackTrackEnergy *= 1.0 - 0.5*G4UniformRand();
      dtaBlackTrackEnergy *= 1.0 - 0.5*G4UniformRand();
    }
//    G4cout << "EvaporationEffects "<<kineticEnergy<<" "
//           <<pnBlackTrackEnergy+dtaBlackTrackEnergy<<endl;
    return (pnBlackTrackEnergy+dtaBlackTrackEnergy)*GeV;
  }

Here is the call graph for this function:

Here is the caller graph for this function:

G4ReactionProductVector * G4Nucleus::Fragmentate

(

)

Definition at line 438 of file G4Nucleus.cc.

  {
    // needs implementation!
    return NULL;
  }

G4int G4Nucleus::GetA_asInt ( void  ) const

inline

Definition at line 109 of file G4Nucleus.hh.

    { return theA; }   

Here is the caller graph for this function:

G4double G4Nucleus::GetAnnihilationDTABlackTrackEnergy ( ) const

inline

Definition at line 159 of file G4Nucleus.hh.

    { return dtaBlackTrackEnergyfromAnnihilation; }

Here is the caller graph for this function:

G4double G4Nucleus::GetAnnihilationPNBlackTrackEnergy ( ) const

inline

Definition at line 156 of file G4Nucleus.hh.

    { return pnBlackTrackEnergyfromAnnihilation; }

Here is the caller graph for this function:

G4ReactionProduct G4Nucleus::GetBiasedThermalNucleus	(	G4double	aMass,
		G4ThreeVector	aVelocity,
		G4double	temp = -1
	)		const

Definition at line 113 of file G4Nucleus.cc.

{
  G4double velMag = aVelocity.mag();
  G4ReactionProduct result;
  G4double value = 0;
  G4double random = 1;
  G4double norm = 3.*std::sqrt(k_Boltzmann*temp*aMass*G4Neutron::Neutron()->GetPDGMass());
  norm /= G4Neutron::Neutron()->GetPDGMass();
  norm *= 5.;
  norm += velMag;
  norm /= velMag;
  const G4int maxNumberOfLoops = 1000000;
  G4int loopCounter = -1;
  while ( (value/norm<random) && ++loopCounter < maxNumberOfLoops )  /* Loop checking, 02.11.2015, A.Ribon */
  {
     result = GetThermalNucleus(aMass, temp);
     G4ThreeVector targetVelocity = 1./result.GetMass()*result.GetMomentum();
     value = (targetVelocity+aVelocity).mag()/velMag;
     random = G4UniformRand();
  }
  if ( loopCounter >= maxNumberOfLoops ) {
    G4ExceptionDescription ed;
    ed << " Failed sampling after maxNumberOfLoops attempts : forced exit! " << G4endl;
    G4Exception( " G4Nucleus::GetBiasedThermalNucleus ", "HAD_NUCLEUS_001", JustWarning, ed );
    result = GetThermalNucleus(aMass, temp);    
  }
  return result;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4Nucleus::GetDTABlackTrackEnergy ( ) const

inline

Definition at line 153 of file G4Nucleus.hh.

    { return dtaBlackTrackEnergy; }

Here is the caller graph for this function:

G4double G4Nucleus::GetEnergyDeposit ( )

inline

Definition at line 184 of file G4Nucleus.hh.

{return excitationEnergy; }

Here is the caller graph for this function:

G4ThreeVector G4Nucleus::GetFermiMomentum

(

)

Definition at line 412 of file G4Nucleus.cc.

  {
    // chv: .. we assume zero temperature!
    
    // momentum is equally distributed in each phasespace volume dpx, dpy, dpz.
    G4double ranflat1=
      G4RandFlat::shoot((G4double)0.,(G4double)fermiMomentum);   
    G4double ranflat2=
      G4RandFlat::shoot((G4double)0.,(G4double)fermiMomentum);   
    G4double ranflat3=
      G4RandFlat::shoot((G4double)0.,(G4double)fermiMomentum);   
    G4double ranmax = (ranflat1>ranflat2? ranflat1: ranflat2);
    ranmax = (ranmax>ranflat3? ranmax : ranflat3);
    
    // Isotropic momentum distribution
    G4double costheta = 2.*G4UniformRand() - 1.0;
    G4double sintheta = std::sqrt(1.0 - costheta*costheta);
    G4double phi = 2.0*pi*G4UniformRand();
    
    G4double pz=costheta*ranmax;
    G4double px=sintheta*std::cos(phi)*ranmax;
    G4double py=sintheta*std::sin(phi)*ranmax;
    G4ThreeVector p(px,py,pz);
    return p;
  }

Here is the call graph for this function:

const G4Isotope* G4Nucleus::GetIsotope ( )

inline

Definition at line 119 of file G4Nucleus.hh.

    { return fIsotope; }

Here is the caller graph for this function:

G4int G4Nucleus::GetN_asInt ( ) const

inline

Definition at line 112 of file G4Nucleus.hh.

    { return theA-theZ; }   

Here is the caller graph for this function:

G4double G4Nucleus::GetPNBlackTrackEnergy ( ) const

inline

Definition at line 150 of file G4Nucleus.hh.

    { return pnBlackTrackEnergy; }

Here is the caller graph for this function:

G4ReactionProduct G4Nucleus::GetThermalNucleus	(	G4double	aMass,
		G4double	temp = -1
	)		const

Definition at line 143 of file G4Nucleus.cc.

{
  G4double currentTemp = temp;
  if (currentTemp < 0) currentTemp = theTemp;
  G4ReactionProduct theTarget;    
  theTarget.SetMass(targetMass*G4Neutron::Neutron()->GetPDGMass());
  G4double px, py, pz;
  px = GetThermalPz(theTarget.GetMass(), currentTemp);
  py = GetThermalPz(theTarget.GetMass(), currentTemp);
  pz = GetThermalPz(theTarget.GetMass(), currentTemp);
  theTarget.SetMomentum(px, py, pz);
  G4double tMom = std::sqrt(px*px+py*py+pz*pz);
  G4double tEtot = std::sqrt((tMom+theTarget.GetMass())*
                             (tMom+theTarget.GetMass())-
                              2.*tMom*theTarget.GetMass());
  //  if(1-tEtot/theTarget.GetMass()>0.001)  this line incorrect (Bug report 1911) 
  if (tEtot/theTarget.GetMass() - 1. > 0.001) {
    // use relativistic energy for higher energies
    theTarget.SetTotalEnergy(tEtot);

  } else {
    // use p**2/2M for lower energies (to preserve precision?) 
    theTarget.SetKineticEnergy(tMom*tMom/(2.*theTarget.GetMass()));
  }    
  return theTarget;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4Nucleus::GetThermalPz	(	const G4double	mass,
		const G4double	temp
	)		const

Definition at line 268 of file G4Nucleus.cc.

  {
    G4double result = G4RandGauss::shoot();
    result *= std::sqrt(k_Boltzmann*temp*mass); // Das ist impuls (Pz),
                                           // nichtrelativistische rechnung
                                           // Maxwell verteilung angenommen
    return result;
  }

Here is the call graph for this function:

Here is the caller graph for this function:

G4int G4Nucleus::GetZ_asInt ( void  ) const

inline

Definition at line 115 of file G4Nucleus.hh.

    { return theZ; }   

G4bool G4Nucleus::operator!= ( const G4Nucleus &  right ) const

inline

Definition at line 89 of file G4Nucleus.hh.

    { return ( this != (G4Nucleus *) &right ); }

G4Nucleus& G4Nucleus::operator= ( const G4Nucleus &  right )

inline

Definition at line 64 of file G4Nucleus.hh.

    {
      if (this != &right) {
        theA=right.theA;
        theZ=right.theZ;
        aEff=right.aEff;
        zEff=right.zEff;
        fIsotope = right.fIsotope;
        pnBlackTrackEnergy=right.pnBlackTrackEnergy; 
        dtaBlackTrackEnergy=right.dtaBlackTrackEnergy;
        pnBlackTrackEnergyfromAnnihilation =
                     right.pnBlackTrackEnergyfromAnnihilation; 
        dtaBlackTrackEnergyfromAnnihilation =
                     right.dtaBlackTrackEnergyfromAnnihilation; 
        theTemp = right.theTemp;
        excitationEnergy = right.excitationEnergy;
        momentum = right.momentum;
        fermiMomentum = right.fermiMomentum;
      }
      return *this;
    }

G4bool G4Nucleus::operator== ( const G4Nucleus &  right ) const

inline

Definition at line 86 of file G4Nucleus.hh.

    { return ( this == (G4Nucleus *) &right ); }

G4DynamicParticle * G4Nucleus::ReturnTargetParticle

(

)

const

Definition at line 241 of file G4Nucleus.cc.

  {
    // choose a proton or a neutron as the target particle
    
    G4DynamicParticle *targetParticle = new G4DynamicParticle;
    if( G4UniformRand() < zEff/aEff )
      targetParticle->SetDefinition( G4Proton::Proton() );
    else
      targetParticle->SetDefinition( G4Neutron::Neutron() );
    return targetParticle;
  }

Here is the call graph for this function:

Here is the caller graph for this function:

void G4Nucleus::SetIsotope ( const G4Isotope *  iso )

inline

Definition at line 122 of file G4Nucleus.hh.

    { 
      fIsotope = iso;
      if(iso) { 
    theZ = iso->GetZ();
        theA = iso->GetN();
        aEff = theA;
        zEff = theZ;
      }
    }

Here is the call graph for this function:

Here is the caller graph for this function:

void G4Nucleus::SetParameters	(	const G4double	A,
		const G4double	Z
	)

Definition at line 212 of file G4Nucleus.cc.

{
  theZ = G4lrint(Z);
  theA = G4lrint(A);   
  if (theA<1 || theZ<0 || theZ>theA) {
    throw G4HadronicException(__FILE__, __LINE__,
            "G4Nucleus::SetParameters called with non-physical parameters");
  }
  aEff = A;  // atomic weight
  zEff = Z;  // atomic number
  fIsotope = 0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void G4Nucleus::SetParameters	(	const G4int	A,
		const G4int	Z
	)

Definition at line 226 of file G4Nucleus.cc.

{
  theZ = Z;
  theA = A;   
  if( theA<1 || theZ<0 || theZ>theA )
    {
      throw G4HadronicException(__FILE__, __LINE__,
                "G4Nucleus::SetParameters called with non-physical parameters");
    }
  aEff = A;  // atomic weight
  zEff = Z;  // atomic number
  fIsotope = 0;
}

Here is the call graph for this function:

---

The documentation for this class was generated from the following files:

• geant4.10.03.p01/source/processes/hadronic/util/include/G4Nucleus.hh
• geant4.10.03.p01/source/processes/hadronic/util/src/G4Nucleus.cc