$\textup{d}\sigma_{\rm \PJpsi\ from\ h_{\rm b}}/\textup{d}y$ as a function of rapidity (left panel) as obtained at midrapidity in this work for $\pt>0$ and by the LHCb collaboration in the forward and backward rapidity regions~ for $0<\pt<14$~GeV/$c$, and $\textup{d}^2\sigma_{\rm \PJpsi\ from\ h_{\rm b}}/\textup{d}y\textup{d}p_{\rm T}$ as a function of \pt (right panel) compared with ATLAS measurements~ (shown up to $\pt = 20$ GeV/$c$). The statistical and systematic uncertainties are shown as vertical error bars and hollow boxes, respectively. The extrapolation uncertainty related to the procedure to go from the visible region ($\pt > 1$~GeV/$c$) to $\pt > 0$ for the ALICE data point in the left panel is indicated with a dashed band For the measurements as a function of \pt, data symbols are placed at the mean value of the \pt distribution within each bin The results are compared to FONLL computations~ with EPPS16~ nPDFs, highlighting the total theoretical uncertainty (empty band) and the contribution from EPPS16 (coloured band). In the right panel, model computations are obtained in the same rapidity range of the ALICE results, namely $-1.37