Helper function to get number of observed events as a function of reco energy given a flux model

[chiarabellenghi]

Helper function that people might want to use.

I think it should allow both single declination and declination bands. It should calculate the events per season. Then one can just sum the output histograms to get the total number of events over multiple seasons.