TDFun.h

Go to the documentation of this file.

#ifndef RAD_TDFun
#define RAD_TDFun
 
#include <cmath>
#include "TRadSpline.h"
 
class TDFun{
protected:
  double fxmin;
  double fxmax;
  double fD0;
  double fb2;
  TRadSpline3 *fd;
public:
  TDFun();
  TDFun(const double&, const double&);
  ~TDFun();
  void Init();
  void Init(const double&, const double&);
 
  inline double EvalSpline(const double &z){return fd->Eval(z);}
 
  inline double Eval(const double &z){  
    double x = 1 - z;
    //    double D = fb2*pow(x,fb2-1)*fD0 - 0.5*fb2*(1+z) - 
    //      fb2*fb2/8*(4*(1+z)*log(x)+(1+3*z*z)/x*log(z)+5+z);
    double D = fb2*pow(x,fb2-1)*fD0 - 0.5*fb2*(1+z) - 
      fb2*fb2/8*(4*(1+z)*log(x)+(1+3*z*z)/x*log1p(-x)+5+z);
    return D;
  }
 
  inline double Eval(const double &z, const double &b2, const double &D0){
    double x = 1 - z;
    //    double D = b2*pow(x,b2-1)*D0 - 0.5*b2*(1+z) - 
    //      b2*b2/8*(4*(1+z)*log(x)+(1+3*z*z)/x*log(z)+5+z);
    double D = b2*pow(x,b2-1)*D0 - 0.5*b2*(1+z) - 
      b2*b2/8*(4*(1+z)*log(x)+(1+3*z*z)/x*log1p(-x)+5+z);
    return D;
  }
 
  inline double EvalSoft(const double &x){
    double z  = 1 - x;
    double lx = log(x);
    double t1 = exp(-fb2*lx);
    double t2 = t1*x;
    //    double D  = fD0 - t2*(0.5*(1+z) + fb2/8*( 4*(1+z)*lx + 4 + (1+z))) - 
    //      fb2/8*t1*(1+3*z*z)*log(z);
    double D  = fD0 - t2*(0.5*(1+z) + fb2/8*( 4*(1+z)*lx + 4 + (1+z))) - 
      fb2/8*t1*(1+3*z*z)*log1p(-x);
    return D;
  }
 
  inline double EvalSoft(const double &x, const double &b2, const double &D0){
    double z  = 1 - x;
    double lx = log(x);
    double t1 = exp(-b2*lx);
    double t2 = t1*x;
    //    double D  = D0 - t2*(0.5*(1+z) + b2/8*( 4*(1+z)*lx + 4 + (1+z))) - 
    //      b2/8*t1*(1+3*z*z)*log(z);
    double D  = D0 - t2*(0.5*(1+z) + b2/8*( 4*(1+z)*lx + 4 + (1+z))) - 
      b2/8*t1*(1+3*z*z)*log1p(-x);
    return D;
  }
 
};
 
#endif //#ifndef RAD_TDFun