...
and if I print out some of the 'bogus' hits from CalXtalRecData, I see
...
If I did this correctly all these baddies are range 1, and I assume maxed out. This may explain the deficit of range 2 crystals. And maybe some effect for getting range 3.
If I exclude hits with energies > 2 GeV, the energy sums look normal.
The plots and printout were done with:
...