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