CLHEP::RandExpZiggurat Class Reference

#include <RandExpZiggurat.h>

Inheritance diagram for CLHEP::RandExpZiggurat:

Public Member Functions
	RandExpZiggurat (HepRandomEngine &anEngine, double mean=1.0)
	RandExpZiggurat (HepRandomEngine *anEngine, double mean=1.0)
virtual	~RandExpZiggurat ()
float	fire ()
float	fire (float mean)
void	fireArray (const int size, float *vect)
void	fireArray (const int size, double *vect)
void	fireArray (const int size, float *vect, float mean)
void	fireArray (const int size, double *vect, double mean)
virtual double	operator() ()
float	operator() (float mean)
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::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)
static void	shootArray (const int size, float *vect, float mean=1.0)
static void	shootArray (const int size, double *vect, double mean=1.0)
static float	shoot (HepRandomEngine *anEngine)
static float	shoot (HepRandomEngine *anEngine, float mean)
static void	shootArray (HepRandomEngine *anEngine, const int size, float *vect, float mean=1.0)
static void	shootArray (HepRandomEngine *anEngine, const int size, double *vect, double mean=1.0)
static std::string	distributionName ()
static bool	ziggurat_init ()
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_REXP (HepRandomEngine *anEngine)
static float	ziggurat_efix (unsigned long jz, HepRandomEngine *anEngine)

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 =RandExpZiggurat::ziggurat_init()
Static Protected Attributes inherited from CLHEP::HepRandom
static const long	seedTable [215][2]

Detailed Description

Author

ATLAS

Definition at line 34 of file RandExpZiggurat.h.

Constructor & Destructor Documentation

CLHEP::RandExpZiggurat::RandExpZiggurat ( HepRandomEngine &  anEngine, double  mean = 1.0  )

inline

Definition at line 156 of file RandExpZiggurat.h.

References ziggurat_init(), and ziggurat_is_init.

                                                                                : localEngine(&anEngine), deleteEngine(false), defaultMean(mean) 
{
  if(!ziggurat_is_init) ziggurat_init();
}

CLHEP::RandExpZiggurat::RandExpZiggurat ( HepRandomEngine *  anEngine, double  mean = 1.0  )

inline

Definition at line 161 of file RandExpZiggurat.h.

References ziggurat_init(), and ziggurat_is_init.

                                                                                : localEngine(anEngine), deleteEngine(true), defaultMean(mean) 
{
  if(!ziggurat_is_init) ziggurat_init();
}

CLHEP::RandExpZiggurat::~RandExpZiggurat ( )

virtual

Definition at line 18 of file RandExpZiggurat.cc.

                                  {
  if ( deleteEngine ) delete localEngine;
}

Member Function Documentation

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

inlinestatic

Definition at line 107 of file RandExpZiggurat.h.

{return "RandExpZiggurat";}  

HepRandomEngine & CLHEP::RandExpZiggurat::engine ( )

virtual

Reimplemented from CLHEP::HepRandom.

Definition at line 16 of file RandExpZiggurat.cc.

{return *localEngine;}

float CLHEP::RandExpZiggurat::fire ( )

inline

Definition at line 83 of file RandExpZiggurat.h.

References fire().

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

{return fire(defaultMean);};

float CLHEP::RandExpZiggurat::fire ( float  mean )

inline

Definition at line 84 of file RandExpZiggurat.h.

References ziggurat_REXP().

{return ziggurat_REXP(localEngine)*mean;};

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

Definition at line 52 of file RandExpZiggurat.cc.

References fire().

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

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

Definition at line 57 of file RandExpZiggurat.cc.

References fire().

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

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

Definition at line 62 of file RandExpZiggurat.cc.

References fire().

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

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

Definition at line 67 of file RandExpZiggurat.cc.

References fire().

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

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

virtual

Reimplemented from CLHEP::HepRandom.

Definition at line 90 of file RandExpZiggurat.cc.

References CLHEP::DoubConv::longs2double(), name(), and CLHEP::possibleKeywordInput().

                                                    {
  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;
  }
  if (possibleKeywordInput(is, "Uvec", defaultMean)) {
    std::vector<unsigned long> t(2);
    is >> defaultMean >> t[0] >> t[1]; defaultMean = DoubConv::longs2double(t); 
    return is;
  }
  // is >> defaultMean encompassed by possibleKeywordInput
  return is;
}

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

virtual

Reimplemented from CLHEP::HepRandom.

Definition at line 14 of file RandExpZiggurat.cc.

Referenced by get(), and put().

{return "RandExpZiggurat";}

double CLHEP::RandExpZiggurat::operator() ( )

virtual

Reimplemented from CLHEP::HepRandom.

Definition at line 27 of file RandExpZiggurat.cc.

References fire().

{
  return fire( defaultMean );
}

float CLHEP::RandExpZiggurat::operator() ( float  mean )

inline

Definition at line 97 of file RandExpZiggurat.h.

References fire().

{return fire( mean );};

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

virtual

Reimplemented from CLHEP::HepRandom.

Definition at line 72 of file RandExpZiggurat.cc.

References CLHEP::DoubConv::dto2longs(), name(), and gammaraytel::pr.

                                                          {
  int pr=os.precision(20);
  std::vector<unsigned long> t(2);
  os << " " << name() << "\n";
  os << "Uvec" << "\n";
  t = DoubConv::dto2longs(defaultMean);
  os << defaultMean << " " << t[0] << " " << t[1] << "\n";
  os.precision(pr);
  return os;
#ifdef REMOVED
  int pr=os.precision(20);
  os << " " << name() << "\n";
  os << defaultMean << "\n";
  os.precision(pr);
  return os;
#endif
}

static float CLHEP::RandExpZiggurat::shoot ( )

inlinestatic

Definition at line 54 of file RandExpZiggurat.h.

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

Referenced by shoot(), and shootArray().

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

static float CLHEP::RandExpZiggurat::shoot ( float  mean )

inlinestatic

Definition at line 55 of file RandExpZiggurat.h.

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

Referenced by shoot().

{return shoot(HepRandom::getTheEngine(),mean);};

static float CLHEP::RandExpZiggurat::shoot ( HepRandomEngine *  anEngine )

inlinestatic

Definition at line 68 of file RandExpZiggurat.h.

References ziggurat_REXP().

{return ziggurat_REXP(anEngine);};

static float CLHEP::RandExpZiggurat::shoot ( HepRandomEngine *  anEngine, float  mean  )

inlinestatic

Definition at line 69 of file RandExpZiggurat.h.

References shoot().

Referenced by shoot().

{return shoot(anEngine)*mean;};

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

static

Definition at line 32 of file RandExpZiggurat.cc.

References shoot().

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

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

static

Definition at line 37 of file RandExpZiggurat.cc.

References shoot().

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

void CLHEP::RandExpZiggurat::shootArray ( HepRandomEngine *  anEngine, const int  size, float *  vect, float  mean = 1.0  )

static

Definition at line 42 of file RandExpZiggurat.cc.

References shoot().

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

void CLHEP::RandExpZiggurat::shootArray ( HepRandomEngine *  anEngine, const int  size, double *  vect, double  mean = 1.0  )

static

Definition at line 47 of file RandExpZiggurat.cc.

References shoot().

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

float CLHEP::RandExpZiggurat::ziggurat_efix ( unsigned long  jz, HepRandomEngine *  anEngine  )

staticprotected

Definition at line 111 of file RandExpZiggurat.cc.

References fe, iz, ke, we, test::x, ziggurat_SHR3(), and ziggurat_UNI().

Referenced by ziggurat_REXP().

{ 
  unsigned long iz=jz&255;
  
  float x;
  for(;;)
  {  
    if(iz==0) return (7.69711-std::log(ziggurat_UNI(anEngine)));          /* iz==0 */
    x=jz*we[iz]; 
    if( fe[iz]+ziggurat_UNI(anEngine)*(fe[iz-1]-fe[iz]) < std::exp(-x) ) return (x);

    /* initiate, try to exit for(;;) loop */
    jz=ziggurat_SHR3(anEngine);
    iz=(jz&255);
    if(jz<ke[iz]) return (jz*we[iz]);
  }
}

bool CLHEP::RandExpZiggurat::ziggurat_init ( )

static

Definition at line 129 of file RandExpZiggurat.cc.

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

Referenced by RandExpZiggurat().

{  
  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;
  return true;
}

static float CLHEP::RandExpZiggurat::ziggurat_REXP ( HepRandomEngine *  anEngine )

inlinestaticprotected

Definition at line 135 of file RandExpZiggurat.h.

References iz, ke, we, ziggurat_efix(), and ziggurat_SHR3().

Referenced by fire(), and shoot().

                                                               {
    unsigned long jz=ziggurat_SHR3(anEngine);
    unsigned long iz=jz&255;
    return (jz<ke[iz]) ? jz*we[iz] : ziggurat_efix(jz,anEngine);
  };

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

inlinestaticprotected

Definition at line 133 of file RandExpZiggurat.h.

Referenced by ziggurat_efix(), and ziggurat_REXP().

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

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

inlinestaticprotected

Definition at line 134 of file RandExpZiggurat.h.

References CLHEP::HepRandomEngine::flat().

Referenced by ziggurat_efix().

{return anEngine->flat();};

Field Documentation

float CLHEP::RandExpZiggurat::fe

staticprotected

Definition at line 129 of file RandExpZiggurat.h.

Referenced by ziggurat_efix(), and ziggurat_init().

float CLHEP::RandExpZiggurat::fn

staticprotected

Definition at line 129 of file RandExpZiggurat.h.

Referenced by ziggurat_init().

unsigned long CLHEP::RandExpZiggurat::ke

staticprotected

Definition at line 128 of file RandExpZiggurat.h.

Referenced by ziggurat_efix(), ziggurat_init(), and ziggurat_REXP().

unsigned long CLHEP::RandExpZiggurat::kn

staticprotected

Definition at line 128 of file RandExpZiggurat.h.

Referenced by ziggurat_init().

float CLHEP::RandExpZiggurat::we

staticprotected

Definition at line 129 of file RandExpZiggurat.h.

Referenced by ziggurat_efix(), ziggurat_init(), and ziggurat_REXP().

float CLHEP::RandExpZiggurat::wn

staticprotected

Definition at line 129 of file RandExpZiggurat.h.

Referenced by ziggurat_init().

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

staticprotected

Definition at line 131 of file RandExpZiggurat.h.

Referenced by RandExpZiggurat(), and ziggurat_init().

---

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

• RandExpZiggurat.h
• RandExpZiggurat.cc