00001 #include "TDFun.h" 00002 //#include "TF1.h" 00003 #include "TRadGlobal.h" 00004 00005 TDFun *djr = NULL; 00006 00007 double tdfun(const double &x){ 00008 return djr->Eval(x); 00009 } 00010 /* 00011 double tdfun(double *x, double *p){ 00012 return djr->Eval(*x); 00013 } 00014 */ 00015 TDFun::TDFun(){ 00016 fd = NULL; 00017 Init(1-gGlobal->Get_XMax(),1-gGlobal->Get_XMin()); 00018 } 00019 00020 TDFun::TDFun(const double &cmin, const double &cmax){ 00021 fd = NULL; 00022 Init(cmin,cmax); 00023 } 00024 00025 void TDFun::Init(){ 00026 Init(1-gGlobal->Get_XMax(),1-gGlobal->Get_XMin()); 00027 } 00028 00029 void TDFun::Init(const double &cmin, const double &cmax){ 00030 fxmin = cmin; 00031 fxmax = cmax; 00032 fb2 = 0.5*gGlobal->Get_Beta(); 00033 fD0 = gGlobal->Get_D0(); 00034 djr = this; 00035 // TF1 fdfun("dfun",tdfun,fxmin,fxmax,0); 00036 if(fd) delete fd; 00037 fd = new TRadSpline3("d", fxmin, fxmax, tdfun, 10001); 00038 } 00039 00040 TDFun::~TDFun(){ 00041 if(fd) delete fd; 00042 }