G4ParticleHPWattSpectrum Class Reference

#include <G4ParticleHPWattSpectrum.hh>

Inheritance diagram for G4ParticleHPWattSpectrum:

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

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

Private Attributes
G4double	expm1
G4ParticleHPVector	theApar
G4ParticleHPVector	theBpar
G4ParticleHPVector	theFractionalProb

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

References expm1, and G4Exp().

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

References G4ParticleHPVector::GetY(), and theFractionalProb.

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

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

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

References eV, G4cout, G4endl, G4Log(), G4UniformRand, G4ParticleHPVector::GetY(), G4INCL::Math::max(), theApar, and theBpar.

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

References G4INCL::KinematicsUtils::energy(), CLHEP::eV, and G4Exp().

Field Documentation

expm1

G4double G4ParticleHPWattSpectrum::expm1

private

Definition at line 80 of file G4ParticleHPWattSpectrum.hh.

Referenced by G4ParticleHPWattSpectrum().

theApar

G4ParticleHPVector G4ParticleHPWattSpectrum::theApar

private

Definition at line 84 of file G4ParticleHPWattSpectrum.hh.

Referenced by Init(), and Sample().

theBpar

G4ParticleHPVector G4ParticleHPWattSpectrum::theBpar

private

Definition at line 85 of file G4ParticleHPWattSpectrum.hh.

Referenced by Init(), and Sample().

theFractionalProb

G4ParticleHPVector G4ParticleHPWattSpectrum::theFractionalProb

private

Definition at line 82 of file G4ParticleHPWattSpectrum.hh.

Referenced by GetFractionalProbability(), and Init().

---

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

• source/processes/hadronic/models/particle_hp/include/G4ParticleHPWattSpectrum.hh
• source/processes/hadronic/models/particle_hp/src/G4ParticleHPWattSpectrum.cc