CLHEP::RandGaussZiggurat Class Reference

#include <RandGaussZiggurat.h>

Inheritance diagram for CLHEP::RandGaussZiggurat:

Public Member Functions
	RandGaussZiggurat (HepRandomEngine &anEngine, double mean=0.0, double stdDev=1.0)
	RandGaussZiggurat (HepRandomEngine *anEngine, double mean=0.0, double stdDev=1.0)
virtual	~RandGaussZiggurat ()
float	fire ()
float	fire (float mean, float stdDev)
void	fireArray (const int size, float *vect)
void	fireArray (const int size, double *vect)
void	fireArray (const int size, float *vect, float mean, float stdDev)
void	fireArray (const int size, double *vect, double mean, double stdDev)
virtual double	operator() ()
virtual double	operator() (double mean, double stdDev)
std::ostream &	put (std::ostream &os) const
std::istream &	get (std::istream &is)
std::string	name () const
HepRandomEngine &	engine ()
Public Member Functions inherited from CLHEP::RandGauss
	RandGauss (HepRandomEngine &anEngine, double mean=0.0, double stdDev=1.0)
	RandGauss (HepRandomEngine *anEngine, double mean=0.0, double stdDev=1.0)
virtual	~RandGauss ()
double	fire ()
double	fire (double mean, double stdDev)
void	fireArray (const int size, double *vect)
void	fireArray (const int size, double *vect, double mean, double stdDev)
bool	getF () const
void	setF (bool val)
Public Member Functions inherited from CLHEP::HepRandom
	HepRandom ()
	HepRandom (long seed)
	HepRandom (HepRandomEngine &algorithm)
	HepRandom (HepRandomEngine *algorithm)
virtual	~HepRandom ()
double	flat ()
void	flatArray (const int size, double *vect)
double	flat (HepRandomEngine *theNewEngine)
void	flatArray (HepRandomEngine *theNewEngine, const int size, double *vect)

Static Public Member Functions
static float	shoot ()
static float	shoot (float mean, float stdDev)
static void	shootArray (const int size, float *vect, float mean=0.0, float stdDev=1.0)
static void	shootArray (const int size, double *vect, double mean=0.0, double stdDev=1.0)
static float	shoot (HepRandomEngine *anotherEngine)
static float	shoot (HepRandomEngine *anotherEngine, float mean, float stdDev)
static void	shootArray (HepRandomEngine *anotherEngine, const int size, float *vect, float mean=0.0, float stdDev=1.0)
static void	shootArray (HepRandomEngine *anotherEngine, const int size, double *vect, double mean=0.0, double stdDev=1.0)
static std::string	distributionName ()
static bool	ziggurat_init ()
Static Public Member Functions inherited from CLHEP::RandGauss
static double	shoot ()
static double	shoot (double mean, double stdDev)
static void	shootArray (const int size, double *vect, double mean=0.0, double stdDev=1.0)
static double	shoot (HepRandomEngine *anEngine)
static double	shoot (HepRandomEngine *anEngine, double mean, double stdDev)
static void	shootArray (HepRandomEngine *anEngine, const int size, double *vect, double mean=0.0, double stdDev=1.0)
static std::string	distributionName ()
static bool	getFlag ()
static void	setFlag (bool val)
static void	saveEngineStatus (const char filename[]="Config.conf")
static void	restoreEngineStatus (const char filename[]="Config.conf")
static std::ostream &	saveFullState (std::ostream &os)
static std::istream &	restoreFullState (std::istream &is)
static std::ostream &	saveDistState (std::ostream &os)
static std::istream &	restoreDistState (std::istream &is)
Static Public Member Functions inherited from CLHEP::HepRandom
static void	setTheSeed (long seed, int lux=3)
static long	getTheSeed ()
static void	setTheSeeds (const long *seeds, int aux=-1)
static const long *	getTheSeeds ()
static void	getTheTableSeeds (long *seeds, int index)
static HepRandom *	getTheGenerator ()
static void	setTheEngine (HepRandomEngine *theNewEngine)
static HepRandomEngine *	getTheEngine ()
static void	saveEngineStatus (const char filename[]="Config.conf")
static void	restoreEngineStatus (const char filename[]="Config.conf")
static std::ostream &	saveFullState (std::ostream &os)
static std::istream &	restoreFullState (std::istream &is)
static std::ostream &	saveDistState (std::ostream &os)
static std::istream &	restoreDistState (std::istream &is)
static std::ostream &	saveStaticRandomStates (std::ostream &os)
static std::istream &	restoreStaticRandomStates (std::istream &is)
static void	showEngineStatus ()
static int	createInstance ()
static std::string	distributionName ()

Static Protected Member Functions
static unsigned long	ziggurat_SHR3 (HepRandomEngine *anEngine)
static float	ziggurat_UNI (HepRandomEngine *anEngine)
static float	ziggurat_RNOR (HepRandomEngine *anEngine)
static float	ziggurat_nfix (long hz, HepRandomEngine *anEngine)
Static Protected Member Functions inherited from CLHEP::RandGauss
static double	getVal ()
static void	setVal (double nextVal)

Static Protected Attributes
static unsigned long	kn [128]
static unsigned long	ke [256]
static float	wn [128]
static float	fn [128]
static float	we [256]
static float	fe [256]
static bool	ziggurat_is_init =RandGaussZiggurat::ziggurat_init()
Static Protected Attributes inherited from CLHEP::HepRandom
static const long	seedTable [215][2]

Additional Inherited Members
Protected Member Functions inherited from CLHEP::RandGauss
double	normal ()
Protected Attributes inherited from CLHEP::RandGauss
double	defaultMean
double	defaultStdDev
shared_ptr< HepRandomEngine >	localEngine

Detailed Description

Author

ATLAS

Definition at line 36 of file RandGaussZiggurat.h.

Constructor & Destructor Documentation

CLHEP::RandGaussZiggurat::RandGaussZiggurat ( HepRandomEngine &  anEngine, double  mean = 0.0, double  stdDev = 1.0  )

inline

Definition at line 135 of file RandGaussZiggurat.h.

References ziggurat_init(), and ziggurat_is_init.

                                                                                          : RandGauss(anEngine, mean, stdDev) 
{
  if(!ziggurat_is_init) ziggurat_init();
}

CLHEP::RandGaussZiggurat::RandGaussZiggurat ( HepRandomEngine *  anEngine, double  mean = 0.0, double  stdDev = 1.0  )

inline

Definition at line 140 of file RandGaussZiggurat.h.

References ziggurat_init(), and ziggurat_is_init.

                                                                                          : RandGauss(anEngine, mean, stdDev) 
{
  if(!ziggurat_is_init) ziggurat_init();
}

CLHEP::RandGaussZiggurat::~RandGaussZiggurat ( )

virtual

Definition at line 15 of file RandGaussZiggurat.cc.

                                      {
}

Member Function Documentation

static std::string CLHEP::RandGaussZiggurat::distributionName ( )

inlinestatic

Definition at line 86 of file RandGaussZiggurat.h.

{return "RandGaussZiggurat";}  

HepRandomEngine & CLHEP::RandGaussZiggurat::engine ( )

virtual

Reimplemented from CLHEP::RandGauss.

Definition at line 13 of file RandGaussZiggurat.cc.

References CLHEP::RandGauss::engine().

{return RandGauss::engine();}

float CLHEP::RandGaussZiggurat::fire ( )

inline

Definition at line 66 of file RandGaussZiggurat.h.

References CLHEP::RandGauss::defaultMean, CLHEP::RandGauss::defaultStdDev, CLHEP::RandGauss::localEngine, and ziggurat_RNOR().

Referenced by fireArray().

{return ziggurat_RNOR(localEngine.get()) * defaultStdDev + defaultMean;};

float CLHEP::RandGaussZiggurat::fire ( float  mean, float  stdDev  )

inline

Definition at line 68 of file RandGaussZiggurat.h.

References CLHEP::RandGauss::localEngine, and ziggurat_RNOR().

{return ziggurat_RNOR(localEngine.get()) * stdDev + mean;};

void CLHEP::RandGaussZiggurat::fireArray	(	const int	size,
		float *	vect
	)

Definition at line 137 of file RandGaussZiggurat.cc.

References CLHEP::RandGauss::defaultMean, CLHEP::RandGauss::defaultStdDev, and fire().

{
   for (int i=0; i<size; ++i) {
     vect[i] = fire( defaultMean, defaultStdDev );
   }
}

void CLHEP::RandGaussZiggurat::fireArray	(	const int	size,
		double *	vect
	)

Definition at line 144 of file RandGaussZiggurat.cc.

References CLHEP::RandGauss::defaultMean, CLHEP::RandGauss::defaultStdDev, and fire().

{
   for (int i=0; i<size; ++i) {
     vect[i] = fire( defaultMean, defaultStdDev );
   }
}

void CLHEP::RandGaussZiggurat::fireArray	(	const int	size,
		float *	vect,
		float	mean,
		float	stdDev
	)

Definition at line 151 of file RandGaussZiggurat.cc.

References fire().

{
   for (int i=0; i<size; ++i) {
     vect[i] = fire( mean, stdDev );
   }
}

void CLHEP::RandGaussZiggurat::fireArray	(	const int	size,
		double *	vect,
		double	mean,
		double	stdDev
	)

Definition at line 158 of file RandGaussZiggurat.cc.

References fire().

{
   for (int i=0; i<size; ++i) {
     vect[i] = fire( mean, stdDev );
   }
}

std::istream & CLHEP::RandGaussZiggurat::get ( std::istream &  is )

virtual

Reimplemented from CLHEP::RandGauss.

Definition at line 173 of file RandGaussZiggurat.cc.

References CLHEP::RandGauss::get(), and name().

                                                      {
  std::string inName;
  is >> inName;
  if (inName != name()) {
    is.clear(std::ios::badbit | is.rdstate());
    std::cerr << "Mismatch when expecting to read state of a "
              << name() << " distribution\n"
          << "Name found was " << inName
          << "\nistream is left in the badbit state\n";
    return is;
  }
  RandGauss::get(is);
  return is;
}

std::string CLHEP::RandGaussZiggurat::name ( ) const

virtual

Reimplemented from CLHEP::RandGauss.

Definition at line 18 of file RandGaussZiggurat.cc.

Referenced by get(), and put().

{
  return "RandGaussZiggurat";
}

double CLHEP::RandGaussZiggurat::operator() ( )

virtual

Reimplemented from CLHEP::RandGauss.

Definition at line 101 of file RandGaussZiggurat.cc.

References CLHEP::RandGauss::defaultMean, CLHEP::RandGauss::defaultStdDev, CLHEP::RandGauss::localEngine, and ziggurat_RNOR().

                                     {
  return ziggurat_RNOR(localEngine.get()) * defaultStdDev + defaultMean;
}

double CLHEP::RandGaussZiggurat::operator() ( double  mean, double  stdDev  )

virtual

Reimplemented from CLHEP::RandGauss.

Definition at line 105 of file RandGaussZiggurat.cc.

References CLHEP::RandGauss::localEngine, and ziggurat_RNOR().

                                                                 {
  return ziggurat_RNOR(localEngine.get()) * stdDev + mean;
}

std::ostream & CLHEP::RandGaussZiggurat::put ( std::ostream &  os ) const

virtual

Reimplemented from CLHEP::RandGauss.

Definition at line 165 of file RandGaussZiggurat.cc.

References name(), gammaraytel::pr, and CLHEP::RandGauss::put().

                                                            {
  int pr=os.precision(20);
  os << " " << name() << "\n";
  RandGauss::put(os);
  os.precision(pr);
  return os;
}

static float CLHEP::RandGaussZiggurat::shoot ( )

inlinestatic

Definition at line 48 of file RandGaussZiggurat.h.

References CLHEP::HepRandom::getTheEngine(), and ziggurat_RNOR().

Referenced by shootArray().

{return ziggurat_RNOR(HepRandom::getTheEngine());};

static float CLHEP::RandGaussZiggurat::shoot ( float  mean, float  stdDev  )

inlinestatic

Definition at line 49 of file RandGaussZiggurat.h.

References shoot().

Referenced by shoot().

{return shoot()*stdDev + mean;};

static float CLHEP::RandGaussZiggurat::shoot ( HepRandomEngine *  anotherEngine )

inlinestatic

Definition at line 57 of file RandGaussZiggurat.h.

References ziggurat_RNOR().

{return ziggurat_RNOR(anotherEngine);};

static float CLHEP::RandGaussZiggurat::shoot ( HepRandomEngine *  anotherEngine, float  mean, float  stdDev  )

inlinestatic

Definition at line 58 of file RandGaussZiggurat.h.

References shoot().

Referenced by shoot().

{return shoot(anotherEngine)*stdDev + mean;};

void CLHEP::RandGaussZiggurat::shootArray ( const int  size, float *  vect, float  mean = 0.0, float  stdDev = 1.0  )

static

Definition at line 109 of file RandGaussZiggurat.cc.

References shoot().

{
   for (int i=0; i<size; ++i) {
     vect[i] = shoot(mean,stdDev);
   }
}

void CLHEP::RandGaussZiggurat::shootArray ( const int  size, double *  vect, double  mean = 0.0, double  stdDev = 1.0  )

static

Definition at line 116 of file RandGaussZiggurat.cc.

References shoot().

{
   for (int i=0; i<size; ++i) {
     vect[i] = shoot(mean,stdDev);
   }
}

void CLHEP::RandGaussZiggurat::shootArray ( HepRandomEngine *  anotherEngine, const int  size, float *  vect, float  mean = 0.0, float  stdDev = 1.0  )

static

Definition at line 123 of file RandGaussZiggurat.cc.

References shoot().

{
   for (int i=0; i<size; ++i) {
     vect[i] = shoot(anEngine,mean,stdDev);
   }
}

void CLHEP::RandGaussZiggurat::shootArray ( HepRandomEngine *  anotherEngine, const int  size, double *  vect, double  mean = 0.0, double  stdDev = 1.0  )

static

Definition at line 130 of file RandGaussZiggurat.cc.

References shoot().

{
   for (int i=0; i<size; ++i) {
     vect[i] = shoot(anEngine,mean,stdDev);
   }
}

bool CLHEP::RandGaussZiggurat::ziggurat_init ( )

static

Definition at line 23 of file RandGaussZiggurat.cc.

References fe, fn, ke, kn, we, wn, and ziggurat_is_init.

Referenced by RandGaussZiggurat().

{  
  const double rzm1 = 2147483648.0, rzm2 = 4294967296.;
  double dn=3.442619855899,tn=dn,vn=9.91256303526217e-3, q;
  double de=7.697117470131487, te=de, ve=3.949659822581572e-3;
  int i;

/* Set up tables for RNOR */
  q=vn/std::exp(-.5*dn*dn);
  kn[0]=(unsigned long)((dn/q)*rzm1);
  kn[1]=0;

  wn[0]=q/rzm1;
  wn[127]=dn/rzm1;

  fn[0]=1.;
  fn[127]=std::exp(-.5*dn*dn);

  for(i=126;i>=1;i--) {
    dn=std::sqrt(-2.*std::log(vn/dn+std::exp(-.5*dn*dn)));
    kn[i+1]=(unsigned long)((dn/tn)*rzm1);
    tn=dn;
    fn[i]=std::exp(-.5*dn*dn);
    wn[i]=dn/rzm1;
  }

/* Set up tables for REXP */
  q = ve/std::exp(-de);
  ke[0]=(unsigned long)((de/q)*rzm2);
  ke[1]=0;

  we[0]=q/rzm2;
  we[255]=de/rzm2;

  fe[0]=1.;
  fe[255]=std::exp(-de);

  for(i=254;i>=1;i--) {
    de=-std::log(ve/de+std::exp(-de));
    ke[i+1]= (unsigned long)((de/te)*rzm2);
    te=de;
    fe[i]=std::exp(-de);
    we[i]=de/rzm2;
  }
  ziggurat_is_init=true;
  
  //std::cout<<"Done RandGaussZiggurat::ziggurat_init()"<<std::endl;
  
  return true;
}

float CLHEP::RandGaussZiggurat::ziggurat_nfix ( long  hz, HepRandomEngine *  anEngine  )

staticprotected

Definition at line 74 of file RandGaussZiggurat.cc.

References fn, iz, kn, wn, test::x, ziggurat_SHR3(), and ziggurat_UNI().

Referenced by ziggurat_RNOR().

{
  const float r = 3.442620f;     /* The start of the right tail */
  float x, y;
  unsigned long iz=hz&127;
  for(;;)
  {  
    x=hz*wn[iz];      /* iz==0, handles the base strip */
    if(iz==0) {
      do { 
        /* change to (1.0 - UNI) as argument to std::log(), because CLHEP generates [0,1), 
           while the original UNI generates (0,1] */
        x=-std::log(1.0 - ziggurat_UNI(anEngine))*0.2904764; /* .2904764 is 1/r */
        y=-std::log(1.0 - ziggurat_UNI(anEngine));
      } while(y+y<x*x);
      return (hz>0)? r+x : -r-x;
    }
    /* iz>0, handle the wedges of other strips */
    if( fn[iz]+(1.0 - ziggurat_UNI(anEngine))*(fn[iz-1]-fn[iz]) < std::exp(-.5*x*x) ) return x;

    /* initiate, try to exit for(;;) for loop*/
    hz=(signed)ziggurat_SHR3(anEngine);
    iz=hz&127;
    if((unsigned long)abs(hz)<kn[iz]) return (hz*wn[iz]);
  }
}

static float CLHEP::RandGaussZiggurat::ziggurat_RNOR ( HepRandomEngine *  anEngine )

inlinestaticprotected

Definition at line 114 of file RandGaussZiggurat.h.

References iz, kn, wn, ziggurat_nfix(), and ziggurat_SHR3().

Referenced by fire(), operator()(), and shoot().

                                                               {
    long hz=(signed)ziggurat_SHR3(anEngine);
    unsigned long iz=hz&127;
    return ((unsigned long)abs(hz)<kn[iz]) ? hz*wn[iz] : ziggurat_nfix(hz,anEngine);
  };

static unsigned long CLHEP::RandGaussZiggurat::ziggurat_SHR3 ( HepRandomEngine *  anEngine )

inlinestaticprotected

Definition at line 112 of file RandGaussZiggurat.h.

Referenced by ziggurat_nfix(), and ziggurat_RNOR().

{return (unsigned int)(*anEngine);};

static float CLHEP::RandGaussZiggurat::ziggurat_UNI ( HepRandomEngine *  anEngine )

inlinestaticprotected

Definition at line 113 of file RandGaussZiggurat.h.

References CLHEP::HepRandomEngine::flat().

Referenced by ziggurat_nfix().

{return anEngine->flat();};

Field Documentation

float CLHEP::RandGaussZiggurat::fe

staticprotected

Definition at line 108 of file RandGaussZiggurat.h.

Referenced by ziggurat_init().

float CLHEP::RandGaussZiggurat::fn

staticprotected

Definition at line 108 of file RandGaussZiggurat.h.

Referenced by ziggurat_init(), and ziggurat_nfix().

unsigned long CLHEP::RandGaussZiggurat::ke

staticprotected

Definition at line 107 of file RandGaussZiggurat.h.

Referenced by ziggurat_init().

unsigned long CLHEP::RandGaussZiggurat::kn

staticprotected

Definition at line 107 of file RandGaussZiggurat.h.

Referenced by ziggurat_init(), ziggurat_nfix(), and ziggurat_RNOR().

float CLHEP::RandGaussZiggurat::we

staticprotected

Definition at line 108 of file RandGaussZiggurat.h.

Referenced by ziggurat_init().

float CLHEP::RandGaussZiggurat::wn

staticprotected

Definition at line 108 of file RandGaussZiggurat.h.

Referenced by ziggurat_init(), ziggurat_nfix(), and ziggurat_RNOR().

bool CLHEP::RandGaussZiggurat::ziggurat_is_init =RandGaussZiggurat::ziggurat_init()

staticprotected

Definition at line 110 of file RandGaussZiggurat.h.

Referenced by RandGaussZiggurat(), and ziggurat_init().

---

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

• RandGaussZiggurat.h
• RandGaussZiggurat.cc