GammaFunctionGenerator.h

Go to the documentation of this file.

#ifndef GammaFunctionGenerator_H
#define GammaFunctionGenerator_H

// FAMOS headers
#include "FastSimulation/Utilities/interface/GammaNumericalGenerator.h"

// CLHEP
#include "CLHEP/GenericFunctions/IncompleteGamma.hh"

//STL
#include <vector>

class RandomEngineAndDistribution;

class GammaFunctionGenerator {
public:
  GammaFunctionGenerator();

  virtual ~GammaFunctionGenerator();

  double shoot(RandomEngineAndDistribution const*) const;

  void setParameters(double a, double b, double xm);

private:
  double gammaFrac(RandomEngineAndDistribution const*) const;
  double gammaInt(RandomEngineAndDistribution const*) const;

private:
  // The integer numerical functions
  std::vector<GammaNumericalGenerator> theGammas;

  // The gamma distribution core coefficients
  std::vector<double> coreCoeff;

  // The gamma distribution core proba
  double coreProba;

  // possibility to store different limits
  std::vector<double> approxLimit;

  // boundaries
  double xmin;
  double xmax;

  // closest lower integer
  unsigned na;
  // alpha-na
  double frac;
  // alpha function parameters
  double alpha, beta;
  //  Incomlete Gamma = Int(0,x)[t^(alpha-1)exp(-t)dt]/Gamma(alpha);
  Genfun::IncompleteGamma myIncompleteGamma;

  // some useful integrals
  std::vector<double> integralToApproxLimit;

  // if xmin>xmax
  bool badRange;
};
#endif