Distributor.cxx

Go to the documentation of this file.

 
#include "eicsmear/smear/Distributor.h"
 
#include <TF1.h>
#include <TRandom.h>
#include <TUUID.h>
 
namespace Smear {
 
Distributor::Distributor()
: mPlus(0.)
, mMinus(0.)
, mDistribution(NULL) {
}
 
Distributor::Distributor(const TString& formula, double lower, double upper,
                         double minimum, double maximum)
: mPlus(0.)
, mMinus(0.)
, mDistribution(new TF1(TUUID().AsString(), formula, minimum, maximum)) {
  mDistribution->SetParameters(0., 1.);
  if (lower > 0.) {
    mMinus = lower;
  }  // if
  if (upper > 0.) {
    mPlus = upper;
  }  // if
}
 
Distributor::~Distributor() {
  if (mDistribution) {
    delete mDistribution;
    mDistribution = NULL;
  }  // if
}
 
double Distributor::Generate(double mean, double sigma) {
  double random(0.);
  if (!mDistribution) {
    random = gRandom->Gaus(mean, sigma);
  } else {
    mDistribution->SetParameters(mean, sigma);
    if (mMinus > 0. || mPlus > 0.) {
      random = mDistribution->GetRandom(mean - mMinus, mean + mPlus);
    } else {
      random = mDistribution->GetRandom();
    }  // if
  }  // if
  return random;
}
 
}  // namespace Smear