G4ParticleHPMadlandNixSpectrum Class Reference

#include <G4ParticleHPMadlandNixSpectrum.hh>

Inheritance diagram for G4ParticleHPMadlandNixSpectrum:

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

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 ( )

inline

Definition at line 60 of file G4ParticleHPMadlandNixSpectrum.hh.

  {
  }

Member Function Documentation

GetFractionalProbability()

G4double G4ParticleHPMadlandNixSpectrum::GetFractionalProbability ( G4double  anEnergy )

inlinevirtual

Implements G4VParticleHPEDis.

Definition at line 74 of file G4ParticleHPMadlandNixSpectrum.hh.

  {
    return theFractionalProb.GetY(anEnergy);
  }

Init()

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

inlinevirtual

Implements G4VParticleHPEDis.

Definition at line 64 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 )

virtual

Implements G4VParticleHPEDis.

Definition at line 74 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