G4ParticleHPWattSpectrum Class Reference

#include <G4ParticleHPWattSpectrum.hh>

Inheritance diagram for G4ParticleHPWattSpectrum:

Public Member Functions
	G4ParticleHPWattSpectrum ()
	~G4ParticleHPWattSpectrum () override=default
void	Init (std::istream &aDataFile) override
G4double	GetFractionalProbability (G4double anEnergy) override
G4double	Sample (G4double anEnergy) override
Public Member Functions inherited from G4VParticleHPEDis
	G4VParticleHPEDis ()=default
virtual	~G4VParticleHPEDis ()=default

Detailed Description

Definition at line 45 of file G4ParticleHPWattSpectrum.hh.

Constructor & Destructor Documentation

G4ParticleHPWattSpectrum()

G4ParticleHPWattSpectrum::G4ParticleHPWattSpectrum ( )

inline

Definition at line 48 of file G4ParticleHPWattSpectrum.hh.

{ expm1 = G4Exp(-1.); }

~G4ParticleHPWattSpectrum()

G4ParticleHPWattSpectrum::~G4ParticleHPWattSpectrum ( )

overridedefault

Member Function Documentation

GetFractionalProbability()

G4double G4ParticleHPWattSpectrum::GetFractionalProbability ( G4double anEnergy )

inlineoverridevirtual

Implements G4VParticleHPEDis.

Definition at line 58 of file G4ParticleHPWattSpectrum.hh.

    {
      return theFractionalProb.GetY(anEnergy);
    }

Init()

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

inlineoverridevirtual

Implements G4VParticleHPEDis.

Definition at line 51 of file G4ParticleHPWattSpectrum.hh.

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

Sample()

G4double G4ParticleHPWattSpectrum::Sample ( G4double anEnergy )

overridevirtual

Implements G4VParticleHPEDis.

Definition at line 36 of file G4ParticleHPWattSpectrum.cc.

{
  G4double a = theApar.GetY(anEnergy) * eV;
  G4double b = theBpar.GetY(anEnergy) / eV;
  G4double result = 0.;
  G4double random, cut, max;
  max = std::sinh(std::sqrt(b * 15. * a));
 
  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;
    }
    random = G4UniformRand();
    result = -a * G4Log(random);
    cut = G4UniformRand();
  } while (cut > std::sinh(std::sqrt(b * result)) / max);  // Loop checking, 11.05.2015, T. Koi
  return result;
}

---

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

• G4ParticleHPWattSpectrum.hh
• G4ParticleHPWattSpectrum.cc