NoiseSaS.cc

Go to the documentation of this file.

#include <utl/MathConstants.h>
#include <utl/Math.h>
#include <utl/NoiseSaS.h>
#include <utl/RandomEngine.h>
#include <CLHEP/Random/RandFlat.h>
#include <utl/Trace.h>
#include <time.h>
#include <utl/MathConstants.h>
#include <utl/RandomEngine.h>
#include <CLHEP/Random/RandExponential.h>
#include <fwk/RandomEngineRegistry.h>
#include <fwk/CentralConfig.h>
#include <utl/config.h>

using namespace std;
using namespace utl;
using namespace fwk;
using CLHEP::RandExponential;


double
NoiseSaS::CreateSaSNoiseInTimeDomain(const double alpha)
{
  // Create an SaS noise sample
  const double v = fRandomGeneratorV();
  const double w = RandExponential::shoot(&fRandomEngine->GetEngine(), 1);
  const int i = random() & 1 ? 1 : -1;
  // calculate the following SaS trace sample
  const double x1 = sin(alpha*v) / pow(cos(v), 1/alpha);
  const double x2 = pow(cos(v-alpha*v)/w, (1 - alpha)/alpha);
  return i*x1*x2;
}