Geant4Random.h

Go to the documentation of this file.

// $Id: $
//==========================================================================
//  AIDA Detector description implementation for LCD
//--------------------------------------------------------------------------
// 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 DD4HEP_GEANT4RANDOM_H
#define DD4HEP_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 Simulation {

    // 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 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);
    };

  }    // End namespace Simulation
}      // End namespace DD4hep
#endif // DD4HEP_GEANT4RANDOM_H