CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-HIN-18-003 ; CERN-EP-2021-028
Study of Drell-Yan dimuon production in proton-lead collisions at ${\sqrt {\smash [b]{s_{_{\mathrm {NN}}}}}} = $ 8.16 TeV
JHEP 05 (2021) 182
Abstract: Differential cross sections for the Drell-Yan process, including Z boson production, using the dimuon decay channel are measured in proton-lead (pPb) collisions at a nucleon-nucleon centre-of-mass energy of 8.16 TeV. A data sample recorded with the CMS detector at the LHC is used, corresponding to an integrated luminosity of 173 nb$^{-1}$. The differential cross section as a function of the dimuon mass is measured in the range 15-600 GeV, for the first time in proton-nucleus collisions. It is also reported as a function of dimuon rapidity over the mass ranges 15-60 GeV and 60-120 GeV, and ratios for the p-going over the Pb-going beam directions are built. In both mass ranges, the differential cross sections as functions of the dimuon transverse momentum ${p_{\mathrm{T}}}$ and of a geometric variable ${\phi^*}$ are measured, where ${\phi^*}$ highly correlates with ${p_{\mathrm{T}}}$ but is determined with higher precision. In the Z mass region, the rapidity dependence of the data indicate a modification of the distribution of partons within a lead nucleus as compared to the proton case. The data are more precise than predictions based upon current models of parton distributions.
Figures & Tables Summary References CMS Publications
Figures

png pdf
Figure 1:
Comparison of the data (black points) with the ${\mathrm{Z} /\gamma ^{*}}$ signal and background expectations (filled histograms, where "EW" includes $ {\mathrm{Z} /\gamma ^{*}} \to \tau^{+} \tau^{-} $ and diboson), estimated as described in the text, as a function of invariant mass (upper) and rapidity in the centre-of-mass frame for 15 $ < {m_{\mu \mu}} < $ 60 GeV (lower left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (lower right). Vertical error bars represent statistical uncertainties. The ratios of data over expectations are shown in the lower panels. The boson ${p_{\mathrm {T}}}$ reweighting described in the text is not applied. The shaded regions show the quadratic sum of the systematic uncertainties (including the integrated luminosity, but excluding acceptance and unfolding uncertainties) and the nPDF uncertainties (CT14+EPPS16).

png pdf
Figure 1-a:
Comparison of the data (black points) with the ${\mathrm{Z} /\gamma ^{*}}$ signal and background expectations (filled histograms, where "EW" includes $ {\mathrm{Z} /\gamma ^{*}} \to \tau^{+} \tau^{-} $ and diboson), estimated as described in the text, as a function of invariant mass (upper) and rapidity in the centre-of-mass frame for 15 $ < {m_{\mu \mu}} < $ 60 GeV (lower left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (lower right). Vertical error bars represent statistical uncertainties. The ratios of data over expectations are shown in the lower panels. The boson ${p_{\mathrm {T}}}$ reweighting described in the text is not applied. The shaded regions show the quadratic sum of the systematic uncertainties (including the integrated luminosity, but excluding acceptance and unfolding uncertainties) and the nPDF uncertainties (CT14+EPPS16).

png pdf
Figure 1-b:
Comparison of the data (black points) with the ${\mathrm{Z} /\gamma ^{*}}$ signal and background expectations (filled histograms, where "EW" includes $ {\mathrm{Z} /\gamma ^{*}} \to \tau^{+} \tau^{-} $ and diboson), estimated as described in the text, as a function of invariant mass (upper) and rapidity in the centre-of-mass frame for 15 $ < {m_{\mu \mu}} < $ 60 GeV (lower left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (lower right). Vertical error bars represent statistical uncertainties. The ratios of data over expectations are shown in the lower panels. The boson ${p_{\mathrm {T}}}$ reweighting described in the text is not applied. The shaded regions show the quadratic sum of the systematic uncertainties (including the integrated luminosity, but excluding acceptance and unfolding uncertainties) and the nPDF uncertainties (CT14+EPPS16).

png pdf
Figure 1-c:
Comparison of the data (black points) with the ${\mathrm{Z} /\gamma ^{*}}$ signal and background expectations (filled histograms, where "EW" includes $ {\mathrm{Z} /\gamma ^{*}} \to \tau^{+} \tau^{-} $ and diboson), estimated as described in the text, as a function of invariant mass (upper) and rapidity in the centre-of-mass frame for 15 $ < {m_{\mu \mu}} < $ 60 GeV (lower left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (lower right). Vertical error bars represent statistical uncertainties. The ratios of data over expectations are shown in the lower panels. The boson ${p_{\mathrm {T}}}$ reweighting described in the text is not applied. The shaded regions show the quadratic sum of the systematic uncertainties (including the integrated luminosity, but excluding acceptance and unfolding uncertainties) and the nPDF uncertainties (CT14+EPPS16).

png pdf
Figure 2:
Comparison of the data (black points) with the ${\mathrm{Z} /\gamma ^{*}}$ signal and background expectations (filled histograms, where "EW" includes $ {\mathrm{Z} /\gamma ^{*}} \to \tau^{+} \tau^{-} $ and diboson), estimated as described in the text, as a function of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right). The first bins of the ${p_{\mathrm {T}}}$ and ${\phi ^*}$ distributions start at 0. Vertical error bars represent statistical uncertainties. The ratios of data over expectations are shown in the lower panels. The boson ${p_{\mathrm {T}}}$ reweighting described in the text is not applied. The shaded regions show the quadratic sum of the systematic uncertainties (including the integrated luminosity, but excluding acceptance and unfolding uncertainties) and the nPDF uncertainties (CT14+EPPS16).

png pdf
Figure 2-a:
Comparison of the data (black points) with the ${\mathrm{Z} /\gamma ^{*}}$ signal and background expectations (filled histograms, where "EW" includes $ {\mathrm{Z} /\gamma ^{*}} \to \tau^{+} \tau^{-} $ and diboson), estimated as described in the text, as a function of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right). The first bins of the ${p_{\mathrm {T}}}$ and ${\phi ^*}$ distributions start at 0. Vertical error bars represent statistical uncertainties. The ratios of data over expectations are shown in the lower panels. The boson ${p_{\mathrm {T}}}$ reweighting described in the text is not applied. The shaded regions show the quadratic sum of the systematic uncertainties (including the integrated luminosity, but excluding acceptance and unfolding uncertainties) and the nPDF uncertainties (CT14+EPPS16).

png pdf
Figure 2-b:
Comparison of the data (black points) with the ${\mathrm{Z} /\gamma ^{*}}$ signal and background expectations (filled histograms, where "EW" includes $ {\mathrm{Z} /\gamma ^{*}} \to \tau^{+} \tau^{-} $ and diboson), estimated as described in the text, as a function of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right). The first bins of the ${p_{\mathrm {T}}}$ and ${\phi ^*}$ distributions start at 0. Vertical error bars represent statistical uncertainties. The ratios of data over expectations are shown in the lower panels. The boson ${p_{\mathrm {T}}}$ reweighting described in the text is not applied. The shaded regions show the quadratic sum of the systematic uncertainties (including the integrated luminosity, but excluding acceptance and unfolding uncertainties) and the nPDF uncertainties (CT14+EPPS16).

png pdf
Figure 2-c:
Comparison of the data (black points) with the ${\mathrm{Z} /\gamma ^{*}}$ signal and background expectations (filled histograms, where "EW" includes $ {\mathrm{Z} /\gamma ^{*}} \to \tau^{+} \tau^{-} $ and diboson), estimated as described in the text, as a function of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right). The first bins of the ${p_{\mathrm {T}}}$ and ${\phi ^*}$ distributions start at 0. Vertical error bars represent statistical uncertainties. The ratios of data over expectations are shown in the lower panels. The boson ${p_{\mathrm {T}}}$ reweighting described in the text is not applied. The shaded regions show the quadratic sum of the systematic uncertainties (including the integrated luminosity, but excluding acceptance and unfolding uncertainties) and the nPDF uncertainties (CT14+EPPS16).

png pdf
Figure 2-d:
Comparison of the data (black points) with the ${\mathrm{Z} /\gamma ^{*}}$ signal and background expectations (filled histograms, where "EW" includes $ {\mathrm{Z} /\gamma ^{*}} \to \tau^{+} \tau^{-} $ and diboson), estimated as described in the text, as a function of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right). The first bins of the ${p_{\mathrm {T}}}$ and ${\phi ^*}$ distributions start at 0. Vertical error bars represent statistical uncertainties. The ratios of data over expectations are shown in the lower panels. The boson ${p_{\mathrm {T}}}$ reweighting described in the text is not applied. The shaded regions show the quadratic sum of the systematic uncertainties (including the integrated luminosity, but excluding acceptance and unfolding uncertainties) and the nPDF uncertainties (CT14+EPPS16).

png pdf
Figure 3:
Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of the dimuon invariant mass (upper) and rapidity in the centre-of-mass frame for 15 $ < {m_{\mu \mu}} < $ 60 GeV (lower left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (lower right).

png pdf
Figure 3-a:
Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of the dimuon invariant mass (upper) and rapidity in the centre-of-mass frame for 15 $ < {m_{\mu \mu}} < $ 60 GeV (lower left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (lower right).

png pdf
Figure 3-b:
Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of the dimuon invariant mass (upper) and rapidity in the centre-of-mass frame for 15 $ < {m_{\mu \mu}} < $ 60 GeV (lower left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (lower right).

png pdf
Figure 3-c:
Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of the dimuon invariant mass (upper) and rapidity in the centre-of-mass frame for 15 $ < {m_{\mu \mu}} < $ 60 GeV (lower left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (lower right).

png pdf
Figure 4:
Correlation matrices for the systematic uncertainties, excluding integrated luminosity, as functions of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right).

png pdf
Figure 4-a:
Correlation matrices for the systematic uncertainties, excluding integrated luminosity, as functions of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right).

png pdf
Figure 4-b:
Correlation matrices for the systematic uncertainties, excluding integrated luminosity, as functions of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right).

png pdf
Figure 4-c:
Correlation matrices for the systematic uncertainties, excluding integrated luminosity, as functions of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right).

png pdf
Figure 4-d:
Correlation matrices for the systematic uncertainties, excluding integrated luminosity, as functions of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right).

png pdf
Figure 5:
Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of the dimuon invariant mass (upper) and rapidity in the centre-of-mass frame for 15 $ < {m_{\mu \mu}} < $ 60 GeV (lower left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (lower right). The error bars on the data represent the quadratic sum of the statistical and systematic uncertainties. Theory predictions from the POWHEG NLO generator are also shown, using CT14 (blue) or CT14+EPPS16 (red). The boxes show the 68% confidence level (n)PDF uncertainty on these predictions. The ratios of predictions over data are shown in the lower panels, where the data and (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 5-a:
Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of the dimuon invariant mass (upper) and rapidity in the centre-of-mass frame for 15 $ < {m_{\mu \mu}} < $ 60 GeV (lower left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (lower right). The error bars on the data represent the quadratic sum of the statistical and systematic uncertainties. Theory predictions from the POWHEG NLO generator are also shown, using CT14 (blue) or CT14+EPPS16 (red). The boxes show the 68% confidence level (n)PDF uncertainty on these predictions. The ratios of predictions over data are shown in the lower panels, where the data and (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 5-b:
Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of the dimuon invariant mass (upper) and rapidity in the centre-of-mass frame for 15 $ < {m_{\mu \mu}} < $ 60 GeV (lower left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (lower right). The error bars on the data represent the quadratic sum of the statistical and systematic uncertainties. Theory predictions from the POWHEG NLO generator are also shown, using CT14 (blue) or CT14+EPPS16 (red). The boxes show the 68% confidence level (n)PDF uncertainty on these predictions. The ratios of predictions over data are shown in the lower panels, where the data and (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 5-c:
Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of the dimuon invariant mass (upper) and rapidity in the centre-of-mass frame for 15 $ < {m_{\mu \mu}} < $ 60 GeV (lower left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (lower right). The error bars on the data represent the quadratic sum of the statistical and systematic uncertainties. Theory predictions from the POWHEG NLO generator are also shown, using CT14 (blue) or CT14+EPPS16 (red). The boxes show the 68% confidence level (n)PDF uncertainty on these predictions. The ratios of predictions over data are shown in the lower panels, where the data and (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 6:
Differential fiducial cross sections (without the acceptance correction) for the DY process measured in the muon channel, as functions of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right). The first bin of the ${p_{\mathrm {T}}}$ and ${\phi ^*}$ measurements starts at 0. The error bars on the data represent the quadratic sum of the statistical and systematic uncertainties. Theory predictions from the POWHEG NLO generator are also shown, using CT14 (blue) or CT14+EPPS16 (red). The boxes show the 68% confidence level (n)PDF uncertainty on these predictions. The ratios of predictions over data are shown in the lower panels, where the data and (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 6-a:
Differential fiducial cross sections (without the acceptance correction) for the DY process measured in the muon channel, as functions of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right). The first bin of the ${p_{\mathrm {T}}}$ and ${\phi ^*}$ measurements starts at 0. The error bars on the data represent the quadratic sum of the statistical and systematic uncertainties. Theory predictions from the POWHEG NLO generator are also shown, using CT14 (blue) or CT14+EPPS16 (red). The boxes show the 68% confidence level (n)PDF uncertainty on these predictions. The ratios of predictions over data are shown in the lower panels, where the data and (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 6-b:
Differential fiducial cross sections (without the acceptance correction) for the DY process measured in the muon channel, as functions of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right). The first bin of the ${p_{\mathrm {T}}}$ and ${\phi ^*}$ measurements starts at 0. The error bars on the data represent the quadratic sum of the statistical and systematic uncertainties. Theory predictions from the POWHEG NLO generator are also shown, using CT14 (blue) or CT14+EPPS16 (red). The boxes show the 68% confidence level (n)PDF uncertainty on these predictions. The ratios of predictions over data are shown in the lower panels, where the data and (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 6-c:
Differential fiducial cross sections (without the acceptance correction) for the DY process measured in the muon channel, as functions of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right). The first bin of the ${p_{\mathrm {T}}}$ and ${\phi ^*}$ measurements starts at 0. The error bars on the data represent the quadratic sum of the statistical and systematic uncertainties. Theory predictions from the POWHEG NLO generator are also shown, using CT14 (blue) or CT14+EPPS16 (red). The boxes show the 68% confidence level (n)PDF uncertainty on these predictions. The ratios of predictions over data are shown in the lower panels, where the data and (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 6-d:
Differential fiducial cross sections (without the acceptance correction) for the DY process measured in the muon channel, as functions of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right). The first bin of the ${p_{\mathrm {T}}}$ and ${\phi ^*}$ measurements starts at 0. The error bars on the data represent the quadratic sum of the statistical and systematic uncertainties. Theory predictions from the POWHEG NLO generator are also shown, using CT14 (blue) or CT14+EPPS16 (red). The boxes show the 68% confidence level (n)PDF uncertainty on these predictions. The ratios of predictions over data are shown in the lower panels, where the data and (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 7:
Differential cross section for the DY process measured in the muon channel, as a function of the dimuon invariant mass (upper) and rapidity in the centre-of-mass frame for 15 $ < {m_{\mu \mu}} < $ 60 GeV (lower left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (lower right). The error bars on the data represent the quadratic sum of the statistical and systematic uncertainties. Theory predictions from the POWHEG NLO generator are also shown, using CT14 (blue) or CT14+EPPS16 (red). The boxes show the 68% confidence level (n)PDF uncertainty on these predictions. The ratios of predictions over data are shown in the lower panels, where the data and (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 7-a:
Differential cross section for the DY process measured in the muon channel, as a function of the dimuon invariant mass (upper) and rapidity in the centre-of-mass frame for 15 $ < {m_{\mu \mu}} < $ 60 GeV (lower left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (lower right). The error bars on the data represent the quadratic sum of the statistical and systematic uncertainties. Theory predictions from the POWHEG NLO generator are also shown, using CT14 (blue) or CT14+EPPS16 (red). The boxes show the 68% confidence level (n)PDF uncertainty on these predictions. The ratios of predictions over data are shown in the lower panels, where the data and (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 7-b:
Differential cross section for the DY process measured in the muon channel, as a function of the dimuon invariant mass (upper) and rapidity in the centre-of-mass frame for 15 $ < {m_{\mu \mu}} < $ 60 GeV (lower left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (lower right). The error bars on the data represent the quadratic sum of the statistical and systematic uncertainties. Theory predictions from the POWHEG NLO generator are also shown, using CT14 (blue) or CT14+EPPS16 (red). The boxes show the 68% confidence level (n)PDF uncertainty on these predictions. The ratios of predictions over data are shown in the lower panels, where the data and (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 7-c:
Differential cross section for the DY process measured in the muon channel, as a function of the dimuon invariant mass (upper) and rapidity in the centre-of-mass frame for 15 $ < {m_{\mu \mu}} < $ 60 GeV (lower left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (lower right). The error bars on the data represent the quadratic sum of the statistical and systematic uncertainties. Theory predictions from the POWHEG NLO generator are also shown, using CT14 (blue) or CT14+EPPS16 (red). The boxes show the 68% confidence level (n)PDF uncertainty on these predictions. The ratios of predictions over data are shown in the lower panels, where the data and (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 8:
Differential cross sections for the DY process measured in the muon channel, as functions of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right). The first bin of the ${p_{\mathrm {T}}}$ and ${\phi ^*}$ measurements starts at 0. The error bars on the data represent the quadratic sum of the statistical and systematic uncertainties. Theory predictions from the POWHEG NLO generator are also shown, using CT14 (blue) or CT14+EPPS16 (red). The boxes show the 68% confidence level (n)PDF uncertainty on these predictions. The ratios of predictions over data are shown in the lower panels, where the data and (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 8-a:
Differential cross sections for the DY process measured in the muon channel, as functions of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right). The first bin of the ${p_{\mathrm {T}}}$ and ${\phi ^*}$ measurements starts at 0. The error bars on the data represent the quadratic sum of the statistical and systematic uncertainties. Theory predictions from the POWHEG NLO generator are also shown, using CT14 (blue) or CT14+EPPS16 (red). The boxes show the 68% confidence level (n)PDF uncertainty on these predictions. The ratios of predictions over data are shown in the lower panels, where the data and (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 8-b:
Differential cross sections for the DY process measured in the muon channel, as functions of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right). The first bin of the ${p_{\mathrm {T}}}$ and ${\phi ^*}$ measurements starts at 0. The error bars on the data represent the quadratic sum of the statistical and systematic uncertainties. Theory predictions from the POWHEG NLO generator are also shown, using CT14 (blue) or CT14+EPPS16 (red). The boxes show the 68% confidence level (n)PDF uncertainty on these predictions. The ratios of predictions over data are shown in the lower panels, where the data and (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 8-c:
Differential cross sections for the DY process measured in the muon channel, as functions of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right). The first bin of the ${p_{\mathrm {T}}}$ and ${\phi ^*}$ measurements starts at 0. The error bars on the data represent the quadratic sum of the statistical and systematic uncertainties. Theory predictions from the POWHEG NLO generator are also shown, using CT14 (blue) or CT14+EPPS16 (red). The boxes show the 68% confidence level (n)PDF uncertainty on these predictions. The ratios of predictions over data are shown in the lower panels, where the data and (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 8-d:
Differential cross sections for the DY process measured in the muon channel, as functions of ${p_{\mathrm {T}}}$ (upper row) and ${\phi ^*}$ (lower row), for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right). The first bin of the ${p_{\mathrm {T}}}$ and ${\phi ^*}$ measurements starts at 0. The error bars on the data represent the quadratic sum of the statistical and systematic uncertainties. Theory predictions from the POWHEG NLO generator are also shown, using CT14 (blue) or CT14+EPPS16 (red). The boxes show the 68% confidence level (n)PDF uncertainty on these predictions. The ratios of predictions over data are shown in the lower panels, where the data and (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 9:
Forward-backward ratios for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right). The error bars on the data points represent the quadratic sum of the statistical and systematic uncertainties. The theory predictions from the POWHEG NLO generator are also shown, using CT14 [51] (blue), CT14+EPPS16 [14] (red), or CT14+nCTEQ15WZ [19] (green) PDF sets. The boxes show the 68% confidence level (n)PDF uncertainty in these predictions. The ratios of predictions over data are shown in the lower panels, where the data and the (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 9-a:
Forward-backward ratios for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right). The error bars on the data points represent the quadratic sum of the statistical and systematic uncertainties. The theory predictions from the POWHEG NLO generator are also shown, using CT14 [51] (blue), CT14+EPPS16 [14] (red), or CT14+nCTEQ15WZ [19] (green) PDF sets. The boxes show the 68% confidence level (n)PDF uncertainty in these predictions. The ratios of predictions over data are shown in the lower panels, where the data and the (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.

png pdf
Figure 9-b:
Forward-backward ratios for 15 $ < {m_{\mu \mu}} < $ 60 GeV (left) and 60 $ < {m_{\mu \mu}} < $ 120 GeV (right). The error bars on the data points represent the quadratic sum of the statistical and systematic uncertainties. The theory predictions from the POWHEG NLO generator are also shown, using CT14 [51] (blue), CT14+EPPS16 [14] (red), or CT14+nCTEQ15WZ [19] (green) PDF sets. The boxes show the 68% confidence level (n)PDF uncertainty in these predictions. The ratios of predictions over data are shown in the lower panels, where the data and the (n)PDF uncertainties are shown separately, as error bars around one and as coloured boxes, respectively.
Tables

png pdf
Table 1:
Range of systematic uncertainties in percentage of the cross section, given separately for 15 $ < {m_{\mu \mu}} < $ 60 and 60 $ < {m_{\mu \mu}} < $ 120 GeV. Systematic uncertainties for the three mass bins above 120 GeV fall in the range given for 15 $ < {m_{\mu \mu}} < $ 60 GeV. For the theoretical component of acceptance and efficiency, the systematic uncertainty related to efficiency alone (for fiducial cross sections) is given between parentheses.

png pdf
Table 2:
$\chi ^2$ values between the data and the POWHEG predictions and associated probability, from the fiducial cross sections, when experimental and theoretical bin-to-bin correlations are taken into account. The integrated luminosity uncertainty is included in the experimental uncertainties.

png pdf
Table 3:
$\chi ^2$ values between the data and the POWHEG predictions and associated probability, from the full phase space cross sections, when experimental and theoretical bin-to-bin correlations are taken into account. The integrated luminosity uncertainty is included in the experimental uncertainties.
Summary
Differential cross section measurements of the Drell-Yan process in the dimuon channel in proton-lead collisions at ${\sqrt {\smash [b]{s_{_{\mathrm {NN}}}}}} = $ 8.16 TeV have been reported, including the transverse momentum (${p_{\mathrm{T}}}$) and rapidity dependencies in the Z boson mass region (60 $ < m_{\mu\mu} < $ 120 GeV). In addition, for the first time in collisions including nuclei, the ${p_{\mathrm{T}}}$ and rapidity dependence for smaller masses 15 $ < m_{\mu\mu} < $ 60 GeV have been measured. The dependence with ${\phi^*}$ (a geometrical variable that highly correlates with dimuon ${p_{\mathrm{T}}}$ but is determined with higher precision) for both 15 $ < m_{\mu\mu} < $ 60 GeV and 60 $ < m_{\mu\mu} < $ 120 GeV and the mass dependence from 15 to 600 GeV have been presented, also for the first time in proton-nucleus collisions. Finally, forward-backward ratios have been built from the rapidity-dependent cross sections for ${y_{\mathrm{CM}}} > $ 0 to ${y_{\mathrm{CM}}} < $ 0 in both mass regions, highlighting the presence of nuclear effects in the parton distribution functions.

Results for 60 $ < m_{\mu\mu} < $ 120 GeV are the most precise to date, featuring smaller uncertainties than the theoretical predictions, and provide novel constraints on the quark and antiquark nuclear parton distribution functions (nPDFs). Measurements in the lower mass range 15 $ < m_{\mu\mu} < $ 60 GeV give access to a new phase space for nPDF studies, extending to lower longitudinal momentum fraction $x$ and lower energy scale $Q^2$. The ${p_{\mathrm{T}}}$- and ${\phi^*}$ -dependent results are also very sensitive to the details of model details, such as soft quantum chromodynamics phenomena, which they may help to better understand in pPb collisions.
References
1 S. D. Drell and T.-M. Yan Massive lepton pair production in hadron-hadron collisions at high-energies PRL 25 (1970) 316, . [Erratum: \DOI10.1103/PhysRevLett.25.902.2]
2 R. Hamberg, W. L. van Neerven, and T. Matsuura A complete calculation of the order $ \alpha^2_s $ correction to the Drell--Yan K-factor NPB 359 (1991) 343, . [Erratum: \DOI10.1016/S0550-3213(02)00814-3]
3 S. Catani et al. Vector boson production at hadron colliders: a fully exclusive QCD calculation at next-to-next-to-leading order PRL 103 (2009) 082001 0903.2120
4 S. Catani and M. Grazzini Next-to-next-to-leading-order subtraction formalism in hadron collisions and its application to Higgs-boson production at the Large Hadron Collider PRL 98 (2007) 222002 hep-ph/0703012
5 K. Melnikov and F. Petriello Electroweak gauge boson production at hadron colliders through O($ \alpha_s^2 $) PRD 74 (2006) 114017 hep-ph/0609070
6 A. D. Martin, W. J. Stirling, and R. G. Roberts Parton distributions of the proton PRD 50 (1994) 6734 hep-ph/9406315
7 V. Kartvelishvili, R. Kvatadze, and R. Shanidze On Z and Z + jet production in heavy ion collisions PLB 356 (1995) 589 hep-ph/9505418
8 R. Vogt Shadowing effects on vector boson production PRC 64 (2001) 044901 hep-ph/0011242
9 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
10 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
11 European Muon Collaboration The ratio of the nucleon structure functions $ F_2^N $ for iron and deuterium PLB 123 (1983) 275
12 D. de Florian, R. Sassot, P. Zurita, and M. Stratmann Global analysis of nuclear parton distributions PRD 85 (2012) 074028 1112.6324
13 H. Khanpour and S. Atashbar Tehrani Global analysis of nuclear parton distribution functions and their uncertainties at next-to-next-to-leading order PRD 93 (2016) 014026 1601.00939
14 K. J. Eskola, P. Paakkinen, H. Paukkunen, and C. A. Salgado EPPS16: nuclear parton distributions with LHC data EPJC 77 (2017) 163 1612.05741
15 K. Kova\vr\'\ik et al. nCTEQ15 - Global analysis of nuclear parton distributions with uncertainties in the CTEQ framework PRD 93 (2016) 085037 1509.00792
16 M. Walt, I. Helenius, and W. Vogelsang Open-source QCD analysis of nuclear parton distribution functions at NLO and NNLO PRD 100 (2019) 096015 1908.03355
17 NNPDF Collaboration Nuclear parton distributions from lepton-nucleus scattering and the impact of an electron-ion collider EPJC 79 (2019) 471 1904.00018
18 R. Abdul Khalek, J. J. Ethier, J. Rojo, and G. van Weelden nNNPDF2.0: quark flavor separation in nuclei from LHC data JHEP 09 (2020) 183 2006.14629
19 A. Kusina et al. Impact of LHC vector boson production in heavy ion collisions on strange PDFs EPJC 80 (2020) 968 2007.09100
20 N. Armesto Nuclear shadowing JPG 32 (2006) R367 hep-ph/0604108
21 ATLAS Collaboration Measurement of the high-mass Drell--Yan differential cross-section in pp collisions at $ \sqrt{s}= $ 7 TeV with the ATLAS detector PLB 725 (2013) 223 1305.4192
22 ATLAS Collaboration Measurement of the low-mass Drell--Yan differential cross section at $ \sqrt{s} = $ 7 TeV using the ATLAS detector JHEP 06 (2014) 112 1404.1212
23 ATLAS Collaboration Measurement of the double-differential high-mass Drell--Yan cross section in pp collisions at $ \sqrt{s}= $ 8 TeV with the ATLAS detector JHEP 08 (2016) 009 1606.01736
24 ATLAS Collaboration Measurement of $ W^{\pm} $-boson and Z-boson production cross-sections in pp collisions at $ \sqrt{s}= $ 2.76 TeV with the ATLAS detector EPJC 79 (2019) 901 1907.03567
25 ATLAS Collaboration Measurement of the transverse momentum distribution of Drell--Yan lepton pairs in proton--proton collisions at $ \sqrt{s}= $ 13 TeV with the ATLAS detector EPJC 80 (2020) 616 1912.02844
26 CMS Collaboration Measurement of the Drell--Yan cross section in pp collisions at $ \sqrt{s}= $ 7 TeV JHEP 10 (2011) 007 CMS-EWK-10-007
1108.0566
27 CMS Collaboration Measurement of the differential and double-differential Drell--Yan cross sections in proton-proton collisions at $ \sqrt{s} = $ 7 TeV JHEP 12 (2013) 030 CMS-SMP-13-003
1310.7291
28 CMS Collaboration Measurements of differential and double-differential Drell--Yan cross sections in proton-proton collisions at 8 TeV EPJC 75 (2015) 147 CMS-SMP-14-003
1412.1115
29 CMS Collaboration Measurement of the differential Drell--Yan cross section in proton-proton collisions at $ \sqrt{\mathrm{s}} = $ 13 TeV JHEP 12 (2019) 059 CMS-SMP-17-001
1812.10529
30 CMS Collaboration Measurements of differential Z boson production cross sections in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JHEP 12 (2019) 061 CMS-SMP-17-010
1909.04133
31 PHENIX Collaboration Measurements of $ \mu\mu $ pairs from open heavy flavor and Drell--Yan in $ p+p $ collisions at $ \sqrt{s}= $ 200 GeV PRD 99 (2019) 072003 1805.02448
32 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
33 ALICE Collaboration Z-boson production in p-Pb collisions at $ \sqrt{s_{\mathrm{NN}}}= $ 8.16 TeV and Pb-Pb collisions at $ \sqrt{s_{\mathrm{NN}}}= $ 5.02 TeV JHEP 09 (2020) 076 2005.11126
34 ATLAS Collaboration $ Z $ boson production in $ p+ $Pb collisions at $ \sqrt{s_{NN}}= $ 5.02 TeV measured with the ATLAS detector PRC 92 (2015) 044915 1507.06232
35 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
36 A. Banfi et al. Optimisation of variables for studying dilepton transverse momentum distributions at hadron colliders EPJC 71 (2011) 1600 1009.1580
37 A. Banfi, M. Dasgupta, S. Marzani, and L. Tomlinson Predictions for Drell--Yan $ \phi^* $ and $ Q_T $ observables at the LHC PLB 715 (2012) 152 1205.4760
38 S. Marzani $ Q_T $ and $ \phi^* $ observables in Drell--Yan processes EPJWeb Conf. 49 (2013) 14007
39 CMS Collaboration CMS luminosity measurement using 2016 proton-nucleus collisions at nucleon-nucleon center-of-mass energy of 8.16 TeV CMS-PAS-LUM-17-002 CMS-PAS-LUM-17-002
40 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004 CMS-00-001
41 CMS Collaboration Performance of the CMS Level-1 trigger in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JINST 15 (2020) P10017 CMS-TRG-17-001
2006.10165
42 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
43 M. Cacciari, G. P. Salam, and G. Soyez The anti-$ {k_{\mathrm{T}}} $ jet clustering algorithm JHEP 04 (2008) 063 0802.1189
44 M. Cacciari, G. P. Salam, and G. Soyez FastJet user manual EPJC 72 (2012) 1896 1111.6097
45 CMS Collaboration Particle-flow reconstruction and global event description with the CMS detector JINST 12 (2017) P10003 CMS-PRF-14-001
1706.04965
46 CMS Collaboration Performance of the CMS muon detector and muon reconstruction with proton-proton collisions at $ \sqrt{s}= $ 13 TeV JINST 13 (2018) P06015 CMS-MUO-16-001
1804.04528
47 P. Nason A new method for combining NLO QCD with shower Monte Carlo algorithms JHEP 11 (2004) 040 hep-ph/0409146
48 S. Frixione, P. Nason, and C. Oleari Matching NLO QCD computations with parton shower simulations: the POWHEG method JHEP 11 (2007) 070 0709.2092
49 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
50 S. Alioli, P. Nason, C. Oleari, and E. Re NLO vector-boson production matched with shower in POWHEG JHEP 07 (2008) 060 0805.4802
51 S. Dulat et al. New parton distribution functions from a global analysis of quantum chromodynamics PRD 93 (2016) 033006 1506.07443
52 T. Sjostrand et al. An introduction to PYTHIA 8.2 CPC 191 (2015) 159 1410.3012
53 CMS Collaboration Event generator tunes obtained from underlying event and multiparton scattering measurements EPJC 76 (2016) 155 CMS-GEN-14-001
1512.00815
54 S. Jadach, J. H. K uhn, and Z. W\c as TAUOLA---a library of Monte Carlo programs to simulate decays of polarized $ \tau $ leptons CPC 64 (1990) 275
55 P. Golonka and Z. W\c as PHOTOS Monte Carlo: a precision tool for QED corrections in Z and W decays EPJC 45 (2006) 97 hep-ph/0506026
56 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
57 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
58 CMS Collaboration Centrality and pseudorapidity dependence of the transverse energy density in pPb collisions at $ \sqrt{s_\mathrm{NN}} = $ 5.02 TeV PRC 100 (2019) 024902 CMS-HIN-14-014
1810.05745
59 GEANT4 Collaboration GEANT4--a simulation toolkit NIMA 506 (2003) 250
60 CMS Collaboration Measurement of the $ {{\mathrm{W}}^{+}}\mathrm{W}^{-} $ cross section in pp collisions at $ \sqrt{s} = $ 8 TeV and limits on anomalous gauge couplings EPJC 76 (2016) 401 CMS-SMP-14-016
1507.03268
61 CMS Collaboration Measurement of the WZ production cross section in pp collisions at $ \sqrt{s} = $ 7 and 8 TeV and search for anomalous triple gauge couplings at $ \sqrt{s} = $ 8 TeV EPJC 77 (2017) 236 CMS-SMP-14-014
1609.05721
62 CMS Collaboration Measurement of the $ {\mathrm{p}}{\mathrm{p}} \to \mathrm{Z}\mathrm{Z} $ production cross section and constraints on anomalous triple gauge couplings in four-lepton final states at $ \sqrt s= $ 8 TeV PLB 740 (2015) 250 CMS-SMP-13-005
1406.0113
63 CMS Collaboration Observation of top quark production in proton-nucleus collisions PRL 119 (2017) 242001 CMS-HIN-17-002
1709.07411
64 CMS Collaboration Observation of nuclear modifications in W$ ^\pm $ boson production in pPb collisions at $ \sqrt{s_\mathrm{NN}} = $ 8.16 TeV PLB 800 (2020) 135048 CMS-HIN-17-007
1905.01486
65 D. Bourilkov Exploring the LHC landscape with dileptons 1609.08994
66 CMS Collaboration Performance of photon reconstruction and identification with the CMS detector in proton-proton collisions at $ \sqrt{s} = $ 8 TeV JINST 10 (2015) P08010 CMS-EGM-14-001
1502.02702
67 A. Bodek et al. Extracting muon momentum scale corrections for hadron collider experiments EPJC 72 (2012) 2194 1208.3710
68 G. D. Cowan Oxford Univ. Press, Oxford
69 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
70 J. Butterworth et al. PDF4LHC recommendations for LHC Run II JPG 43 (2016) 023001 1510.03865
71 A. Buckley et al. LHAPDF6: parton density access in the LHC precision era EPJC 75 (2015) 132 1412.7420
72 I. P. Lokhtin and A. M. Snigirev A model of jet quenching in ultrarelativistic heavy ion collisions and high-$ p_{\textrm{T}} $ hadron spectra at RHIC EPJC 45 (2006) 211 hep-ph/0506189
73 H.-L. Lai et al. New parton distributions for collider physics PRD 82 (2010) 074024 1007.2241
74 F. Arleo and S. Peign\'e Disentangling shadowing from coherent energy loss using the Drell--Yan process PRD 95 (2017) 011502 1512.01794
Compact Muon Solenoid
LHC, CERN