G4NeutronHPFissionSpectrum Class Reference

#include <G4NeutronHPFissionSpectrum.hh>

Inheritance diagram for G4NeutronHPFissionSpectrum:

Public Member Functions
	G4NeutronHPFissionSpectrum ()
	~G4NeutronHPFissionSpectrum ()
void	Init (std::ifstream &aDataFile)
G4double	GetFractionalProbability (G4double anEnergy)
G4double	Sample (G4double anEnergy)
Public Member Functions inherited from G4VNeutronHPEDis
	G4VNeutronHPEDis ()
virtual	~G4VNeutronHPEDis ()
virtual void	Init (std::ifstream &theData)=0
virtual G4double	GetFractionalProbability (G4double anEnergy)=0
virtual G4double	Sample (G4double anEnergy)=0

Detailed Description

Definition at line 43 of file G4NeutronHPFissionSpectrum.hh.

Constructor & Destructor Documentation

G4NeutronHPFissionSpectrum()

G4NeutronHPFissionSpectrum::G4NeutronHPFissionSpectrum ( )

inline

Definition at line 46 of file G4NeutronHPFissionSpectrum.hh.

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

~G4NeutronHPFissionSpectrum()

G4NeutronHPFissionSpectrum::~G4NeutronHPFissionSpectrum ( )

inline

Definition at line 50 of file G4NeutronHPFissionSpectrum.hh.

  {
  }

Member Function Documentation

GetFractionalProbability()

G4double G4NeutronHPFissionSpectrum::GetFractionalProbability ( G4double  anEnergy )

inlinevirtual

Implements G4VNeutronHPEDis.

Definition at line 60 of file G4NeutronHPFissionSpectrum.hh.

  {
    return theFractionalProb.GetY(anEnergy);
  }

Init()

void G4NeutronHPFissionSpectrum::Init ( std::ifstream &  aDataFile )

inlinevirtual

Implements G4VNeutronHPEDis.

Definition at line 54 of file G4NeutronHPFissionSpectrum.hh.

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

Sample()

G4double G4NeutronHPFissionSpectrum::Sample ( G4double  anEnergy )

inlinevirtual

Implements G4VNeutronHPEDis.

Definition at line 65 of file G4NeutronHPFissionSpectrum.hh.

  {
    G4double theta = theThetaDist.GetY(anEnergy);
    // here we need to sample Maxwells distribution, if 
    // need be.
    G4double result, cut;
    G4double range =50*CLHEP::MeV;
    G4double max = Maxwell((theta*CLHEP::eV)/2., theta);
    G4double value;
    do
    {
      result = range*G4UniformRand();
      value = Maxwell(result, theta);
      cut = G4UniformRand();
    }
    while(cut > value/max);
    return result;
  }

---

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

• G4NeutronHPFissionSpectrum.hh