The transform function problem is fixed.
Plot of the simulated shower results using the original transform function:
Plot of the simulated shower results using LDF6 and a new transform function that allows theta out to 60 degrees:
There were two problems:
1. Of the 3 transform functions included in Brant's original code, only
one had the property of an even cutoff in angle. Using one of the others
gives a limit which changes with phi.
2. Even the correct transform function was giving me odd results, and I believe
this is because it got compiled with incorrect limits on x and y.
The following unweildy PDF has some plots: Transform Function