remove the gsl dependence: use TGraph interpolation methods instead