ROOT 6.18/05 Reference Guide |
This is the base class for the ROOT Random number generators.
This class defines the ROOT Random number interface and it should not be instantiated directly but used via its derived classes. The generator provided in TRandom itself is a LCG (Linear Congruential Generator), the BSD rand
generator, that it should not be used because its period is only 2**31, i.e. approximatly 2 billion events, that can be generated in just few seconds.
To generate random numbers, one should use the derived class, which are :
Using the template TRandomGen class (template on the contained Engine type), it is possible to add any generator based on the standard C++ random library (see the C++ random documentation.) or different variants of the MIXMAX generator using the ROOT::Math::MixMaxEngine. Some of the listed generator above (e.g. TRandomMixMax256 or TRandomMT64) are convenient typedef's of generator built using the template TRandomGen class.
Please note also that this class (TRandom) implements also a very simple generator (linear congruential) with period = \(10^9\), known to have defects (the lower random bits are correlated) and it is failing the majority of the random number generator tests. Therefore it should NOT be used in any statistical study.
The following table shows some timings (in nanoseconds/call) for the random numbers obtained using a macbookpro 2.6 GHz Intel Core i7 CPU:
The following methods are provided to generate random numbers disctributed according to some basic distributions:
::Exp(tau)
::Integer(imax)
::Gaus(mean,sigma)
::Rndm()
::Uniform(x1)
::Landau(mpv,sigma)
::Poisson(mean)
::Binomial(ntot,prob)
Random numbers distributed according to 1-d, 2-d or 3-d distributions contained in TF1, TF2 or TF3 objects can also be generated. For example, to get a random number distributed following abs(sin(x)/x)*sqrt(x) you can do :
or you can use the UNURAN package. You need in this case to initialize UNURAN to the function you would like to generate.
The techniques of using directly a TF1,2 or 3 function is powerful and can be used to generate numbers in the defined range of the function. Getting a number from a TF1,2,3 function is also quite fast. UNURAN is a powerful and flexible tool which containes various methods for generate random numbers for continuous distributions of one and multi-dimension. It requires some set-up (initialization) phase and can be very fast when the distribution parameters are not changed for every call.
The following table shows some timings (in nanosecond/call) for basic functions, TF1 functions and using UNURAN obtained running the tutorial math/testrandom.C Numbers have been obtained on an Intel Xeon Quad-core Harpertown (E5410) 2.33 GHz running Linux SLC4 64 bit and compiled with gcc 3.4
Note that the time to generate a number from an arbitrary TF1 function using TF1::GetRandom or using TUnuran is independent of the complexity of the function.
TH1::FillRandom(TH1 *) or TH1::FillRandom(const char *tf1name) can be used to fill an histogram (1-d, 2-d, 3-d from an existing histogram or from an existing function.
Note this interesting feature when working with objects. You can use several TRandom objects, each with their "independent" random sequence. For example, one can imagine
eventGenerator
can be used to generate the event kinematics. tracking can be used to track the generated particles with random numbers independent from eventGenerator. This very interesting feature gives the possibility to work with simple and very fast random number generators without worrying about random number periodicity as it was the case with Fortran. One can use TRandom::SetSeed to modify the seed of one generator.
A TRandom object may be written to a Root file
gRandom->Write("Random")
) ; Public Member Functions | |
TRandom (UInt_t seed=65539) | |
Default constructor. For seed see SetSeed(). More... | |
virtual | ~TRandom () |
Default destructor. More... | |
virtual Int_t | Binomial (Int_t ntot, Double_t prob) |
Generates a random integer N according to the binomial law. More... | |
virtual Double_t | BreitWigner (Double_t mean=0, Double_t gamma=1) |
Return a number distributed following a BreitWigner function with mean and gamma. More... | |
virtual void | Circle (Double_t &x, Double_t &y, Double_t r) |
Generates random vectors, uniformly distributed over a circle of given radius. More... | |
virtual Double_t | Exp (Double_t tau) |
Returns an exponential deviate. More... | |
virtual Double_t | Gaus (Double_t mean=0, Double_t sigma=1) |
Samples a random number from the standard Normal (Gaussian) Distribution with the given mean and sigma. More... | |
virtual UInt_t | GetSeed () const |
virtual UInt_t | Integer (UInt_t imax) |
Returns a random integer uniformly distributed on the interval [ 0, imax-1 ]. More... | |
virtual Double_t | Landau (Double_t mean=0, Double_t sigma=1) |
Generate a random number following a Landau distribution with location parameter mu and scale parameter sigma: Landau( (x-mu)/sigma ) Note that mu is not the mpv(most probable value) of the Landa distribution and sigma is not the standard deviation of the distribution which is not defined. More... | |
virtual Int_t | Poisson (Double_t mean) |
Generates a random integer N according to a Poisson law. More... | |
virtual Double_t | PoissonD (Double_t mean) |
Generates a random number according to a Poisson law. More... | |
virtual void | Rannor (Double_t &a, Double_t &b) |
Return 2 numbers distributed following a gaussian with mean=0 and sigma=1. More... | |
virtual void | Rannor (Float_t &a, Float_t &b) |
Return 2 numbers distributed following a gaussian with mean=0 and sigma=1. More... | |
virtual void | ReadRandom (const char *filename) |
Reads saved random generator status from filename. More... | |
virtual Double_t | Rndm () |
Machine independent random number generator. More... | |
virtual Double_t | Rndm (Int_t) |
virtual void | RndmArray (Int_t n, Double_t *array) |
Return an array of n random numbers uniformly distributed in ]0,1]. More... | |
virtual void | RndmArray (Int_t n, Float_t *array) |
Return an array of n random numbers uniformly distributed in ]0,1]. More... | |
virtual void | SetSeed (ULong_t seed=0) |
Set the random generator seed. More... | |
virtual void | Sphere (Double_t &x, Double_t &y, Double_t &z, Double_t r) |
Generates random vectors, uniformly distributed over the surface of a sphere of given radius. More... | |
virtual Double_t | Uniform (Double_t x1, Double_t x2) |
Returns a uniform deviate on the interval (x1, x2). More... | |
virtual Double_t | Uniform (Double_t x1=1) |
Returns a uniform deviate on the interval (0, x1). More... | |
virtual void | WriteRandom (const char *filename) const |
Writes random generator status to filename. More... | |
Public Member Functions inherited from TNamed | |
TNamed () | |
TNamed (const char *name, const char *title) | |
TNamed (const TNamed &named) | |
TNamed copy ctor. More... | |
TNamed (const TString &name, const TString &title) | |
virtual | ~TNamed () |
TNamed destructor. More... | |
virtual void | Clear (Option_t *option="") |
Set name and title to empty strings (""). More... | |
virtual TObject * | Clone (const char *newname="") const |
Make a clone of an object using the Streamer facility. More... | |
virtual Int_t | Compare (const TObject *obj) const |
Compare two TNamed objects. More... | |
virtual void | Copy (TObject &named) const |
Copy this to obj. More... | |
virtual void | FillBuffer (char *&buffer) |
Encode TNamed into output buffer. More... | |
virtual const char * | GetName () const |
Returns name of object. More... | |
virtual const char * | GetTitle () const |
Returns title of object. More... | |
virtual ULong_t | Hash () const |
Return hash value for this object. More... | |
virtual Bool_t | IsSortable () const |
virtual void | ls (Option_t *option="") const |
List TNamed name and title. More... | |
TNamed & | operator= (const TNamed &rhs) |
TNamed assignment operator. More... | |
virtual void | Print (Option_t *option="") const |
Print TNamed name and title. More... | |
virtual void | SetName (const char *name) |
Set the name of the TNamed. More... | |
virtual void | SetNameTitle (const char *name, const char *title) |
Set all the TNamed parameters (name and title). More... | |
virtual void | SetTitle (const char *title="") |
Set the title of the TNamed. More... | |
virtual Int_t | Sizeof () const |
Return size of the TNamed part of the TObject. More... | |
Public Member Functions inherited from TObject | |
TObject () | |
TObject constructor. More... | |
TObject (const TObject &object) | |
TObject copy ctor. More... | |
virtual | ~TObject () |
TObject destructor. More... | |
void | AbstractMethod (const char *method) const |
Use this method to implement an "abstract" method that you don't want to leave purely abstract. More... | |
virtual void | AppendPad (Option_t *option="") |
Append graphics object to current pad. More... | |
virtual void | Browse (TBrowser *b) |
Browse object. May be overridden for another default action. More... | |
ULong_t | CheckedHash () |
Check and record whether this class has a consistent Hash/RecursiveRemove setup (*) and then return the regular Hash value for this object. More... | |
virtual const char * | ClassName () const |
Returns name of class to which the object belongs. More... | |
virtual void | Clear (Option_t *="") |
virtual TObject * | Clone (const char *newname="") const |
Make a clone of an object using the Streamer facility. More... | |
virtual Int_t | Compare (const TObject *obj) const |
Compare abstract method. More... | |
virtual void | Copy (TObject &object) const |
Copy this to obj. More... | |
virtual void | Delete (Option_t *option="") |
Delete this object. More... | |
virtual Int_t | DistancetoPrimitive (Int_t px, Int_t py) |
Computes distance from point (px,py) to the object. More... | |
virtual void | Draw (Option_t *option="") |
Default Draw method for all objects. More... | |
virtual void | DrawClass () const |
Draw class inheritance tree of the class to which this object belongs. More... | |
virtual TObject * | DrawClone (Option_t *option="") const |
Draw a clone of this object in the current selected pad for instance with: gROOT->SetSelectedPad(gPad) . More... | |
virtual void | Dump () const |
Dump contents of object on stdout. More... | |
virtual void | Error (const char *method, const char *msgfmt,...) const |
Issue error message. More... | |
virtual void | Execute (const char *method, const char *params, Int_t *error=0) |
Execute method on this object with the given parameter string, e.g. More... | |
virtual void | Execute (TMethod *method, TObjArray *params, Int_t *error=0) |
Execute method on this object with parameters stored in the TObjArray. More... | |
virtual void | ExecuteEvent (Int_t event, Int_t px, Int_t py) |
Execute action corresponding to an event at (px,py). More... | |
virtual void | Fatal (const char *method, const char *msgfmt,...) const |
Issue fatal error message. More... | |
virtual TObject * | FindObject (const char *name) const |
Must be redefined in derived classes. More... | |
virtual TObject * | FindObject (const TObject *obj) const |
Must be redefined in derived classes. More... | |
virtual Option_t * | GetDrawOption () const |
Get option used by the graphics system to draw this object. More... | |
virtual const char * | GetIconName () const |
Returns mime type name of object. More... | |
virtual const char * | GetName () const |
Returns name of object. More... | |
virtual char * | GetObjectInfo (Int_t px, Int_t py) const |
Returns string containing info about the object at position (px,py). More... | |
virtual Option_t * | GetOption () const |
virtual const char * | GetTitle () const |
Returns title of object. More... | |
virtual UInt_t | GetUniqueID () const |
Return the unique object id. More... | |
virtual Bool_t | HandleTimer (TTimer *timer) |
Execute action in response of a timer timing out. More... | |
virtual ULong_t | Hash () const |
Return hash value for this object. More... | |
Bool_t | HasInconsistentHash () const |
Return true is the type of this object is known to have an inconsistent setup for Hash and RecursiveRemove (i.e. More... | |
virtual void | Info (const char *method, const char *msgfmt,...) const |
Issue info message. More... | |
virtual Bool_t | InheritsFrom (const char *classname) const |
Returns kTRUE if object inherits from class "classname". More... | |
virtual Bool_t | InheritsFrom (const TClass *cl) const |
Returns kTRUE if object inherits from TClass cl. More... | |
virtual void | Inspect () const |
Dump contents of this object in a graphics canvas. More... | |
void | InvertBit (UInt_t f) |
virtual Bool_t | IsEqual (const TObject *obj) const |
Default equal comparison (objects are equal if they have the same address in memory). More... | |
virtual Bool_t | IsFolder () const |
Returns kTRUE in case object contains browsable objects (like containers or lists of other objects). More... | |
R__ALWAYS_INLINE Bool_t | IsOnHeap () const |
virtual Bool_t | IsSortable () const |
R__ALWAYS_INLINE Bool_t | IsZombie () const |
virtual void | ls (Option_t *option="") const |
The ls function lists the contents of a class on stdout. More... | |
void | MayNotUse (const char *method) const |
Use this method to signal that a method (defined in a base class) may not be called in a derived class (in principle against good design since a child class should not provide less functionality than its parent, however, sometimes it is necessary). More... | |
virtual Bool_t | Notify () |
This method must be overridden to handle object notification. More... | |
void | Obsolete (const char *method, const char *asOfVers, const char *removedFromVers) const |
Use this method to declare a method obsolete. More... | |
void | operator delete (void *ptr) |
Operator delete. More... | |
void | operator delete[] (void *ptr) |
Operator delete []. More... | |
void * | operator new (size_t sz) |
void * | operator new (size_t sz, void *vp) |
void * | operator new[] (size_t sz) |
void * | operator new[] (size_t sz, void *vp) |
TObject & | operator= (const TObject &rhs) |
TObject assignment operator. More... | |
virtual void | Paint (Option_t *option="") |
This method must be overridden if a class wants to paint itself. More... | |
virtual void | Pop () |
Pop on object drawn in a pad to the top of the display list. More... | |
virtual void | Print (Option_t *option="") const |
This method must be overridden when a class wants to print itself. More... | |
virtual Int_t | Read (const char *name) |
Read contents of object with specified name from the current directory. More... | |
virtual void | RecursiveRemove (TObject *obj) |
Recursively remove this object from a list. More... | |
void | ResetBit (UInt_t f) |
virtual void | SaveAs (const char *filename="", Option_t *option="") const |
Save this object in the file specified by filename. More... | |
virtual void | SavePrimitive (std::ostream &out, Option_t *option="") |
Save a primitive as a C++ statement(s) on output stream "out". More... | |
void | SetBit (UInt_t f) |
void | SetBit (UInt_t f, Bool_t set) |
Set or unset the user status bits as specified in f. More... | |
virtual void | SetDrawOption (Option_t *option="") |
Set drawing option for object. More... | |
virtual void | SetUniqueID (UInt_t uid) |
Set the unique object id. More... | |
virtual void | SysError (const char *method, const char *msgfmt,...) const |
Issue system error message. More... | |
R__ALWAYS_INLINE Bool_t | TestBit (UInt_t f) const |
Int_t | TestBits (UInt_t f) const |
virtual void | UseCurrentStyle () |
Set current style settings in this object This function is called when either TCanvas::UseCurrentStyle or TROOT::ForceStyle have been invoked. More... | |
virtual void | Warning (const char *method, const char *msgfmt,...) const |
Issue warning message. More... | |
virtual Int_t | Write (const char *name=0, Int_t option=0, Int_t bufsize=0) |
Write this object to the current directory. More... | |
virtual Int_t | Write (const char *name=0, Int_t option=0, Int_t bufsize=0) const |
Write this object to the current directory. More... | |
Public Member Functions inherited from ROOT::Math::TRandomEngine | |
virtual | ~TRandomEngine () |
virtual double | Rndm ()=0 |
Protected Attributes | |
UInt_t | fSeed |
Protected Attributes inherited from TNamed | |
TString | fName |
TString | fTitle |
Additional Inherited Members | |
Public Types inherited from TObject | |
enum | { kIsOnHeap = 0x01000000 , kNotDeleted = 0x02000000 , kZombie = 0x04000000 , kInconsistent = 0x08000000 , kBitMask = 0x00ffffff } |
enum | { kSingleKey = BIT(0) , kOverwrite = BIT(1) , kWriteDelete = BIT(2) } |
enum | EDeprecatedStatusBits { kObjInCanvas = BIT(3) } |
enum | EStatusBits { kCanDelete = BIT(0) , kMustCleanup = BIT(3) , kIsReferenced = BIT(4) , kHasUUID = BIT(5) , kCannotPick = BIT(6) , kNoContextMenu = BIT(8) , kInvalidObject = BIT(13) } |
Static Public Member Functions inherited from TObject | |
static Long_t | GetDtorOnly () |
Return destructor only flag. More... | |
static Bool_t | GetObjectStat () |
Get status of object stat flag. More... | |
static void | SetDtorOnly (void *obj) |
Set destructor only flag. More... | |
static void | SetObjectStat (Bool_t stat) |
Turn on/off tracking of objects in the TObjectTable. More... | |
Protected Member Functions inherited from TObject | |
virtual void | DoError (int level, const char *location, const char *fmt, va_list va) const |
Interface to ErrorHandler (protected). More... | |
void | MakeZombie () |
#include <TRandom.h>
TRandom::TRandom | ( | UInt_t | seed = 65539 | ) |
Default constructor. For seed see SetSeed().
Definition at line 175 of file TRandom.cxx.
|
virtual |
Default destructor.
Can reset gRandom to 0 if gRandom points to this generator.
Definition at line 184 of file TRandom.cxx.
Generates a random integer N according to the binomial law.
Coded from Los Alamos report LA-5061-MS.
N is binomially distributed between 0 and ntot inclusive with mean prob*ntot and prob is between 0 and 1.
Note: This function should not be used when ntot is large (say >100). The normal approximation is then recommended instead (with mean =*ntot+0.5 and standard deviation sqrt(ntot*prob*(1-prob)).
Definition at line 200 of file TRandom.cxx.
Return a number distributed following a BreitWigner function with mean and gamma.
Definition at line 214 of file TRandom.cxx.
Generates random vectors, uniformly distributed over a circle of given radius.
Input : r = circle radius Output: x,y a random 2-d vector of length r
Definition at line 228 of file TRandom.cxx.
Samples a random number from the standard Normal (Gaussian) Distribution with the given mean and sigma.
Uses the Acceptance-complement ratio from W. Hoermann and G. Derflinger This is one of the fastest existing method for generating normal random variables. It is a factor 2/3 faster than the polar (Box-Muller) method used in the previous version of TRandom::Gaus. The speed is comparable to the Ziggurat method (from Marsaglia) implemented for example in GSL and available in the MathMore library.
REFERENCE: - W. Hoermann and G. Derflinger (1990): The ACR Method for generating normal random variables, OR Spektrum 12 (1990), 181-185.
Implementation taken from UNURAN (c) 2000 W. Hoermann & J. Leydold, Institut f. Statistik, WU Wien
Definition at line 263 of file TRandom.cxx.
|
inlinevirtual |
Returns a random integer uniformly distributed on the interval [ 0, imax-1 ].
Note that the interfal contains the values of 0 and imax-1 but not imax.
Definition at line 349 of file TRandom.cxx.
Generate a random number following a Landau distribution with location parameter mu and scale parameter sigma: Landau( (x-mu)/sigma ) Note that mu is not the mpv(most probable value) of the Landa distribution and sigma is not the standard deviation of the distribution which is not defined.
For mu =0 and sigma=1, the mpv = -0.22278
The Landau random number generation is implemented using the function landau_quantile(x,sigma), which provides the inverse of the landau cumulative distribution. landau_quantile has been converted from CERNLIB ranlan(G110).
Definition at line 369 of file TRandom.cxx.
Generates a random integer N according to a Poisson law.
Prob(N) = exp(-mean)*mean^N/Factorial(N)
Use a different procedure according to the mean value. The algorithm is the same used by CLHEP. For lower value (mean < 25) use the rejection method based on the exponential. For higher values use a rejection method comparing with a Lorentzian distribution, as suggested by several authors. This routine since is returning 32 bits integer will not work for values larger than 2*10**9. One should then use the Trandom::PoissonD for such large values.
Definition at line 391 of file TRandom.cxx.
Generates a random number according to a Poisson law.
Prob(N) = exp(-mean)*mean^N/Factorial(N)
This function is a variant of TRandom::Poisson returning a double instead of an integer.
Definition at line 443 of file TRandom.cxx.
Return 2 numbers distributed following a gaussian with mean=0 and sigma=1.
Definition at line 504 of file TRandom.cxx.
Return 2 numbers distributed following a gaussian with mean=0 and sigma=1.
Definition at line 489 of file TRandom.cxx.
|
virtual |
Reads saved random generator status from filename.
Definition at line 519 of file TRandom.cxx.
|
virtual |
Machine independent random number generator.
Based on the BSD Unix (Rand) Linear congrential generator. Produces uniformly-distributed floating points between 0 and 1. Identical sequence on all machines of >= 32 bits. Periodicity = 2**31, generates a number in (0,1). Note that this is a generator which is known to have defects (the lower random bits are correlated) and therefore should NOT be used in any statistical study).
Implements ROOT::Math::TRandomEngine.
Reimplemented in TRandom1, TRandom1, TRandom2, TRandom2, TRandom3, TRandom3, TRandomGen< Engine >, and TRandomGen< Engine >.
Definition at line 541 of file TRandom.cxx.
Return an array of n random numbers uniformly distributed in ]0,1].
Reimplemented in TRandom2, TRandom3, TRandomGen< Engine >, and TRandom1.
Definition at line 564 of file TRandom.cxx.
Return an array of n random numbers uniformly distributed in ]0,1].
Reimplemented in TRandom2, TRandom3, TRandomGen< Engine >, and TRandom1.
Definition at line 577 of file TRandom.cxx.
Set the random generator seed.
Note that default value is zero, which is different than the default value used when constructing the class. If the seed is zero the seed is set to a random value which in case of TRandom depends on the lowest 4 bytes of TUUID The UUID will be identical if SetSeed(0) is called with time smaller than 100 ns Instead if a different generator implementation is used (TRandom1, 2 or 3) the seed is generated using a 128 bit UUID. This results in different seeds and then random sequence for every SetSeed(0) call.
Reimplemented in TRandom1, TRandom2, TRandom3, and TRandomGen< Engine >.
Definition at line 597 of file TRandom.cxx.
Generates random vectors, uniformly distributed over the surface of a sphere of given radius.
Input : r = sphere radius Output: x,y,z a random 3-d vector of length r Method: (based on algorithm suggested by Knuth and attributed to Robert E Knop) which uses less random numbers than the CERNLIB RN23DIM algorithm
Definition at line 617 of file TRandom.cxx.
Returns a uniform deviate on the interval (x1, x2).
Definition at line 644 of file TRandom.cxx.
Returns a uniform deviate on the interval (0, x1).
Definition at line 635 of file TRandom.cxx.
|
virtual |
Writes random generator status to filename.
Definition at line 653 of file TRandom.cxx.