Geant4Random.h

Go to the documentation of this file.

//==========================================================================
//  AIDA Detector description implementation 
//--------------------------------------------------------------------------
// Copyright (C) Organisation europeenne pour la Recherche nucleaire (CERN)
// All rights reserved.
//
// For the licensing terms see $DD4hepINSTALL/LICENSE.
// For the list of contributors see $DD4hepINSTALL/doc/CREDITS.
//
// Author     : M.Frank
//
//==========================================================================
#ifndef DDG4_GEANT4RANDOM_H
#define DDG4_GEANT4RANDOM_H
 
// Framework include files
#include <DDG4/Geant4Action.h>
 
// C/C++ include files
#include <string>
 
// Forward declarations
class TRandom;
 
namespace CLHEP   {  class HepRandomEngine;   }
 
namespace dd4hep {
 
  namespace sim {
 
    // Forward declarations
    class Geant4Exec;
 
 
    class Geant4Random : public Geant4Action   {
      friend class Geant4Exec;
 
    protected:
      std::string  m_file;
      std::string  m_engineType;
      long         m_seed, m_luxury;
      bool         m_replace;
      
      CLHEP::HepRandomEngine* m_engine;
 
      TRandom *m_rootRandom, *m_rootOLD;
      bool m_inited;
      
    public:
 
      Geant4Random(Geant4Context* context, const std::string& name);
      virtual ~Geant4Random();
 
      static Geant4Random* instance(bool throw_exception=true);
 
      static Geant4Random* setMainInstance(Geant4Random* ptr);
 
 
      void initialize();
      
      CLHEP::HepRandomEngine* engine()      {  return m_engine;   }
      
      virtual void setSeed(long seed);
 
      virtual void setSeeds(const long * seeds, int size);
      virtual void saveStatus( const char filename[] = "Config.conf") const;
      virtual void restoreStatus( const char filename[] = "Config.conf" );
      virtual void showStatus() const;
 
 
      double rndm_clhep();
 
      double rndm(int i=0);
      void   rndmArray(int n, float *array);
      void   rndmArray(int n, double *array);
      double uniform(double x1=1);
      double uniform(double x1, double x2);
 
      double exp(double tau);
      double gauss(double mean=0, double sigma=1);
      double gamma(double k, double lambda);
      double landau(double mean=0, double sigma=1);
      void   circle(double &x, double &y, double r);
      void   sphere(double &x, double &y, double &z, double r);
      double poisson(double mean=1e0 );
      double breit_wigner(double mean=0e0, double gamma=1e0);
    };
 
  }    // End namespace sim
}      // End namespace dd4hep
#endif // DDG4_GEANT4RANDOM_H