Figure 4

$\textup{d}\sigma_{\rm prompt\ \JJpsi}/\textup{d}y$ as a function of the rapidity in the centre-of-mass frame (left panel) as obtained in this work at midrapidity and by the LHCb collaboration in the forward and backward rapidity regions , and $\textup{d}^2\sigma_{\rm prompt\ \JJpsi}/\textup{d}y\textup{d}p_{\rm T}$ as a function of \pt (right panel) compared with ATLAS results  (reported up to $\pt = 20$ GeV/$c$). The results for inclusive \JJpsi are also shown superimposed in the right hand panel. The vertical error bars and the boxes on top of each point represent the statistical and systematic uncertainties. For the measurements as a function of \pt, the data symbols are placed within each bin at the mean of the \pt distribution determined from MC simulations. The results of a model  including nuclear shadowing based on the EPPS16  and nCTEQ15  nPDFs are shown superimposed on both panels (see text for details). In the right panel, the computations refer only to the ALICE rapidity range.