Nuclear modification factor of inclusive $\Upsilon$(1S) in p-Pb collisions at $\snn=5.02$ TeV as a function of rapidity. The vertical error bars represent the statistical uncertainties and the open boxes the uncorrelated systematic uncertainties. The full boxes around $R_{\rm pPb}=1$ show the size of the correlated uncertainties. Also shown are several model calculations:
(top) parton energy loss with and without EPS09 shadowing at NLO and CEM with EPS09 shadowing at NLO;
(bottom) CGC based and CSM with EPS09 shadowing at LO. For the latter the effect of variation in the shadowing and EMC curves is highlighted as described in the text.