The correlation matrix corresponding to the statistical uncertainties on the differential cross-section ($d\sigma/dp_T$) fit results for $W^-$. To combine with $W^+$ results (Fig8a), use the rows and columns ordered as $W^+$ and then $W^-$. Assume no correlation in the statistical uncertainties between $W^+$ and $W^-$ (zero entries in the off-diagonal blocks).

The correlation matrix corresponding to all of the systematic uncertainties on the differential cross-section ($d\sigma/dp_T$) fit results for $W^+$ and $W^-$ combined together, with rows and columns ordered as $W^+$ and then $W^-$.

Differential cross section times dimuon branching fraction of Y(1S) as a function of $y^{Y}_{CM}$ in pPb collisions. The global uncertainty arises from the integrated luminosity uncertainty in pPb collisions.

Differential cross section times dimuon branching fraction of Y(2S) as a function of $y^{Y}_{CM}$ in pPb collisions. The global uncertainty arises from the integrated luminosity uncertainty in pPb collisions.

Differential cross section times dimuon branching fraction of Y(3S) as a function of $y^{Y}_{CM}$ in pPb collisions. The global uncertainty arises from the integrated luminosity uncertainty in pPb collisions.

Differential cross section times dimuon branching fraction of Y(1S) as a function of pT in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(2S) as a function of pT in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(3S) as a function of pT in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(1S) as a function of $|y^{Y}_{CM}|$ in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(2S) as a function of $|y^{Y}_{CM}|$ in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(3S) as a function of $|y^{Y}_{CM}|$ in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Nuclear modification factor of Y(1S) as a function of pT. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) as a function of pT. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) as a function of pT. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(1S) as a function of $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) as a function of $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) as a function of $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(1S) at forward and backward $y^{Y}_{CM}$ for pT < 6 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) at forward and backward $y^{Y}_{CM}$ for pT < 6 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) at forward and backward $y^{Y}_{CM}$ for pT < 6 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(1S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

RFB of Y(1S) versus $N^{|\eta_{lab}|<2.4}_{tracks}$.

RFB of Y(2S) versus $N^{|\eta_{lab}|<2.4}_{tracks}$.

RFB of Y(3S) versus $N^{|\eta_{lab}|<2.4}_{tracks}$.

RFB of Y(1S) versus $E^{|\eta_{lab}|>4}_{T}$.

RFB of Y(2S) versus $E^{|\eta_{lab}|>4}_{T}$.

RFB of Y(3S) versus $E^{|\eta_{lab}|>4}_{T}$.

Nuclear modification factor of Y(1S), Y(2S), and Y(3S) integrated over pT and $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

$D^+ D^0$~mass distributions for selected candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^+ D^+$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^+ D^0 \pi^+$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Background-subtracted distributions for the multiplicity of tracks reconstructed in the vertex detector for $T_{cc}^+\to D^0 D^0 \pi^+$ signal, low-mass $D^0\bar{D}^0$ and $D^0D^0$ pairs. The binning scheme is chosen to have an approximately uniform distribution for $D^0\bar{D}^0$ pairs. The distributions for the $D^0\bar{D}^0$ and $D^0D^0$ pairs are normalised to the same yields as the $T_{cc}^+\to D^0 D^0 \pi^+$ signal. For better visualisation, the points are slightly displaced from the bin centres. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Background-subtracted transverse momentum spectra for $T_{cc}^+\to D^0 D^0 \pi^+$ signal, low-mass $D^0\bar{D}^0$ and $D^0D^0$ pairs. The binning scheme is chosen to have an approximately uniform distribution for $D^0\bar{D}^0$ pairs. The distributions for the $D^0\bar{D}^0$ and $D^0D^0$ pairs are normalised to the same yields as the $T_{cc}^+\to D^0 D^0 \pi^+$ signal. For better visualisation, the points are slightly displaced from the bin centres. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^0 \bar{D}^0$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^0 D^0$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $\bar{D}^0D^0\pi^+$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^0D^0\pi^+$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^0D^-$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^0D^+$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Distribution of $D^0 D^0 \pi^+$ mass where the contribution of the non-$D^0$ background has been statistically subtracted by assigning the a weight to every candidate.

Mass distribution for $D^0 \pi^+$ pairs from selected $D^0 D^0 \pi^+$ candidates with a mass below the $D^{*+}D^0$ mass threshold with non-$D^0$ background subtracted by assigning the a weight to every candidate.

$D^0 D^0$~mass distributions for selected candidates with the $D^0$ background subtracted by assigning the a weight to every candidate.

$D^+ D^0$~mass distributions for selected candidates with the $D^0$ background subtracted by assigning the a weight to every candidate.

Mass distributions for selected $D^+ D^+$ candidates with the $D^0$ background subtracted by assigning the a weight to every candidate.

Mass distributions for selected $D^+ D^0 \pi^+$ candidates with the $D^0$ background subtracted by assigning the a weight to every candidate.

Mean $p_T$ values of the $\Upsilon(1$S$)$, $\Upsilon(2$S$)$, and $\Upsilon(3$S$)$ states with $p_T\,>\,0\,GeV$ and $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $0<p_T<5\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $5<p_T<7\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $7<p_T<9\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $9<p_T<11\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $11<p_T<15\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $15<p_T<20\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $20<p_T<50\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $0<p_T<5\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $5<p_T<7\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $7<p_T<9\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $9<p_T<11\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $11<p_T<15\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $15<p_T<20\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $20<p_T<50\,$GeV range

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of "forward" track multiplicity $N_{track}^{\Delta\phi}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$. Forward tracks are those with momentum direction in $\Delta\phi\,<\,\pi/3$ w.r.t. the $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of "transverse" track multiplicity $N_{track}^{\Delta\phi}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$. Transverse tracks are those with momentum direction in $\pi/3\,<\,\Delta\phi\,<\,2\pi/3$ w.r.t. the $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of "backward" track multiplicity $N_{track}^{\Delta\phi}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$. Backward tracks are those with momentum direction in $\Delta\phi\,>\,2\pi/3$ w.r.t. the $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$ and $N^{\Delta R}_{track}\,=\,0$ charged particles produced in $\Delta R\,<\,0.5$ cone around $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$ and $N^{\Delta R}_{track}\,=\,1$ charged particles produced in $\Delta R\,<\,0.5$ cone around $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$ and $N^{\Delta R}_{track}\,=\,2$ charged particles produced in $\Delta R\,<\,0.5$ cone around $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$ and $N^{\Delta R}_{track}\,>\,2$ charged particles produced in $\Delta R\,<\,0.5$ cone around $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ in transverse sphericity bin $0.00\,<\,S_T\,<\,0.55$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ in transverse sphericity bin $0.55\,<\,S_T\,<\,0.70$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ in transverse sphericity bin $0.70\,<\,S_T\,<\,0.85$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ in transverse sphericity bin $0.85\,<\,S_T\,<\,1.00$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$.

Efficiency-corrected multiplicity bins used in the $\Upsilon(n$S$)$ ratio analysis and the corresponding mean number of charged particle tracks.

The fiducial cross section measured in bins of the azimuthal angle difference between the mesons for events in the fiducial region with 2 Y(1S) with absolute rapidity less than 2.0.

The expected and observed 95% CL limits on the product of the cross section and branching fraction for a tetraquark. The branching fraction is the product of the branching fraction for the decay of the resonance to a Y(1S) meson and two muons, and the Y(1S) meson dimuon branching fraction.

The expected and observed 95% CL limits on the product of the cross section and branching fraction for a scalar state. The branching fraction is the product of the branching fraction for the decay of the resonance to a Y(1S) meson and two muons, and the Y(1S) meson dimuon branching fraction.

The expected and observed 95% CL limits on the product of the cross section and branching fraction for a pseudoscalar state. The branching fraction is the product of the branching fraction for the decay of the resonance to a Y(1S) meson and two muons, and the Y(1S) meson dimuon branching fraction.

The expected and observed 95% CL limits on the product of the cross section and branching fraction for a spin-2 state. The branching fraction is the product of the branching fraction for the decay of the resonance to a Y(1S) meson and two muons, and the Y(1S) meson dimuon branching fraction.

Kinematic properties and charges of the four muons selected in data in the signal region for the measurement of the YY production cross section. MUON1 and MUON2 form the dimuon pair associated to the dimuon object required at trigger-level. The invariant mass of both muon pairs is between 8.6 and 11 GeV. All events pass the trigger requirements and have four muons with a total charge of 0 and with vertex probability above 5%. The muons have pT above 3.5 GeV if they have an absolute pseudorapidity below 0.9 and 2.5 GeV otherwise. The muons pass the medium muon identification criterion.

The covariance matrix containing the experimental uncertainties of Figure 2a is shown.

The covariance matrix containing the model uncertainties of Figure 2a is shown.

The covariance matrix containing the statistical uncertainties of Figure 2b is shown.

The covariance matrix containing the experimental uncertainties of Figure 2b is shown.

The covariance matrix containing the model uncertainties of Figure 2b is shown.

Values of the predicted signal yields in each of the 15 signal regions (for $ x=0.75 $).

Pre-fit background-only covariance among different search bins.

Expected $95\%$ CL upper limit on the production cross section of top squarks (for $ x=0.25 $).

Top squark - neutralino mass points lying on the expected exclusion contour (for $ x=0.25 $).

Top squark - neutralino mass points lying on the expected-$1\sigma_\text{expected}$ exclusion contour (for $ x=0.25 $).

Top squark - neutralino mass points lying on the expected+$1\sigma_\text{expected}$ exclusion contour (for $ x=0.25 $).

Observed $95\%$ CL upper limit on the production cross section of top squarks (for $ x=0.25 $).

Top squark - neutralino mass points lying on the observed exclusion contour (for $ x=0.25 $).

Top squark - neutralino mass points lying on the observed-$1\sigma_\text{observed}$ exclusion contour (for $ x=0.25 $).

Top squark - neutralino mass points lying on the observed+$1\sigma_\text{observed}$ exclusion contour (for $ x=0.25 $).

Expected $95\%$ CL upper limit on the production cross section of top squarks (for $ x=0.5 $).

Top squark - neutralino mass points lying on the expected exclusion contour (for $ x=0.5 $).

Top squark - neutralino mass points lying on the expected-$1\sigma_\text{expected}$ exclusion contour (for $ x=0.5 $).

Top squark - neutralino mass points lying on the expected+$1\sigma_\text{expected}$ exclusion contour (for $ x=0.5 $).

Observed $95\%$ CL upper limit on the production cross section of top squarks (for $ x=0.5 $).

Top squark - neutralino mass points lying on the observed exclusion contour (for $ x=0.5 $).

Top squark - neutralino mass points lying on the observed-$1\sigma_\text{observed}$ exclusion contour (for $ x=0.5 $).

Top squark - neutralino mass points lying on the observed+$1\sigma_\text{observed}$ exclusion contour (for $ x=0.5 $).

Expected $95\%$ CL upper limit on the production cross section of top squarks (for $ x=0.75 $).

Top squark - neutralino mass points lying on the expected exclusion contour (for $ x=0.75 $).

Top squark - neutralino mass points lying on the expected-$1\sigma_\text{expected}$ exclusion contour (for $ x=0.75 $).

Top squark - neutralino mass points lying on the expected+$1\sigma_\text{expected}$ exclusion contour (for $ x=0.75 $).

Observed $95\%$ CL upper limit on the production cross section of top squarks (for $ x=0.75 $).

Top squark - neutralino mass points lying on the observed exclusion contour (for $ x=0.75 $).

Top squark - neutralino mass points lying on the observed-$1\sigma_\text{observed}$ exclusion contour (for $ x=0.75 $).

Top squark - neutralino mass points lying on the observed+$1\sigma_\text{observed}$ exclusion contour (for $ x=0.75 $).