G4Poisson.hh File Reference

#include <CLHEP/Units/PhysicalConstants.h>
#include "G4Types.hh"
#include "G4Exp.hh"
#include "Randomize.hh"

Include dependency graph for G4Poisson.hh:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions
G4long	G4Poisson (G4double mean)

Function Documentation

G4long G4Poisson ( G4double  mean )

inline

Definition at line 51 of file G4Poisson.hh.

{
  G4long number = 0;
  const G4int border = 16;
  const G4double limit = 2e9;

  if(mean <= border)
  {
    G4double position = G4UniformRand();
    G4double poissonValue = G4Exp(-mean);
    G4double poissonSum = poissonValue;

    while(poissonSum <= position)
    {
      ++number;
      poissonValue *= mean/number;
      poissonSum += poissonValue;
    }
    return number;
  } // the case of mean <= 16

  G4double t = std::sqrt(-2.*std::log(G4UniformRand()))*
               std::cos(2.*CLHEP::pi*G4UniformRand());
  G4double value = mean + t*std::sqrt(mean) + 0.5;
  if(value < 0.) {return 0;}
  return (value >= limit) ? G4long(limit) : G4long(value);
}

Here is the call graph for this function:

Here is the caller graph for this function: