/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_13_patch3/src/FastSimulation/Utilities/interface/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 RandomEngine;

class GammaFunctionGenerator
{
 public:

  GammaFunctionGenerator(const RandomEngine* engine);

  virtual ~GammaFunctionGenerator();
  
  double shoot() const;
  
    void setParameters(double a,double b, double xm);

 private:
  double gammaFrac() const;
  double gammaInt() 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;

  const RandomEngine* random;

};
#endif