Class that stores isotopic abundances for a given element. More...

#include <G4INCLNaturalIsotopicDistributions.hh>

Public Member Functions
	IsotopicDistribution (IsotopeVector const &aVector)
	Constructor.

G4int	drawRandomIsotope () const
	Draw a random isotope based on the abundance vector.

IsotopeVector const &	getIsotopes () const
	Get the isotope vector.

Detailed Description

Class that stores isotopic abundances for a given element.

Definition at line 64 of file G4INCLNaturalIsotopicDistributions.hh.

Constructor & Destructor Documentation

◆ IsotopicDistribution()

G4INCL::IsotopicDistribution::IsotopicDistribution ( IsotopeVector const & aVector )

Constructor.

Definition at line 59 of file G4INCLNaturalIsotopicDistributions.cc.

                                                                         :
    theIsotopes(aVector)
  {
    G4double previousAbundance = 0.;
    // Cumulate the abundances
    for(IsotopeIter i=theIsotopes.begin(), e=theIsotopes.end(); i!=e; ++i) {
      i->theAbundance += previousAbundance;
      previousAbundance = i->theAbundance;
    }
    // Normalize the abundances to 1
    const G4double normalisation = 1./theIsotopes.back().theAbundance;
    for(IsotopeIter i=theIsotopes.begin(), e=theIsotopes.end(); i!=e; ++i)
      i->theAbundance *= normalisation;
  }

Member Function Documentation

◆ drawRandomIsotope()

G4int G4INCL::IsotopicDistribution::drawRandomIsotope ( ) const

Draw a random isotope based on the abundance vector.

Definition at line 74 of file G4INCLNaturalIsotopicDistributions.cc.

                                                      {
    const G4double r = Random::shoot();
    for(unsigned int i=0; i<theIsotopes.size()-1; ++i) {
      if(r<=theIsotopes.at(i).theAbundance)
        return theIsotopes.at(i).theA;
    }
    return theIsotopes.back().theA;
  }

Referenced by G4INCL::NaturalIsotopicDistributions::drawRandomIsotope(), and G4INCL::ParticleTable::drawRandomNaturalIsotope().

◆ getIsotopes()

IsotopeVector const & G4INCL::IsotopicDistribution::getIsotopes ( ) const

Get the isotope vector.

Definition at line 83 of file G4INCLNaturalIsotopicDistributions.cc.

                                                               {
    return theIsotopes;
  }

Referenced by G4INCL::INCL::initUniverseRadiusForAntiprotonAtRest().