StarRandom.h

#ifndef __StarRandom_h__
#define __StarRandom_h__

#include "TObject.h"
#include "Math/GSLRndmEngines.h"
#include "TVector2.h"
#include "TRandom.h"
#include "TArrayC.h"

class StarRandom : public TObject
{
  StarRandom();
 public:
  ~StarRandom();

  static StarRandom &Instance();

  static void set( ROOT::Math::GSLRandomEngine *engine );

  static void capture();

  static void seed( UInt_t s );

  static void seed( UShort_t seed1, UShort_t seed2 );

  Double_t operator()() const;
  Double_t operator()( const Double_t mn, const Double_t mx ) const;

  Double_t operator()( Int_t *idummy ){ return (*this)(); }


  Double_t flat() const;
  Double_t flat( const Double_t mn, const Double_t mx ) const;

  Double_t expo( const Double_t mu ) const;
  Double_t landau() const;
  Double_t gauss( const Double_t sigma ) const ;
  TVector2 gauss2d( const Double_t sx, const Double_t sy, const Double_t rho ) const;

  UInt_t poisson( const Double_t mu ) const;

 private:
 protected:

  static StarRandom *sInstance;                       
  static ROOT::Math::GSLRandomEngine *mEngine;        

  UInt_t   mSeed;                                     // Random number seed
  TArrayC  mState;                                    // Random number state
  
  ClassDef(StarRandom,1);

};

#endif