WCSimRandomParameters.hh

Go to the documentation of this file.

#ifndef WCSimRandomParameters_h
#define WCSimRandomParameters_h 1
#include "WCSimRandomMessenger.hh"
#include "WCSimRootOptions.hh"
#include "CLHEP/Random/Random.h"
#include "CLHEP/Random/RanluxEngine.h"
#include "CLHEP/Random/JamesRandom.h"
#include "CLHEP/Random/RanecuEngine.h"

class WCSimRandomParameters
{

public:
  WCSimRandomParameters() 
    {
      seed=31415; 
      generator=RANDOM_E_HEPJAMES;
      RandomMessenger = new WCSimRandomMessenger(this);
    }
  ~WCSimRandomParameters() {delete RandomMessenger;}

  WCSimRandomGenerator_t GetGenerator() {return generator; }
  void SetGenerator(WCSimRandomGenerator_t rng)
  {
    switch (rng)
      {
        case RANDOM_E_RANLUX:
          {
            printf("Setting the random number generator to RANLUX\n");
            CLHEP::RanluxEngine *newluxengine = new CLHEP::RanluxEngine(31415,4);  // highest luxury level
            CLHEP::HepRandom::setTheEngine(newluxengine);
          }
          break;
      case RANDOM_E_RANECU:
          {
            printf("Setting the random number generator to RANECU\n");
            CLHEP::RanecuEngine *newecuengine = new CLHEP::RanecuEngine();
            CLHEP::HepRandom::setTheEngine(newecuengine);
          }
          break;
          
        case RANDOM_E_HEPJAMES:
          {
            printf("Setting the random number generator to HEPJAMES\n");
            CLHEP::HepJamesRandom *newjamesengine = new CLHEP::HepJamesRandom();
            CLHEP::HepRandom::setTheEngine(newjamesengine);
          }  
          break;
      default:
        {
          printf("Random number generator type not understood: %d\n",rng);
          exit(0);
        }
      }
  };
  int GetSeed() {return CLHEP::HepRandom::getTheSeed();}
  void SetSeed(int iseed) 
    { 
      CLHEP::HepRandom::setTheSeed(iseed);
      printf("Setting the Random Seed to: %d\n",iseed); 
      seed = iseed;
    }

  void SaveOptionsToOutput(WCSimRootOptions * wcopt)
  {
    wcopt->SetRandomSeed(seed);
    wcopt->SetRandomGenerator(generator);
  }

private:
  WCSimRandomGenerator_t generator;
  int seed;
  WCSimRandomMessenger *RandomMessenger;
};

#endif