Comparisons between ALICE $\jpsi$ and $\psi(2S)$ data and summed NRQCD and FONLL model calculations from [36] and [37]. The first five panels correspond to inclusive $\jpsi$ production cross sections as a function of $y$ in pp collisions at $\sqrts=13$, $8$ and $7$, $5.02$ and $2.76$ TeV (red), while the next three panels to inclusive $\psi(2S)$ cross sections as a function of $y$ at $\sqrts=13$, $8$ and $7$ TeV (blue).