G4ParticleHPMadlandNixSpectrum Class Reference

#include <G4ParticleHPMadlandNixSpectrum.hh>

Inheritance diagram for G4ParticleHPMadlandNixSpectrum:

Public Member Functions
	G4ParticleHPMadlandNixSpectrum ()
	~G4ParticleHPMadlandNixSpectrum () override=default
void	Init (std::istream &aDataFile) override
G4double	GetFractionalProbability (G4double anEnergy) override
G4double	Sample (G4double anEnergy) override
Public Member Functions inherited from G4VParticleHPEDis
	G4VParticleHPEDis ()=default
virtual	~G4VParticleHPEDis ()=default

Detailed Description

Definition at line 51 of file G4ParticleHPMadlandNixSpectrum.hh.

Constructor & Destructor Documentation

G4ParticleHPMadlandNixSpectrum()

G4ParticleHPMadlandNixSpectrum::G4ParticleHPMadlandNixSpectrum ( )

inline

Definition at line 54 of file G4ParticleHPMadlandNixSpectrum.hh.

    {
      expm1 = G4Exp(-1.);
      theAvarageKineticPerNucleonForLightFragments = 0.0;
      theAvarageKineticPerNucleonForHeavyFragments = 0.0;
    }

~G4ParticleHPMadlandNixSpectrum()

G4ParticleHPMadlandNixSpectrum::~G4ParticleHPMadlandNixSpectrum ( )

overridedefault

Member Function Documentation

GetFractionalProbability()

G4double G4ParticleHPMadlandNixSpectrum::GetFractionalProbability ( G4double anEnergy )

inlineoverridevirtual

Implements G4VParticleHPEDis.

Definition at line 72 of file G4ParticleHPMadlandNixSpectrum.hh.

    {
      return theFractionalProb.GetY(anEnergy);
    }

Init()

void G4ParticleHPMadlandNixSpectrum::Init ( std::istream & aDataFile )

inlineoverridevirtual

Implements G4VParticleHPEDis.

Definition at line 62 of file G4ParticleHPMadlandNixSpectrum.hh.

    {
      theFractionalProb.Init(aDataFile);
      aDataFile >> theAvarageKineticPerNucleonForLightFragments;
      theAvarageKineticPerNucleonForLightFragments *= CLHEP::eV;
      aDataFile >> theAvarageKineticPerNucleonForHeavyFragments;
      theAvarageKineticPerNucleonForHeavyFragments *= CLHEP::eV;
      theMaxTemp.Init(aDataFile);
    }

Sample()

G4double G4ParticleHPMadlandNixSpectrum::Sample ( G4double anEnergy )

overridevirtual

Implements G4VParticleHPEDis.

Definition at line 73 of file G4ParticleHPMadlandNixSpectrum.cc.

{
  G4double tm = theMaxTemp.GetY(anEnergy);
  G4double last = 0, buff, current = 100 * MeV;
  G4double precision = 0.001;
  G4double newValue = 0., oldValue = 0.;
  G4double random = G4UniformRand();
 
  G4int icounter = 0;
  G4int icounter_max = 1024;
  do {
    icounter++;
    if (icounter > icounter_max) {
      G4cout << "Loop-counter exceeded the threshold value at " << __LINE__ << "th line of "
             << __FILE__ << "." << G4endl;
      break;
    }
    oldValue = newValue;
    newValue = FissionIntegral(tm, current);
    if (newValue < random) {
      buff = current;
      current += std::abs(current - last) / 2.;
      last = buff;
      if (current > 190 * MeV)
        throw G4HadronicException(__FILE__, __LINE__,
                                  "Madland-Nix Spectrum has not converged in sampling");
    }
    else {
      buff = current;
      current -= std::abs(current - last) / 2.;
      last = buff;
    }
  } while (std::abs(oldValue - newValue)
           > precision * newValue);  // Loop checking, 11.05.2015, T. Koi
  return current;
}

---

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

• G4ParticleHPMadlandNixSpectrum.hh
• G4ParticleHPMadlandNixSpectrum.cc