TDFun.h

Go to the documentation of this file.

#ifndef RAD_TDFun
#define RAD_TDFun
 
#include "TRadSpline.h"
#include <cmath>
 
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