We introduce a global analysis of collinearly factorized nuclear parton distribution functions (PDFs) including, for the first time, data constraints from LHC proton–lead collisions. In comparison to our previous analysis, EPS09, where data only from charged-lepton–nucleus deep inelastic scattering (DIS), Drell–Yan (DY) dilepton production in proton–nucleus collisions and inclusive pion production in deuteron–nucleus collisions were the input, we now increase the variety of data constraints to cover also neutrino–nucleus DIS and low-mass DY production in pion–nucleus collisions. The new LHC data significantly extend the kinematic reach of the data constraints. We now allow much more freedom for the flavor dependence of nuclear effects than in other currently available analyses. As a result, especially the uncertainty estimates are more objective flavor by flavor. The neutrino DIS plays a pivotal role in obtaining a mutually consistent behavior for both up and down valence quarks, and the LHC dijet data clearly constrain gluons at large momentum fraction. Mainly for insufficient statistics, the pion–nucleus DY and heavy-gauge-boson production in proton–lead collisions impose less visible constraints. The outcome – a new set of next-to-leading order nuclear PDFs called EPPS16 – is made available for applications in high-energy nuclear collisions.

We scrutinize the recent LHCb data for D 0-meson production in p+Pb collisions within a next-to-leading order QCD framework. Our calculations are performed in the SACOT-m T variant of the general-mass variable-flavour-number scheme (GM-VFNS), which has previously been shown to provide a realistic description of the LHC p+p data. Using the EPPS16 and nCTEQ15 nuclear parton distribution functions (PDFs) we show that a very good agreement is obtained also in the p+Pb case both for cross sections and nuclear modification ratios in the wide rapidity range covered by the LHCb data. Encouraged by the good correspondence, we quantify the impact of these data on the nuclear PDFs by the Hessian reweighting technique. We find compelling direct evidence of gluon shadowing at small momentum fractions x, with no signs of parton dynamics beyond the collinear factorization. We also compare our theoretical framework to a fixed-order calculation supplemented with a parton shower. While the two frameworks differ in the absolute cross sections, these differences largely cancel in the nuclear modification ratios. Thus, the constraints for nuclear PDFs appear solid.