CLHEP::RandGeneral Class Reference

#include <RandGeneral.h>

Inheritance diagram for CLHEP::RandGeneral:

Public Member Functions
HepRandomEngine &	engine ()
double	fire ()
void	fireArray (const int size, double *vect)
double	flat ()
double	flat (HepRandomEngine *theNewEngine)
void	flatArray (const int size, double *vect)
void	flatArray (HepRandomEngine *theNewEngine, const int size, double *vect)
std::istream &	get (std::istream &is)
std::string	name () const
double	operator() ()
std::ostream &	put (std::ostream &os) const
	RandGeneral (const double *aProbFunc, int theProbSize, int IntType=0)
	RandGeneral (HepRandomEngine &anEngine, const double *aProbFunc, int theProbSize, int IntType=0)
	RandGeneral (HepRandomEngine *anEngine, const double *aProbFunc, int theProbSize, int IntType=0)
double	shoot ()
double	shoot (HepRandomEngine *anEngine)
void	shootArray (const int size, double *vect)
void	shootArray (HepRandomEngine *anEngine, const int size, double *vect)
virtual	~RandGeneral ()

Static Public Member Functions
static int	createInstance ()
static std::string	distributionName ()
static HepRandomEngine *	getTheEngine ()
static HepRandom *	getTheGenerator ()
static long	getTheSeed ()
static const long *	getTheSeeds ()
static void	getTheTableSeeds (long *seeds, int index)
static std::istream &	restoreDistState (std::istream &is)
static void	restoreEngineStatus (const char filename[]="Config.conf")
static std::istream &	restoreFullState (std::istream &is)
static std::istream &	restoreStaticRandomStates (std::istream &is)
static std::ostream &	saveDistState (std::ostream &os)
static void	saveEngineStatus (const char filename[]="Config.conf")
static std::ostream &	saveFullState (std::ostream &os)
static std::ostream &	saveStaticRandomStates (std::ostream &os)
static void	setTheEngine (HepRandomEngine *theNewEngine)
static void	setTheSeed (long seed, int lxr=3)
static void	setTheSeeds (const long *seeds, int aux=-1)
static void	showEngineStatus ()

Static Protected Attributes
static const long	seedTable [215][2]

Private Member Functions
double	mapRandom (double rand) const
void	prepareTable (const double *aProbFunc)
void	useFlatDistribution ()

Private Attributes
int	InterpolationType
std::shared_ptr< HepRandomEngine >	localEngine
int	nBins
double	oneOverNbins
std::vector< double >	theIntegralPdf

Detailed Description

Author

Definition at line 38 of file RandGeneral.h.

Constructor & Destructor Documentation

RandGeneral() [1/3]

CLHEP::RandGeneral::RandGeneral	(	const double *	aProbFunc,
		int	theProbSize,
		int	IntType = 0
	)

Definition at line 66 of file RandGeneral.cc.

  : HepRandom(),
    localEngine(HepRandom::getTheEngine(), do_nothing_deleter()),
    nBins(theProbSize), 
    InterpolationType(IntType)
{
  prepareTable(aProbFunc);
}

References prepareTable().

RandGeneral() [2/3]

CLHEP::RandGeneral::RandGeneral	(	HepRandomEngine &	anEngine,
		const double *	aProbFunc,
		int	theProbSize,
		int	IntType = 0
	)

Definition at line 77 of file RandGeneral.cc.

: HepRandom(),
  localEngine(&anEngine, do_nothing_deleter()), 
  nBins(theProbSize),
  InterpolationType(IntType)
{
  prepareTable(aProbFunc);
}

References prepareTable().

RandGeneral() [3/3]

CLHEP::RandGeneral::RandGeneral	(	HepRandomEngine *	anEngine,
		const double *	aProbFunc,
		int	theProbSize,
		int	IntType = 0
	)

Definition at line 89 of file RandGeneral.cc.

: HepRandom(),
  localEngine(anEngine), 
  nBins(theProbSize),
  InterpolationType(IntType)
{
  prepareTable(aProbFunc);
}

References prepareTable().

~RandGeneral()

CLHEP::RandGeneral::~RandGeneral ( )

virtual

Definition at line 174 of file RandGeneral.cc.

                          {
}

Member Function Documentation

createInstance()

int CLHEP::HepRandom::createInstance ( )

staticinherited

Definition at line 311 of file Random.cc.

{
  return static_cast<int>( theDefaults().ensureInitialized() );
}

References CLHEP::anonymous_namespace{Random.cc}::defaults::ensureInitialized(), and CLHEP::anonymous_namespace{Random.cc}::theDefaults().

distributionName()

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

inlinestatic

Definition at line 124 of file RandGeneral.h.

{return "RandGeneral";}  

engine()

HepRandomEngine & CLHEP::RandGeneral::engine ( )

virtual

Reimplemented from CLHEP::HepRandom.

Definition at line 59 of file RandGeneral.cc.

{return *localEngine;}

References localEngine.

fire()

double CLHEP::RandGeneral::fire

(

)

Referenced by fireArray().

fireArray()

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

Definition at line 240 of file RandGeneral.cc.

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

References fire().

flat() [1/2]

double CLHEP::HepRandom::flat ( )

inherited

Definition at line 199 of file Random.cc.

{
  return theDefaults().theEngine->flat();
}

References CLHEP::anonymous_namespace{Random.cc}::theDefaults(), and CLHEP::anonymous_namespace{Random.cc}::defaults::theEngine.

Referenced by CLHEP::RandExponential::shoot().

flat() [2/2]

double CLHEP::HepRandom::flat ( HepRandomEngine *  theNewEngine )

inlineinherited

flatArray() [1/2]

void CLHEP::HepRandom::flatArray ( const int  size, double *  vect  )

inherited

Definition at line 204 of file Random.cc.

{
  theDefaults().theEngine->flatArray(size,vect);
}

References CLHEP::anonymous_namespace{Random.cc}::theDefaults(), and CLHEP::anonymous_namespace{Random.cc}::defaults::theEngine.

flatArray() [2/2]

void CLHEP::HepRandom::flatArray ( HepRandomEngine *  theNewEngine, const int  size, double *  vect  )

inlineinherited

get()

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

virtual

Reimplemented from CLHEP::HepRandom.

Definition at line 266 of file RandGeneral.cc.

                                                {
  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", nBins)) {
    std::vector<unsigned long> t(2);
    is >> nBins >> oneOverNbins >> InterpolationType;
    is >> t[0] >> t[1]; oneOverNbins = DoubConv::longs2double(t); 
    theIntegralPdf.resize(nBins+1);
    for (unsigned int i=0; i<theIntegralPdf.size(); ++i) {
      is >> theIntegralPdf[i] >> t[0] >> t[1];
      theIntegralPdf[i] = DoubConv::longs2double(t); 
    }
    return is;
  }
  // is >> nBins encompassed by possibleKeywordInput
  is >> oneOverNbins >> InterpolationType;
  theIntegralPdf.resize(nBins+1);
  for (unsigned int i=0; i<theIntegralPdf.size(); ++i) is >> theIntegralPdf[i];
  return is;
}

References InterpolationType, CLHEP::DoubConv::longs2double(), name(), nBins, oneOverNbins, CLHEP::possibleKeywordInput(), and theIntegralPdf.

getTheEngine()

HepRandomEngine * CLHEP::HepRandom::getTheEngine ( )

staticinherited

Definition at line 268 of file Random.cc.

{
  return theDefaults().theEngine.get();
}

References CLHEP::anonymous_namespace{Random.cc}::theDefaults(), and CLHEP::anonymous_namespace{Random.cc}::defaults::theEngine.

Referenced by CLHEP::StaticRandomStates::restore(), CLHEP::RandFlat::restoreEngineStatus(), CLHEP::RandGauss::restoreEngineStatus(), CLHEP::RandFlat::saveEngineStatus(), CLHEP::RandGauss::saveEngineStatus(), CLHEP::RandExponential::shoot(), CLHEP::RandExpZiggurat::shoot(), CLHEP::RandFlat::shoot(), CLHEP::RandGauss::shoot(), CLHEP::RandGaussZiggurat::shoot(), CLHEP::RandChiSquare::shoot(), CLHEP::RandStudentT::shoot(), CLHEP::RandBreitWigner::shoot(), CLHEP::RandGamma::shoot(), CLHEP::RandPoisson::shoot(), CLHEP::RandPoissonQ::shoot(), CLHEP::RandBinomial::shoot(), and CLHEP::RandFlat::shootArray().

getTheGenerator()

HepRandom * CLHEP::HepRandom::getTheGenerator ( )

staticinherited

Definition at line 263 of file Random.cc.

{
  return theDefaults().theGenerator.get();
}

References CLHEP::anonymous_namespace{Random.cc}::theDefaults(), and CLHEP::anonymous_namespace{Random.cc}::defaults::theGenerator.

getTheSeed()

long CLHEP::HepRandom::getTheSeed ( )

staticinherited

Definition at line 239 of file Random.cc.

{
  return theDefaults().theEngine->getSeed();
}

References CLHEP::anonymous_namespace{Random.cc}::theDefaults(), and CLHEP::anonymous_namespace{Random.cc}::defaults::theEngine.

getTheSeeds()

const long * CLHEP::HepRandom::getTheSeeds ( )

staticinherited

Definition at line 249 of file Random.cc.

{
  return theDefaults().theEngine->getSeeds();
}

References CLHEP::anonymous_namespace{Random.cc}::theDefaults(), and CLHEP::anonymous_namespace{Random.cc}::defaults::theEngine.

getTheTableSeeds()

void CLHEP::HepRandom::getTheTableSeeds ( long *  seeds, int  index  )

staticinherited

Definition at line 254 of file Random.cc.

{
  if ((index >= 0) && (index < 215)) {
    seeds[0] = seedTable[index][0];
    seeds[1] = seedTable[index][1];
  }
  else seeds = NULL;
}

Referenced by CLHEP::HepJamesRandom::HepJamesRandom(), CLHEP::MTwistEngine::MTwistEngine(), CLHEP::RanecuEngine::RanecuEngine(), CLHEP::Ranlux64Engine::Ranlux64Engine(), CLHEP::RanluxEngine::RanluxEngine(), and CLHEP::RanecuEngine::setSeed().

mapRandom()

double CLHEP::RandGeneral::mapRandom ( double  rand ) const

private

Definition at line 182 of file RandGeneral.cc.

                                               {
//
// Private method to take the random (however it is created) and map it
// according to the distribution.
//
 
  int nbelow = 0;     // largest k such that I[k] is known to be <= rand
  int nabove = nBins;     // largest k such that I[k] is known to be >  rand
  int middle;
  
  while (nabove > nbelow+1) {
    middle = (nabove + nbelow+1)>>1;
    if (rand >= theIntegralPdf[middle]) {
      nbelow = middle;
    } else {
      nabove = middle;
    }
  } // after this loop, nabove is always nbelow+1 and they straddle rad:
    assert ( nabove == nbelow+1 );
    assert ( theIntegralPdf[nbelow] <= rand );
    assert ( theIntegralPdf[nabove] >= rand );  
        // If a defective engine produces rand=1, that will 
        // still give sensible results so we relax the > rand assertion
 
  if ( InterpolationType == 1 ) {
 
    return nbelow * oneOverNbins;
 
  } else {
 
    double binMeasure = theIntegralPdf[nabove] - theIntegralPdf[nbelow];
    // binMeasure is always aProbFunc[nbelow], 
    // but we don't have aProbFunc any more so we subtract.
 
    if ( binMeasure == 0 ) { 
    // rand lies right in a bin of measure 0.  Simply return the center
    // of the range of that bin.  (Any value between k/N and (k+1)/N is 
    // equally good, in this rare case.)
        return (nbelow + .5) * oneOverNbins;
    }
 
    double binFraction = (rand - theIntegralPdf[nbelow]) / binMeasure;
 
    return (nbelow + binFraction) * oneOverNbins;
  }
 
} // mapRandom(rand)

References InterpolationType, nBins, oneOverNbins, and theIntegralPdf.

name()

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

virtual

Reimplemented from CLHEP::HepRandom.

Definition at line 58 of file RandGeneral.cc.

{return "RandGeneral";}

Referenced by source.g4viscp.G4Scene::create_scene(), get(), mcscore.MCParticle::printout(), put(), and source.g4viscp.G4Scene::update_scene().

operator()()

double CLHEP::RandGeneral::operator() ( )

virtual

Reimplemented from CLHEP::HepRandom.

prepareTable()

void CLHEP::RandGeneral::prepareTable ( const double *  aProbFunc )

private

Definition at line 101 of file RandGeneral.cc.

                                                      {
//
// Private method called only by constructors.  Prepares theIntegralPdf.
//
  if (nBins < 1) {
    std::cerr << 
    "RandGeneral constructed with no bins - will use flat distribution\n";
    useFlatDistribution();
    return;
  }
 
  theIntegralPdf.resize(nBins+1);
  theIntegralPdf[0] = 0;
  int ptn;
  double weight;
 
  for ( ptn = 0; ptn<nBins; ++ptn ) {
    weight = aProbFunc[ptn];
    if ( weight < 0 ) {
    // We can't stomach negative bin contents, they invalidate the 
    // search algorithm when the distribution is fired.
      std::cerr << 
    "RandGeneral constructed with negative-weight bin " << ptn <<
    " = " << weight << " \n   -- will substitute 0 weight \n";
      weight = 0;
    }
    // std::cout << ptn << "  " << weight << "  " << theIntegralPdf[ptn] << "\n";
    theIntegralPdf[ptn+1] = theIntegralPdf[ptn] + weight;
  } 
 
  if ( theIntegralPdf[nBins] <= 0 ) {
    std::cerr << 
      "RandGeneral constructed nothing in bins - will use flat distribution\n";
    useFlatDistribution();
    return;
  }
 
  for ( ptn = 0; ptn < nBins+1; ++ptn ) {
    theIntegralPdf[ptn] /= theIntegralPdf[nBins];
    // std::cout << ptn << "  " << theIntegralPdf[ptn] << "\n";
  }
 
  // And another useful variable is ...
  oneOverNbins = 1.0 / nBins;
 
  // One last chore:
 
  if ( (InterpolationType != 0) && (InterpolationType != 1) ) {
    std::cerr << 
      "RandGeneral does not recognize IntType " << InterpolationType 
      << "\n Will use type 0 (continuous linear interpolation \n";
    InterpolationType = 0;
  }
 
} // prepareTable()

References InterpolationType, nBins, oneOverNbins, theIntegralPdf, and useFlatDistribution().

Referenced by RandGeneral().

put()

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

virtual

Reimplemented from CLHEP::HepRandom.

Definition at line 249 of file RandGeneral.cc.

                                                      {
  int pr=os.precision(20);
  std::vector<unsigned long> t(2);
  os << " " << name() << "\n";
  os << "Uvec" << "\n";
  os << nBins << " " << oneOverNbins << " " << InterpolationType << "\n";    
  t = DoubConv::dto2longs(oneOverNbins);
  os << t[0] << " " << t[1] << "\n";
  assert (static_cast<int>(theIntegralPdf.size())==nBins+1);
  for (unsigned int i=0; i<theIntegralPdf.size(); ++i) {
    t = DoubConv::dto2longs(theIntegralPdf[i]);
    os << theIntegralPdf[i] << " " << t[0] << " " << t[1] << "\n";
  }
  os.precision(pr);
  return os;
}

References CLHEP::DoubConv::dto2longs(), InterpolationType, name(), nBins, oneOverNbins, and theIntegralPdf.

restoreDistState()

static std::istream & CLHEP::HepRandom::restoreDistState ( std::istream &  is )

inlinestaticinherited

Definition at line 136 of file Random.h.

{return is;}

restoreEngineStatus()

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

staticinherited

Definition at line 283 of file Random.cc.

{
  theDefaults().theEngine->restoreStatus( filename );
}

References CLHEP::anonymous_namespace{Random.cc}::theDefaults(), and CLHEP::anonymous_namespace{Random.cc}::defaults::theEngine.

restoreFullState()

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

staticinherited

Definition at line 293 of file Random.cc.

                                                          {
  is >> *getTheEngine();
  return is;
}

Referenced by CLHEP::RandFlat::restoreFullState(), and CLHEP::RandGauss::restoreFullState().

restoreStaticRandomStates()

std::istream & CLHEP::HepRandom::restoreStaticRandomStates ( std::istream &  is )

staticinherited

Definition at line 302 of file Random.cc.

                                                                   {
  return StaticRandomStates::restore(is);
}

References CLHEP::StaticRandomStates::restore().

saveDistState()

static std::ostream & CLHEP::HepRandom::saveDistState ( std::ostream &  os )

inlinestaticinherited

Definition at line 133 of file Random.h.

{return os;}

saveEngineStatus()

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

staticinherited

Definition at line 278 of file Random.cc.

{
  theDefaults().theEngine->saveStatus( filename );
}

References CLHEP::anonymous_namespace{Random.cc}::theDefaults(), and CLHEP::anonymous_namespace{Random.cc}::defaults::theEngine.

Referenced by G4CascadeInterface::ApplyYourself(), G4HadronicProcess::PostStepDoIt(), and G4CascadeInterface::Propagate().

saveFullState()

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

staticinherited

Definition at line 288 of file Random.cc.

                                                       {
  os << *getTheEngine();
  return os;
}

Referenced by G4EventManager::DoProcessing(), G4EventManager::ProcessOneEvent(), CLHEP::RandFlat::saveFullState(), and CLHEP::RandGauss::saveFullState().

saveStaticRandomStates()

std::ostream & CLHEP::HepRandom::saveStaticRandomStates ( std::ostream &  os )

staticinherited

Definition at line 298 of file Random.cc.

                                                                {
  return StaticRandomStates::save(os);
}

References CLHEP::StaticRandomStates::save().

setTheEngine()

void CLHEP::HepRandom::setTheEngine ( HepRandomEngine *  theNewEngine )

staticinherited

Definition at line 273 of file Random.cc.

{
  theDefaults().theEngine.reset( theNewEngine, do_nothing_deleter() );
}

References CLHEP::anonymous_namespace{Random.cc}::theDefaults(), and CLHEP::anonymous_namespace{Random.cc}::defaults::theEngine.

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

setTheSeed()

void CLHEP::HepRandom::setTheSeed ( long  seed, int  lxr = 3  )

staticinherited

Definition at line 234 of file Random.cc.

{
  theDefaults().theEngine->setSeed(seed,lux);
}

References CLHEP::lux, CLHEP::anonymous_namespace{Random.cc}::theDefaults(), and CLHEP::anonymous_namespace{Random.cc}::defaults::theEngine.

setTheSeeds()

void CLHEP::HepRandom::setTheSeeds ( const long *  seeds, int  aux = -1  )

staticinherited

Definition at line 244 of file Random.cc.

{
  theDefaults().theEngine->setSeeds(seeds,aux);
}

References CLHEP::anonymous_namespace{Random.cc}::theDefaults(), and CLHEP::anonymous_namespace{Random.cc}::defaults::theEngine.

shoot() [1/2]

double CLHEP::RandGeneral::shoot ( )

inline

Referenced by G4JAEAPolarizedElasticScatteringModel::SampleSecondaries(), and shootArray().

shoot() [2/2]

double CLHEP::RandGeneral::shoot

(

HepRandomEngine *

anEngine

)

shootArray() [1/2]

void CLHEP::RandGeneral::shootArray ( const int  size, double *  vect  )

inline

shootArray() [2/2]

void CLHEP::RandGeneral::shootArray	(	HepRandomEngine *	anEngine,
		const int	size,
		double *	vect
	)

Definition at line 230 of file RandGeneral.cc.

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

References shoot().

showEngineStatus()

void CLHEP::HepRandom::showEngineStatus ( )

staticinherited

Definition at line 306 of file Random.cc.

{
  theDefaults().theEngine->showStatus();
}

References CLHEP::anonymous_namespace{Random.cc}::theDefaults(), and CLHEP::anonymous_namespace{Random.cc}::defaults::theEngine.

useFlatDistribution()

void CLHEP::RandGeneral::useFlatDistribution ( )

private

Definition at line 157 of file RandGeneral.cc.

                                      {
//
// Private method called only by prepareTables in case of user error. 
//
    nBins = 1;
    theIntegralPdf.resize(2);
    theIntegralPdf[0] = 0;
    theIntegralPdf[1] = 1;
    oneOverNbins = 1.0;
    return;
 
} // UseFlatDistribution()

References nBins, oneOverNbins, and theIntegralPdf.

Referenced by prepareTable().

Field Documentation

InterpolationType

int CLHEP::RandGeneral::InterpolationType

private

Definition at line 134 of file RandGeneral.h.

Referenced by get(), mapRandom(), prepareTable(), and put().

localEngine

std::shared_ptr<HepRandomEngine> CLHEP::RandGeneral::localEngine

private

Definition at line 130 of file RandGeneral.h.

Referenced by engine().

nBins

int CLHEP::RandGeneral::nBins

private

Definition at line 132 of file RandGeneral.h.

Referenced by get(), mapRandom(), prepareTable(), put(), and useFlatDistribution().

oneOverNbins

double CLHEP::RandGeneral::oneOverNbins

private

Definition at line 133 of file RandGeneral.h.

Referenced by get(), mapRandom(), prepareTable(), put(), and useFlatDistribution().

seedTable

const long CLHEP::HepRandom::seedTable

staticprotectedinherited

Definition at line 156 of file Random.h.

theIntegralPdf

std::vector<double> CLHEP::RandGeneral::theIntegralPdf

private

Definition at line 131 of file RandGeneral.h.

Referenced by get(), mapRandom(), prepareTable(), put(), and useFlatDistribution().

---

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

• source/externals/clhep/include/CLHEP/Random/RandGeneral.h
• source/externals/clhep/src/RandGeneral.cc