// AvdS/AP function, 5 parameters

Double_t e_mich(Double_t *x, Double_t *par) {
  Double_t result = 0; 

// par[0]=3.80361e+01 par[1]=-1.48152e+00
// par[2]=4.60782e+01 par[3]=5.02199e+01 par[4]=4.07456e+00

 //without gaussian edges (gives errors)
 result=fabs(TMath::Erf(x[0])*(x[0]-par[1])**2*(3-2*x[0]/par[2])*(1-TMath::Erf((x[0]-par[3])/par[4]))); 

 //with Gaussians at edges
 //if (x[0]<=2) 
 //  result=exp(x[0]);
 //if (x[0]>2 && x[0]<=par[3])
 //  result=fabs((x[0]-par[1])**2*(3-2*x[0]/par[2])*(1-TMath::Erf((x[0]-par[3])/par[4])));
 //if (x[0]>par[3]) 
 //  result=1000*exp(par[3]-x[0]);

 //original e_mich
 //if (x[0]<=2)
 //  result=exp(x[0])
 //if ((x[0]>2 && x[0]<=par[3])
 //  result=abs((x[0]-par[1])**2*(3-2*x[0]/par[2])*(1-TMath::Erf((x[0]-par[3])/par[4])));
 //if (x[0]>par[3])
 //  result=1000*exp(par[3]-x[0]);

 result=par[0]*result;

 return result;
}