Update:Plots all include a cut on theta <= 45.
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:

