UnivTimeKGGamma.h

Go to the documentation of this file.

#ifndef __UnivTimeKG_GammaTimeModel_h_
#define __UnivTimeKG_GammaTimeModel_h_

#include <tls/UnivTimeKG.h>
#include <utl/Polynomial.h>

#include <vector>


namespace UnivTimeKG {

  class GammaTimeModel : public TimeModel {

  public:
    GammaTimeModel();
    GammaTimeModel(const int i);
    virtual ~GammaTimeModel() { }

    void CalculateModel(const double DX, const double r, const double psi,
                        const double theta, const double lgE, std::vector<double>& output);

    double getShapeParameter(const unsigned int ipar, const std::vector<double>& pars, const double DX);

    void setShapeParameters(const double DX, const double r, const double psi, const double theta, const double lgE);
    void setShapeParametersDirectly(const double mm, const double ss, const double ll);

    void setParameterOffsets(const double m, const double s);

    double pdf(const double);
    double cdf(const double);
    double cdf_analytical(const double* const x, const double* const p);
    double invcdf(const double);

    std::vector< std::vector<double> > fRpoints;
    std::vector< std::vector<utl::Polynomial> > fPolys;

    double m, s, l;
    double moff, soff;

  };

}


#endif