G4ParticleHPWattSpectrum Class Reference

#include <G4ParticleHPWattSpectrum.hh>

Inheritance 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 ()
virtual void	Init (std::istream &theData)=0
virtual G4double	GetFractionalProbability (G4double anEnergy)=0
virtual G4double	Sample (G4double anEnergy)=0

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

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

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

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

---

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

• G4ParticleHPWattSpectrum.hh
• G4ParticleHPWattSpectrum.cc