24#ifndef RandGaussZiggurat_h
25#define RandGaussZiggurat_h 1
28#include "CLHEP/Random/defs.h"
29#include "CLHEP/Random/RandGauss.h"
37class RandGaussZiggurat :
public RandGauss {
50 static inline float shoot(
float mean,
float stdDev ) {
return shoot()*stdDev + mean;};
52 static void shootArray (
const int size,
float* vect,
float mean=0.0,
float stdDev=1.0 );
53 static void shootArray (
const int size,
double* vect,
double mean=0.0,
double stdDev=1.0 );
73 void fireArray (
const int size,
float* vect,
float mean,
float stdDev );
74 void fireArray (
const int size,
double* vect,
double mean,
double stdDev );
81 std::ostream &
put ( std::ostream & os )
const;
82 std::istream &
get ( std::istream & is );
108 static unsigned long kn[128],
ke[256];
109 static float wn[128],
fn[128],
we[256],
fe[256];
117 unsigned long iz=hz&127;
134#ifdef ENABLE_BACKWARDS_COMPATIBILITY
136using namespace CLHEP;
static HepRandomEngine * getTheEngine()
static float shoot(float mean, float stdDev)
void fireArray(const int size, double *vect, double mean, double stdDev)
std::istream & get(std::istream &is)
static bool ziggurat_init()
void fireArray(const int size, float *vect)
RandGaussZiggurat(HepRandomEngine *anEngine, double mean=0.0, double stdDev=1.0)
static void shootArray(HepRandomEngine *anotherEngine, const int size, double *vect, double mean=0.0, double stdDev=1.0)
static float ziggurat_nfix(long hz, HepRandomEngine *anEngine)
static float ziggurat_UNI(HepRandomEngine *anEngine)
static bool ziggurat_init()
void fireArray(const int size, double *vect)
std::ostream & put(std::ostream &os) const
static void shootArray(HepRandomEngine *anotherEngine, const int size, float *vect, float mean=0.0, float stdDev=1.0)
static float ziggurat_RNOR(HepRandomEngine *anEngine)
virtual double operator()(double mean, double stdDev)
static unsigned long ke[256]
float fire(float mean, float stdDev)
void fireArray(const int size, float *vect, float mean, float stdDev)
static unsigned long kn[128]
static bool ziggurat_is_init
virtual ~RandGaussZiggurat()
HepRandomEngine & engine()
static float shoot(HepRandomEngine *anotherEngine, float mean, float stdDev)
RandGaussZiggurat(HepRandomEngine &anEngine, double mean=0.0, double stdDev=1.0)
static float shoot(HepRandomEngine *anotherEngine)
static unsigned long ziggurat_SHR3(HepRandomEngine *anEngine)
virtual double operator()()
static void shootArray(const int size, double *vect, double mean=0.0, double stdDev=1.0)
static void shootArray(const int size, float *vect, float mean=0.0, float stdDev=1.0)
static std::string distributionName()
static float ziggurat_nfix(long hz, HepRandomEngine *anEngine)
shared_ptr< HepRandomEngine > localEngine