G4NeutronHPWattSpectrum Class Reference

#include <G4NeutronHPWattSpectrum.hh>

Inheritance diagram for G4NeutronHPWattSpectrum:

Public Member Functions
	G4NeutronHPWattSpectrum ()
	~G4NeutronHPWattSpectrum ()
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 42 of file G4NeutronHPWattSpectrum.hh.

Constructor & Destructor Documentation

G4NeutronHPWattSpectrum::G4NeutronHPWattSpectrum ( )

inline

Definition at line 45 of file G4NeutronHPWattSpectrum.hh.

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

G4NeutronHPWattSpectrum::~G4NeutronHPWattSpectrum ( )

inline

Definition at line 49 of file G4NeutronHPWattSpectrum.hh.

  {
  }

Member Function Documentation

G4double G4NeutronHPWattSpectrum::GetFractionalProbability ( G4double  anEnergy )

inlinevirtual

Implements G4VNeutronHPEDis.

Definition at line 60 of file G4NeutronHPWattSpectrum.hh.

References G4NeutronHPVector::GetY().

  {
    return theFractionalProb.GetY(anEnergy);
  }

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

inlinevirtual

Implements G4VNeutronHPEDis.

Definition at line 53 of file G4NeutronHPWattSpectrum.hh.

References G4NeutronHPVector::Init().

  {
    theFractionalProb.Init(aDataFile, CLHEP::eV);
    theApar.Init(aDataFile, CLHEP::eV);
    theBpar.Init(aDataFile, CLHEP::eV);
  }

G4double G4NeutronHPWattSpectrum::Sample ( G4double  anEnergy )

virtual

Implements G4VNeutronHPEDis.

Definition at line 33 of file G4NeutronHPWattSpectrum.cc.

References test::a, test::b, python.hepunit::eV, G4UniformRand, G4NeutronHPVector::GetY(), and G4INCL::Math::max().

  {
    G4double a = theApar.GetY(anEnergy)*eV;
    G4double b = theBpar.GetY(anEnergy)/eV;
    G4double result;
    G4double random, cut, max;
    max = std::sinh(std::sqrt(b*15.*a));
    do
    {
      random = G4UniformRand();
      result = -a*std::log(random);
      cut = G4UniformRand();
    }
    while(cut>std::sinh(std::sqrt(b*result))/max);
    return result;
  }

---

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

• G4NeutronHPWattSpectrum.hh
• G4NeutronHPWattSpectrum.cc