CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-PAS-HIN-17-007
Constraints on nuclear parton distributions from W boson production in pPb collisions at ${\sqrt {s_{_{\mathrm {NN}}}}} = $ 8.16 TeV
Abstract: The production of W bosons has been measured in proton-lead collisions at a nucleon-nucleon center-of-mass energy of ${\sqrt {s_{_{\mathrm {NN}}}}} = $ 8.16 TeV. The measurement of W bosons is performed in the $\mathrm{W}\rightarrow\mu\nu_{\mu}$ channel, using a data sample with integrated luminosity of 173.4 $\pm$ 8.7 nb$^{-1}$ collected by the CMS experiment at the CERN LHC. The number of positive and negative W bosons is determined in the muon pseudorapidity region $|\eta_{\mathrm{lab}}| < $ 2.4 and transverse momentum $p_{\mathrm{T}}^{\mu} > $ 25 GeV/$c$. The W boson differential cross sections, muon charge asymmetry, and the ratios of W boson yields on the p-going over the Pb-going directions are reported as a function of the muon pseudorapidity in the center-of-mass frame. The measurements are compared to the predictions from theoretical calculations based on parton distribution functions (PDF). The results of this analysis favor PDF calculations including nuclear modifications and provide constraints on the nuclear PDF global fits.
Figures & Tables Summary Additional Figures References CMS Publications
Figures

png pdf
Figure 1:
The missing transverse momentum $ {{p_{\mathrm {T}}} ^\text {miss}} $ distribution for $ {{\mathrm {W^+}}\to {\mu ^+} {\nu _ {\mu}}} $ events within the $-0.2 < {\eta ^{\mu}_{\textrm {CM}}} < 0.0$ (left) and $ 1.8 < {\eta ^{\mu}_{\textrm {CM}}} < 1.93$ (right) ranges. Unbinned fits to the data (black points) are performed with six contributions, stacked from bottom to top: ${{\mathrm {t}\overline {\mathrm {t}}}}$ (orange), ${{\mathrm {Z}/\gamma *} \to {{\tau} {\overline {\tau}}}}$ (dark blue), ${{\mathrm {W^+}}\to {\overline {\tau}} {\nu _ {\tau}}}$ (red), ${{\mathrm {Z}/\gamma *} \to {\mu ^+\mu ^-}}$ (green), QCD multijet (light blue) and $ {{\mathrm {W^+}}\to {\mu ^+} {\nu _ {\mu}}} $ (yellow). The $ {\eta ^{\mu}_{\textrm {CM}}} $ regions are defined such that the proton is moving towards positive pseudorapidity. Error bars represent statistical uncertainties. The lower panels display the data divided by the result of the fit, for each $ {{p_{\mathrm {T}}} ^\text {miss}} $ bin.

png pdf
Figure 1-a:
The missing transverse momentum $ {{p_{\mathrm {T}}} ^\text {miss}} $ distribution for $ {{\mathrm {W^+}}\to {\mu ^+} {\nu _ {\mu}}} $ events within the $-0.2 < {\eta ^{\mu}_{\textrm {CM}}} < 0.0$ range. Unbinned fits to the data (black points) are performed with six contributions, stacked from bottom to top: ${{\mathrm {t}\overline {\mathrm {t}}}}$ (orange), ${{\mathrm {Z}/\gamma *} \to {{\tau} {\overline {\tau}}}}$ (dark blue), ${{\mathrm {W^+}}\to {\overline {\tau}} {\nu _ {\tau}}}$ (red), ${{\mathrm {Z}/\gamma *} \to {\mu ^+\mu ^-}}$ (green), QCD multijet (light blue) and $ {{\mathrm {W^+}}\to {\mu ^+} {\nu _ {\mu}}} $ (yellow). The $ {\eta ^{\mu}_{\textrm {CM}}} $ regions are defined such that the proton is moving towards positive pseudorapidity. Error bars represent statistical uncertainties. The lower panel displays the data divided by the result of the fit, for each $ {{p_{\mathrm {T}}} ^\text {miss}} $ bin.

png pdf
Figure 1-b:
The missing transverse momentum $ {{p_{\mathrm {T}}} ^\text {miss}} $ distribution for $ {{\mathrm {W^+}}\to {\mu ^+} {\nu _ {\mu}}} $ events within the $ 1.8 < {\eta ^{\mu}_{\textrm {CM}}} < 1.93$ range. Unbinned fits to the data (black points) are performed with six contributions, stacked from bottom to top: ${{\mathrm {t}\overline {\mathrm {t}}}}$ (orange), ${{\mathrm {Z}/\gamma *} \to {{\tau} {\overline {\tau}}}}$ (dark blue), ${{\mathrm {W^+}}\to {\overline {\tau}} {\nu _ {\tau}}}$ (red), ${{\mathrm {Z}/\gamma *} \to {\mu ^+\mu ^-}}$ (green), QCD multijet (light blue) and $ {{\mathrm {W^+}}\to {\mu ^+} {\nu _ {\mu}}} $ (yellow). The $ {\eta ^{\mu}_{\textrm {CM}}} $ regions are defined such that the proton is moving towards positive pseudorapidity. Error bars represent statistical uncertainties. The lower panel displays the data divided by the result of the fit, for each $ {{p_{\mathrm {T}}} ^\text {miss}} $ bin.

png pdf
Figure 2:
Production cross sections for $ {{\mathrm {W^+}}\to {\mu ^+} {\nu _ {\mu}}} $ (left) and $ {{\mathrm {W^-}}\to {\mu ^-} {\overline {\nu}_ {\mu}}} $ (right), as a function of the muon pseudorapidity in the center-of-mass frame. The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. The global luminosity uncertainty of $\pm $5.0% is not shown. Calculations using with CT14 PDF (red line), CT14+EPPS16 nPDF (green line) and CT14+nCTEQ15 nPDF (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval. The bottom panels show the ratio of data (black points), CT14+EPPS16 (green line) and CT14+nCTEQ15 (brown line) with respect to CT14.

png pdf
Figure 2-a:
Production cross sections for $ {{\mathrm {W^+}}\to {\mu ^+} {\nu _ {\mu}}} $, as a function of the muon pseudorapidity in the center-of-mass frame. The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. The global luminosity uncertainty of $\pm $5.0% is not shown. Calculations using with CT14 PDF (red line), CT14+EPPS16 nPDF (green line) and CT14+nCTEQ15 nPDF (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval. The bottom panel shows the ratio of data (black points), CT14+EPPS16 (green line) and CT14+nCTEQ15 (brown line) with respect to CT14.

png pdf
Figure 2-b:
Production cross sections for $ {{\mathrm {W^-}}\to {\mu ^-} {\overline {\nu}_ {\mu}}} $, as a function of the muon pseudorapidity in the center-of-mass frame. The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. The global luminosity uncertainty of $\pm $5.0% is not shown. Calculations using with CT14 PDF (red line), CT14+EPPS16 nPDF (green line) and CT14+nCTEQ15 nPDF (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval. The bottom panel shows the ratio of data (black points), CT14+EPPS16 (green line) and CT14+nCTEQ15 (brown line) with respect to CT14.

png pdf
Figure 3:
Forward-backward ratios, $N_{{{\mu}}}(+ {\eta ^{\mu}_{\textrm {CM}}})/N_{{{\mu}}}(- {\eta ^{\mu}_{\textrm {CM}}})$, for the positive (left) and negative (right) muons. The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. Calculations using with CT14 PDF (red line), CT14+EPPS16 nPDF (green line) and CT14+nCTEQ15 nPDF (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval.

png pdf
Figure 3-a:
Forward-backward ratios, $N_{{{\mu}}}(+ {\eta ^{\mu}_{\textrm {CM}}})/N_{{{\mu}}}(- {\eta ^{\mu}_{\textrm {CM}}})$, for the positive muons. The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. Calculations using with CT14 PDF (red line), CT14+EPPS16 nPDF (green line) and CT14+nCTEQ15 nPDF (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval.

png pdf
Figure 3-b:
Forward-backward ratios, $N_{{{\mu}}}(+ {\eta ^{\mu}_{\textrm {CM}}})/N_{{{\mu}}}(- {\eta ^{\mu}_{\textrm {CM}}})$, for the negative muons. The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. Calculations using with CT14 PDF (red line), CT14+EPPS16 nPDF (green line) and CT14+nCTEQ15 nPDF (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval.

png pdf
Figure 4:
The forward-backward ratio of all muons, $N_{{\mu}}(+ {\eta ^{\mu}_{\textrm {CM}}})/N_{{\mu}}(- {\eta ^{\mu}_{\textrm {CM}}})$, as a function of $ {\eta ^{\mu}_{\textrm {CM}}} $. The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. Calculations using with CT14 PDF (red line), CT14+EPPS16 nPDF (green line) and CT14+nCTEQ15 nPDF (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval.

png pdf
Figure 5:
Muon charge asymmetry, $(N_{{{\mu}}}^{+} - N_{{{\mu}}}^{-})/(N_{{{\mu}}}^{+} + N_{{{\mu}}}^{-})$, as a function of the muon pseudorapidity in the center-of-mass frame. The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. Calculations using with CT14 PDF (red line), CT14+EPPS16 nPDF (green line) and CT14+nCTEQ15 nPDF (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval.

png pdf
Figure 6:
Comparison of the muon charge asymmetry measured at 8.16 TeV (black points) and at 5.02 TeV [12] (blue squares). The muon pseudorapidity has been shifted according to ${{\eta ^{\mu}_{\textrm {ref}}} = {\eta ^{\mu}_{\textrm {CM}}} \pm \textrm {ln} (8.16~\mathrm{TeV} / {\sqrt {s_{_{\text {NN}}}}} )}$. The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. Calculations using with CT14+EPPS16 nPDF at 8.16 TeV (green line) and at 5.02 TeV (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval.

png pdf
Figure 7:
Distribution of the $\chi ^2$/ndf values from the comparison of data (absolute cross sections) and theoretical calculations, for the CT14, nCTEQ15 and EPPS16 individual sets. The vertical dashed lines represent the prediction corresponding to the central set of CT14, nCTEQ15, and EPPS16.
Tables

png pdf
Table 1:
Maximum uncertainty of the measured observables as a function of $ {\eta ^{\mu}_{\textrm {CM}}} $ determined for each category. The uncertainties of the cross sections are relative while for the asymmetries are absolute. The global luminosity uncertainty of $\pm $5.0% is not included in the total systematic uncertainty of the cross sections.

png pdf
Table 2:
Results of the $\chi ^{2}$ statistical test between the measurements and the theory calculations from the CT14 PDF, CT14+EPPS16 nPDF and CT14+nCTEQ15 nPDF calculations. The value of the $\chi ^{2}$, the number of degrees of freedom (ndf) and the $\chi ^{2}$ probability (Prob.), are presented for the W differential cross sections, the muon charge asymmetry, the charged muon forward-backward ratios, and the forward-backward ratio of all muons, respectively.
Summary
A measurement of W boson production in pPb collisions at ${\sqrt {s_{_{\mathrm {NN}}}}} = $ 8.16 TeV is reported, using the muon decay channel for muons with $ {p_{\mathrm{T}}} $ greater than 25 GeV/$c$ and $|{{\eta_{\textrm{lab}}} }| < $ 2.4. The absolute differential production cross sections for positive and negative $ \mathrm{W} \to\mu\nu_{\mu} $ decays, the muon charge asymmetry and the muon forward-backward ratios, are measured as a function of the muon pseudorapidity in the center-of-mass frame, in the range ${\eta^{\mu}_{\textrm{CM}}} \in [-2.86, 1.93]$.

The measurements are compared to theoretical predictions assuming both proton PDF (CT14) and nuclear PDF (CT14+EPPS16, CT14+nCTEQ15) sets. The absolute cross sections and the forward-backward asymmetries exhibit significant deviations from the CT14 case, revealing nuclear modifications of the PDFs unambiguously for the first time in the production of weak bosons in nuclear collisions. On the contrary, both the CT14+EPPS16 and the CT14+nCTEQ15 calculations show a good overall agreement with data, though the data favors the former global fit. In the latter case, nevertheless, only the individual sets which exhibit the smallest nuclear PDF modifications at small values of $x$ (in the shadowing region) turn out to be compatible with experimental measurements. The small experimental uncertainties reported in the measurements will enable a significant reduction of the current uncertainties of the quark and antiquark nPDFs, in the range $10^{-3}$-$10^{-1}$.
Additional Figures

png pdf
Additional Figure 1:
Comparison of the production cross section for $\mathrm {W}^{-}\rightarrow \mu ^{-}\bar{\nu}_{\mu}$ (left) and $\mathrm {W}^{+}\rightarrow \mu ^{+}\nu _{\mu}$ (right) measured at 8.16 TeV (black points) and at 5.02 TeV [Phys. Lett. B 750 (2015) 565] (blue squares). The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. The global luminosity uncertainty of $\pm $5.0% for the data at 8.16 TeV and of $\pm $3.5% for the data at 5.02 TeV are not shown. Calculations using with CT14+EPPS16 nPDF [Eur. Phys. J. C 77, 163 (2017)] at 8.16 TeV (green line) and at 5.02 TeV (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval.

png pdf
Additional Figure 1-a:
Comparison of the production cross section for $\mathrm {W}^{-}\rightarrow \mu ^{-}\bar{\nu}_{\mu}$ measured at 8.16 TeV (black points) and at 5.02 TeV [Phys. Lett. B 750 (2015) 565] (blue squares). The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. The global luminosity uncertainty of $\pm $5.0% for the data at 8.16 TeV and of $\pm $3.5% for the data at 5.02 TeV are not shown. Calculations using with CT14+EPPS16 nPDF [Eur. Phys. J. C 77, 163 (2017)] at 8.16 TeV (green line) and at 5.02 TeV (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval.

png pdf
Additional Figure 1-b:
Comparison of the production cross section for $\mathrm {W}^{+}\rightarrow \mu ^{+}\nu _{\mu}$ measured at 8.16 TeV (black points) and at 5.02 TeV [Phys. Lett. B 750 (2015) 565] (blue squares). The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. The global luminosity uncertainty of $\pm $5.0% for the data at 8.16 TeV and of $\pm $3.5% for the data at 5.02 TeV are not shown. Calculations using with CT14+EPPS16 nPDF [Eur. Phys. J. C 77, 163 (2017)] at 8.16 TeV (green line) and at 5.02 TeV (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval.

png pdf
Additional Figure 2:
Comparison of the muon charge asymmetry measured at 8.16 TeV (black points) and at 5.02 TeV [Phys. Lett. B 750 (2015) 565] (blue squares). The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. Calculations using with CT14+EPPS16 nPDF [Eur. Phys. J. C 77, 163 (2017)] at 8.16 TeV (green line) and at 5.02 TeV (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval.

png pdf
Additional Figure 3:
Comparison of the production cross section for $\mathrm {W}^{-}\rightarrow \mu ^{-}\bar{\nu}_{\mu}$ measured at 8.16 TeV (black points) and at 5.02 TeV [Phys. Lett. B 750 (2015) 565] (blue squares). The comparison is presented in two muon pseudorapidity regions: $\eta _{\mathrm {CM}} < $ 0 (left plot) and $\eta _{\mathrm {CM}} > $ 0 (right plot). The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. The global luminosity uncertainty of $\pm $5.0% for the data at 8.16 TeV and of $\pm $3.5% for the data at 5.02 TeV are not shown.

png pdf
Additional Figure 3-a:
Comparison of the production cross section for $\mathrm {W}^{-}\rightarrow \mu ^{-}\bar{\nu}_{\mu}$ measured at 8.16 TeV (black points) and at 5.02 TeV [Phys. Lett. B 750 (2015) 565] (blue squares). The comparison is presented in pseudorapidity regions $\eta _{\mathrm {CM}} < $ 0. The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. The global luminosity uncertainty of $\pm $5.0% for the data at 8.16 TeV and of $\pm $3.5% for the data at 5.02 TeV are not shown.

png pdf
Additional Figure 3-b:
Comparison of the production cross section for $\mathrm {W}^{-}\rightarrow \mu ^{-}\bar{\nu}_{\mu}$ measured at 8.16 TeV (black points) and at 5.02 TeV [Phys. Lett. B 750 (2015) 565] (blue squares). The comparison is presented in pseudorapidity regions $\eta _{\mathrm {CM}} > $ 0. The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. The global luminosity uncertainty of $\pm $5.0% for the data at 8.16 TeV and of $\pm $3.5% for the data at 5.02 TeV are not shown.

png pdf
Additional Figure 4:
Comparison of the production cross section for $\mathrm {W}^{+}\rightarrow \mu ^{+}\nu _{\mu}$ measured at 8.16 TeV (black points) and at 5.02 TeV [Phys. Lett. B 750 (2015) 565] (blue squares). The comparison is presented in two muon pseudorapidity regions: $\eta _{\mathrm {CM}} < $ 0 (left plot) and $\eta _{\mathrm {CM}} > $ 0 (right plot). The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. The global luminosity uncertainty of $\pm $5.0% for the data at 8.16 TeV and of $\pm $3.5% for the data at 5.02 TeV are not shown.

png pdf
Additional Figure 4-a:
Comparison of the production cross section for $\mathrm {W}^{+}\rightarrow \mu ^{+}\nu _{\mu}$ measured at 8.16 TeV (black points) and at 5.02 TeV [Phys. Lett. B 750 (2015) 565] (blue squares). The comparison is presented in muon pseudorapidity region $\eta _{\mathrm {CM}} < $ 0. The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. The global luminosity uncertainty of $\pm $5.0% for the data at 8.16 TeV and of $\pm $3.5% for the data at 5.02 TeV are not shown.

png pdf
Additional Figure 4-b:
Comparison of the production cross section for $\mathrm {W}^{+}\rightarrow \mu ^{+}\nu _{\mu}$ measured at 8.16 TeV (black points) and at 5.02 TeV [Phys. Lett. B 750 (2015) 565] (blue squares). The comparison is presented in muon pseudorapidity region $\eta _{\mathrm {CM}} > $ 0. The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. The global luminosity uncertainty of $\pm $5.0% for the data at 8.16 TeV and of $\pm $3.5% for the data at 5.02 TeV are not shown.

png pdf
Additional Figure 5:
Comparison of the muon charge asymmetry measured at 8.16 TeV (black points) and at 5.02 TeV [Phys. Lett. B 750 (2015) 565] (blue squares). The muon pseudorapidity has been shifted according to ${{\eta ^{\mu}_{\textrm {ref}}} = {\eta ^{\mu}_{\textrm {CM}}} \pm \textrm {ln}\left (8.16~\text{TeV}/ {\sqrt {s_{_{\text {NN}}}}} \right)}$. The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only.

png pdf
Additional Figure 6:
Comparison of the production cross section for $\mathrm {W}^{-}\rightarrow \mu ^{-}\bar{\nu}_{\mu}$ measured at 8.16 TeV (black points) and at 5.02 TeV [Phys. Lett. B 750 (2015) 565] (blue squares). The comparison is presented in two muon pseudorapidity regions: $\eta _{\mathrm {CM}} < $ 0 (left plot) and $\eta _{\mathrm {CM}} > $ 0 (right plot). The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. The global luminosity uncertainty of $\pm $5.0% for the data at 8.16 TeV and of $\pm $3.5% for the data at 5.02 TeV are not shown. Calculations using with CT14+EPPS16 nPDF [Eur. Phys. J. C 77, 163 (2017)] at 8.16 TeV (green line) and at 5.02 TeV (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval.

png pdf
Additional Figure 6-a:
Comparison of the production cross section for $\mathrm {W}^{-}\rightarrow \mu ^{-}\bar{\nu}_{\mu}$ measured at 8.16 TeV (black points) and at 5.02 TeV [Phys. Lett. B 750 (2015) 565] (blue squares). The comparison is presented in muon pseudorapidity region $\eta _{\mathrm {CM}} < $ 0. The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. The global luminosity uncertainty of $\pm $5.0% for the data at 8.16 TeV and of $\pm $3.5% for the data at 5.02 TeV are not shown. Calculations using with CT14+EPPS16 nPDF [Eur. Phys. J. C 77, 163 (2017)] at 8.16 TeV (green line) and at 5.02 TeV (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval.

png pdf
Additional Figure 6-b:
Comparison of the production cross section for $\mathrm {W}^{-}\rightarrow \mu ^{-}\bar{\nu}_{\mu}$ measured at 8.16 TeV (black points) and at 5.02 TeV [Phys. Lett. B 750 (2015) 565] (blue squares). The comparison is presented in muon pseudorapidity region $\eta _{\mathrm {CM}} > $ 0. The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. The global luminosity uncertainty of $\pm $5.0% for the data at 8.16 TeV and of $\pm $3.5% for the data at 5.02 TeV are not shown. Calculations using with CT14+EPPS16 nPDF [Eur. Phys. J. C 77, 163 (2017)] at 8.16 TeV (green line) and at 5.02 TeV (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval.

png pdf
Additional Figure 7:
Comparison of the production cross section for $\mathrm {W}^{+}\rightarrow \mu ^{+}\nu _{\mu}$ measured at 8.16 TeV (black points) and at 5.02 TeV [Phys. Lett. B 750 (2015) 565] (blue squares). The comparison is presented in two muon pseudorapidity regions: $\eta _{\mathrm {CM}} < $ 0 (left plot) and $\eta _{\mathrm {CM}} > $ 0 (right plot). The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. The global luminosity uncertainty of $\pm $5.0% for the data at 8.16 TeV and of $\pm $3.5% for the data at 5.02 TeV are not shown. Calculations using with CT14+EPPS16 nPDF [Eur. Phys. J. C 77, 163 (2017)] at 8.16 TeV (green line) and at 5.02 TeV (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval.

png pdf
Additional Figure 7-a:
Comparison of the production cross section for $\mathrm {W}^{+}\rightarrow \mu ^{+}\nu _{\mu}$ measured at 8.16 TeV (black points) and at 5.02 TeV [Phys. Lett. B 750 (2015) 565] (blue squares). The comparison is presented in muon pseudorapidity region $\eta _{\mathrm {CM}} < $ 0. The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. The global luminosity uncertainty of $\pm $5.0% for the data at 8.16 TeV and of $\pm $3.5% for the data at 5.02 TeV are not shown. Calculations using with CT14+EPPS16 nPDF [Eur. Phys. J. C 77, 163 (2017)] at 8.16 TeV (green line) and at 5.02 TeV (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval.

png pdf
Additional Figure 7-b:
Comparison of the production cross section for $\mathrm {W}^{+}\rightarrow \mu ^{+}\nu _{\mu}$ measured at 8.16 TeV (black points) and at 5.02 TeV [Phys. Lett. B 750 (2015) 565] (blue squares). The comparison is presented in muon pseudorapidity region $\eta _{\mathrm {CM}} > $ 0. The brackets represent the statistical and systematic uncertainties summed in quadrature, while the error bars show the statistical uncertainties only. The global luminosity uncertainty of $\pm $5.0% for the data at 8.16 TeV and of $\pm $3.5% for the data at 5.02 TeV are not shown. Calculations using with CT14+EPPS16 nPDF [Eur. Phys. J. C 77, 163 (2017)] at 8.16 TeV (green line) and at 5.02 TeV (brown line), are also displayed, including their PDF uncertainty bands at 68% confidence interval.

png pdf
Additional Figure 8:
Distribution of the $\chi ^2$/ndf values from the comparison of data (muon charge asymmetry) and theoretical calculations, for the CT14 [Phys. Rev. D 93, 033006 (2016)], nCTEQ15 [Phys. Rev. D 93, 085037 (2016)] and EPPS16 [Eur. Phys. J. C 77, 163 (2017)] individual sets. The vertical dashed lines represent the prediction corresponding to the central set of CT14, nCTEQ15, and EPPS16.

png pdf
Additional Figure 9:
Distribution of the $\chi ^2$/ndf values from the comparison of data (muon forward-backward ratio) and theoretical calculations, for the CT14 [Phys. Rev. D 93, 033006 (2016)], nCTEQ15 [Phys. Rev. D 93, 085037 (2016)] and EPPS16 [Eur. Phys. J. C 77, 163 (2017)] individual sets. The vertical dashed lines represent the prediction corresponding to the central set of CT14, nCTEQ15, and EPPS16.

png pdf
Additional Figure 10:
Distribution of the $\chi ^2$/ndf values from the comparison of data (forward-backward ratio of all muons) and theoretical calculations, for the CT14 [Phys. Rev. D 93, 033006 (2016)], nCTEQ15 [Phys. Rev. D 93, 085037 (2016)] and EPPS16 [Eur. Phys. J. C 77, 163 (2017)] individual sets. The vertical dashed lines represent the prediction corresponding to the central set of CT14, nCTEQ15, and EPPS16.

png pdf
Additional Figure 11:
Distribution of the simulated Bjorken x in the Pb ion as a function of the muon pseudorapidity for $\mathrm {W}^{-}\rightarrow \mu ^{-}\bar{\nu}_{\mu}$ (left) and $\mathrm {W}^{+}\rightarrow \mu ^{+}\nu _{\mu}$ (right) generated at 8.16 TeV. The W boson production have been simulated using the NLO generator POWHEG v2 with the POWHEG BOX W.ew.BMNNP using the CT14+EPPS16 nPDF [Eur. Phys. J. C 77, 163 (2017)].

png pdf
Additional Figure 11-a:
Distribution of the simulated Bjorken x in the Pb ion as a function of the muon pseudorapidity for $\mathrm {W}^{-}\rightarrow \mu ^{-}\bar{\nu}_{\mu}$ generated at 8.16 TeV. The W boson production have been simulated using the NLO generator POWHEG v2 with the POWHEG BOX W.ew.BMNNP using the CT14+EPPS16 nPDF [Eur. Phys. J. C 77, 163 (2017)].

png pdf
Additional Figure 11-b:
Distribution of the simulated Bjorken x in the Pb ion as a function of the muon pseudorapidity for $\mathrm {W}^{+}\rightarrow \mu ^{+}\nu _{\mu}$ generated at 8.16 TeV. The W boson production have been simulated using the NLO generator POWHEG v2 with the POWHEG BOX W.ew.BMNNP using the CT14+EPPS16 nPDF [Eur. Phys. J. C 77, 163 (2017)].

png pdf
Additional Figure 12:
Correlation matrix of the W boson production cross section measurements. The statistical and systematic uncertainties have been included except for the luminosity uncertainty. The black lines are used to distinguish the different bins of muon charge which are ordered from top to bottom as: Minus-Minus, Minus-Plus, Plus-Minus and Plus-Plus.
References
1 J. Butterworth et al. PDF4LHC recommendations for LHC Run II Journal of Physics G: Nuclear and Particle Physics 43 (2016) 023001 1510.03865
2 S. Dulat et al. New parton distribution functions from a global analysis of quantum chromodynamics PRD 93 (2016) 033006 1506.07443
3 NNPDF Collaboration Parton distributions for the LHC Run II JHEP 04 (2015) 040 1410.8849
4 L. A. Harland-Lang, A. D. Martin, P. Motylinski, and R. S. Thorne Parton distributions in the LHC era: MMHT 2014 PDFs EPJC 75 (2015) 204 1412.3989
5 R. Vogt Shadowing effects on vector boson production PRC 64 (2001) 044901 hep-ph/0011242
6 X.-F. Zhang and G. I. Fai $ Z^0 $ production as a test of nuclear effects at the LHC PLB 545 (2002) 91 hep-ph/0205155
7 H. Paukkunen and C. A. Salgado Constraints for the nuclear parton distributions from Z and W production at the LHC JHEP 03 (2011) 071 1010.5392
8 N. Armesto Nuclear shadowing JPG 32 (2006) R367 hep-ph/0604108
9 CDF Collaboration Direct Measurement of the $ W $ Production Charge Asymmetry in $ {\rm p}\overline{\rm p} $ Collisions at $ \sqrt{s} = $ 1.96 TeV PRL 102 (2009) 181801 0901.2169
10 CMS Collaboration Measurement of the muon charge asymmetry in inclusive pp~$ \to W+X $ production at $ \sqrt s = $ 7 TeV and an improved determination of light parton distribution functions PRD 90 (2014) 032004 CMS-SMP-12-021
1312.6283
11 CMS Collaboration Measurement of the electron charge asymmetry in inclusive $ W $ production in pp collisions at $ \sqrt{s}= $ 7 ~TeV PRL 109 (2012) 111806 CMS-SMP-12-001
1206.2598
12 CMS Collaboration Study of W boson production in pPb collisions at $ \sqrt{s_{\mathrm{NN}}} = $ 5.02 TeV PLB 750 (2015) 565 CMS-HIN-13-007
1503.05825
13 Z. Conesa del Valle et al. Effect of heavy-quark energy loss on the muon differential production cross-section in Pb-Pb collisions at $ \sqrt{s_{\rm NN}} = $ 5.5 ~TeV PLB 663 (2008) 202 0712.0051
14 Z. Conesa del Valle Vector bosons in heavy-ion collisions at the LHC EPJC 61 (2009) 729 0903.1432
15 ATLAS Collaboration Measurement of the centrality dependence of J/$ \psi $ yields and observation of Z production in lead-lead collisions with the ATLAS detector at the LHC PLB 697 (2011) 294 1012.5419
16 ATLAS Collaboration Measurement of $ Z $ boson Production in PbPb Collisions at $ \sqrt{s_{\rm NN}}= $ 2.76 TeV with the ATLAS Detector PRL 110 (2013) 022301 1210.6486
17 ATLAS Collaboration Measurement of the production and lepton charge asymmetry of $ W $ bosons in Pb+Pb collisions at $ \sqrt{{s}_{\mathrm {NN}}}= $ 2.76 TeV with the ATLAS detector EPJC 75 (2015) 23 1408.4674
18 CMS Collaboration Study of Z boson production in PbPb collisions at nucleon-nucleon centre of mass energy = 2.76 TeV PRL 106 (2011) 212301 CMS-HIN-10-003
1102.5435
19 CMS Collaboration Study of Z production in PbPb and pp collisions at $ \sqrt{s_{\mathrm{NN}}}= $ 2.76 TeV in the dimuon and dielectron decay channels JHEP 03 (2015) 022 CMS-HIN-13-004
1410.4825
20 CMS Collaboration Study of $ W $ boson production in PbPb and $ pp $ collisions at $ \sqrt{s_{NN}}= $ 2.76 TeV PLB 715 (2012) 66 CMS-HIN-11-008
1205.6334
21 ATLAS Collaboration Measurement of $ W\rightarrow\mu\nu $ production in $ p $+Pb collision at $ \sqrt{s_{_\text{NN}}}= $ 5.02 TeV with ATLAS detector at the LHC ATLAS-CONF-2015-056
22 ALICE Collaboration W and Z boson production in p-Pb collisions at $ \sqrt{s_{\rm NN}} = $ 5.02 TeV JHEP 02 (2017) 077 1611.03002
23 CMS Collaboration Study of Z boson production in pPb collisions at $ \sqrt {s_{NN}} = $ 5.02 TeV PLB 759 (2016) 36 CMS-HIN-15-002
1512.06461
24 K. J. Eskola, P. Paakkinen, H. Paukkunen, and C. A. Salgado EPPS16: Nuclear parton distributions with LHC data EPJC 77 (2017) 163 1612.05741
25 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004 CMS-00-001
26 E. Todesco and J. Wenninger Large Hadron Collider momentum calibration and accuracy PRAccel. Beams 20 (2017) 081003
27 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
28 CMS Collaboration Performance of CMS muon reconstruction in $ pp $ collision events at $ \sqrt{s} = $ 7 TeV JINST 7 (2012) P10002 CMS-MUO-10-004
1206.4071
29 CMS Collaboration Observation of top quark production in proton-nucleus collisions PRL 119 (2017) 242001 CMS-HIN-17-002
1709.07411
30 CMS Collaboration Particle-flow reconstruction and global event description with the CMS detector JINST 12 (2017) P10003 CMS-PRF-14-001
1706.04965
31 S. Frixione, P. Nason, and C. Oleari Matching NLO QCD computations with parton shower simulations: the POWHEG method JHEP 11 (2007) 070 0709.2092
32 P. Nason A new method for combining NLO QCD with shower Monte Carlo algorithms JHEP 11 (2004) 040 hep-ph/0409146
33 S. Alioli, P. Nason, C. Oleari, and E. Re A general framework for implementing NLO calculations in shower Monte Carlo programs: the POWHEG BOX JHEP 06 (2010) 043 1002.2581
34 L. Barze et al. Implementation of electroweak corrections in the POWHEG BOX: single W production JHEP 04 (2012) 037 1202.0465
35 L. Barze et al. Neutral-current Drell--Yan with combined QCD and electroweak corrections in the POWHEG BOX EPJC 73 (2013) 2474 1302.4606
36 S. Frixione, P. Nason, and G. Ridolfi A positive-weight next-to-leading-order Monte Carlo for heavy flavour hadroproduction JHEP 09 (2007) 126 0707.3088
37 T. Sjostrand et al. An Introduction to PYTHIA 8.2 CPC 191 (2015) 159 1410.3012
38 CMS Collaboration Event generator tunes obtained from underlying event and multiparton scattering measurements EPJC 76 (2016) 155 CMS-GEN-14-001
1512.00815
39 CMS Collaboration Study of the underlying event at forward rapidity in pp collisions at $ \sqrt{s} = $ 0.9, 2.76, and 7 TeV JHEP 04 (2013) 072 CMS-FWD-11-003
1302.2394
40 T. Pierog et al. EPOS LHC: Test of collective hadronization with data measured at the CERN Large Hadron Collider PRC 92 (2015) 034906 1306.0121
41 CMS Collaboration Pseudorapidity distributions of charged hadrons in proton-lead collisions at $ \sqrt{s_{_\mathrm{NN}}} = $ 5.02 and 8.16 TeV JHEP 01 (2018) 045 CMS-HIN-16-021
1710.09355
42 GEANT4 Collaboration GEANT4: A Simulation toolkit NIMA 506 (2003) 250
43 CMS Collaboration Missing transverse energy performance of the CMS detector JINST 6 (2011) P09001 CMS-JME-10-009
1106.5048
44 CMS Collaboration Measurements of Inclusive $ W $ and $ Z $ Cross Sections in $ pp $ Collisions at $ \sqrt{s}= $ 7 TeV JHEP 01 (2011) 080 CMS-EWK-10-002
1012.2466
45 Particle Data Group Collaboration Review of Particle Physics CPC 40 (2016) 100001
46 S. Alioli, P. Nason, C. Oleari, and E. Re NLO vector-boson production matched with shower in POWHEG JHEP 07 (2008) 060 0805.4802
47 R. Boughezal et al. Color singlet production at NNLO in MCFM EPJC 77 (2017) 7 1605.08011
48 K. Kova\ifmmode \checkr\else \vr\fi\'\ik et al. nCTEQ15: Global analysis of nuclear parton distributions with uncertainties in the CTEQ framework PRD 93 (Apr, 2016) 085037 1509.00792
49 F. Arleo, \'E. Chapon, and H. Paukkunen Scaling properties of inclusive W$ ^\pm $ production at hadron colliders EPJC 76 (2016) 214 1509.03993
50 J. Pumplin et al. Uncertainties of predictions from parton distribution functions. 2. The Hessian method PRD 65 (2001) 014013 hep-ph/0101032
Compact Muon Solenoid
LHC, CERN