Bayesian upper exclusion limits on the ratio $g_{W'}/g_{W}$ for an SSM-like W$\prime$ boson are shown. The unity coupling ratio (blue dotted curve) corresponds to the SSM common benchmark. The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian lower exclusion limits on the NUGIM G(221) mixing angle $\cot(\theta_{E})$ are shown as a function of the W$\prime$ boson mass. The theoretically excluded region is shaded in grey. The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian upper exclusion limits at 95% CL on the product of the production cross section and branching fraction of a QBH in an associated $ au$ lepton and neutrino final state. Minimum threshold masses $m_{th}$ of up to 6.6 TeV are excluded at 95% CL. The observed limit (solid line) is obtained from the intersection with the LO QBH cross section (blue dotted curve). The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian upper limits at 95% CL on the cross section of the process $pp\rightarrow\tau\nu$ mediated via LQ exchange in the t-channel. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The predicted LQ cross section at LO in the three coupling benchmark scenarios is depicted in different colors for $g_{U}=1$. The uncertainty bandy correspond to the sum in quadrature of PDF and scale variations. The first benchmark scenario considers only couplings to left-handed SM fermions (i.e. $\beta_{\text{R}}^{ij}=0$) and is referred to as "best fit LH". The second benchmark, referred to as "best fit LH+RH", considers $|\beta_{\text{R}}^{\text{b}\tau}|=1$ and all other $\beta_{\text{R}}^{ij}=0$. In the third "democratic" benchmark, equal couplings only to LH fermions are assumed, i.e. $\beta_{\text{L}}^{ij}=1$ and $\beta_{\text{R}}^{ij}=0$ for all $i$ and $j$.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the LH scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the LH+RH scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the democratic scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Bayesian upper exclusion limits at 95% CL on each of the Wilson coefficients described by the EFT model based on 2016-2018 data. The three different coupling types represent a left-handed vector coupling ($\epsilon^{cb}_{L}$), tensor-like coupling ($\epsilon^{cb}_{T}$), and scalar-tensor-like coupling ($\epsilon^{cb}_{S_{L}}$). The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Summary of exclusion limits (expected and observed) calculated at 95% CL for full Run-2 CMS data.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "year_binnumber", following the binning from Figure 4.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "year_binnumber", following the binning used in Figure 4.

Predicted signal yields for the 2017 data-taking period, corresponding to 41.3 fb$^{-1}$, after the application of each search requirement (cumulative) for various signal hypothesis. The requirements listed are presented as total efficiencies w.r.t. the previous selection step.

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

Experimental Xi values and FJF predictions for the four NRQCD terms using BK LDME parameters

Experimental Xi values and FJF predictions for the four NRQCD terms using BCKL LDME parameters

Experimental Xi values and FJF predictions for the four NRQCD terms using BK LDME parameters

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of Z pt on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of Z pt on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of |y(Z)| on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of |y(Z)| on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of $\phi^{\scriptscriptstyle *}_\eta$ on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of $\phi^{\scriptscriptstyle *}_\eta$ on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of Z pt on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of Z pt on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of |y(Z)| on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of |y(Z)| on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of $\phi^{\scriptscriptstyle *}_\eta$ on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of $\phi^{\scriptscriptstyle *}_\eta$ on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

The measured absolute cross sections in bins of |y(Z)|, using dressed level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of Z pt, using dressed level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of $\phi^{\scriptscriptstyle *}_\eta$, using dressed level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of Z pt different |y(Z)| bins, using dressed level leptons. The cross sections are normalized by the bin width.

The measured normalized cross sections in bins of Z pt, using dressed level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of $\phi^{\scriptscriptstyle *}_\eta$, using dressed level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of |y(Z)|, using dressed level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of Z pt in |y(Z)| bins, using dressed level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of |y(Z)| in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of |y(Z)| in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of |y(Z)| in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of |y(Z)| in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of |y(Z)| in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of |y(Z)| in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of |y(Z)| in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of |y(Z)| in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

The measured absolute cross sections in bins of |y(Z)|, using born level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of Z pt, using born level leptons. The cross sections are normalized by the bin width. The first bin (0-1 GeV) is not shown as large differences were observed in aMC@NLO and POWHEG predictions at Born level.

The measured absolute cross sections in bins of $\phi^{\scriptscriptstyle *}_\eta$, using born level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of Z pt different |y(Z)| bins, using born level leptons. The cross sections are normalized by the bin width. The first bin (0-1 GeV) is not shown as large differences were observed in aMC@NLO and POWHEG predictions at Born level.

The measured normalized cross sections in bins of Z pt, using born level leptons. The first bin (0-1 GeV) is not shown as large differences were observed in aMC@NLO and POWHEG predictions at Born level. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of $\phi^{\scriptscriptstyle *}_\eta$, using born level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of |y(Z)|, using born level leptons.

The measured normalized cross sections (left) in bins of Z pt in |y(Z)| bins, using born level leptons. The first bin (0-1 GeV) is not shown as large differences were observed in aMC@NLO and POWHEG predictions at Born level. The cross sections are divided by the bin width and then the distributions are normalized to 1.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of |y(Z)| in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of |y(Z)| in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of |y(Z)| in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of |y(Z)| in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dilepton final state.

Elliptic-flow coefficients $v_2$ based on the six-particle correlations technique, as functions of transverse momentum and in bins of centrality. The results correspond to the range $|\eta| < 2.4$.

Elliptic-flow coefficients $v_2$ based on the eight-particle correlations technique, as functions of transverse momentum and in bins of centrality. The results correspond to the range $|\eta| < 2.4$.

Triangular-flow coefficients $v_3$ based on the two-particle correlations technique, as functions of transverse momentum and in bins of centrality. The results correspond to the range $|\eta| < 2.4$.

Triangular-flow coefficients $v_3$ based on the scalar-product technique, as functions of transverse momentum and in bins of centrality. The results correspond to the range $|\eta| < 0.8$.

Triangular-flow coefficients $v_3$ based on the four-particle correlations technique, as functions of transverse momentum and in bins of centrality. The results correspond to the range $|\eta| < 2.4$.

The $v_4$ coefficients based on the two-particle correlations technique, as functions of transverse momentum and in bins of centrality. The results correspond to the range $|\eta| < 2.4$.

The $v_4$ coefficients based on the scalar-product technique, as functions of transverse momentum and in bins of centrality. The results correspond to the range $|\eta| < 0.8$.

Centrality dependence of the spectrum-weighted $v_2$ flow harmonics with $0.3 < p_{\mathrm{T}} < 3.0~\mathrm{GeV}/c$. The $v_2$ results are shown for two-, four-, six-, and eight-particle correlations.

Centrality dependence of the spectrum-weighted $v_3$ flow harmonics with $0.3 < p_{\mathrm{T}} < 3.0~\mathrm{GeV}/c$. The results are shown for two- and four-particle correlations.

Centrality dependence of the spectrum-weighted $v_4$ flow harmonics with $0.3 < p_{\mathrm{T}} < 3.0~\mathrm{GeV}/c$. The results are shown for two-particle correlations.

Centrality dependence of $v_2\{4\}/v_2\{2\}$ ratios.

Centrality dependence of $v_2\{6\}/v_2\{4\}$ ratios.

Centrality dependence of $v_3\{4\}/v_3\{2\}$ ratios.

The $v_2$ results measured with two-particle correlations from PbPb collisions at $5.02~$TeV, shown as a function of $p_{\mathrm{T}}$ in eleven centrality bins.

The $v_3$ results measured with two-particle correlations from PbPb collisions at $5.02~$TeV, shown as a function of $p_{\mathrm{T}}$ in eleven centrality bins.

The $v_4$ results measured with two-particle correlations from PbPb collisions at $5.02~$TeV, shown as a function of $p_{\mathrm{T}}$ in eleven centrality bins.

Ratios of the $v_2$ harmonic coefficients from two-particle correlations in XeXe and PbPb collisions as functions of $p_{\mathrm{T}}$ in 11 centrality bins.

Ratios of the $v_3$ harmonic coefficients from two-particle correlations in XeXe and PbPb collisions as functions of $p_{\mathrm{T}}$ in 11 centrality bins.

Ratios of the $v_4$ harmonic coefficients from two-particle correlations in XeXe and PbPb collisions as functions of $p_{\mathrm{T}}$ in 11 centrality bins.

Centrality dependence of the spectrum-weighted $v_2$, $v_3$, and $v_4$ harmonic coefficients from two-particle correlations method for $0.3 < p_{\mathrm{T}} < 3.0 \mathrm{GeV}/c$ for PbPb collisions at $5.02$~TeV.

Ratios of the $v_2$, $v_3$, and $v_4$ harmonic coefficients from two-particle correlations in XeXe and PbPb collisions as functions or $0.3 < p_{\mathrm{T}} < 3.0~\mathrm{GeV}/c$ as a function of centrality.

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.

Charged particle yield distributions $Y(\Delta r)$ of b jets with $1 < p_{\text{T}}^{\text{trk}} < 12$ GeV are presented as functions of $\Delta r$.b jets with $p_T > 120$ GeV and charged particles with $1 < p^{\text{trk}}_{\text{T}} < 12$ GeV are used to construct the distributions as functions of $\Delta r$ differential $p_{\text{T}}^{\text{trk}}$ bins.

Charged particle yield differences of b to inclusive jets for $1 < p_{\text{T}}^{\text{trk}} < 12$ GeV are presented as functions of $\Delta r$.b jets with $p_T > 120$ GeV and charged particles with $1 < p^{\text{trk}}_{\text{T}} < 12$ GeV are used to construct the distributions as functions of $\Delta r$ differential $p_{\text{T}}^{\text{trk}}$ bins.

The jet shapes $\rho(\Delta r)$ of inclusive jets with $p_{\text{T}} > 120$ GeV and $p^{trk}_T > 1$ GeV are presented as functions of $\Delta r$ for differential $p_{\text{T}}^{\text{trk}}$ bin.

The jet shape distributions $\rho(\Delta r)$ of b jets with $p_{\text{T}} > 120$ GeV and $1< p^{\text{trk}}_{\text{T}} < 12$ GeV are presented as functions of $\Delta r$ for differential $p_{\text{T}}^{\text{trk}}$ bin.

The jet shapes $\rho (\Delta r)$ of inclusive jets with $p_{\text{T}} > 120$ GeV and $p^{\text{trk}}_{\text{T}} > 1$ GeV are presented as functions of $\Delta r.

The jet shapes $\rho (\Delta r)$ of b jets with $p_{\text{T}} > 120$ GeV and $p^{\text{trk}}_{\text{T}} > 1$ GeV are presented as functions of $\Delta r.

The b-to-inclusive jet shape ratios for $p_{\text{T}} > 120$ GeV and $p^{\text{trk}}_{\text{T}} > 1$ GeV are presented as functions of $\Delta r.