#include <G4ParticleHPFissionSpectrum.hh>

Inheritance diagram for G4ParticleHPFissionSpectrum:

Public Member Functions
	G4ParticleHPFissionSpectrum ()

G4double	GetFractionalProbability (G4double anEnergy)

void	Init (std::istream &aDataFile)

G4double	Sample (G4double anEnergy)

	~G4ParticleHPFissionSpectrum ()

Private Member Functions
G4double	Maxwell (G4double anEnergy, G4double theta)

Private Attributes
G4double	expm1

G4ParticleHPVector	theFractionalProb

G4ParticleHPVector	theThetaDist

Detailed Description

Definition at line 44 of file G4ParticleHPFissionSpectrum.hh.

Constructor & Destructor Documentation

◆ G4ParticleHPFissionSpectrum()

G4ParticleHPFissionSpectrum::G4ParticleHPFissionSpectrum ( )

inline

Definition at line 47 of file G4ParticleHPFissionSpectrum.hh.

  {
    expm1 = G4Exp(-1.);
  }

References expm1, and G4Exp().

◆ ~G4ParticleHPFissionSpectrum()

G4ParticleHPFissionSpectrum::~G4ParticleHPFissionSpectrum ( )

inline

Definition at line 51 of file G4ParticleHPFissionSpectrum.hh.

52 {

53 }

Member Function Documentation

◆ GetFractionalProbability()

G4double G4ParticleHPFissionSpectrum::GetFractionalProbability ( G4double anEnergy )

inlinevirtual

Implements G4VParticleHPEDis.

Definition at line 61 of file G4ParticleHPFissionSpectrum.hh.

  {
    return theFractionalProb.GetY(anEnergy);
  }

References G4ParticleHPVector::GetY(), and theFractionalProb.

◆ Init()

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

inlinevirtual

Implements G4VParticleHPEDis.

Definition at line 55 of file G4ParticleHPFissionSpectrum.hh.

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

References CLHEP::eV, G4ParticleHPVector::Init(), theFractionalProb, and theThetaDist.

◆ Maxwell()

G4double G4ParticleHPFissionSpectrum::Maxwell	(	G4double	anEnergy,
		G4double	theta
	)

inlineprivate

Definition at line 95 of file G4ParticleHPFissionSpectrum.hh.

  {
    G4double result = std::sqrt(anEnergy/CLHEP::eV)*G4Exp(-anEnergy/CLHEP::eV/theta);
    return result;
  }

References CLHEP::eV, and G4Exp().

Referenced by Sample().

◆ Sample()

G4double G4ParticleHPFissionSpectrum::Sample ( G4double anEnergy )

inlinevirtual

Implements G4VParticleHPEDis.

Definition at line 66 of file G4ParticleHPFissionSpectrum.hh.

  {
    G4double theta = theThetaDist.GetY(anEnergy);
    // here we need to sample Maxwells distribution, if 
    // need be.
    G4double result=0., cut;
    G4double range =50*CLHEP::MeV;
    G4double max = Maxwell((theta*CLHEP::eV)/2., theta);
    G4double value;
    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;
      }
      result = range*G4UniformRand();
      value = Maxwell(result, theta);
      cut = G4UniformRand();
    }
    while(cut > value/max); // Loop checking, 11.05.2015, T. Koi
    return result;
  }