// Polinomial fit

Double_t r_dif_pol5(Double_t *x, Double_t *par) {

  // pol5

  return par[0]+par[1]*x[0]+par[2]*x[0]*x[0]+
         par[3]*x[0]*x[0]*x[0]+par[4]*x[0]*x[0]*x[0]*x[0];

  }