// GLG4TimeGen.cc
// Contact person: Phil Jones <p.g.jones@qmul.ac.uk>
// See GLG4TimeGen.hh for more details
//———————————————————————//

#include <RAT/GLG4TimeGen.hh>
#include <Randomize.hh>
#include <RAT/GLG4StringUtil.hh>
#include "CLHEP/Units/SystemOfUnits.h"

using namespace CLHEP;

void GLG4TimeGen_Uniform::SetState(G4String state)
{
  state = util_strip_default(state);
  if (state.length() == 0) {
    // print help and current state
    G4cout << "Current state of this GLG4TimeGen:\n"
	   << " \"" << GetState() << "\"\n" << G4endl;
    G4cout << "Format of argument to GLG4TimeGen::SetState: \n"
      " \"rate\"\n"
      " where rate is in 1/sec."
	   << G4endl;
    return;
  }

  fRate = util_to_double(state) / s;
}

G4String GLG4TimeGen_Uniform::GetState() const
{
  return util_dformat("%lf", fRate);
}

// *****************************************************

double GLG4TimeGen_Poisson::GenerateEventTime(double offset)
{
  return -log(1.0-G4UniformRand())/fRate + offset;
}