...
The visible deviations of fit from the spectrum is explained by low energy tails increasing with theta angle, as confirmed by simulation:
Using peak position for theta < 10 deg as normalization I found the deviation for linearity as a function of crystal energy (not corrected for angle):
...
The high energy tail of GCR peak is explained by GCR energy spectrum and rather well reproduced by simulation: