Nuclear modification factor for quarkonium production in high energy proton-nucleus collisions within the parton saturation approach

In this work, we analyze nuclear effects in the inclusive production of quarkonium in proton-nucleus collisions at high energy regime. A theoretical framework that includes initial state effects like nuclear shadowing and gluon density saturation is considered. Numerical results for nuclear modification factor, $R_{pA}(y)$, as a function of meson rapidity in proton-nucleus collisions are presented. The parameter-free predictions are compared to the available data from the Large Hadron Collider. Discussion on the main theoretical uncertainties is made, with emphasis on the phenomenological models for the nuclear saturation scale.