Error on integral method + solar mod not used for primary flux

[GallegoSav]

Hello,

The quad integrate method from scipy on CreateBackgroundSpectrumMEGAlib.py lead to several errors. Using the method trapz from numpy solved the issue.
I just had to replace the line

LEOBackground/CreateBackgroundSpectrumMEGAlib.py

Line 64 in b8898ce

IntSpectrum, err = quad(Megalibfunc[i], 10**Elow, 10**Ehigh)

per IntSpectrum = np.trapz(Megalibfunc[i](Energies),Energies)

Also for the computation of the flux for the primary particles, the solar mod factor is computed but not implemented :

LEOBackground/LEOBackgroundGenerator.py

Line 502 in b8898ce

""" Solar modulation factor from Gleeson & Axford 1968"""

the output should look like that : f(E)*redfac*solmodfac ?