Fig. 11 right: The $p_\mathrm{T}$ differential production cross section of charged-particle anti-$k_{T}$ $R=0.4$ b jets in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}~=~5.02$ TeV using the IP method.

Fig. 12 left: The combined differential production cross section of charged-particle anti-$k_{T}$ $R=0.4$ b jets in pp collisions at $\sqrt{s}~=~5.02$ TeV.

Fig. 12 right: The combined differential production cross section of charged-particle anti-$k_{T}$ $R=0.4$ b jets in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}~=~5.02$ TeV.

Fig. 13 left: The b-jet fraction in pp collisions at $\sqrt{s}~=~5.02$ TeV.

Fig. 13 right: The b-jet fraction in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}~=~5.02$ TeV.

Fig. 14 left: The nuclear modification factor $R_{\mathrm{pPb}}$ of b jets as a function of $P_T$ from the SV method.

Fig. 14 left: The nuclear modification factor $R_{\mathrm{pPb}}$ of b jets as a function of $P_T$ from the IP method.

Fig. 14 right: The combined nuclear modification factor $R_{\mathrm{pPb}}$ of b jets as a function of $P_T$.

95% CL expected and observed frequentist upper limit for the LbL cross section within the $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$ search region.

95% expected and observed one-dimensional limits on $\zeta_1$ and $\zeta_2$ anomalous LbyL production parameters, when the other parameter is set to zero. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

95% expected and observed limits on $\zeta_1$ and $\zeta_2$ anomalous LbyL production parameters. The parametric elliptic form is assumed: $a_0\zeta_1^2+a_1\zeta_1\zeta_2+a_2\zeta_2^2=1$. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1400-1500 MeV.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1500-1600 MeV.

$M_{p\eta}$ invariant mass distribution under the conditions: $E_\gamma$ = 1400-1500 MeV; 1120 MeV $\leq M_{p\pi^0} \leq$ 1220 MeV; $\theta_p^{lab} \leq$ 25°; 25° $\leq \theta_\gamma^{lab} \leq$ 155°.

$M_{p\eta}$ invariant mass distribution under the conditions: $E_\gamma$ = 1400-1500 MeV; 1120 MeV $\leq M_{p\pi^0} \leq$ 1220 MeV; $\theta_p^{lab} \leq$ 25°; 1° $\leq \theta_\gamma^{lab} \leq$ 155°.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1400-1500 MeV and $M_{p\pi^0} \leq$ 1190 MeV, allowing for all proton and photon laboratory angles.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1170 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1180 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1200 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1210 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1220 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1400-1420 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1460 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1460-1500 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1500-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1540-1600 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{p\eta}$.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{\pi^0\eta}$.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{p\pi^0}$, with a cut on the structure in the $M^2_{p\eta}$ invariant mass distribution.

Excess yield relative to the PWA distributions in the $p-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $p-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $\pi^0-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $\pi^0-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Differential cross section $d\sigma/dM_{\pi^0\eta}$ for $E_\gamma$ = 1540-1600 MeV.

$M_{p\eta}$ invariant mass distribution, taken as difference to the PWA calculations (s. Fig. 11).

$M_{\pi^0\eta}$ invariant mass distribution, taken as difference to the PWA calculations (s. Fig. 11).

$\theta_{p\eta}$ opening angle distribution, taken as difference to the PWA calculations (s. Fig. 12).

$\theta_{\pi^0\eta}$ opening angle distribution, taken as difference to the PWA calculations (s. Fig. 12).

Two-dimensional distributions of the leading electron vs the leading muon $|d_{0}|$, for the events that pass the e$\mu$ preselection. Data events and $\tilde{t} \to b\ell$ signal events with a $\tilde{t}$ mass of 700 GeV and a proper decay length of 10 mm that correspond to 2018 data-taking conditions are shown. In each $|d_{0}|$-$|d_{0}|$ bin, the number of events divided by the bin area is plotted. The inclusive signal region covers the region between 100 $\mu$m and 10 cm in each $|d_{0}|$ variable shown.

Two-dimensional distributions of the leading vs the subleading electron $|d_0|$, for the events that pass the ee preselection. Data events and $\tilde{t} \to b\ell$ signal events with a $\tilde{t}$ mass of 700 GeV and a proper decay length of 10 mm that correspond to 2018 data-taking conditions are shown. In each $|d_{0}|$-$|d_{0}|$ bin, the number of events divided by the bin area is plotted. The inclusive signal region covers the region between 100 $\mu$m and 10 cm in each $|d_{0}|$ variable shown.

Two-dimensional distributions of the leading vs the subleading muon $|d_{0}|$, for the events that pass the $\mu\mu$ preselection. Data events and $\tilde{t} \to b\ell$ signal events with a $\tilde{t}$ mass of 700 GeV and a proper decay length of 10 mm that correspond to 2018 data-taking conditions are shown. In each $|d_{0}|$-$|d_{0}|$ bin, the number of events divided by the bin area is plotted. The inclusive signal region covers the region between 100 $\mu$m and 10 cm in each $|d_{0}|$ variable shown.

The number of observed and estimated background events in each channel and SR, with a representative signal overlaid. For each background estimate and signal yield, the total uncertainty (statistical plus systematic) is given. The distributions shown correspond to the events before the final maximum likelihood fit to the data.

The number of observed and estimated background events in each channel and SR in 2016, with a representative signal overlaid. For each background estimate and signal yield, the total uncertainty is given. The distributions shown correspond to the events before the final maximum likelihood fit to the data.

The number of observed and estimated background events in each channel and SR in 2017 and 2018, with a representative signal overlaid. For each background estimate and signal yield, the total uncertainty is given. The distributions shown correspond to the events before the final maximum likelihood fit to the data.

The observed 95% CL upper limits on the long-lived top squark production cross section, in the $c\tau_0$ -- mass plane, for the three channels combined. The $\tilde{t} \to b\ell $ process is shown. These limits assume $\mathcal{B}(\tilde{t} \to b\ell)$ is 100%, and each lepton has an equal probability of being an electron, a muon, or a tau lepton. The area to the left of the black curve represents the observed exclusion region, and the dashed red lines indicate the expected limits and their 68% confidence intervals.

The observed 95% CL upper limits on the long-lived top squark production cross section, in the $c\tau_0$ -- mass plane, for the three channels combined. The $\tilde{t} \to d\ell$ process is shown. These limits assume $\mathcal{B}(\tilde{t} \to d\ell)$ is 100%, and each lepton has an equal probability of being an electron, a muon, or a tau lepton. The area to the left of the black curve represents the observed exclusion region, and the dashed red lines indicate the expected limits and their 68% confidence intervals.

The observed 95% CL upper limits on the long-lived top squark production cross section, in the $c\tau_0$ -- mass plane, for the three channels combined. The $\tilde{t} \to b\ell $ process is shown. These limits assume $\mathcal{B}(\tilde{t} \to b\ell)$ is 100%, and each lepton has an equal probability of being an electron, a muon, or a tau lepton. The area to the left of the black curve represents the observed exclusion region, and the dashed red lines indicate the expected limits and their 68% confidence intervals.

The observed 95% CL upper limits on the long-lived top squark production cross section, in the $c\tau_0$ -- mass plane, for the three channels combined. The $\tilde{t} \to d\ell$ process is shown. These limits assume $\mathcal{B}(\tilde{t} \to d\ell)$ is 100%, and each lepton has an equal probability of being an electron, a muon, or a tau lepton. The area to the left of the black curve represents the observed exclusion region, and the dashed red lines indicate the expected limits and their 68% confidence intervals.

The 95% CL upper limits on the long-lived top squark proper decay length ($c\tau_0$) as a function of its mass, for the e$\mu$, ee, and $\mu\mu$ channels, and their combination. The $\tilde{t} \to b\ell $ process is shown. These limits assume $\mathcal{B}(\tilde{t} \to b\ell)$ is 100%, and each lepton has an equal probability of being an electron, a muon, or a tau lepton.

The 95% CL upper limits on the long-lived top squark proper decay length ($c\tau_0$) as a function of its mass, for the e$\mu$, ee, and $\mu\mu$ channels, and their combination. The $\tilde{t} \to d\ell$ process is shown. These limits assume $\mathcal{B}(\tilde{t} \to d\ell)$ is 100%, and each lepton has an equal probability of being an electron, a muon, or a tau lepton.

The 95% CL constraints on the long-lived slepton $c\tau_{0}$ and mass. The $\tilde{\tau}$ and co-NLSP limits are shown for the three channels combined, while the $\tilde{e}$ and $\tilde{\mu}$ NLSP limits are shown for the ee and $\mu\mu$ channels, respectively. These limits assume that the superpartners of the left- and right-handed leptons are degenerate in mass and $\mathcal{B}(\tilde{\ell} \to \ell\tilde{G})$ is 100%. The area to the left of the solid curves represents the observed exclusion region, and the dashed lines indicate the expected limits.

The observed 95% CL upper limits on the long-lived slepton production cross section, in the $c\tau_0$ -- mass plane. The co-NLSP limits are shown for the three channels combined. These limits assume that the superpartners of the left- and right-handed leptons are degenerate in mass and $\mathcal{B}(\tilde{\ell} \to \ell\tilde{G})$ is 100%. The area to the left of the black curve represents the observed exclusion region, and the dashed red lines indicate the expected limits and their 68% confidence intervals.

The observed 95% CL upper limits on the long-lived slepton production cross section, in the $c\tau_0$ -- mass plane. The co-NLSP limits are shown for the three channels combined. These limits assume that the superpartners of the left- and right-handed leptons are degenerate in mass and $\mathcal{B}(\tilde{\ell} \to \ell\tilde{G})$ is 100%. The area to the left of the black curve represents the observed exclusion region, and the dashed red lines indicate the expected limits and their 68% confidence intervals.

The 95% CL upper limits on the $\mathrm{H} \to \mathrm{S}\mathrm{S}, \mathrm{S} \to \ell^{+}\ell^{-}$ branching fraction as a function of $c\tau_0$, for a Higgs boson with a mass of 125 GeV and a long-lived scalar with a mass of 30 GeV or 50 GeV, for the three channels combined. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two electrons or two muons. The area above the solid (dashed) curve represents the observed (expected) exclusion region.

The 95% CL upper limits on the product of the cross section and branching fraction $\mathrm{H} \to \mathrm{S}\mathrm{S}, \mathrm{S} \to \ell^{+}\ell^{-}$ as a function of $c\tau_0$, for a 300 GeV Higgs boson and several long-lived scalar masses, for the three channels combined. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two electrons or two muons. The area above the solid (dashed) curves represents the observed (expected) exclusion region.

The 95% CL upper limits on the product of the cross section and branching fraction $\mathrm{H} \to \mathrm{S}\mathrm{S}, \mathrm{S} \to \ell^{+}\ell^{-}$ as a function of $c\tau_0$, for a 400 GeV Higgs boson and several long-lived scalar masses, for the three channels combined. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two electrons or two muons. The area above the solid (dashed) curves represents the observed (expected) exclusion region.

The 95% CL upper limits on the product of the cross section and branching fraction $\mathrm{H} \to \mathrm{S}\mathrm{S}, \mathrm{S} \to \ell^{+}\ell^{-}$ as a function of $c\tau_0$, for a 600 GeV Higgs boson and several long-lived scalar masses, for the three channels combined. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two electrons or two muons. The area above the solid (dashed) curves represents the observed (expected) exclusion region.

The 95% CL upper limits on the product of the cross section and branching fraction $\mathrm{H} \to \mathrm{S}\mathrm{S}, \mathrm{S} \to \ell^{+}\ell^{-}$ as a function of $c\tau_0$, for a 800 GeV Higgs boson and several long-lived scalar masses, for the three channels combined. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two electrons or two muons. The area above the solid (dashed) curves represents the observed (expected) exclusion region.

The 95% CL upper limits on the product of the cross section and branching fraction $\mathrm{H} \to \mathrm{S}\mathrm{S}, \mathrm{S} \to \ell^{+}\ell^{-}$ as a function of $c\tau_0$, for a 1000 GeV Higgs boson and several long-lived scalar masses, for the three channels combined. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two electrons or two muons. The area above the solid (dashed) curves represents the observed (expected) exclusion region.

The 95% CL upper limits on the $\mathrm{H} \to \mathrm{S}\mathrm{S} \to 4\ell$ branching fraction as a function of $c\tau_0$, for a Higgs boson with a mass of 125 GeV and a long-lived scalar with a mass of 30 GeV or 50 GeV, for the ee channel only. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two electrons. The area above the solid (dashed) curve represents the observed (expected) exclusion region. The curves are not smooth because very few ee channel events pass the preselection for the Higgs boson and scalar masses shown in this figure.

The 95% CL upper limits on the product of the cross section and branching fraction $\mathrm{H} \to \mathrm{S}\mathrm{S} \to 4\ell$ as a function of $c\tau_0$, for a 600 GeV Higgs boson and several long-lived scalar masses, for the ee channel only. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two electrons. The area above the solid (dashed) curves represents the observed (expected) exclusion region.

The 95% CL upper limits on the product of the cross section and branching fraction $\mathrm{H} \to \mathrm{S}\mathrm{S} \to 4\ell$ as a function of $c\tau_0$, for a 1000 GeV Higgs boson and several long-lived scalar masses, for the ee channel only. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two electrons. The area above the solid (dashed) curves represents the observed (expected) exclusion region.

The 95% CL upper limits on the $\mathrm{H} \to \mathrm{S}\mathrm{S} \to 4\ell$ branching fraction as a function of $c\tau_0$, for a Higgs boson with a mass of 125 GeV and a long-lived scalar with a mass of 30 GeV or 50 GeV, for the $\mu\mu$ channel only. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two muons. The area above the solid (dashed) curve represents the observed (expected) exclusion region.

The 95% CL upper limits on the product of the cross section and branching fraction $\mathrm{H} \to \mathrm{S}\mathrm{S} \to 4\ell$ as a function of $c\tau_0$, for a 600 GeV Higgs boson and several long-lived scalar masses, for the $\mu\mu$ channel only. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two muons. The area above the solid (dashed) curves represents the observed (expected) exclusion region.

The 95% CL upper limits on the product of the cross section and branching fraction $\mathrm{H} \to \mathrm{S}\mathrm{S} \to 4\ell$ as a function of $c\tau_0$, for a 1000 GeV Higgs boson and several long-lived scalar masses, for the $\mu\mu$ channel only. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two muons. The area above the solid (dashed) curves represents the observed (expected) exclusion region.

Observed likelihood scan in the ($\kappa_\tau$, $\tilde{\kappa}_\tau$) plane.

Observed likelihood scan in the ($\alpha^{\mathrm{H}\tau\tau}$, $\mu_{\text{ggH}}$) plane. The $\mu_{\text{qqH}}$ signal strength modifier is profiled.

Observed likelihood scan in the ($\alpha^{\mathrm{H}\tau\tau}$, $\mu_{\text{qqH}}$) plane. The $\mu_{\text{ggH}}$ signal strength modifier is profiled.

Observed likelihood scan in the ($\mu_{\text{qqH}}$, $\mu_{\text{ggH}}$) plane. The value of $\alpha^{\mathrm{H}\tau\tau}$ is profiled.

Upper limits on $\mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ at 95% CL including the BDT cut as a function of the signal mass.

Upper limits on $\mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ at 95% CL without the BDT cut as a function of the signal mass.

The number of events when applying each cut in the analysis for $m_a$ = 16 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 20 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 30 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 40 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 50 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 60 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

Data/MC ratio as a function of reconstructed $\pi^0$ momentum for the 2$\gamma$0p selection

Energy spectra for the 1$\gamma$1p selected events. The figure shows the unconstrained background prediction and breakdowns as a function of reconstructed shower energy.

Energy spectra for the 1$\gamma$1p selected events. This figure shows the total background prediction with systematic uncertainty both before and after the 2$\gamma$ constraint.

Energy spectra for the 1$\gamma$0p selected events. The figure shows the unconstrained background prediction and breakdowns as a function of reconstructed shower energy.

Energy spectra for the 1$\gamma$0p selected events. This figure shows the total background prediction with systematic uncertainty both before and after the 2$\gamma$ constraint.

The observed event rate for the 1$\gamma$1p event sample, and comparisons to unconstrained (left) and constrained (right) background and LEE model predictions.

The observed event rate for the 1$\gamma$0p event sample, and comparisons to unconstrained (left) and constrained (right) background and LEE model predictions.

The fractional systematic covariance matrix for the final selected 1$\gamma$1p, 1$\gamma$0p, 2$\gamma$1p, and 2$\gamma$0p event samples, in bins of reconstructed shower energy for the 1$\gamma$ distributions, and reconstructed $\pi^{0}$ momentum for the 2$\gamma$ distributions. The bin boundaries are defined in Sec. II in supplementary material. This covariance matrix is evaluated at central value prediction and all variances have been rounded to three significant digits.

The fractional systematic covariance matrix for the final selected 1$\gamma$1p, 1$\gamma$0p, 2$\gamma$1p, and 2$\gamma$0p event samples, in bins of reconstructed shower energy for the 1$\gamma$1p signal, 1$\gamma$1p background, 1$\gamma$0p signal and 1$\gamma$0p background distributions, and reconstructed $\pi^{0}$ momentum for the 2$\gamma$ distributions. The bin boundaries are defined in Sec. II in supplementary material. This covariance matrix is evaluated at central value prediction and variances have been rounded to three significant digits.

Rapidity Differential production cross sections of J/$\psi$ in the range 2 $\leq$ p$_{\rm T}$ < 4 GeV/c

Rapidity Differential production cross sections of J/$\psi$ in the range 4 $\leq$ p$_{\rm T}$ < 6 GeV/c

Rapidity Differential production cross sections of J/$\psi$ in the range 6 $\leq$ p$_{\rm T}$ < 12 GeV/c

Differential production cross sections of $\psi$(2S) as a function of $p_{\rm T}$.

Differential production cross sections of $\psi$(2S) as a function of rapidity.

$\psi$(2S) over J/$\psi$ ratio as a function of $p_{\rm T}$.

$\psi$(2S) over J/$\psi$ ratio as a function of rapidity.

Differential production cross sections of $\Upsilon(1S)$ as a function of $p_{\rm T}$.

Differential production cross sections of $\Upsilon(1S)$ as a function of rapidity.

J/$\psi$ 5-to-13 TeV cross section ratio as a function of $p_{\rm T}$.

J/$\psi$ 7-to-13 TeV cross section ratio as a function of $p_{\rm T}$.

J/$\psi$ 8-to-13 TeV cross section ratio as a function of $p_{\rm T}$.

J/$\psi$ 5-to-13 TeV cross section ratio as a function of rapidity.

J/$\psi$ 7-to-13 TeV cross section ratio as a function of rapidity.

J/$\psi$ 8-to-13 TeV cross section ratio as a function of rapidity.

$\psi$(2S) 5-to-13 TeV cross section ratio as a function of rapidity.

$\psi$(2S) 7-to-13 TeV cross section ratio as a function of rapidity.

$\psi$(2S) 8-to-13 TeV cross section ratio as a function of rapidity.

$\psi$(2S) 5-to-13 TeV cross section ratio as a function of $p_{\rm T}$.

$\psi$(2S) 7-to-13 TeV cross section ratio as a function of $p_{\rm T}$.

$\psi$(2S) 8-to-13 TeV cross section ratio as a function of $p_{\rm T}$.

$\psi$(2S)-to-J/$\psi$ cross section ratio as a function of $\sqrt{s}$.

Integrated production cross section of J/$\psi$.

Integrated production cross section of $\psi$(2S).

Integrated production cross section of $\Upsilon$(1S).

Integrated production cross section of $\Upsilon$(2S).

Integrated production cross section of $\Upsilon$(3S).

$\Upsilon$(2S)-to-$\Upsilon$(1S) cross section ratio.

$\Upsilon$(3S)-to-$\Upsilon$(1S) cross section ratio.

Comparison of the pT spectrum for the fourth-leading from data to different PYTHIA8 (P8),HERWIG++ (H++),and HERWIG7 (H7) tunes.

Comparison of the eta spectrum for the leading jet from data to different PYTHIA8 (P8),HERWIG++ (H++),and HERWIG7 (H7) tunes.

Comparison of the eta spectrum for the sub-leading from data to different PYTHIA8 (P8),HERWIG++ (H++),and HERWIG7 (H7) tunes.

Comparison of the eta spectrum for the third-leading from data to different PYTHIA8 (P8),HERWIG++ (H++),and HERWIG7 (H7) tunes.

Comparison of the eta spectrum for the fourth-leading from data to different PYTHIA8 (P8),HERWIG++ (H++),and HERWIG7 (H7) tunes.

Comparison of the DeltaPhiSoft distribution from data to different PYTHIA8 (P8),HERWIG++ (H++), and HERWIG7 (H7) tunes. All distributions have been normalized to regions where a reduced DPS contribution is expected.

Comparison of the DeltaPhiMin distribution from data to different PYTHIA8 (P8),HERWIG++ (H++), and HERWIG7 (H7) tunes. All distributions have been normalized to regions where a reduced DPS contribution is expected.

Comparison of the DeltaY distribution from data to different PYTHIA8 (P8),HERWIG++ (H++), and HERWIG7 (H7) tunes. All distributions have been normalized to regions where a reduced DPS contribution is expected.

Comparison of the phi_ij distribution from data to different PYTHIA8 (P8),HERWIG++ (H++), and HERWIG7 (H7) tunes. All distributions have been normalized to regions where a reduced DPS contribution is expected.

Comparison of the Deltap_{T,Soft} distribution from data to different PYTHIA 8 (P8), HERWIG ++ (H++), and HERWIG 7 (H7) tunes. All distributions have been normalized to regions where a reduced DPS contribution is expected.

Comparison of the DeltaS distribution from data to different PYTHIA 8 (P8), HERWIG ++ (H++), and HERWIG 7 (H7) tunes. All distributions have been normalized to regions where a reduced DPS contribution is expected.

Comparison of the pT spectrum for the leading jet from data with different KATIE(KT), MADGRAPH5aMC@NLO(MG5), and POWHEG(PW) models

Comparison of the pT spectrum for the sub-leading from data with different KATIE(KT), MADGRAPH5aMC@NLO(MG5), and POWHEG(PW) models

Comparison of the pT spectrum for the third-leading from data to different PYTHIA8 (P8),HERWIG++ (H++),and HERWIG7 (H7) tunes.

Comparison of the pT spectrum for the fourth-leading from data with different KATIE(KT), MADGRAPH5aMC@NLO(MG5), and POWHEG(PW) models

Comparison of the eta spectrum for the leading jet from data with different KATIE(KT), MADGRAPH5aMC@NLO(MG5), and POWHEG(PW) models

Comparison of the eta spectrum for the sub-leading from data with different KATIE(KT), MADGRAPH5aMC@NLO(MG5), and POWHEG(PW) models

Comparison of the eta spectrum for the third-leading from data with different KATIE(KT), MADGRAPH5aMC@NLO(MG5), and POWHEG(PW) models

Comparison of the eta spectrum for the fourth-leading from data with different KATIE(KT), MADGRAPH5aMC@NLO(MG5), and POWHEG(PW) models

Comparison of the DeltaPhiSoft distribution from data to differentKATIE(KT), MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations. All distributions have been normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the DeltaPhiMin, distribution from data to differentKATIE(KT), MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations. All distributions have been normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the DeltaY and distribution from data to differentKATIE(KT), MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations. All distributions have been normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the phi_ij distribution from data to differentKATIE(KT), MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations. All distributions have been normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the Deltap_{T,Soft} distribution from data to different KATIE(KT),MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations. All distributions havebeen normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the DeltaS distribution from data to different KATIE(KT),MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations. All distributions havebeen normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the pT spectrum for the leading jet from data with different SPS+DPS KATIE( KT) and PYTHIA8 (P8) models.

Comparison of the pT spectrum for the sub-leading from data with different SPS+DPS KATIE( KT) and PYTHIA8 (P8) models.

Comparison of the pT spectrum for the third-leading from data with different SPS+DPS KATIE( KT) and PYTHIA8 (P8) models.

Comparison of the pT spectrum for the fourth-leading from data with different SPS+DPS KATIE( KT) and PYTHIA8 (P8) models.

Comparison of the eta spectrum for the leading jet from data with different SPS+DPS KATIE( KT) and PYTHIA8 (P8) models.

Comparison of the eta spectrum for the sub-leading from data with different SPS+DPS KATIE( KT) and PYTHIA8 (P8) models.

Comparison of the eta spectrum for the third-leading from data with different SPS+DPS KATIE( KT) and PYTHIA8 (P8) models.

Comparison of the eta spectrum for the fourth-leading from data with different SPS+DPS KATIE( KT) and PYTHIA8 (P8) models.

Comparison of the DeltaPhiSoft distribution from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models. All distributions have been normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the DeltaPhiMin distribution from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models. All distributions have been normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the DeltaY distribution from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models. All distributions have been normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the phi_ij distribution from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models. All distributions have been normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the Deltap_{T,Soft} distributions from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models. All distributions havebeen normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the DeltaS distributions from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models. All distributions have been normalized to regions where a reduced DPS contribution is expected.

The DeltaS distribution obtained from the mixed data sample compared to predictions from the pure DPS sample in PYTHIA 8 (P8) and KATIE (KT). The distributions are normalized to unity.

The results of the template fit for the POWHEG (PW) NLO 2 -> 2 model without the hard MPI removed. As the DeltaS distribution obtained from the mixed data sample carries a statistical and systematic uncertainty, so does the total fitted sample.

Comparison of the DeltaPhiSoft distribution from data to different PYTHIA8 (P8),HERWIG++ (H++), and HERWIG7 (H7) tunes.

Comparison of the DeltaPhiMin distribution from data to different PYTHIA8 (P8),HERWIG++ (H++), and HERWIG7 (H7) tunes.

Comparison of the DeltaY distribution from data to different PYTHIA8 (P8),HERWIG++ (H++), and HERWIG7 (H7) tunes.

Comparison of the phi_ij distribution from data to different PYTHIA8 (P8),HERWIG++ (H++), and HERWIG7 (H7) tunes.

Comparison of the Deltap_{T,Soft} distribution from data to different PYTHIA 8 (P8), HERWIG ++ (H++), and HERWIG 7 (H7) tunes.

Comparison of the DeltaS distribution from data to different PYTHIA 8 (P8), HERWIG ++ (H++), and HERWIG 7 (H7) tunes.

Comparison of the DeltaPhiSoft distribution from data to different KATIE(KT),MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations.

Comparison of the DeltaPhiMin distribution from data to different KATIE(KT),MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations.

Comparison of the DeltaY distribution from data to different KATIE(KT),MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations.

Comparison of the phi_ij distribution from data to different KATIE(KT),MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations.

Comparison of the Deltap_{T,Soft} distribution from data to different KATIE(KT),MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations.

Comparison of the DeltaS distribution from data to different KATIE(KT),MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations.

Comparison of the DeltaPhiSoft distribution from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models.

Comparison of the DeltaPhiMin distribution from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models.

Comparison of the DeltaY distribution from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models.

Comparison of the phi_ij distribution from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models.

Comparison of the Deltap_{T,Soft} distribution from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models.

Comparison of the DeltaS distributions from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models.

(Anti)proton spectrum in HM III V0M multiplicity class

(Anti)deuteron spectrum in HM I V0M multiplicity class

(Anti)deuteron spectrum in HM II V0M multiplicity class

(Anti)deuteron spectrum in HM III V0M multiplicity class

(Anti)triton spectrum in HM V0M multiplicity class

(Anti)triton spectrum in MB V0M multiplicity class

(Anti)helion spectrum in HM V0M multiplicity class

(Anti)helion spectrum in MB V0M multiplicity class

(Anti)helion spectrum in MB I V0M multiplicity class

(Anti)helion spectrum in MB II V0M multiplicity class

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in HM I V0M multiplicity class

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in HM II V0M multiplicity class

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in HM III V0M multiplicity class

Coalescence parameter $B_3$ as a function of $p_{\mathrm{T}}$ in HM V0M multiplicity class

Coalescence parameter $B_3$ as a function of $p_{\mathrm{T}}$ in MB V0M multiplicity class

Coalescence parameter $B_3$ as a function of $p_{\mathrm{T}}$ in MB I V0M multiplicity class

Coalescence parameter $B_3$ as a function of $p_{\mathrm{T}}$ in MB II V0M multiplicity class

Deuteron over proton ratio in high-multiplicity pp collisions as a function of the charged particle multiplicity density at mid-rapidity

Helion over proton ratio in minimum-bias pp collisions as a function of the charged particle multiplicity density at mid-rapidity

Helion over proton ratio in high-multiplicity pp collisions as a function of the charged particle multiplicity density at mid-rapidity

Triton over helion ratio in high-multiplicity pp collisions as a function of the transverse momentum