CLHEP::RandGauss Class Reference

#include <RandGauss.h>

Inheritance diagram for CLHEP::RandGauss:

Public Member Functions
	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)
virtual double	operator() ()
virtual double	operator() (double mean, double stdDev)
std::string	name () const
HepRandomEngine &	engine ()
std::ostream &	put (std::ostream &os) const
std::istream &	get (std::istream &is)
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 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 ()

Protected Member Functions
double	normal ()

Static Protected Member Functions
static double	getVal ()
static void	setVal (double nextVal)

Protected Attributes
double	defaultMean
double	defaultStdDev
shared_ptr< HepRandomEngine >	localEngine

Additional Inherited Members
Static Protected Attributes inherited from CLHEP::HepRandom
static const long	seedTable [215][2]

Detailed Description

Author

Definition at line 41 of file RandGauss.h.

Constructor & Destructor Documentation

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

inline

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

inline

CLHEP::RandGauss::~RandGauss ( )

virtual

Definition at line 49 of file RandGauss.cc.

                      {
}

Member Function Documentation

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

inlinestatic

Definition at line 99 of file RandGauss.h.

Referenced by restoreDistState(), and saveDistState().

{return "RandGauss";}  

HepRandomEngine & CLHEP::RandGauss::engine ( )

virtual

Reimplemented from CLHEP::HepRandom.

Reimplemented in CLHEP::RandGaussQ, and CLHEP::RandGaussZiggurat.

Definition at line 43 of file RandGauss.cc.

References localEngine.

Referenced by CLHEP::RandGaussZiggurat::engine(), and CLHEP::RandGaussQ::engine().

{return *localEngine;}

double CLHEP::RandGauss::fire

(

)

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

double CLHEP::RandGauss::fire ( double  mean, double  stdDev  )

inline

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

Definition at line 156 of file RandGauss.cc.

References defaultMean, defaultStdDev, fire(), and test::v.

{
  for( double* v = vect; v != vect + size; ++v )
    *v = fire( defaultMean, defaultStdDev );
}

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

Definition at line 162 of file RandGauss.cc.

References fire(), and test::v.

{
  for( double* v = vect; v != vect + size; ++v )
    *v = fire( mean, stdDev );
}

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

virtual

Reimplemented from CLHEP::HepRandom.

Reimplemented in CLHEP::RandGaussQ, and CLHEP::RandGaussZiggurat.

Definition at line 255 of file RandGauss.cc.

References plottest35::c1, defaultMean, defaultStdDev, CLHEP::DoubConv::longs2double(), name(), and CLHEP::possibleKeywordInput().

Referenced by CLHEP::RandGaussZiggurat::get(), and CLHEP::RandGaussQ::get().

                                              {
  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;
  }
  std::string c1;
  std::string c2;
  if (possibleKeywordInput(is, "Uvec", c1)) {
    std::vector<unsigned long> t(2);
    is >> defaultMean >> t[0] >> t[1]; defaultMean = DoubConv::longs2double(t); 
    is >> defaultStdDev>>t[0]>>t[1]; defaultStdDev = DoubConv::longs2double(t); 
    std::string ng;
    is >> ng;
    set = false;
    if (ng == "nextGauss") {
      is >> nextGauss >> t[0] >> t[1]; nextGauss = DoubConv::longs2double(t);
      set = true;
    }
    return is;
  }
  // is >> c1 encompassed by possibleKeywordInput
  is >> defaultMean >> c2 >> defaultStdDev;
  if ( (!is) || (c1 != "Mean:") || (c2 != "Sigma:") ) {
    std::cerr << "i/o problem while expecting to read state of a "
              << name() << " distribution\n"
          << "default mean and/or sigma could not be read\n";
    return is;
  }
  is >> c1 >> c2 >> nextGauss;
  if ( (!is) || (c1 != "RANDGAUSS") ) {
    is.clear(std::ios::badbit | is.rdstate());
    std::cerr << "Failure when reading caching state of RandGauss\n";
    return is;
  }
  if (c2 == "CACHED_GAUSSIAN:") {
    set = true;
  } else if (c2 == "NO_CACHED_GAUSSIAN:") {
    set = false;  
  } else {
    is.clear(std::ios::badbit | is.rdstate());
    std::cerr << "Unexpected caching state keyword of RandGauss:" << c2
          << "\nistream is left in the badbit state\n";
  } 
  return is;
} // get

bool CLHEP::RandGauss::getF ( ) const

inline

Definition at line 115 of file RandGauss.h.

{return set;}

static bool CLHEP::RandGauss::getFlag ( )

inlinestatic

Definition at line 111 of file RandGauss.h.

Referenced by saveDistState(), saveEngineStatus(), and shoot().

{return set_st;}

static double CLHEP::RandGauss::getVal ( )

inlinestaticprotected

Definition at line 145 of file RandGauss.h.

Referenced by saveDistState(), saveEngineStatus(), and shoot().

{return nextGauss_st;}

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

virtual

Reimplemented from CLHEP::HepRandom.

Reimplemented in CLHEP::RandGaussQ, and CLHEP::RandGaussZiggurat.

Definition at line 42 of file RandGauss.cc.

Referenced by get(), and put().

{return "RandGauss";}

double CLHEP::RandGauss::normal ( )

protected

Definition at line 130 of file RandGauss.cc.

References localEngine.

{
  // Gaussian random numbers are generated two at the time, so every other
  // time this is called we just return a number generated the time before.

  if ( set ) {
    set = false;
    return nextGauss;
  }

  double r;
  double v1,v2,fac,val;

  do {
    v1 = 2.0 * localEngine->flat() - 1.0;
    v2 = 2.0 * localEngine->flat() - 1.0;
    r = v1*v1 + v2*v2;
  } while ( r > 1.0 );

  fac = std::sqrt(-2.0*std::log(r)/r);
  val = v1*fac;
  nextGauss = val;
  set = true;
  return v2*fac;
}

double CLHEP::RandGauss::operator() ( )

virtual

Reimplemented from CLHEP::HepRandom.

Reimplemented in CLHEP::RandGaussQ, and CLHEP::RandGaussZiggurat.

Definition at line 52 of file RandGauss.cc.

References defaultMean, defaultStdDev, and fire().

                             {
  return fire( defaultMean, defaultStdDev );
}

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

virtual

Reimplemented in CLHEP::RandGaussQ, and CLHEP::RandGaussZiggurat.

Definition at line 56 of file RandGauss.cc.

References fire().

                                                         {
  return fire( mean, stdDev );
}

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

virtual

Reimplemented from CLHEP::HepRandom.

Reimplemented in CLHEP::RandGaussQ, and CLHEP::RandGaussZiggurat.

Definition at line 236 of file RandGauss.cc.

References defaultMean, defaultStdDev, CLHEP::DoubConv::dto2longs(), and name().

Referenced by CLHEP::RandGaussZiggurat::put(), and CLHEP::RandGaussQ::put().

                                                    {
  os << name() << "\n";
  int prec = os.precision(20);
  std::vector<unsigned long> t(2);
  os << "Uvec\n";
  t = DoubConv::dto2longs(defaultMean);
  os << defaultMean << " " << t[0] << " " << t[1] << "\n";
  t = DoubConv::dto2longs(defaultStdDev);
  os << defaultStdDev << " " << t[0] << " " << t[1] << "\n";
  if ( set ) {
    t = DoubConv::dto2longs(nextGauss);
    os << "nextGauss " << nextGauss << " " << t[0] << " " << t[1] << "\n";
  } else {
    os << "no_cached_nextGauss \n";
  }
  os.precision(prec);
  return os;
} // put   

std::istream & CLHEP::RandGauss::restoreDistState ( std::istream &  is )

static

Definition at line 324 of file RandGauss.cc.

References plottest35::c1, distributionName(), CLHEP::DoubConv::longs2double(), CLHEP::possibleKeywordInput(), and setFlag().

Referenced by CLHEP::StaticRandomStates::restore(), and restoreFullState().

                                                           {
  std::string inName;
  is >> inName;
  if (inName != distributionName()) {
    is.clear(std::ios::badbit | is.rdstate());
    std::cerr << "Mismatch when expecting to read static state of a "
              << distributionName() << " distribution\n"
          << "Name found was " << inName
          << "\nistream is left in the badbit state\n";
    return is;
  }
  std::string c1;
  std::string c2;
  if (possibleKeywordInput(is, "Uvec", c1)) {
    std::vector<unsigned long> t(2);
    std::string ng;
    is >> ng;
    setFlag (false);
    if (ng == "nextGauss_st") {
      is >> nextGauss_st >> t[0] >> t[1]; 
      nextGauss_st = DoubConv::longs2double(t);
      setFlag (true);
    }
    return is;
  }
  // is >> c1 encompassed by possibleKeywordInput
  is >> c2 >> nextGauss_st;
  if ( (!is) || (c1 != "RANDGAUSS") ) {
    is.clear(std::ios::badbit | is.rdstate());
    std::cerr << "Failure when reading caching state of static RandGauss\n";
    return is;
  }
  if (c2 == "CACHED_GAUSSIAN:") {
    setFlag(true);
  } else if (c2 == "NO_CACHED_GAUSSIAN:") {
    setFlag(false);  
  } else {
    is.clear(std::ios::badbit | is.rdstate());
    std::cerr << "Unexpected caching state keyword of static RandGauss:" << c2
          << "\nistream is left in the badbit state\n";
  } 
  return is;
} 

void CLHEP::RandGauss::restoreEngineStatus ( const char  filename[] = "Config.conf" )

static

Definition at line 189 of file RandGauss.cc.

References CLHEP::HepRandom::getTheEngine(), CLHEP::DoubConv::longs2double(), CLHEP::possibleKeywordInput(), CLHEP::HepRandomEngine::restoreStatus(), and setFlag().

                                                           {

  // First restore the engine status just like the base class would do:
  getTheEngine()->restoreStatus( filename );

  // Now find the line describing the cached variate:

  std::ifstream infile ( filename, std::ios::in );
  if (!infile) return;

  char inputword[] = "NO_KEYWORD    "; // leaves room for 14 characters plus \0
  while (true) {
    infile.width(13);
    infile >> inputword;
    if (strcmp(inputword,"RANDGAUSS")==0) break;
    if (infile.eof()) break;
    // If the file ends without the RANDGAUSS line, that means this
    // was a file produced by an earlier version of RandGauss.  We will
    // replicated the old behavior in that case:  set_st is cleared.
  }

  // Then read and use the caching info:

  if (strcmp(inputword,"RANDGAUSS")==0) {
    char setword[40];   // the longest, staticFirstUnusedBit: has length 21
    infile.width(39);
    infile >> setword;  // setword should be CACHED_GAUSSIAN:
    if (strcmp(setword,"CACHED_GAUSSIAN:") ==0) {
      if (possibleKeywordInput(infile, "Uvec", nextGauss_st)) {
        std::vector<unsigned long> t(2);
        infile >> nextGauss_st >> t[0] >> t[1]; 
        nextGauss_st = DoubConv::longs2double(t); 
      }
      // is >> nextGauss_st encompassed by possibleKeywordInput
      setFlag(true);
    } else {
      setFlag(false);
      infile >> nextGauss_st; // because a 0 will have been output
    }
  } else {
    setFlag(false);
  }

} // restoreEngineStatus

std::istream & CLHEP::RandGauss::restoreFullState ( std::istream &  is )

static

Definition at line 374 of file RandGauss.cc.

References restoreDistState(), and CLHEP::HepRandom::restoreFullState().

                                                           {
  HepRandom::restoreFullState(is);
  restoreDistState(is);
  return is;
}

std::ostream & CLHEP::RandGauss::saveDistState ( std::ostream &  os )

static

Definition at line 309 of file RandGauss.cc.

References distributionName(), CLHEP::DoubConv::dto2longs(), getFlag(), and getVal().

Referenced by saveFullState().

                                                        {
  int prec = os.precision(20);
  std::vector<unsigned long> t(2);
  os << distributionName() << "\n";
  os << "Uvec\n";
  if ( getFlag() ) {
    t = DoubConv::dto2longs(getVal());
    os << "nextGauss_st " << getVal() << " " << t[0] << " " << t[1] << "\n";
  } else {
    os << "no_cached_nextGauss_st \n";
  }
  os.precision(prec);
  return os;
}    

void CLHEP::RandGauss::saveEngineStatus ( const char  filename[] = "Config.conf" )

static

Definition at line 169 of file RandGauss.cc.

References testem0::app, CLHEP::DoubConv::dto2longs(), getFlag(), CLHEP::HepRandom::getTheEngine(), getVal(), and CLHEP::HepRandomEngine::saveStatus().

                                                         {

  // First save the engine status just like the base class would do:
  getTheEngine()->saveStatus( filename );

  // Now append the cached variate, if any:

  std::ofstream outfile ( filename, std::ios::app );

  if ( getFlag() ) {
    std::vector<unsigned long> t(2);
    t = DoubConv::dto2longs(getVal());
    outfile << "RANDGAUSS CACHED_GAUSSIAN: Uvec " 
        << getVal() << " " << t[0] << " " << t[1] << "\n";
  } else {
    outfile << "RANDGAUSS NO_CACHED_GAUSSIAN: 0 \n" ;
  }

} // saveEngineStatus

std::ostream & CLHEP::RandGauss::saveFullState ( std::ostream &  os )

static

Definition at line 368 of file RandGauss.cc.

References saveDistState(), and CLHEP::HepRandom::saveFullState().

Referenced by CLHEP::StaticRandomStates::save().

                                                        {
  HepRandom::saveFullState(os);
  saveDistState(os);
  return os;
}

void CLHEP::RandGauss::setF ( bool  val )

inline

Definition at line 117 of file RandGauss.h.

{set = val;}

static void CLHEP::RandGauss::setFlag ( bool  val )

inlinestatic

Definition at line 113 of file RandGauss.h.

Referenced by restoreDistState(), restoreEngineStatus(), and shoot().

{set_st = val;}

static void CLHEP::RandGauss::setVal ( double  nextVal )

inlinestaticprotected

Definition at line 147 of file RandGauss.h.

Referenced by shoot().

{nextGauss_st = nextVal;}

double CLHEP::RandGauss::shoot ( )

static

Definition at line 60 of file RandGauss.cc.

References CLHEP::HepRandomEngine::flat(), getFlag(), CLHEP::HepRandom::getTheEngine(), getVal(), setFlag(), setVal(), and test::x.

Referenced by shootArray().

{
  // Gaussian random numbers are generated two at the time, so every other
  // time this is called we just return a number generated the time before.

  if ( getFlag() ) {
    setFlag(false);
    double x = getVal();
    return x; 
    // return getVal();
  } 

  double r;
  double v1,v2,fac,val;
  HepRandomEngine* anEngine = HepRandom::getTheEngine();

  do {
    v1 = 2.0 * anEngine->flat() - 1.0;
    v2 = 2.0 * anEngine->flat() - 1.0;
    r = v1*v1 + v2*v2;
  } while ( r > 1.0 );

  fac = std::sqrt(-2.0*std::log(r)/r);
  val = v1*fac;
  setVal(val);
  setFlag(true);
  return v2*fac;
}

static double CLHEP::RandGauss::shoot ( double  mean, double  stdDev  )

inlinestatic

double CLHEP::RandGauss::shoot ( HepRandomEngine *  anEngine )

static

Definition at line 96 of file RandGauss.cc.

References CLHEP::HepRandomEngine::flat(), getFlag(), getVal(), setFlag(), and setVal().

{
  // Gaussian random numbers are generated two at the time, so every other
  // time this is called we just return a number generated the time before.

  if ( getFlag() ) {
    setFlag(false);
    return getVal();
  }

  double r;
  double v1,v2,fac,val;

  do {
    v1 = 2.0 * anEngine->flat() - 1.0;
    v2 = 2.0 * anEngine->flat() - 1.0;
    r = v1*v1 + v2*v2;
  } while ( r > 1.0 );

  fac = std::sqrt( -2.0*std::log(r)/r);
  val = v1*fac;
  setVal(val);
  setFlag(true);
  return v2*fac;
}

static double CLHEP::RandGauss::shoot ( HepRandomEngine *  anEngine, double  mean, double  stdDev  )

inlinestatic

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

static

Definition at line 89 of file RandGauss.cc.

References shoot(), and test::v.

{
  for( double* v = vect; v != vect + size; ++v )
    *v = shoot(mean,stdDev);
}

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

static

Definition at line 122 of file RandGauss.cc.

References shoot(), and test::v.

{
  for( double* v = vect; v != vect + size; ++v )
    *v = shoot(anEngine,mean,stdDev);
}

Field Documentation

double CLHEP::RandGauss::defaultMean

protected

Definition at line 151 of file RandGauss.h.

Referenced by CLHEP::RandGaussZiggurat::fire(), CLHEP::RandGaussZiggurat::fireArray(), CLHEP::RandGaussQ::fireArray(), fireArray(), get(), CLHEP::RandGaussZiggurat::operator()(), CLHEP::RandGaussQ::operator()(), operator()(), and put().

double CLHEP::RandGauss::defaultStdDev

protected

Definition at line 152 of file RandGauss.h.

Referenced by CLHEP::RandGaussZiggurat::fire(), CLHEP::RandGaussZiggurat::fireArray(), CLHEP::RandGaussQ::fireArray(), fireArray(), get(), CLHEP::RandGaussZiggurat::operator()(), CLHEP::RandGaussQ::operator()(), operator()(), and put().

shared_ptr<HepRandomEngine> CLHEP::RandGauss::localEngine

protected

Definition at line 154 of file RandGauss.h.

Referenced by engine(), CLHEP::RandGaussZiggurat::fire(), normal(), CLHEP::RandGaussZiggurat::operator()(), and CLHEP::RandGaussQ::operator()().

---

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

• RandGauss.h
• RandGauss.cc