BabayagaNLORandom.cxx

Go to the documentation of this file.

//--------------------------------------------------------------------------
//
// Environment:
//      This software is part of the EvtGen package developed jointly
//      for the BaBar and CLEO collaborations.  If you use all or part
//      of it, please give an appropriate acknowledgement.
//
// Copyright Information: See EvtGen/COPYRIGHT
//      Copyright (C) 1998      Caltech, UCSB
//
// Module: EvtRandom.cc
//
// Description: routines to get random numbers from
//              random number generator.
//
// Modification history:
//
//    DJL/RYD   September 25, 1996           Module created
//    Wang Boqun September 06, 2010          Add fortran interface
//
//------------------------------------------------------------------------
//
 
#include "BabayagaNLORandom.h"
 
#include <iostream>
using namespace std;
using namespace CLHEP;
 
HepRandomEngine* BabayagaNLORandom::_randomEngine = 0;
 
void BabayagaNLORandom::setRandomEngine( CLHEP::HepRandomEngine* randomEngine ) {
  _randomEngine = randomEngine;
}
 
double BabayagaNLORandom::random() {
 
  if ( _randomEngine == 0 )
  {
    cerr << "No random engine available in "
         << "BabayagaNLORandom::random()." << endl;
  }
 
  return _randomEngine->flat();
}
 
double BabayagaNLORandom::Flat( double min, double max ) {
 
  if ( min > max )
  {
    cerr << "min>max in BabayagaNLORandom::Flat(" << min << "," << max << ")" << endl;
    double tmp = min;
    min        = max;
    max        = tmp;
  }
 
  return BabayagaNLORandom::random() * ( max - min ) + min;
}
 
double BabayagaNLORandom::Flat( double max ) { return max * BabayagaNLORandom::random(); }
 
double BabayagaNLORandom::Flat() { return BabayagaNLORandom::random(); }
 
void BabayagaNLORandom::FlatArray( double* vect, const int size ) {
  if ( _randomEngine == 0 )
    cout << "Can not get randomEngine pointer in BabayagaNLORandom::FlatArray" << endl;
  else { _randomEngine->flatArray( size, vect ); }
}
 
void BabayagaNLORandom::FlatArrayF( float* vect, const int size ) {
 
  if ( _randomEngine == 0 )
  {
    cout << "Can not get randomEngine pointer in BabayagaNLORandom::FlatArray" << endl;
    return;
  }
  else
  {
    for ( int ir = 0; ir < size; ir++ ) vect[ir] = _randomEngine->flat();
  }
}
 
extern "C" {
 
void babayaga_ranlux_( float* vect, const int& size ) {
  BabayagaNLORandom::FlatArrayF( vect, size );
}
 
void babayaga_ranluxd_( double* vect, const int& size ) {
  BabayagaNLORandom::FlatArray( vect, size );
}
}