G4NeutronHPMadlandNixSpectrum Class Reference

#include <G4NeutronHPMadlandNixSpectrum.hh>

Inheritance diagram for G4NeutronHPMadlandNixSpectrum:

Public Member Functions
	G4NeutronHPMadlandNixSpectrum ()
	~G4NeutronHPMadlandNixSpectrum ()
void	Init (std::istream &aDataFile)
G4double	GetFractionalProbability (G4double anEnergy)
G4double	Sample (G4double anEnergy)
Public Member Functions inherited from G4VNeutronHPEDis
	G4VNeutronHPEDis ()
virtual	~G4VNeutronHPEDis ()

Detailed Description

Definition at line 47 of file G4NeutronHPMadlandNixSpectrum.hh.

Constructor & Destructor Documentation

G4NeutronHPMadlandNixSpectrum::G4NeutronHPMadlandNixSpectrum ( )

inline

Definition at line 50 of file G4NeutronHPMadlandNixSpectrum.hh.

  {
    expm1 = std::exp(-1.);
  }

G4NeutronHPMadlandNixSpectrum::~G4NeutronHPMadlandNixSpectrum ( )

inline

Definition at line 54 of file G4NeutronHPMadlandNixSpectrum.hh.

  {
  }

Member Function Documentation

G4double G4NeutronHPMadlandNixSpectrum::GetFractionalProbability ( G4double  anEnergy )

inlinevirtual

Implements G4VNeutronHPEDis.

Definition at line 68 of file G4NeutronHPMadlandNixSpectrum.hh.

References G4NeutronHPVector::GetY().

  {
    return theFractionalProb.GetY(anEnergy);
  }

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

inlinevirtual

Implements G4VNeutronHPEDis.

Definition at line 58 of file G4NeutronHPMadlandNixSpectrum.hh.

References G4NeutronHPVector::Init().

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

G4double G4NeutronHPMadlandNixSpectrum::Sample ( G4double  anEnergy )

virtual

Implements G4VNeutronHPEDis.

Definition at line 72 of file G4NeutronHPMadlandNixSpectrum.cc.

References G4UniformRand, G4NeutronHPVector::GetY(), and python.hepunit::MeV.

  {
    G4double tm = theMaxTemp.GetY(anEnergy);
    G4double last=0, buff, current = 100*MeV;
    G4double precision = 0.001;
    G4double newValue = 0., oldValue=0.;
    G4double random = G4UniformRand();
    
    do
    {
      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);
    return current;
  }

---

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

• G4NeutronHPMadlandNixSpectrum.hh
• G4NeutronHPMadlandNixSpectrum.cc