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