Best-fit model (curves) and unfolded IceCube energy spectrum (points with error bars) for showers (upper panel) and tracks (lower panel) using our method described in Appendix A, where the best fit for tracks in Table 2 has been used. Black dots (and error bars) refer to HESE data and pink/light dots (and error bars) to TGM, which contains the 4.5 PeV reconstructed energy track. The energy scale refers to reconstructed neutrino energy (data) and incident energy (model). The different contributions (and the total) for the multi-component model are shown, as indicated in the plot legends. Compared to the IceCube analysis, the residual atmospheric background is shown explicitly both in data and model.