10 const Double_t dxbin = (0.17-0.13)/40;
11 const Double_t sigma = 0.0012;
14 Double_t fdm5(Double_t *xx, Double_t *par)
17 if (x <= 0.13957)
return 0;
18 Double_t xp3 = (x-par[3])*(x-par[3]);
19 Double_t res = dxbin*(par[0]*TMath::Power(x-0.13957, par[1])
20 + par[2] / 2.5066/par[4]*TMath::Exp(-xp3/2/par[4]/par[4]));
25 Double_t fdm2(Double_t *xx, Double_t *par)
28 if (x <= 0.13957)
return 0;
29 Double_t xp3 = (x-0.1454)*(x-0.1454);
30 Double_t res = dxbin*(par[0]*TMath::Power(x-0.13957, 0.25)
31 + par[1] / 2.5066/sigma*TMath::Exp(-xp3/2/sigma/sigma));