Below in Figure 1 are the results of the inelastic proton dilution factor as a function of detector number after rebinning simulated FASTBUS inelastic+elastic protons in the fashion of the French TOF spectra. The inelastic+elastic protons are fitted with two gaussians. The results of the fit are numerically integrated to find the contribution of the elastic and inelastic protons in the elastic cut window (taken to be +/-2sigma).
The top plot in Figure 1 is a comparison of the inelastic dilution factor by the method of Lars(magenta crosses),by the numerical integration of the fit of the simulated data(blue triangles), by use of the Erf function on the fit parameters(red square)of the simulated data, and the "true" dilution factor and found by getting the bin contents of the inelastic and elastic histograms(black circle) as a function of detector number. The values of the numerical integration(and the use of the Erf function) of the fit is systematically smaller than the Lars results, this is due to the fact that the simulated FASTBUS data was generated by fit parameters found by Benoit's fitting program for the French detectors in run 16096.
The bottom plot in Figure 1 is a comparison of the "true" inelastic proton dilution factor(that is, the dilution factor found by getting the bin contents of the individual elastic and inelastic histograms within the elastic cut window) to the inelastic proton dilution factor found by using the two gaussian fit and then using numerical integration or the Erf function to find the elastic and inelastic contributions within the +/-2sigma elastic cut window.
Figure 1