G4BetaSpectrumSampler Class Reference

#include <G4BetaSpectrumSampler.hh>

Static Public Member Functions
static G4double	shoot (const G4int npoints, const G4double *aCDF, const G4double estep)

Detailed Description

Definition at line 42 of file G4BetaSpectrumSampler.hh.

Member Function Documentation

shoot()

G4double G4BetaSpectrumSampler::shoot ( const G4int npoints, const G4double * aCDF, const G4double estep )

static

Definition at line 40 of file G4BetaSpectrumSampler.cc.

{
  G4double prob = aCDF[npoints - 1]*G4UniformRand();
  G4int i = 0;
  for (; i<npoints; ++i) { if (prob <= aCDF[i]) { break; } }
  const G4double p1 = (i > 0) ? aCDF[i - 1] : aCDF[0];
  const G4double p2 = aCDF[i];
  const G4double delta = p2 - p1;
  const G4double x = (delta > 0.0) ? estep*i - estep*(p2 - prob)/delta : estep*i;
  return x;
}

Referenced by G4BetaMinusDecay::DecayIt(), and G4BetaPlusDecay::DecayIt().

---

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

• G4BetaSpectrumSampler.hh
• G4BetaSpectrumSampler.cc