// Gaussian fit Double_t z_p2e_gaus(Double_t *x, Double_t *par) { // par[0]=norm, par[2]=peak, par[2]=rms // par[0]=7127, par[1]=1.505, par[2]=1.501 Double_t t = (x[0]-par[1])/par[2]); return par[0]*TMath::Exp(-t*t); }