Geant4.10
•All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions | Static Public Member Functions
CLHEP::RandStudentT Class Reference

#include <RandStudentT.h>

Inheritance diagram for CLHEP::RandStudentT:
CLHEP::HepRandom

Public Member Functions

 RandStudentT (HepRandomEngine &anEngine, double a=1.0)
 
 RandStudentT (HepRandomEngine *anEngine, double a=1.0)
 
virtual ~RandStudentT ()
 
std::ostream & put (std::ostream &os) const
 
std::istream & get (std::istream &is)
 
double fire ()
 
double fire (double a)
 
void fireArray (const int size, double *vect)
 
void fireArray (const int size, double *vect, double a)
 
double operator() ()
 
double operator() (double a)
 
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 double shoot ()
 
static double shoot (double a)
 
static void shootArray (const int size, double *vect, double a=1.0)
 
static double shoot (HepRandomEngine *anEngine)
 
static double shoot (HepRandomEngine *anEngine, double a)
 
static void shootArray (HepRandomEngine *anEngine, const int size, double *vect, double a=1.0)
 
static std::string distributionName ()
 
- 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 ()
 

Additional Inherited Members

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

Detailed Description

Author

Definition at line 42 of file RandStudentT.h.

Constructor & Destructor Documentation

CLHEP::RandStudentT::RandStudentT ( HepRandomEngine &  anEngine,
double  a = 1.0 
)
inline
CLHEP::RandStudentT::RandStudentT ( HepRandomEngine *  anEngine,
double  a = 1.0 
)
inline
CLHEP::RandStudentT::~RandStudentT ( )
virtual

Definition at line 30 of file RandStudentT.cc.

31 {
32 }

Member Function Documentation

static std::string CLHEP::RandStudentT::distributionName ( )
inlinestatic

Definition at line 100 of file RandStudentT.h.

100 {return "RandStudentT";}
HepRandomEngine & CLHEP::RandStudentT::engine ( )
virtual

Reimplemented from CLHEP::HepRandom.

Definition at line 28 of file RandStudentT.cc.

28 {return *localEngine;}
double CLHEP::RandStudentT::fire ( )
inline

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

double CLHEP::RandStudentT::fire ( double  a)

Definition at line 98 of file RandStudentT.cc.

References test::v.

98  {
99 
100  double u,v,w;
101 
102  do
103  {
104  u = 2.0 * localEngine->flat() - 1.0;
105  v = 2.0 * localEngine->flat() - 1.0;
106  }
107  while ((w = u * u + v * v) > 1.0);
108 
109  return(u * std::sqrt( a * ( std::exp(- 2.0 / a * std::log(w)) - 1.0) / w));
110 }
void CLHEP::RandStudentT::fireArray ( const int  size,
double *  vect 
)

Definition at line 112 of file RandStudentT.cc.

References fire(), and test::v.

113 {
114  for( double* v = vect; v != vect + size; ++v )
115  *v = fire(defaultA);
116 }
void CLHEP::RandStudentT::fireArray ( const int  size,
double *  vect,
double  a 
)

Definition at line 118 of file RandStudentT.cc.

References fire(), and test::v.

120 {
121  for( double* v = vect; v != vect + size; ++v )
122  *v = fire(a);
123 }
std::istream & CLHEP::RandStudentT::get ( std::istream &  is)
virtual

Reimplemented from CLHEP::HepRandom.

Definition at line 150 of file RandStudentT.cc.

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

150  {
151  std::string inName;
152  is >> inName;
153  if (inName != name()) {
154  is.clear(std::ios::badbit | is.rdstate());
155  std::cerr << "Mismatch when expecting to read state of a "
156  << name() << " distribution\n"
157  << "Name found was " << inName
158  << "\nistream is left in the badbit state\n";
159  return is;
160  }
161  if (possibleKeywordInput(is, "Uvec", defaultA)) {
162  std::vector<unsigned long> t(2);
163  is >> defaultA >> t[0] >> t[1]; defaultA = DoubConv::longs2double(t);
164  return is;
165  }
166  // is >> defaultA encompassed by possibleKeywordInput
167  return is;
168 }
CLHEP::RandStudentT::name
std::string name() const
Definition: RandStudentT.cc:27
CLHEP::possibleKeywordInput
bool possibleKeywordInput(IS &is, const std::string &key, T &t)
Definition: RandomEngine.h:167
CLHEP::DoubConv::longs2double
static double longs2double(const std::vector< unsigned long > &v)
Definition: DoubConv.cc:114
std::string CLHEP::RandStudentT::name ( ) const
virtual

Reimplemented from CLHEP::HepRandom.

Definition at line 27 of file RandStudentT.cc.

Referenced by get(), and put().

27 {return "RandStudentT";}
double CLHEP::RandStudentT::operator() ( )
virtual

Reimplemented from CLHEP::HepRandom.

Definition at line 34 of file RandStudentT.cc.

References fire().

34  {
35  return fire( defaultA );
36 }
double CLHEP::RandStudentT::operator() ( double  a)

Definition at line 38 of file RandStudentT.cc.

References fire().

38  {
39  return fire( a );
40 }
std::ostream & CLHEP::RandStudentT::put ( std::ostream &  os) const
virtual

Reimplemented from CLHEP::HepRandom.

Definition at line 139 of file RandStudentT.cc.

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

139  {
140  int pr=os.precision(20);
141  std::vector<unsigned long> t(2);
142  os << " " << name() << "\n";
143  os << "Uvec" << "\n";
144  t = DoubConv::dto2longs(defaultA);
145  os << defaultA << " " << t[0] << " " << t[1] << "\n";
146  os.precision(pr);
147  return os;
148 }
CLHEP::RandStudentT::name
std::string name() const
Definition: RandStudentT.cc:27
CLHEP::DoubConv::dto2longs
static std::vector< unsigned long > dto2longs(double d)
Definition: DoubConv.cc:98
static double CLHEP::RandStudentT::shoot ( )
inlinestatic

Referenced by shootArray().

double CLHEP::RandStudentT::shoot ( double  a)
static

Definition at line 42 of file RandStudentT.cc.

References DBL_MAX, CLHEP::HepRandomEngine::flat(), CLHEP::HepRandom::getTheEngine(), and test::v.

42  {
43 /******************************************************************
44  * *
45  * Student-t Distribution - Polar Method *
46  * *
47  ******************************************************************
48  * *
49  * The polar method of Box/Muller for generating Normal variates *
50  * is adapted to the Student-t distribution. The two generated *
51  * variates are not independent and the expected no. of uniforms *
52  * per variate is 2.5464. *
53  * *
54  ******************************************************************
55  * *
56  * FUNCTION : - tpol samples a random number from the *
57  * Student-t distribution with parameter a > 0. *
58  * REFERENCE : - R.W. Bailey (1994): Polar generation of random *
59  * variates with the t-distribution, Mathematics *
60  * of Computation 62, 779-781. *
61  * SUBPROGRAM : - ... (0,1)-Uniform generator *
62  * *
63  * *
64  * Implemented by F. Niederl, 1994 *
65  ******************************************************************/
66 
67  double u,v,w;
68 
69 // Check for valid input value
70 
71  if( a < 0.0) return (DBL_MAX);
72 
73  do
74  {
75  u = 2.0 * HepRandom::getTheEngine()->flat() - 1.0;
76  v = 2.0 * HepRandom::getTheEngine()->flat() - 1.0;
77  }
78  while ((w = u * u + v * v) > 1.0);
79 
80  return(u * std::sqrt( a * ( std::exp(- 2.0 / a * std::log(w)) - 1.0) / w));
81 }
CLHEP::HepRandomEngine::flat
virtual double flat()=0
CLHEP::HepRandom::getTheEngine
static HepRandomEngine * getTheEngine()
Definition: Random.cc:165
DBL_MAX
#define DBL_MAX
Definition: templates.hh:83
static double CLHEP::RandStudentT::shoot ( HepRandomEngine *  anEngine)
inlinestatic
double CLHEP::RandStudentT::shoot ( HepRandomEngine *  anEngine,
double  a 
)
static

Definition at line 125 of file RandStudentT.cc.

References CLHEP::HepRandomEngine::flat(), and test::v.

125  {
126 
127  double u,v,w;
128 
129  do
130  {
131  u = 2.0 * anEngine->flat() - 1.0;
132  v = 2.0 * anEngine->flat() - 1.0;
133  }
134  while ((w = u * u + v * v) > 1.0);
135 
136  return(u * std::sqrt( a * ( std::exp(- 2.0 / a * std::log(w)) - 1.0) / w));
137 }
void CLHEP::RandStudentT::shootArray ( const int  size,
double *  vect,
double  a = 1.0 
)
static

Definition at line 83 of file RandStudentT.cc.

References shoot(), and test::v.

85 {
86  for( double* v = vect; v != vect + size; ++v )
87  *v = shoot(a);
88 }
CLHEP::RandStudentT::shoot
static double shoot()
void CLHEP::RandStudentT::shootArray ( HepRandomEngine *  anEngine,
const int  size,
double *  vect,
double  a = 1.0 
)
static

Definition at line 90 of file RandStudentT.cc.

References shoot(), and test::v.

93 {
94  for( double* v = vect; v != vect + size; ++v )
95  *v = shoot(anEngine,a);
96 }
CLHEP::RandStudentT::shoot
static double shoot()
The documentation for this class was generated from the following files:
Generated on Wed Apr 30 2014 15:57:37 for Geant4.10 by   doxygen 1.8.7