I ran the simShowers program to produce 18000 simulated showers. (1000 per interval of 0.2 in Log(E), from 10^16 to 10^19.6 eV.)
Output of simShowers2: sim1000.txt
Summary of input parameters: inputdata1000.txt
Fit results: fits1000.dat
I used the readshowers macro to graph the simulated data:





Individual parameters:












Distribution of reconstructed energies for each input energy bin:

















