Double_t ffit(Double_t *x,Double_t *par)
{
  Double_t x1 = x[0]-par[2];
  Double_t gauss5 = 5.0*par[1]*TMath::Exp(-0.5*x1^2/par[3]^2)/sqrt(2*3.1415926)*par[3]+par[4]+par[5]*x1;
  return gauss5;
}