...
I plan to do it and hope it will improve data/MC agreement and will allow to use big angle protons for calibration.
GCR peak shape depends on momentum of incident ion
p = 3.5 GeV/ nuclon
p=5.6 GeV/nuclon
p=11 GeV/nuclon
p=18 GeV/nuclon
Peaks are perfectly fitted with lognormal function, but its parameters significantly depend on ion momentum.
Explaination: while for p=3.5 GeV/nucl the maximum energy of secondary electron is 12 MeV, for p=18 GeV/nucl this energy becomes 320 MeV
and high energy tail increases.
To get proper peak shape for fitting data we have to find peak shape versus momentum and then fold it with GCR spectrum.