The spin transfer $D_{LL}$ to (a) $\Lambda$ and (b) $\bar{\Lambda}$ hyperons produced at positive pseudorapidity with respect to the polarized proton beam from $MB$, $JP$, and $HT$ data versus hyperon transverse momenta $p_{T}$. The sizes of the statistical and systematic uncertainties are indicated by the vertical bars and bands, respectively. For clarity, the HT data points have been shifted slightly in $p_{T}$. The dotted vertical lines indicate the $p_{T}$ intervals in the analysis of HT and JP data.

Comparison of $\Lambda$ and $\bar{\Lambda}$ spin transfer $D_{LL}$ in polarized proton-proton collisions at $\sqrt{s} = 200 GeV$ for (a) positive and (b) negative $\eta$ versus $p_{T}$. The vertical bars and bands indicate the sizes of the statistical and systematic uncertainties, respectively. The $\bar{\Lambda}$ data points have been shifted slightly in $p_{T}$ for clarity. The dotted vertical lines indicate the $p_{T}$ intervals in the analysis of HT and JP data. The horizontal lines show model predictions evaluated at $\eta$ and largest $p_{T}$ of the data.

Comparison of $\Lambda$ and $\bar{\Lambda}$ spin transfer $D_{LL}$ in polarized proton-proton collisions at $\sqrt{s} = 200$ GeV for (a) positive and (b) negative $\eta$ versus $p_{T}$. The vertical bars and bands indicate the sizes of the statistical and systematic uncertainties, respectively. The $\bar{\Lambda}$ data points have been shifted slightly in $p_{T}$ for clarity. The dotted vertical lines indicate the $p_{T}$ intervals in the analysis of HT and JP data. The horizontal lines show model predictions evaluated at $\eta$ and largest $p_{T}$ of the data.

The spin transfer $D_{LL}$ to (a) $\Lambda$ and (b) $\bar{\Lambda}$ hyperons, calculated using the new $\alpha_\Lambda=0.732$ updated on PDG2020, produced at positive pseudorapidity with respect to the polarized proton beam from $MB$, $JP$, and $HT$ data versus hyperon transverse momenta $p_{T}$. The sizes of the statistical and systematic uncertainties are indicated by the vertical bars and bands, respectively. For clarity, the HT data points have been shifted slightly in $p_{T}$. The dotted vertical lines indicate the $p_{T}$ intervals in the analysis of HT and JP data.

The spin transfer $D_{LL}$ to (a) $\Lambda$ and (b) $\bar{\Lambda}$ hyperons, calculated using the new $\alpha_\Lambda=0.732$ updated on PDG2020, produced at positive pseudorapidity with respect to the polarized proton beam from $MB$, $JP$, and $HT$ data versus hyperon transverse momenta $p_{T}$. The sizes of the statistical and systematic uncertainties are indicated by the vertical bars and bands, respectively. For clarity, the HT data points have been shifted slightly in $p_{T}$. The dotted vertical lines indicate the $p_{T}$ intervals in the analysis of HT and JP data.

Comparison of $\Lambda$ and $\bar{\Lambda}$ spin transfer $D_{LL}$, calculated using the new $\alpha_\Lambda=0.732$ updated on PDG2020, in polarized proton-proton collisions at $\sqrt{s} = 200$ GeV for (a) positive and (b) negative $\eta$ versus $p_{T}$. The vertical bars and bands indicate the sizes of the statistical and systematic uncertainties, respectively. The $\bar{\Lambda}$ data points have been shifted slightly in $p_{T}$ for clarity. The dotted vertical lines indicate the $p_{T}$ intervals in the analysis of HT and JP data. The horizontal lines show model predictions evaluated at $\eta$ and largest $p_{T}$ of the data.

Observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance. The acceptance $\mathcal{A}$ is defined as $\mathcal{A} = N$(events with $m_{X}^{GEN} > 85\% m_{X}^{input}$) / $N$(events generated in the full phase space defined by the CMS default generator settings). Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Mass exclusion ranges of the benchmark signal scenarios are depicted with hatched areas inside the black contours.

Observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance. The acceptance $\mathcal{A}$ is defined as $\mathcal{A} = N$(events with $m_{X}^{GEN} > 85\% m_{X}^{input}$) / $N$(events generated in the full phase space defined by the CMS default generator settings). Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Mass exclusion ranges of the benchmark signal scenarios are depicted with hatched areas inside the black contours.

Efficiencies of the selection requirements on the benchmark signal processes: $Z_{R} \to ggg$ with nominal width ($\Gamma_{X}/m_{X}\sim 3\%$). A value of -1 means that the corresponding efficiency is not calculated for this year.

Efficiencies of the selection requirements on the benchmark signal processes: $Z_{R} \to ggg$ with narrow width ($\Gamma_{X}/m_{X}\sim 0.01\%$). A value of -1 means that the corresponding efficiency is not calculated for this year.

Efficiencies of the selection requirements on the benchmark signal processes: $G_{KK} \to {\varphi}(gg)g$, where $\varphi$ is the radion. A value of -1 means that the corresponding efficiency is not calculated for this year.

Efficiencies of the selection requirements on the benchmark signal processes: $q^{*} \to V(qq)q$, where $V$ is a beyond-the-SM vector boson. A value of -1 means that the corresponding efficiency is not calculated for this year.

Acceptance of the signal selection requirement $m_{X}^{\text{GEN}}/m_{X}^{\text{input}} > 85\%$ on the benchmark signal process $Z_{R} \to ggg$. The acceptance is defined as $\mathcal{A} = N$(events with $m_{X}^{\text{GEN}}/m_{X}^{\text{input}} > 85\%$)/$N$(events generated in the full phase space defined by the CMS default generator settings).

Acceptance of the signal selection requirement $m_{X}^{\text{GEN}}/m_{X}^{\text{input}} > 85\%$ on the benchmark signal process $G_{KK} \to {\varphi}(gg)g$.

Acceptance of the signal selection requirement $m_{X}^{\text{GEN}}/m_{X}^{\text{input}} > 85\%$ on the benchmark signal process $q^{*} \to V(qq)q$.

Observed local significance for a 3-body decay $ggg$ resonance, shown for resonances with nominal width (blue solid line) and narrow width (red dashed line).

Observed local significance for a cascade decay $ggg$ resonance.

Observed local significance for a cascade decay $qqq$ resonance.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.2$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.3$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.4$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.5$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.6$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.7$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.8$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance $m_{Y} / m_{X} = 0.2$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.3$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.4$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.5$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.6$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.7$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.8$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal

Cut flow table of the selection requirments for the $Z_{R}$ model scenarios. Values shown are absolute efficiencies, e.g., values shown in the column of 'Efficiency ($\Delta_{R}^{max} < 3.0$)' represent the cumulative efficiencies achieved through all selection requirements applied up to and including the current selection criterion.

Cut flow table of the selection requirments for the $G_{KK}$ model scenarios. Values shown are absolute efficiencies, e.g., values shown in the column of 'Efficiency ($\Delta_{R}^{max} < 3.0$)' represent the cumulative efficiencies achieved through all selection requirements applied up to and including the current selection criterion.

Cut flow table of the selection requirments for the $q^{*}$ model scenarios. Values shown are absolute efficiencies, e.g., values shown in the column of 'Efficiency ($\Delta_{R}^{max} < 3.0$)' represent the cumulative efficiencies achieved through all selection requirements applied up to and including the current selection criterion.

Predicted production cross section of the benchmark signal process $pp \to Z_{R} \to ggg$.

Predicted production cross section of the benchmark signal process $pp \to G_{KK} \to \varphi(gg)g$.

Predicted production cross section of the benchmark signal process $pp \to q^{*} \to V(qq)q$.

Correlation matrix for the differential branching fraction extraction between different $q^2$ bins in the simultaneous fit.

The product of acceptance and efficiency ($A\epsilon$) of the $\mathcal{B}(\mathrm{B}^+\to\mathrm{K}^+\mu^+\mu^-)$ channel, as a function of the muon pair $q^2$, as measured in simulated signal events, after all the corrections applied. In the figure, regions corresponding to resonances are displayed with red markers.

Same as Table 5, but for the bins containing the J$\psi$ and $\psi$(2S) resonances.

Integrated branching fraction $\mathcal{B}(\mathrm{B}^\pm \to \mathrm{K}^\pm\mu^+\mu^-)$ in the low-$q^2$ region.

Ratio of branching fractions $\mathcal{B}(\text{B}^\pm \to \text{K}^\pm\mu^+\mu^-)$ to $\mathcal{B}(\text{B}^\pm \to \text{K}^\pm\text{e}^+\text{e}^-)$, determined from the measured double ratio $R$(K) of these decays to the respective branching fractions of the decay chains $\text{B}^\pm\to\text{J/}\psi\text{K}^\pm$ with $\text{J/}\psi\to\mu^+\mu^-$ and $\text{e}^+\text{e}^-$.

Negative log likelihood function from the fit profiled as a function of the inverse ratio $R(\mathrm{K})^{-1}$.

Postfit score of the $\text{t}\bar{\text{t}}\text{Z}$ output node from the multiclass classifier in $\text{SR}_\text{3l,3j}$ for events with at least 2 b jets

Postfit score of the tWZ output node from the binary classifier in $\text{SR}_\text{3l,2j}$

$\text{p}_\text{T}$ of the leading jet in the $\text{SR}_\text{3l,2j}$

$\text{p}_\text{T}$ of the leading lepton in the $\text{SR}_\text{3l,2j}$

Max $\text{p}_\text{T}$ of the lepton-jet in the $\text{SR}_\text{3l,3j}$

Invariant mass of the system of leptons, jets, and $\vec{\text{p}}^\text{miss}_\text{T}$ in the $\text{SR}_\text{3l,3j}$

Two-dimensional likelihood scan of the signal strenghts of tWZ ($\mu_\text{tWZ}$) vs. $\text{t}\bar{\text{t}}\text{Z}$ ($\mu_{\text{t}\bar{\text{t}}\text{Z}}$). The blue lines show the 68%, 95%, and 99% confidence level contours. The black cross represents the best-fit value, while the red diamond the SM expectation

Two-dimensional likelihood scan of the signal strenghts of tWZ ($\mu_\text{tWZ}$) vs. $\text{t}\bar{\text{t}}\text{Z}$ ($\mu_{\text{t}\bar{\text{t}}\text{Z}}$). The blue lines show the 68%, 95%, and 99% confidence level contours. The black cross represents the best-fit value, while the red diamond the SM expectation

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{1})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}^{\mathrm{add.}}_{2})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{2})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{add.}})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\Delta\mathrm{R}_{\mathrm{b}\mathrm{b}}^{\mathrm{avg}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}_{3})|$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}_{3})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}_{3})$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}_{3})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}_{4})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}_{4})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{b}}_{\mathrm{T}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{b}}_{\mathrm{T}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $|\eta(\mathrm{b}^{\mathrm{extra}}_{1})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{1})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}^{\mathrm{extra}}_{2})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{2})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{extra}})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $H^{\mathrm{j}}_{\mathrm{T}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{j}}_{\mathrm{T}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\mathrm{m}_{\mathrm{b}\mathrm{b}}^{\mathrm{max}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{light}}_{\mathrm{T}}$ in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $H^{\mathrm{light}}_{\mathrm{T}}$ in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $N_{\mathrm{jets}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $N_{\mathrm{jets}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $N_{\mathrm{b}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{add.}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\Delta\mathrm{R}_{\mathrm{b}\mathrm{b}}^{\mathrm{avg}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}_{4})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}_{4})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{extra}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\mathrm{m}_{\mathrm{b}\mathrm{b}}^{\mathrm{max}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $N_{\mathrm{jets}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $N_{\mathrm{jets}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $N_{\mathrm{b}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{add.}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\Delta\mathrm{R}_{\mathrm{b}\mathrm{b}}^{\mathrm{avg}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}_{4})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}_{4})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{extra}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\mathrm{m}_{\mathrm{b}\mathrm{b}}^{\mathrm{max}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $N_{\mathrm{jets}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $N_{\mathrm{jets}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $N_{\mathrm{b}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Fiducial cross sections from the measurements of all observables, compared to predictions from different ttbb simulation approaches. For each of the normalized differential measurements the fiducial cross section in the respective phase space is also determined. In the paper only one representative observable is quoted for each fiducial phase space, while here the measured cross section with the uncertainties from the fit to the respective observable is summarized.

Correlation matrix associated to the statistical covariance matrix of the data and MC for the primary Lund jet plane density for AK8 jets in a one-dimensional representation with bin indices. The mapping between the bin indices and the physical binning can be imported from the XML file attached to this HepData record using the TUnfoldBinningXML class of ROOT (qualitatively, it corresponds to slicing the Lund plane horizontally from low kT to high kT).

Observed profiled likelihood on $f_{\Lambda1}$ using Approach 1. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{a3}$ using Approach 1. The signal strength modifiers are treated as free parameters.

Observed profiled likelihood on $f_{a3}$ using Approach 1. The signal strength modifiers are treated as free parameters.

Expected profiled likelihood on $f_{\Lambda1}^{Z\gamma}$ using Approach 1. The signal strength modifiers are treated as free parameters.

Observed profiled likelihood on $f_{\Lambda1}^{Z\gamma}$ using Approach 1. The signal strength modifiers are treated as free parameters.

Expected profiled likelihood on $f_{a2}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Observed profiled likelihood on $f_{a2}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{\Lambda1}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Observed profiled likelihood on $f_{\Lambda1}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{a3}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters.

Observed profiled likelihood on $f_{a3}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters.

Expected profiled likelihood on $f_{a2}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Observed profiled likelihood on $f_{a2}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{\Lambda1}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Observed profiled likelihood on $f_{\Lambda1}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{a3}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters.

Observed profiled likelihood on $f_{a3}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters.

Expected profiled likelihood on the $\delta c_{Z}$ coupling of the SMEFT Higgs basis.

Observed profiled likelihood on the $\delta c_{Z}$ coupling of the SMEFT Higgs basis.

Expected profiled likelihood on the $c_{Z\Box}$ coupling of the SMEFT Higgs basis.

Observed profiled likelihood on the $c_{Z\Box}$ coupling of the SMEFT Higgs basis.

Expected profiled likelihood on the $c_{ZZ}$ coupling of the SMEFT Higgs basis.

Observed profiled likelihood on the $c_{ZZ}$ coupling of the SMEFT Higgs basis.

Expected profiled likelihood on the $\tilde{c}_{ZZ}$ coupling of the SMEFT Higgs basis.

Observed profiled likelihood on the $\tilde{c}_{ZZ}$ coupling of the SMEFT Higgs basis.

Expected profiled likelihood on $f_{a3}^{ggH}$. The signal strength modifiers and the CP-odd HVV anomalous coupling cross section fraction are treated as free parameters.

Observed profiled likelihood on $f_{a3}^{ggH}$. The signal strength modifiers and the CP-odd HVV anomalous coupling cross section fraction are treated as free parameters.

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.33, 0.33, 0.33) and in the Majorana scenario.

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.0, 1.0, 0.0) and in the Dirac-like scenario.

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.0, 0.5, 0.5) and in the Dirac-like scenario.

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.5, 0.5, 0.0) and in the Dirac-like scenario.

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.33, 0.33, 0.33) and in the Dirac-like scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 1.0 GeV in the Majorana scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 1.5 GeV in the Majorana scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 2.0 GeV in the Majorana scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 1.0 GeV in the Dirac-like scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 1.5 GeV in the Dirac-like scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 2.0 GeV in the Dirac-like scenario.

Distributions of $S_{\mathrm{ML}}$ for data, simulated background and signal events with $n_{\mathrm{track}}$ of 4. The distributions are shown for split-SUSY signals with a gluino mass of 2000 GeV and neutralino mass of 1800 GeV. Different gluino proper decay lengths are shown as $c\tau$ in the legend. All distributions are normalized to unity.

Distributions of $S_{\mathrm{ML}}$ for data, simulated background and signal events with $n_{\mathrm{track}}$ of $\geq$ 5. The distributions are shown for split-SUSY signals with a gluino mass of 2000 GeV and neutralino mass of 1900 GeV. Different gluino proper decay lengths are shown as $c\tau$ in the legend. All distributions are normalized to unity.

Distributions of $S_{\mathrm{ML}}$ for data, simulated background and signal events with $n_{\mathrm{track}}$ of $\geq$ 5. The distributions are shown for split-SUSY signals with a gluino mass of 2000 GeV and neutralino mass of 1800 GeV. Different gluino proper decay lengths are shown as $c\tau$ in the legend. All distributions are normalized to unity.

The distribution of $n_{\mathrm{track}}$ in different $S_{\mathrm{ML}}$ regions for simulated background events. Events with 0 $ < S_{\mathrm{ML}} < $ 0.2 (blue), 0.2 $ < S_{\mathrm{ML}} < $ 0.6 (red), and 0.6 $ < S_{\mathrm{ML}} < $ 1.0 (green) are compared. All distributions are normalized to unity. The similar $n_{\mathrm{track}}$ distributions demonstrate that $n_{\mathrm{track}}$ and $S_{\mathrm{ML}}$ are decorrelated.

The distribution of $d_{\mathrm{BV}}$ in $K_{\mathrm{S}}^{\mathrm{0}}$ vertices in data (black) and simulation (purple). The lower panel shows the ratio between data and simulation.

The vertex reconstruction efficiency for artificially displaced vertices in data (black) and simulation (red). In this example, the artificially displaced vertices are corrected to mimic split-SUSY signal events with gluino mass of 2000 GeV and neutralino mass of 1800 GeV.

The ML tagging efficiency for artificially displaced vertices in data (black) and simulation (red). In this example, the artificially displaced vertices are corrected to mimic split-SUSY signal events with gluino mass of 2000 GeV and neutralino mass of 1800 GeV.

Number of predicted and observed events in the control, validation, and search regions. Predictions are calculated using Eqs. (2) and (3) and fitting the data under the background-only hypothesis. Regions are organized by $S_{\mathrm{ML}}$ and $n_{\mathrm{track}}$ values, and region names corresponding with Fig. 7 are given in parentheses. The predicted number of events that pass the $S_{\mathrm{ML}}$ selection and the observed number of events that pass or fail the $S_{\mathrm{ML}}$ selection are shown in seperate rows.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the split-SUSY signal model with a mass splitting of 100 GeV, shown as a function of gluino mass and $c\tau$. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the split-SUSY signal model with a mass splitting of 100 GeV, shown as a function of gluino mass and $c\tau$. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the split-SUSY model with a $c\tau$ of 10 mm, shown as a function of gluino mass and mass splitting. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the split-SUSY model with a $c\tau$ of 10 mm, shown as a function of gluino mass and mass splitting. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the GMSB SUSY signal model, shown as a function of gluino mass and $c\tau$. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the GMSB SUSY signal model, shown as a function of gluino mass and $c\tau$. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

Observed upper limits at 95% CL on B($\text{H} \rightarrow \text{a}_{1}\text{a}_{1}$) in percent, obtained from the combination of the $\mu\mu$bb and $\tau\tau$bb channels. The results are obtained as functions of $m_{\text{a}_{1}}$ for 2HDM+S Type I (independent of tan$\beta$), Type II (tan$\beta$=2.0), Type III (tan$\beta$=2.0), and Type IV (tan$\beta$=0.6), respectively.

Unfolded E2C distributions in data compared to MC predictions.

Unfolded E2C distributions in data compared to MC predictions.

Unfolded E2C distributions in data compared to MC predictions.

Unfolded E2C distributions in data compared to MC predictions.

Unfolded E2C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

The fitted slopes of the E3C/E2C data distributions as a function of jet pt are used to illustrate the dependency of alphas on jet pt.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

The chi2 scan result using eq.3 for different alphas.

The correlation matrix is composed by 10 jet pt region, each region is represented by a block in the plot. Inside each block, there are 22 xL bins same as the E2C, E3C and E3C/E2C distributions. Therefore, the x and y bins of the correlation matrix is given by, binNumber = pT_index * 22 + xL_index.

The correlation matrix is composed by 10 jet pt region, each region is represented by a block in the plot. Inside each block, there are 22 xL bins same as the E2C, E3C and E3C/E2C distributions. Therefore, the x and y bins of the correlation matrix is given by, binNumber = pT_index * 22 + xL_index.

The correlation matrix is composed by 10 jet pt region, each region is represented by a block in the plot. Inside each block, there are 22 xL bins same as the E2C, E3C and E3C/E2C distributions. Therefore, the x and y bins of the correlation matrix is given by, binNumber = pT_index * 22 + xL_index.

The energy weight of E2C as defined in eq.1 in the paper is used in the unfolding, the binning is listed in this table.

The energy weight of E3C as defined in eq.1 in the paper is used in the unfolding, the binning is listed in this table.

Upper 95% CL limits on HH signal cross section scanned over the $\kappa_{\lambda}$ parameter while fixing the $\kappa_{2V}$ and $\kappa_{V}$ to their SM-predicted values.

Upper 95% CL limits on VHH signal cross section scanned over the $\kappa_{2V}$ parameter while fixing the $\kappa_{\lambda}$ and $\kappa_{V}$ to their SM-predicted values.

Upper 95% CL limits on HH signal cross section scanned over the $\kappa_{2V}$ parameter while fixing the $\kappa_{\lambda}$ and $\kappa_{V}$ to their SM-predicted values.

Upper 95% CL limits on VHH signal cross section scanned over the $\kappa_{V}$ parameter while fixing the $\kappa_{2V}$ and $\kappa_{\lambda}$ to their SM-predicted values.

Upper 95% CL limits on HH signal cross section scanned over the $\kappa_{V}$ parameter while fixing the $\kappa_{2V}$ and $\kappa_{\lambda}$ to their SM-predicted values.

Differential cross sections normalized to the fiducial cross section as a function of the $p_T$ of the second-highest-$p_T$ jet

Differential cross sections normalized to the fiducial cross section as a function of the $|\eta|$ of the highest-$p_T$ jet

Differential cross sections normalized to the fiducial cross section as a function of the $|\eta|$ of the second-highest-$p_T$ jet

Differential cross sections normalized to the fiducial cross section as a function of the pseudorapidity difference between the two highest-$p_T$ jets in the event

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the two highest-$p_T$ jets in the event

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system, in the full $m_{4l}$ range

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with 0 jet, in the on-shell ZZ region

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with 1 jet, in the on-shell ZZ region

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with 2 jets, in the on-shell ZZ region

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with at least 3 jets, in the on-shell ZZ region

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with 0 jet, in the full $m_{4l}$ range

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with 1 jet, in the full $m_{4l}$ range

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with 2 jets, in the full $m_{4l}$ range

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with at least 3 jets, in the full $m_{4l}$ range

Invariant mass distributions of the $t\to jj^{'}b$ candidates, $m_{top}$, in the 0 b category after all selections. The error bars indicate the statistical uncertainties.

Invariant mass distributions of the $t\to jj^{'}b$ candidates, $m_{top}$, in the 1 b category after all selections. The error bars indicate the statistical uncertainties.

Invariant mass distributions of the $t\to jj^{'}b$ candidates, $m_{top}$, in the $\geq$2 b category after all selections. The error bars indicate the statistical uncertainties.

The measured cross section in the $\mu$+jets, $\rm{e}$+jets and $\ell$+jets channels. The first uncertainty is the statistical, the second is the systematic.

The observed distribution of the dijet invariant mass in the $\ell$+jets decay channel in the 0 b tagged jet category. The vertical bars on the data points represent the statistical uncertainties. The integral in the $\rm{e}$($\mu$)+jets channel (including overflow events) is 1375 (2406).

The observed distribution of the dijet invariant mass in the $\ell$+jets decay channel in the 1 b tagged jet category. The vertical bars on the data points represent the statistical uncertainties. The integral in the $\rm{e}$($\mu$)+jets channel (including overflow events) is 61 (129).

The observed distribution of the dijet invariant mass in the $\ell$+jets decay channel in the $\geq$2 b tagged jet category. The vertical bars on the data points represent the statistical uncertainties. The integral in the $\rm{e}$($\mu$)+jets channel (including overflow events) is 19 (33).

The observed distribution of the the minimum angular distance $\Delta R$ between all pairs of jets $j$ and $j'$ in the $\ell$+jets decay channel in the 0 b tagged jet category. The vertical bars on the data points represent the statistical uncertainties. The integral in the $\rm{e}$($\mu$)+jets channel (including overflow events) is 1375 (2406).

The observed distribution of the the minimum angular distance $\Delta R$ between all pairs of jets $j$ and $j'$ in the $\ell$+jets decay channel in the 1 b tagged jet category. The vertical bars on the data points represent the statistical uncertainties. The integral in the $\rm{e}$($\mu$)+jets channel (including overflow events) is 61 (129).

The observed distribution of the the minimum angular distance $\Delta R$ between all pairs of jets $j$ and $j'$ in the $\ell$+jets decay channel in the $\geq$2 b tagged jet category. The vertical bars on the data points represent the statistical uncertainties. The integral in the $\rm{e}$($\mu$)+jets channel (including overflow events) is 19 (33).

The observed distribution of the jet multiplicity in the $\rm{e}^\pm \mu^\mp$ decay channel for events passing the dilepton criteria. The error bars indicate the statistical uncertainties. The last bin of the distributions contains the overflow events.

The observed distribution of the scalar $p_{\rm{T}}$ sum of all jets in the $\rm{e}^\pm \mu^\mp$ decay channel for events passing the dilepton criteria. The error bars indicate the statistical uncertainties. The last bin of the distributions contains the overflow events.

The observed distribution of the invariant mass of the lepton pair in the $\rm{e}^\pm \mu^\mp$ decay channel after requiring at least two jets. The error bars indicate the statistical uncertainties. The integral is 24. The last bin of the distributions contains the overflow events.

The observed distribution of the $p_{\rm{T}}$ of the lepton pair in the $\rm{e}^\pm \mu^\mp$ decay channel after requiring at least two jets. The error bars indicate the statistical uncertainties. The integral is 24. The last bin of the distributions contains the overflow events.

The observed distribution of the $p_{\rm{T}}^{\rm{miss}}$ in the $\mu^\pm \mu^\mp$ decay channel for events passing the dilepton criteria and Z boson veto. The error bars indicate the statistical uncertainties. The last bin of the distributions contains the overflow events.

The observed distribution of the invariant mass of the lepton pair in the $\mu^\pm \mu^\mp$ decay channel for events passing the dilepton criteria, Z boson veto, and $p_{\rm{T}}^{\rm{miss}}>35$ GeV requirement. The error bars indicate the statistical uncertainties. The last bin of the distributions contains the overflow events.

x_J distribution of incluisve dijets in pp collisions

x_J distribution of b-quark dijets in pp collisions

x_J distribution of incluisve dijets in PbPb collisions (0-10%)

x_J distribution of b-quark dijets in PbPb collisions (0-10%)

x_J distribution of incluisve dijets in PbPb collisions (10-30%)

x_J distribution of b-quark dijets in PbPb collisions (10-30%)

x_J distribution of incluisve dijets in PbPb collisions (30-100%)

x_J distribution of b-quark dijets in PbPb collisions (30-100%)

Mean x_J of incluisve dijets in pp and PbPb collisions vs. number of participants

Mean x_J of b-quark dijets in pp and PbPb collisions vs. numberof participants

Mean x_J of inclusive jets and b-quark dijets in PbPb collisions minus mean x_J in pp collisions vs. numberof participants

Distribution of the transverse momentum ratio between the Z boson and the jet in PbPb collisions. The distribution is normalized by the number of Z events.

The mean value of the X_JZ distribution in pp collisions. The resolutions of the measured jet energy and azimuthal angle are smeared to match those of the PbPb sample.

The mean value of the X_JZ distribution in PbPb collisions.

The average number of jet partners per Z boson in pp collisions. The resolutions of the measured jet energy and azimuthal angle are smeared to match those of the PbPb sample.

The average number of jet partners per Z boson in PbPb collisions.

Distribution of the transverse momentum ratio between the Z boson and the jet in pp collisions. The distribution is normalized by the number of Z events. No resolution smearing is applied to the measured jet energy and azimuthal angle.

Total excess correlated yield observed in the PbPb data with respect to the reference measured in pp collisions, shown as a function of track pT in in centrality interval 10-30% for both leading jets with pT >120GeV and subleading jets with pT > 50 GeV.

Total excess correlated yield observed in the PbPb data with respect to the reference measured in pp collisions, shown as a function of track pT in in centrality interval 0-10% for both leading jets with pT >120GeV and subleading jets with pT > 50 GeV.

Self-normalized zg distributions in PbPb and smeared pp collisions in the 30-50 centrality event class for jets with PTJET 160-180 GeV.

Self-normalized zg distributions in PbPb and smeared pp collisions in the 10-30 centrality event class for jets with PTJET 160-180 GeV.

Self-normalized zg distributions in PbPb and smeared pp collisions in the 0-10 centrality event class for jets with PTJET 160-180 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 50-80 centrality event class for jets with PTJET 160-180 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 30-50 centrality event class for jets with PTJET 160-180 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 10-30 centrality event class for jets with PTJET 160-180 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 0-10 centrality event class for jets with PTJET 160-180 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 10% most central events for jets with PTJET 140-160 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 10% most central events for jets with PTJET 160-180 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 10% most central events for jets with PTJET 180-200 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 10% most central events for jets with PTJET 200-250 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 10% most central events for jets with PTJET 250-300 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 10% most central events for jets with PTJET 300-500 GeV.

p-T differential invariant yield of Lambda + anti-Lambda in pp collisions with center-of-mass energy/nucleon = 7 TeV.

p-T differential invariant yield of Xi+ + Xi- in pp collisions with center-of-mass energy/nucleon = 7 TeV.

p-T differential invariant yield of Xi+ + Xi- in pp collisions with center-of-mass energy/nucleon = 7 TeV.

p-T differential invariant yield of K0s in pPb collisions with center-of-mass energy/nucleon = 5.02 TeV.

p-T differential invariant yield of K0s in pPb collisions with center-of-mass energy/nucleon = 5.02 TeV.

p-T differential invariant yield of Lambda + anti-Lambda in pPb collisions with center-of-mass energy/nucleon = 5.02 TeV.

p-T differential invariant yield of Lambda + anti-Lambda in pPb collisions with center-of-mass energy/nucleon = 5.02 TeV.

p-T differential invariant yield of Xi+ + Xi- in pPb collisions with center-of-mass energy/nucleon = 5.02 TeV.

p-T differential invariant yield of Xi+ + Xi- in pPb collisions with center-of-mass energy/nucleon = 5.02 TeV.

p-T differential invariant yield of K0s in PbPb collisions with center-of-mass energy/nucleon = 2.76 TeV.

p-T differential invariant yield of K0s in PbPb collisions with center-of-mass energy/nucleon = 2.76 TeV.

p-T differential invariant yield of Lambda + anti-Lambda in PbPb collisions with center-of-mass energy/nucleon = 2.76 TeV.

p-T differential invariant yield of Lambda + anti-Lambda in PbPb collisions with center-of-mass energy/nucleon = 2.76 TeV.

p-T differential invariant yield of Xi+ + Xi- in PbPb collisions with center-of-mass energy/nucleon = 2.76 TeV.

p-T differential invariant yield of Xi+ + Xi- in PbPb collisions with center-of-mass energy/nucleon = 2.76 TeV.

pT-dependent Lambda+anti-Lambda to K0S production ratios in pp collisions with centre-of-mass energy = 7 TeV.

pT-dependent Lambda+anti-Lambda to K0S production ratios in pPb collisions with centre-of-mass energy/nucleon = 5.02 TeV.

pT-dependent Lambda+anti-Lambda to K0S production ratios in PbPb collisions with centre-of-mass energy/nucleon = 2.76 TeV.

pT-dependent Xi+ + Xi- to Lambda+anti-Lambda production ratios in pp collisions with centre-of-mass energy = 7 TeV.

pT-dependent Xi+ + Xi- to Lambda+anti-Lambda production ratios in pPb collisions with centre-of-mass energy/nucleon = 5.02 TeV.

pT-dependent Xi+ + Xi- to Lambda+anti-Lambda production ratios in PbPb collisions with centre-of-mass energy/nucleon = 2.76 TeV.

Simultenous blast-wave fits parameters in pp collisions with center-of-mass energy/nucleon = 7 TeV.

Simultenous blast-wave fits parameters in pPb collisions with center-of-mass energy/nucleon = 5.02 TeV.

Simultenous blast-wave fits parameters in PbPb collisions with center-of-mass energy/nucleon = 2.76 TeV.

Mean transverse kinetic energy MEAN(KE_T) of K0s in pp collisions with centre-of-mass energy/nucleon=7 TeV.

Mean transverse kinetic energy MEAN(KE_T) of Lambda+anti-Lambda in pp collisions with centre-of-mass energy/nucleon=7 TeV.

Mean transverse kinetic energy MEAN(KE_T) of K0s in pPb collisions with centre-of-mass energy/nucleon=5.02 TeV.

Mean transverse kinetic energy MEAN(KE_T) of Lambda+anti-Lambda in pPb collisions with centre-of-mass energy/nucleon=5.02 TeV.

Mean transverse kinetic energy MEAN(KE_T) of Xi+ + Xi- in pPb collisions with centre-of-mass energy/nucleon=5.02 TeV.

Mean transverse kinetic energy MEAN(KE_T) of K0s in PbPb collisions with centre-of-mass energy/nucleon=2.76 TeV.

Mean transverse kinetic energy MEAN(KE_T) of Lambda + anti-Lambda in PbPb collisions with centre-of-mass energy/nucleon=2.76 TeV.

p-T differential invariant yield of K0s in pPb collisions with center-of-mass energy/nucleon = 5.02 TeV.

p-T differential invariant yield of Lambda + anti-Lambda in pPb collisions with center-of-mass energy/nucleon = 5.02 TeV.

p-T differential invariant yield of K0s in pPb collisions with center-of-mass energy/nucleon = 5.02 TeV.

p-T differential invariant yield of Lambda + anti-Lambda in pPb collisions with center-of-mass energy/nucleon = 5.02 TeV.

p-T differential invariant yield of K0s in pPb collisions with center-of-mass energy/nucleon = 5.02 TeV.

p-T differential invariant yield of Lambda + anti-Lambda in pPb collisions with center-of-mass energy/nucleon = 5.02 TeV.

pT-dependent Lambda+Anti-Lambda to K0S production ratios in pPb collisions with centre-of-mass energy/nucleon = 5.02 TeV.

Mean transverse kinetic energy MEAN(KE_T) of K0s in pPb collisions with centre-of-mass energy/nucleon=5.02 TeV.

Mean transverse kinetic energy MEAN(KE_T) of Lambda + anti-Lambda in pPb collisions with centre-of-mass energy/nucleon=5.02 TeV.

The Symmetric cumulant $SC(n,m)$ results as a function of multiplicity ($N_{trk}^{offline}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

The Symmetric cumulant $SC(n,m)$ results as a function of multiplicity ($N_{trk}^{offline}$) in pPb collisions at $\sqrt{s_{NN}}$ = 8.16 TeV.

The Symmetric cumulant $SC(n,m)$ results as a function of multiplicity ($N_{trk}^{offline}$) in pp collisions at $\sqrt{s}$ = 13.0 TeV.

The Normalized Symmetric cumulant $NSC(n,m)$ results as a function of multiplicity ($N_{trk}^{offline}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. The results are normalized by $\left<(v_2^{\text{sub}})^{2}\right> \, \left<(v_3^{\text{sub}})^{2}\right>$ and $\left<(v_2^{\text{sub}})^{2}\right> \, \left<(v_4^{\text{sub}})^{2}\right>$ from two-particle correlations.

The Normalized Symmetric cumulant $NSC(n,m)$ results as a function of multiplicity ($N_{trk}^{offline}$) in pPb collisions at $\sqrt{s_{NN}}$ = 8.16 TeV. The results are normalized by $\left<(v_2^{\text{sub}})^{2}\right> \, \left<(v_3^{\text{sub}})^{2}\right>$ and $\left<(v_2^{\text{sub}})^{2}\right> \, \left<(v_4^{\text{sub}})^{2}\right>$ from two-particle correlations.

The Normalized Symmetric cumulant $NSC(n,m)$ results as a function of multiplicity ($N_{trk}^{offline}$) in pp collisions at $\sqrt{s}$ = 13.0 TeV. The results are normalized by $\left<(v_2^{\text{sub}})^{2}\right> \, \left<(v_3^{\text{sub}})^{2}\right>$ and $\left<(v_2^{\text{sub}})^{2}\right> \, \left<(v_4^{\text{sub}})^{2}\right>$ from two-particle correlations.

95% CL limits on the Lagrangian approach (a.k.a. LEP parametrization) parameters in Wgamma events.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) elliptic flow, $v^{(\alpha)}_2$, as a function of $p_T$ in 10-20% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) elliptic flow, $v^{(\alpha)}_2$, as a function of $p_T$ in 20-30% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) elliptic flow, $v^{(\alpha)}_2$, as a function of $p_T$ in 30-40% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) elliptic flow, $v^{(\alpha)}_2$, as a function of $p_T$ in 40-50% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) elliptic flow, $v^{(\alpha)}_2$, as a function of $p_T$ in 50-60% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ in 0-0.2% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ in 0-5% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ in 0-10% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ in 10-20% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ in 20-30% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ in 30-40% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ in 40-50% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ in 50-60% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) elliptic flow, $v^{(\alpha)}_2$, as a function of $p_T$ for multiplicity class $220 \leq N^{offline}_{trk}<260$ in pPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) elliptic flow, $v^{(\alpha)}_2$, as a function of $p_T$ for multiplicity class $185 \leq N^{offline}_{trk}<220$ in pPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) elliptic flow, $v^{(\alpha)}_2$, as a function of $p_T$ for multiplicity class $150 \leq N^{offline}_{trk}<185$ in pPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) elliptic flow, $v^{(\alpha)}_2$, as a function of $p_T$ $120 \leq N^{offline}_{trk}<150$ in pPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ for multiplicity class $220 \leq N^{offline}_{trk}<260$ in pPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ for multiplicity class $185 \leq N^{offline}_{trk}<220$ in pPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ for multiplicity class $150 \leq N^{offline}_{trk}<185$ in pPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ for multiplicity class $120 \leq N^{offline}_{trk}<150$ in pPb collisions.

Pearson correlation coefficients $r_2$ and $r_3$ reconstructed using the leading and subleading flows as a function of $p^{a}_{T} - p^{b}_{T}$ in bin $p^{a}_{T}$ for centrality class 0-0.2% in PbPb collisions.

Pearson correlation coefficients $r_2$ and $r_3$ reconstructed using the leading and subleading flows as a function of $p^{a}_{T} - p^{b}_{T}$ in bin $p^{a}_{T}$ for centrality class 0-5% in PbPb collisions.

Pearson correlation coefficients $r_2$ and $r_3$ reconstructed using the leading and subleading flows as a function of $p^{a}_{T} - p^{b}_{T}$ in bin $p^{a}_{T}$ for centrality class 10-20% in PbPb collisions.

Pearson correlation coefficients $r_2$ and $r_3$ reconstructed using the leading and subleading flows as a function of $p^{a}_{T} - p^{b}_{T}$ in bin $p^{a}_{T}$ for centrality class 20-30% in PbPb collisions.

Pearson correlation coefficients $r_2$ and $r_3$ reconstructed using the leading and subleading flows as a function of $p^{a}_{T} - p^{b}_{T}$ in bin $p^{a}_{T}$ for centrality class 30-40% in PbPb collisions.

Pearson correlation coefficients $r_2$ and $r_3$ reconstructed using the leading and subleading flows as a function of $p^{a}_{T} - p^{b}_{T}$ in bin $p^{a}_{T}$ for centrality class 40-50% in PbPb collisions.

Elliptic and triangular ratios of subleading and leading flow as a function of multiplicity for the highest bin 2.5 < $p_T$ < 3.0 GeV in PbPb collisions.

Elliptic and triangular ratios of subleading and leading flow as a function of multiplicity for the highest bin 2.5 < $p_T$ < 3.0 GeV in pPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) multiplicity modes, $v^{(\alpha)}_0$, as a function of $p_T$ in 0-0.2% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) multiplicity modes, $v^{(\alpha)}_0$, as a function of $p_T$ in 0-5% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) multiplicity modes, $v^{(\alpha)}_0$, as a function of $p_T$ in 0-10% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) multiplicity modes, $v^{(\alpha)}_0$, as a function of $p_T$ in 10-20% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) multiplicity modes, $v^{(\alpha)}_0$, as a function of $p_T$ in 20-30% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) multiplicity modes, $v^{(\alpha)}_0$, as a function of $p_T$ in 30-40% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) multiplicity modes, $v^{(\alpha)}_0$, as a function of $p_T$ in 40-50% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) multiplicity modes, $v^{(\alpha)}_0$, as a function of $p_T$ in 50-60% centrality PbPb collisions.

The ratio of the production cross sections of inclusive $\Lambda_{c}^{+}$ to prompt $D^{0}$ versus $p_{T}$ from $0-100\%$ centrality PbPb collisions. The $31\%$ normalization uncertainty in PbPb collisions, is not included in the statistical and systematic uncertainties for each data point.

Differential cross section for Y(2S) states as a function of their transverse momentum and per unit of rapidity in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(3S) states as a function of their transverse momentum and per unit of rapidity in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(3S) states as a function of their transverse momentum and per unit of rapidity in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(1S) states as a function of their transverse momentum and per unit of rapidity in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(1S) states as a function of their transverse momentum and per unit of rapidity in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(2S) states as a function of their transverse momentum and per unit of rapidity in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(2S) states as a function of their transverse momentum and per unit of rapidity in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(1S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(1S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(2S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(2S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(3S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(3S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(1S) states as a function of their rapidity and integrated over transverse momentum in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(1S) states as a function of their rapidity and integrated over transverse momentum in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(2S) states as a function of their rapidity and integrated over transverse momentum in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(2S) states as a function of their rapidity and integrated over transverse momentum in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Nuclear modification factor for Υ(1S) states in PbPb collisions as a function of pT. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(1S) states in PbPb collisions as a function of pT. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(2S) states in PbPb collisions as a function of pT. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(2S) states in PbPb collisions as a function of pT. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(1S) states in PbPb collisions as a function of |y|. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(1S) states in PbPb collisions as a function of |y|. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(2S) states in PbPb collisions as a function of |y|. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(2S) states in PbPb collisions as a function of |y|. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factors for Υ(1S) meson production in PbPb collisions, as a function of centrality, displayed as the average number of participating nucleons. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainties from the PbPb data (3.2%) or from the pp reference (6.3%) are displayed at unity as empty and filled red black boxes, respectively.

Nuclear modification factors for Υ(1S) meson production in PbPb collisions, as a function of centrality, displayed as the average number of participating nucleons. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainties from the PbPb data (3.2%) or from the pp reference (6.3%) are displayed at unity as empty and filled red black boxes, respectively.

Nuclear modification factors for Υ(2S) meson production in PbPb collisions, as a function of centrality, displayed as the average number of participating nucleons. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainties from the PbPb data (3.2%) or from the pp reference (6.9%) are displayed at unity as empty and filled black black boxes, respectively.

Nuclear modification factors for Υ(2S) meson production in PbPb collisions, as a function of centrality, displayed as the average number of participating nucleons. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainties from the PbPb data (3.2%) or from the pp reference (6.9%) are displayed at unity as empty and filled black black boxes, respectively.

Nuclear modification factors for Υ(1S), Υ(2S) and Υ(3S) meson production in PbPb collisions, integrated over centrality. For the Υ(3S), the upper limit derived on the nuclear modification factor is given.

Nuclear modification factors for Υ(1S), Υ(2S) and Υ(3S) meson production in PbPb collisions, integrated over centrality. For the Υ(3S), the upper limit derived on the nuclear modification factor is given.

Three-particle correlation with respect to the 2nd order event plane from Pb-going side in pPb collisions as a function of multiplicity with individual track pT between 0.3 to 3.0 GeV/c.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions as a function of multiplicity with individual track pT between 0.3 to 3.0 GeV/c.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions as a function of centrality with individual track pT between 0.3 to 3.0 GeV/c, where the average multiplicity value is plotted on figure in order to compare with pPb collisions.

Three-particle correlation with respect to the 2nd order event plane from Pb-going side, p-going side in pPb collisions as a function of multiplicity with individual track pT between 0.3 to 3.0 GeV/c.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions as a function of multiplicity with individual track pT between 0.3 to 3.0 GeV/c.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions as a function of centrality with individual track pT between 0.3 to 3.0 GeV/c, where the average multiplicity value is plotted on figure in order to compare with pPb collisions.

Three-particle correlation with respect to the 2nd order event plane from Pb-going and p-going direction in pPb collisions and in PbPb collisions as a function DeltaEta for multiplicity [185,220) with individual track pT between 0.3 to 3.0 GeV/c. Data points are plotted at the bin center.

Three-particle correlation with respect to the 3rd order event plane from Pb-going side in pPb collisions.

Three-particle correlation with respect to the 3rd order event plane from p-going side in pPb collisions.

Three-particle correlation with respect to the 3rd order event plane in PbPb collisions.

Two-particle correlation in pPb collisions.

Two-particle correlation in PbPb collisions.

Three-particle correlation with respect to the 2nd order event plane from Pb-going side in pPb collisions.

Three-particle correlation with respect to the 2nd order event plane from p-going side in pPb collisions.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions.

Three-particle correlation with respect to the 3rd order event plane from Pb-going side in pPb collisions.

Three-particle correlation with respect to the 3rd order event plane from p-going side in pPb collisions.

Three-particle correlation with respect to the 3rd order event plane in PbPb collisions.

Two-particle correlation in pPb collisions.

Two-particle correlation in PbPb collisions.

Three-particle correlation with respect to the 2nd order event plane from Pb-going side in pPb collisions.

Three-particle correlation with respect to the 2nd order event plane from p-going side in pPb collisions.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions.

Three-particle correlation with respect to the 3rd order event plane from Pb-going side in pPb collisions.

Three-particle correlation with respect to the 3rd order event plane from p-going side in pPb collisions.

Three-particle correlation with respect to the 3rd order event plane in PbPb collisions.

Two-particle correlation in pPb collisions.

Two-particle correlation in PbPb collisions.

Three-particle correlation with respect to the 2nd order event plane from Pb-going side in pPb collisions.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions.

Three-particle correlation with respect to the 2nd order event plane from Pb-going side in pPb collisions.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions.

Three-particle correlation with respect to the 3rd order event plane from Pb-going side in pPb collisions.

Three-particle correlation with respect to the 3rd order event plane in PbPb collisions.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions.

Two-particle correlation in pPb collisions.

Two-particle correlation in PbPb collisions.

Two-particle correlation in PbPb collisions.

Three-particle correlation with respect to the 2nd order event plane from Pb-going side in pPb collisions.

Three-particle correlation with respect to the 2nd order event plane from p-going side in pPb collisions.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions.

Three-particle correlation with respect to the 2nd order event plane from Pb-going side in pPb collisions.

Three-particle correlation with respect to the 2nd order event plane from p-going side in pPb collisions.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions.

Three-particle correlation with respect to the 2nd order event plane from Pb-going side in pPb collisions.

Three-particle correlation with respect to the 2nd order event plane from p-going side in pPb collisions.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions.

Three-particle correlation with respect to the 3rd order event plane from Pb-going side in pPb collisions.

Three-particle correlation with respect to the 3rd order event plane from p-going side in pPb collisions.

Three-particle correlation with respect to the 3rd order event plane in PbPb collisions.

Three-particle correlation with respect to the 3rd order event plane from Pb-going side in pPb collisions.

Three-particle correlation with respect to the 3rd order event plane from p-going side in pPb collisions.

Three-particle correlation with respect to the 3rd order event plane in PbPb collisions.

Three-particle correlation with respect to the 3rd order event plane from Pb-going side in pPb collisions.

Three-particle correlation with respect to the 3rd order event plane from p-going side in pPb collisions.

Three-particle correlation with respect to the 3rd order event plane in PbPb collisions.

Two-particle correlation in pPb collisions.

Two-particle correlation in PbPb collisions.

Two-particle correlation in pPb collisions.

Two-particle correlation in PbPb collisions.

Two-particle correlation in pPb collisions.

Two-particle correlation in PbPb collisions.

Three-particle correlation with respect to the 2nd order event plane in Pb-going direction in pPb collisions.

Three-particle correlation with respect to the 2nd order event plane in p-going direction in pPb collisions.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions.

Three-particle correlation with respect to the 3rd order event plane in Pb-going direction in pPb collisions.

Three-particle correlation with respect to the 3rd order event plane in PbPb collisions.

Three-particle correlation with respect to the 3rd order event plane in PbPb collisions.

Two-particle correlation in pPb collisions.

Two-particle correlation in PbPb collisions.

Two-particle correlation in PbPb collisions.

Three-particle to two-particle correlation ratio from Pb-going side in pPb collisions.

Three-particle to two-particle correlation ratio in PbPb collisions.

Three-particle to two-particle correlation ratio in PbPb collisions.

Three-particle to two-particle correlation ratio in pPb (Pb- and p-going direction) collisions.

Three-particle to two-particle correlation ratio in pPb (Pb- and p-going direction) collisions.

Three-particle to two-particle correlation ratio in pPb (Pb- and p-going direction) collisions.

Three-particle to two-particle correlation ratio in PbPb collisions.

Three-particle to two-particle correlation ratio in PbPb collisions.

Three-particle to two-particle correlation ratio in PbPb collisions.

Three-particle correlation with respect to the 2nd order event plane from Pb-going side as a function of v2 in pPb collisions.

Three-particle correlation with respect to the 2nd order event plane from p-going side as a function of v2 in pPb collisions.

Three-particle correlation with respect to the 2nd order event plane as a function of v2 in PbPb collisions.

Three-particle correlation (OS-SS) with respect to the 2nd order event plane as a function of v2 in pPb collisions.

Three-particle correlation (OS-SS) with respect to the 2nd order event plane as a function of v2 in PbPb collisions.

Three-particle correlation (OS-SS) with respect to the 2nd order event plane as a function of v2 in PbPb collisions.

Two-particle correlation (OS-SS) as a function of v2 in pPb collisions.

Two-particle correlation (OS-SS) as a function of v2 in PbPb collisions.

Ratio of three-particle and two-particle correlation (OS-SS) as a function of v2 in pPb collisions.

Ratio of three-particle and two-particle correlation (OS-SS) as a function of v2 in PbPb collisions.

Extracted intercept parameter b_norm in PbPb collisions.

Extracted intercept parameter b_norm in pPb collisions.

Upper limit on v2-independent gamma_112 correlator component in PbPb collisions.

Upper limit on v2-independent gamma_112 correlator component in pPb collisions.

The $\langle x_{\mathrm{j}\gamma} \rangle$ and $R_{\mathrm{j}\gamma}$, the number of associated jets per photon, in 0-30% and 30-100% centrality PbPb collisions. The smeared pp data are added for comparison. A value of 999 in the bin range indicates that there is no upper limit.

The $I_{\mathrm{AA}}^{\mathrm{jet}}$ vs. $p_{\mathrm{T}}^{\mathrm{jet}}$ for 0-30% and 30-100% centrality PbPb collisions. Bins for which there are no data points are denoted as 'empty'.

Detector-level photon $\phi$ distribution

Detector-level dipoton invariant mass distribution

Detector-level diphoton rapidity distribution

Detector-level dipoton p$_{T}$ distribution

Observed exclusion limits at 95% CL in the ALP-photon coupling g$_{a \gamma}$ versus ALP mass $m_{a}$ plane.

Differential cross section of J/psi mesons from b hadrons as a function of dimuon pT in pp and PbPb collisions. The PbPb cross sections are normalised by TAA for direct comparison. Global uncertainties arise from the integrated luminosity uncertainty in pp collisions, and the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross section of prompt J/psi mesons as a function of dimuon rapidty in pp and PbPb collisions. The PbPb cross sections are normalised by TAA for direct comparison. Global uncertainties arise from the integrated luminosity uncertainty in pp collisions, and the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross section of J/psi mesons from b hadrons as a function of dimuon rapidty in pp and PbPb collisions. The PbPb cross sections are normalised by TAA for direct comparison. Global uncertainties arise from the integrated luminosity uncertainty in pp collisions, and the number of minimum bias events and TAA uncertainties for PbPb collisions.

Nuclear modification factor of prompt J/psi mesons as a function of dimuon pT in the forward rapidity range. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.