Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and a $\tilde{\chi}_1^\pm$ that decays to $W^\pm\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction to $q_i\bar{q}_j W\pm\tilde{\chi}_1^0$.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to bottom quarks ($\tilde{g}\to b\bar{b}\tilde{\chi}_1^0$). Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to top quarks ($\tilde{g}\to t\bar{t}\tilde{\chi}_1^0$). Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for light-flavor squark pair production, where the squarks decay to $q\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay and 1-fold degeneracy in the light-flavor squarks (corresponding to the inner set of curves in the limit plot). To get the theory cross section for other N-fold degeneracy assumptions (e.g. 8-fold for the outer curves in the limit plot), just multiply by N.

Exclusion limits at 95% CL for bottom squark pair production, where the squarks decay to $b\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay to $t\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay to $b\tilde{\chi}_1^\pm$ and the $\tilde{\chi}_1^0$ decay to $W^\pm\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay either to $b\tilde{\chi}_1^\pm\to bW^\pm\tilde{\chi}_1^0$ or to $t\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay to $c\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for the mono-$\phi$ model, in which a resonantly-produced colored scalar decays to a massive Dirac fermion and a quark. Signal cross sections are calculated at leading order in $\alpha_S$, assuming unity branching fraction for the given decay.

Cross section limits for $\mathrm{LQ}\to\mathrm{q}\nu$, where $q=u,\,d,\,s,\,\mathrm{or}\,c$. Limits are at the 95% confidence level. Theory cross sections are LO for vector LQ, and NLO for scalar LQ. Branching ratio is assumed to be 100% to $\mathrm{q}\nu$.

Cross section limits for $\mathrm{LQ}\to\mathrm{b}\nu$. Limits are at the 95% confidence level. Theory cross sections are LO for vector LQ, and NLO for scalar LQ. Branching ratio is assumed to be 100% to $\mathrm{b}\nu$.

Cross section limits for $\mathrm{LQ}\to\mathrm{t}\nu$. Limits are at the 95% confidence level. Theory cross sections are LO for vector LQ, and NLO for scalar LQ. Branching ratios are assumed to be $\mathcal{B}(\mathrm{LQ}\to\mathrm{t}\nu)=1-\beta$, and $\mathcal{B}(\mathrm{LQ}\to\mathrm{b}\tau)=\beta$.

Predictions and observations for monojet signal regions

Predictions and observations for signal regions with $250 \leq H_\mathrm{T} < 450$ GeV

Predictions and observations for signal regions with $450 \leq H_\mathrm{T} < 575$ GeV and $N_\mathrm{j}<7$

Predictions and observations for signal regions with $450 \leq H_\mathrm{T} < 575$ GeV and $N_\mathrm{j}\geq7$

Predictions and observations for signal regions with $575 \leq H_\mathrm{T} < 1200$ GeV and $N_\mathrm{j}^\mathrm{hi}<4$

Predictions and observations for signal regions with $575 \leq H_\mathrm{T} < 1200$ GeV and $4\leq N_\mathrm{j}^\mathrm{hi}<7$

Predictions and observations for signal regions with $575 \leq H_\mathrm{T} < 1200$ GeV and $N_\mathrm{j}\geq7$

Predictions and observations for signal regions with $1200 \leq H_\mathrm{T} < 1500$ GeV and $N_\mathrm{j}^\mathrm{hi}<4$

Predictions and observations for signal regions with $1200 \leq H_\mathrm{T} < 1500$ GeV and $4\leq N_\mathrm{j}^\mathrm{hi}<7$

Predictions and observations for signal regions with $1200 \leq H_\mathrm{T} < 1500$ GeV and $N_\mathrm{j}\geq7$

Predictions and observations for signal regions with $H_\mathrm{T} \geq 1500$ GeV and $N_\mathrm{j}<7$

Predictions and observations for signal regions with $H_\mathrm{T} \geq 1500$ GeV and $N_\mathrm{j}\geq7$

Covariance matrix for the 282 signal regions of the inclusive $M_\mathrm{T2}$ search

Correlation matrix for the 282 signal regions of the inclusive $M_\mathrm{T2}$ search

Bin number definitions for the $M_\mathrm{T2}$ covariance and correlation matrices

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 10$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 50$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 200$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct light squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 10$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark.

Exclusion limits at 95% CL for direct light squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 50$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark.

Exclusion limits at 95% CL for direct light squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 200$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino is long-lived with $c\tau_0 = 10$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino is long-lived with $c\tau_0 = 50$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino is long-lived with $c\tau_0 = 200$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

The maximum chargino mass excluded at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the gluino mass is fixed to 1900 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

The maximum chargino mass excluded at 95% CL for direct squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the squark mass is fixed to 900 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

The maximum chargino mass excluded at 95% CL for direct squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the squark mass is fixed to 1500 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, and the eight light squarks' masses are assumed to be degenerate. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the stop mass is fixed to 1000 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 5$ to 1000 cm while the gluino mass is fixed to 1600 GeV and the neutralino's mass is fixed to 1575 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 5$ to 1000 cm while the squark mass is fixed to 2000 GeV and the neutralino's mass is fixed to 1000 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, and the eight light squarks' masses are assumed to be degenerate.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 5$ to 1000 cm while the stop mass is fixed to 1100 GeV and the neutralino's mass is fixed to 1000 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Predictions and observations for 2016 disappearing track signal regions

Predictions and observations for 2017-2018 pixel track signal regions

Predictions and observations for 2017-2018 medium (M) and long (L) length track signal regions

Covariance matrix for the 68 signal regions of the disappearing tracks $M_\mathrm{T2}$ search

Correlation matrix for the 68 signal regions of the disappearing tracks $M_\mathrm{T2}$ search

Measured and predicted fiducial cross-sections in both the lllljj and ll$\nu\nu$jj channels for the inclusive ZZjj processes. Uncertainties due to different sources are presented.

Measured and predicted fiducial cross-sections in both the lllljj and ll$\nu\nu$jj channels for the inclusive ZZjj processes. Uncertainties due to different sources are presented.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ QCD CR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ QCD CR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\nu\nu jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\nu\nu jj$ SR.

Measured unfolded differential cross-section as a function of the dilepton transverse momentum $p^{ll}_{\mathrm{T}}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the the four-body transverse momentum $p^{ll\gamma\gamma}_{\mathrm{T}}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the diphoton invariant mass $m_{\gamma\gamma}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the four-body invariant mass $m_{ll\gamma\gamma}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Expected and observed $95\%$ confidence intervals for the coupling parameters $f_{T,j}/\Lambda^{4}$ of transverse dimension-8 operators. All parameter values outside of the stated range are excluded at the chosen confidence level. No unitarity constraints are applied.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,8}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,0}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,1}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,2}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,5}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,6}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,7}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,9}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Ranking of the systematic uncertainties on the measured cross-section, normalised to the predicted value, in the inclusive fit to data. The impact of each nuisance parameter, $\Delta \sigma_{\text{inc}}/\sigma^{\text{pred.}}_{\text{inc}}$, is computed by comparing the nominal best-fit value of $\sigma_{\text{inc}}/\sigma^{\text{pred}}_{\text{inc}}$ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\theta$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta$ ($\pm \Delta \hat{\theta}$). The figure shows the effect of the ten most significant uncertainties.

Ranking of the systematic uncertainties on the measured cross-section, normalised to the predicted value, in the fiducial fit to data. The impact of each nuisance parameter, $\Delta \sigma_{\text{fid}}/\sigma^{\text{pred.}}_{\text{fid}}$, is computed by comparing the nominal best-fit value of $\sigma_{\text{fid}}/\sigma^{\text{pred}}_{\text{fid}}$ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\theta$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta$ ($\pm \Delta \hat{\theta}$). The figure shows the effect of the ten most significant uncertainties.

Ranking of the systematic uncertainties on the measured cross-section, normalised to the predicted value, in the fiducial fit to data. The impact of each nuisance parameter, $\Delta \sigma_{\text{fid}}/\sigma^{\text{pred.}}_{\text{fid}}$, is computed by comparing the nominal best-fit value of $\sigma_{\text{fid}}/\sigma^{\text{pred}}_{\text{fid}}$ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\theta$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta$ ($\pm \Delta \hat{\theta}$). The figure shows the effect of the ten most significant uncertainties.

Impact of different categories of systematic uncertainties on the fiducial and inclusive measurements. The quoted values are obtained by repeating the fit, fixing a set of nuisance parameters of the sources corresponding to the considered category, and subtracting in quadrature the resulting uncertainty from the total uncertainty of the nominal fit. The total uncertainty is different from the sum in quadrature of the different components due to correlations between nuisance parameters built by the fit.

Impact of different categories of systematic uncertainties on the fiducial and inclusive measurements. The quoted values are obtained by repeating the fit, fixing a set of nuisance parameters of the sources corresponding to the considered category, and subtracting in quadrature the resulting uncertainty from the total uncertainty of the nominal fit. The total uncertainty is different from the sum in quadrature of the different components due to correlations between nuisance parameters built by the fit.

Fiducial region definition

Fiducial region definition

Azimuthal correlation distribution between heavy-flavour decay electrons and charged hadrons, scaled by the number of electrons in EMCal triggered events in the electron transverse momentum range 4.5-6 GeV/c.

Relative beauty contribution to the heavy-flavour electron yield obtained with the method based on the track impact parameter.

Relative beauty contribution to the heavy-flavour electron yield obtained with the method based on the track impact parameter.

Relative beauty contribution to the heavy-flavour electron yield obtained from the method based on the azimuthal correlation of heavy-flavour decay electrons and charged tracks.

Relative beauty contribution to the heavy-flavour electron yield obtained from the method based on the azimuthal correlation of heavy-flavour decay electrons and charged tracks.

The $p_{\rm T}$-differential inclusive cross section of electrons (multiplied by $10^3$ for readability) from beauty hadron decays using the electron-hadron azimuthal correlation method.

The pT-differential inclusive cross section of electrons from beauty hadron decays using the electron-hadron azimuthal correlation method.

The $p_{\rm T}$-differential inclusive cross section of electrons (multiplied by $10^3$ for readability) from beauty hadron decays. A 1.9% normalization uncertainty not shown.

The pT-differential inclusive cross section of electrons from beauty hadron decays. A 1.9% normalization uncertainty not shown.

Inclusive bbbar production cross section per rapidity unit measured in mid-rapidity.

Inclusive bbbar production cross section per rapidity unit measured.

Total bbbar production cross section.

Total bbbar production cross section.

Nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of the average number of participants $\langle N_{\mathrm{part}} \rangle$ or as a function of the collision centrality.

Nuclear modification factor of $\Upsilon(2\mathrm{S})$ as a function of the average number of participants $\langle N_{\mathrm{part}} \rangle$ or as a function of the collision centrality.

Ratio of $\Upsilon(2\mathrm{S})$ and $\Upsilon(1\mathrm{S})$ yields (cf. equation 2 in the article for the definition) as a function of the average number of participants $\langle N_{\mathrm{part}} \rangle$ or as a function of the collision centrality. The global uncertainty is the quadratic sum of the branching ratio uncertainties.

Relative nuclear modification factor or double yield ratio between $\Upsilon(2\mathrm{S})$ and $\Upsilon(1\mathrm{S})$ as a function of the average number of participants $\langle N_{\mathrm{part}} \rangle$ or as a function of the collision centrality. The global uncertainty corresponds to the systematic uncertainty on the cross-section ratio in proton–proton collisions.

Nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of transverse momentum for the 0–90% centrality interval.

Nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of rapidity for the 0–90% centrality interval.

Nuclear modification factor of $\Upsilon(2\mathrm{S})$ as a function of rapidity for the 0–90% centrality interval.

Interpolated production cross sections and cross-section ratio for $\Upsilon(1\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ and $\Upsilon(2\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ in proton–proton collisions, integrated over the rapidity interval 2.5–4.0 and $p_{\mathrm{T}} <$ 15 GeV/$c$. Results used for the determination of the integrated nuclear modification factors and of the nuclear modification factors as a function of centrality (Figures 4 and 5 (right)).

Interpolated $p_{\mathrm{T}}$-differential cross section for $\Upsilon(1\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ in proton–proton collisions. Results used for the determination of the nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of its transverse momentum (Figure 6).

Interpolated rapidity-differential cross section for $\Upsilon(1\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ in proton–proton collisions. Results used for the determination of the nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of rapidity (Figure 7).

Interpolated rapidity-differential cross section for $\Upsilon(2\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ in proton–proton collisions. Results used for the determination of the nuclear modification factor of $\Upsilon(2\mathrm{S})$ as a function of rapidity (Figure 7).