UnivTimeKG.h

Go to the documentation of this file.

#ifndef __UnivTimeKG_UnivTimeKG_h_
#define __UnivTimeKG_UnivTimeKG_h_

#include <vector>
#include <string>

namespace UnivTimeKG {

  unsigned int FindIndex(const std::vector<double>& array, const double value);


  double int1_2p(const double x0, const double y0, const double x1, const double y1, const double x);


  double int1(const std::vector<double>& xs, const std::vector<double>& ys, const double x);


  double int1(const double xs[], const double ys[], const double x, const int n);


  double int2(const std::vector<double>& x1s, const std::vector<double>& x2s,
              const std::vector<std::vector<double> >& Ys, const double x1, const double x2);


  struct IntTable {
    std::vector<double> rs;
    std::vector<double> psis;
    std::vector<std::vector<double> > ys;
  };


  class TimeModel {

  public:
    TimeModel();
    virtual ~TimeModel();

    void addInterpolationTable(const std::vector<double>& xs,
                               const std::vector<double>& ys, const std::vector<std::vector<double> >& data);
    void addInterpolationTable(const std::string& filename);
    void clearInterpolationTables();
    virtual double getShapeParameter(const unsigned int ipar, const std::vector<double>& pars, const double DX) = 0;
    void interpolateParameters(const double DX, const double r, const double psi, std::vector<double>& output);

    virtual void setShapeParameters(const double DX, const double r, const double psi, const double theta, const double lgE) = 0;

    virtual void setParameterOffsets(const double m, const double s) = 0;

    virtual double pdf(const double) = 0;
    virtual double cdf(const double) = 0;
    double firstParticlePdf(const double, const double);
    double firstParticlePdfSmeared(const double, const double);
    virtual double invcdf(const double) = 0;
    double getRisetime();
    double getFalltime();

    std::vector<IntTable> tables;
    static const int interpMode = 1;
    unsigned int nParams;

    unsigned int icomp;

  };

}


#endif