G4ParticleHPWattSpectrum Class Reference

#include <G4ParticleHPWattSpectrum.hh>

Inheritance diagram for G4ParticleHPWattSpectrum:

Collaboration diagram for G4ParticleHPWattSpectrum:

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

Private Member Functions
G4double	Watt (G4double anEnergy, G4double a, G4double b)

Private Attributes
G4double	expm1
G4ParticleHPVector	theFractionalProb
G4ParticleHPVector	theApar
G4ParticleHPVector	theBpar

Detailed Description

Definition at line 44 of file G4ParticleHPWattSpectrum.hh.

Constructor & Destructor Documentation

G4ParticleHPWattSpectrum()

G4ParticleHPWattSpectrum::G4ParticleHPWattSpectrum ( )

inline

Definition at line 47 of file G4ParticleHPWattSpectrum.hh.

  {
    expm1 = G4Exp(-1.);
  }

Here is the call graph for this function:

~G4ParticleHPWattSpectrum()

G4ParticleHPWattSpectrum::~G4ParticleHPWattSpectrum ( )

inline

Definition at line 51 of file G4ParticleHPWattSpectrum.hh.

  {
  }

Member Function Documentation

GetFractionalProbability()

G4double G4ParticleHPWattSpectrum::GetFractionalProbability ( G4double  anEnergy )

inlinevirtual

Implements G4VParticleHPEDis.

Definition at line 62 of file G4ParticleHPWattSpectrum.hh.

  {
    return theFractionalProb.GetY(anEnergy);
  }

Here is the call graph for this function:

Init()

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

inlinevirtual

Implements G4VParticleHPEDis.

Definition at line 55 of file G4ParticleHPWattSpectrum.hh.

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

Here is the call graph for this function:

Sample()

G4double G4ParticleHPWattSpectrum::Sample ( G4double  anEnergy )

virtual

Implements G4VParticleHPEDis.

Definition at line 35 of file G4ParticleHPWattSpectrum.cc.

  {
    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));

    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;
  }

Here is the call graph for this function:

Here is the caller graph for this function:

Watt()

G4double G4ParticleHPWattSpectrum::Watt ( G4double  anEnergy, G4double  a, G4double  b  )

inlineprivate

Definition at line 71 of file G4ParticleHPWattSpectrum.hh.

  {
    G4double energy = anEnergy/CLHEP::eV;
    G4double result = G4Exp(-energy/a)*std::sinh(std::sqrt(b*energy));
    return result;
  }

Here is the call graph for this function:

Member Data Documentation

expm1

G4double G4ParticleHPWattSpectrum::expm1

private

Definition at line 80 of file G4ParticleHPWattSpectrum.hh.

theApar

G4ParticleHPVector G4ParticleHPWattSpectrum::theApar

private

Definition at line 84 of file G4ParticleHPWattSpectrum.hh.

theBpar

G4ParticleHPVector G4ParticleHPWattSpectrum::theBpar

private

Definition at line 85 of file G4ParticleHPWattSpectrum.hh.

theFractionalProb

G4ParticleHPVector G4ParticleHPWattSpectrum::theFractionalProb

private

Definition at line 82 of file G4ParticleHPWattSpectrum.hh.

---

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

• G4ParticleHPWattSpectrum.hh
• G4ParticleHPWattSpectrum.cc