Breakdown of relative systematic uncertainties in the measured ${t\bar{t}}$ cross-section. The quoted uncertainties are obtained by repeating the fit with a group of nuisance parameters fixed to their fitted values and subtracting in quadrature the resulting total uncertainty from the uncertainty of the complete fit. However, the total uncertainty is not the quadratic sum of the grouped impacts, as this approach neglects the correlation among the different groups.

The figure shows the measurement of the top quark pair production cross-section, $\sigma_{t\bar{t}}$, in lead-lead (Pb+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV, based on 1.9 nb$^{-1}$ of data.

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the resolved qqbb signal region.

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the resolved lvbb signal region.

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the merged lvbb signal region

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the merged qqbb signal region

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 5jex3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 5jex4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 6jin3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 6jin4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 5jex3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 5jex4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 6jin3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 6jin4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 5jex3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 5jex4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 6jin3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 6jin4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

The mWh distributions in the Low-NN-score SRs of the merged qqbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Low-NN-score SRs of the merged qqbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged qqbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged qqbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Low-NN-score SRs of the merged lvbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Low-NN-score SRs of the merged lvbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Medium-NN-score SRs of the merged lvbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Medium-NN-score SRs of the merged lvbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged lvbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged lvbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

Cutflow table for a few selected signal samples after applying the selection requirements of the resolved analysis. At the initial cut stage, all events from the tbH ± → W ± h(→ b b) production process (including those from the fully hadronic decay channel) are considered.

Cutflow table for a few selected signal samples after applying the selection requirements of the merged analysis. At the initial cut stage, all events from the tbH ± → W ± h(→ bb) production process (including those from the fully hadronic decay channel) are considered.

Distributions in $m_{Z'}^{rec}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the muon channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Distributions in $m_{Z'}^{T}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the invisible channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Distributions in $m_{Z'}^{T}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the invisible channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Distributions in $m_{Z'}^{rec}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the electron channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Distributions in $m_{Z'}^{rec}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the electron channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Expected upper limits at 95% CL on the product of the production cross section and the branching fraction as a function of the Z' boson mass.Each line corresponds to the three different final states and their combined result.

Expected upper limits at 95% CL on the product of the production cross section and the branching fraction as a function of the Z' boson mass.Each line corresponds to the three different final states and their combined result.

Expected and observed upper limits at 95% CL on the product of the production cross section and the branching fraction as functions of the Z' boson mass from the combination of all final states. The limits are compared with the predictions of the HVT model and the expected limits (shown by red curves) from a previous analysis.

Expected and observed upper limits at 95% CL on the product of the production cross section and the branching fraction as functions of the Z' boson mass from the combination of all final states. The limits are compared with the predictions of the HVT model and the expected limits (shown by red curves) from a previous analysis.

Prediction from the HVT model.

Prediction from the HVT model.

Observed upper limits at 95% CL on gF for different Z' boson masses as functions of the product of $g_{H}$ with the sign of $g_{F}$. The two benchmark scenarios of the HVT model are shown by the black markers.

Observed upper limits at 95% CL on gF for different Z' boson masses as functions of the product of $g_{H}$ with the sign of $g_{F}$. The two benchmark scenarios of the HVT model are shown by the black markers.

Distributions in $S_T$ in the SR for the electron channel, after a background-only fit to the data. The signal distributions are scaled to the cross section predicted by the theory. The hatched bands show the post-fit uncertainty band, combining all sources of uncertainty. The ratio of data to the background predictions is shown in the panels below the distributions.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM1.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM1.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM1.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM1.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM1.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM1.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM2.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM2.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM2.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM2.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM2.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM2.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM3.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM3.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM3.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM3.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM3.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM3.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp1.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp1.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp1.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp1.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp1.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp1.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp2.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp2.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp2.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp2.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp2.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp2.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp3.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp3.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp3.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp3.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp3.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp3.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp-1c.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp-1c.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to top or charm quarks and 1 GeV neutralinos for different branching fraction assumptions.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to top or charm quarks and 1 GeV neutralinos for different branching fraction assumptions.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to top or charm quarks and 1 GeV neutralinos for different branching fraction assumptions.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to top or charm quarks and 1 GeV neutralinos for different branching fraction assumptions.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to top or charm quarks and 1 GeV neutralinos for different branching fraction assumptions.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to top or charm quarks and 1 GeV neutralinos for different branching fraction assumptions.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to top or charm quarks and 200 GeV neutralinos for different branching fraction assumptions.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to top or charm quarks and 200 GeV neutralinos for different branching fraction assumptions.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to top or charm quarks and 200 GeV neutralinos for different branching fraction assumptions.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to top or charm quarks and 200 GeV neutralinos for different branching fraction assumptions.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to top or charm quarks and 200 GeV neutralinos for different branching fraction assumptions.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to top or charm quarks and 200 GeV neutralinos for different branching fraction assumptions.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of up-type scalar LQs coupled to first generation leptons and second generation quarks.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of up-type scalar LQs coupled to first generation leptons and second generation quarks.

Expected exclusion limit at 95% CL for pair production of up-type scalar LQs coupled to first generation leptons and second generation quarks.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of up-type scalar LQs coupled to first generation leptons and second generation quarks.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of up-type scalar LQs coupled to first generation leptons and second generation quarks.

Observed exclusion limit at 95% CL for pair production of up-type scalar LQs coupled to first generation leptons and second generation quarks.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of up-type scalar LQs coupled to second generation leptons and quarks.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of up-type scalar LQs coupled to second generation leptons and quarks.

Expected exclusion limit at 95% CL for pair production of up-type scalar LQs coupled to second generation leptons and quarks.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of up-type scalar LQs coupled to second generation leptons and quarks.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of up-type scalar LQs coupled to second generation leptons and quarks.

Observed exclusion limit at 95% CL for pair production of up-type scalar LQs coupled to second generation leptons and quarks.

Observed upper limit on the cross-section at 95% CL for pair production of top squarks decaying into charm quarks and neutralinos for the fit with the best expected CL$_\mathrm{S}$.

Observed upper limit on the cross-section for pair production of up-type scalar LQs coupled to first generation leptons and second generation quarks.

Observed upper limit on the cross-section for pair production of up-type scalar LQs coupled to second generation leptons and quarks.

Expected upper limit on the cross-section $(-1\sigma)$ for $U(1)$ vector LQ models in the minimal coupling scenario with $\beta_{23}=1$, corresponding to a branching fraction $\mathcal{B}(\mathrm{vLQ}^\mathrm{u,min}_{23}\rightarrow c\nu) = 0.5$.

Expected upper limit on the cross-section $(+1\sigma)$ for $U(1)$ vector LQ models in the minimal coupling scenario with $\beta_{23}=1$, corresponding to a branching fraction $\mathcal{B}(\mathrm{vLQ}^\mathrm{u,min}_{23}\rightarrow c\nu) = 0.5$.

Expected upper limit on the cross-section for $U(1)$ vector LQ models in the minimal coupling scenario with $\beta_{23}=1$, corresponding to a branching fraction $\mathcal{B}(\mathrm{vLQ}^\mathrm{u,min}_{23}\rightarrow c\nu) = 0.5$.

Observed upper limit on the cross-section for $U(1)$ vector LQ models in the minimal coupling scenario with $\beta_{23}=1$, corresponding to a branching fraction $\mathcal{B}(\mathrm{vLQ}^\mathrm{u,min}_{23}\rightarrow c\nu) = 0.5$.

Expected upper limit on the cross-section $(-1\sigma)$ for $U(1)$ vector LQ models in the Yang-Mills coupling scenario with $\beta_{23}=1$, corresponding to a branching fraction $\mathcal{B}(\mathrm{vLQ}^\mathrm{u,YM}_{23}\rightarrow c\nu) = 0.5$.

Expected upper limit on the cross-section $(+1\sigma)$ for $U(1)$ vector LQ models in the Yang-Mills coupling scenario with $\beta_{23}=1$, corresponding to a branching fraction $\mathcal{B}(\mathrm{vLQ}^\mathrm{u,YM}_{23}\rightarrow c\nu) = 0.5$.

Expected upper limit on the cross-section for $U(1)$ vector LQ models in the Yang-Mills coupling scenario with $\beta_{23}=1$, corresponding to a branching fraction $\mathcal{B}(\mathrm{vLQ}^\mathrm{u,YM}_{23}\rightarrow c\nu) = 0.5$.

Observed upper limit on the cross-section for $U(1)$ vector LQ models in the Yang-Mills coupling scenario with $\beta_{23}=1$, corresponding to a branching fraction $\mathcal{B}(\mathrm{vLQ}^\mathrm{u,YM}_{23}\rightarrow c\nu) = 0.5$.

Distribution of $R_{\mathrm{ISR}}$ in SR-Comp-1c showing the data and the expected backgrounds, after simultaneously fitting all the control regions. The pre-fit yields for a representative $m(\tilde{t}_1,\tilde{\chi}_1^0)=(450,430)$ GeV signal is included. Processes with top quarks and multiple vector bosons are included in "Other". The right-most bin does not include the overflow entry but the $x$-axis range is chosen to include all observed data.

Distribution of $E_\mathrm{T}^\mathrm{miss}\mathrm{ Sig.}$ in SR-HM1 showing the data and the expected backgrounds, after simultaneously fitting all the control regions. The pre-fit yields for a representative $m(\tilde{t}_1,\tilde{\chi}_1^0)=(1000,1)$ GeV signal is included. Processes with top quarks and multiple vector bosons are included in "Other". The right-most bin does not include the overflow entry but the $x$-axis range is chosen to include all observed data.

Distribution of $m_{\mathrm{T}}^{\mathrm{min}}(c)$ in SR-HM2 showing the data and the expected backgrounds, after simultaneously fitting all the control regions. The pre-fit yields for a representative $m(\tilde{t}_1,\tilde{\chi}_1^0)=(750,450)$ GeV signal is included. Processes with top quarks and multiple vector bosons are included in "Other". The right-most bin does not include the overflow entry but the $x$-axis range is chosen to include all observed data.

Distribution of $R_{\mathrm{ISR}}$ in SR-Comp1 showing the data and the expected backgrounds, after simultaneously fitting all the control regions. The pre-fit yields for a representative $m(\tilde{t}_1,\tilde{\chi}_1^0)=(600,550)$ GeV signal is included. Processes with top quarks and multiple vector bosons are included in "Other". The right-most bin does not include the overflow entry but the $x$-axis range is chosen to include all observed data.

Distribution of $R_{\mathrm{ISR}}$ in SR-Comp2 showing the data and the expected backgrounds, after simultaneously fitting all the control regions. The pre-fit yields for a representative $m(\tilde{t}_1,\tilde{\chi}_1^0)=(550,470)$ GeV signal is included. Processes with top quarks and multiple vector bosons are included in "Other". The right-most bin does not include the overflow entry but the $x$-axis range is chosen to include all observed data.

Distribution of $R_{\mathrm{ISR}}$ in SR-Comp3 showing the data and the expected backgrounds, after simultaneously fitting all the control regions. The pre-fit yields for a representative $m(\tilde{t}_1,\tilde{\chi}_1^0)=(550,375)$ GeV signal is included. Processes with top quarks and multiple vector bosons are included in "Other". The right-most bin does not include the overflow entry but the $x$-axis range is chosen to include all observed data.

Acceptance across the SUSY signal grid for SR-Comp-1c.

Acceptance across the SUSY signal grid for SR-HM1.

Acceptance across the SUSY signal grid for SR-HM2.

Acceptance across the SUSY signal grid for SR-HM3.

Acceptance across the SUSY signal grid for SR-HM-Disc.

Acceptance across the SUSY signal grid for SR-Comp1.

Acceptance across the SUSY signal grid for SR-Comp2.

Acceptance across the SUSY signal grid for SR-Comp3.

Acceptance across the $\mathrm{LQ}^\mathrm{u}_{21}$ signal grid for SR-HM1.

Acceptance across the $\mathrm{LQ}^\mathrm{u}_{21}$ signal grid for SR-HM2.

Acceptance across the $\mathrm{LQ}^\mathrm{u}_{21}$ signal grid for SR-HM3.

Acceptance across the $\mathrm{LQ}^\mathrm{u}_{21}$ signal grid for SR-HM-Disc.

Acceptance across the $\mathrm{LQ}^\mathrm{u}_{22}$ signal grid for SR-HM1.

Acceptance across the $\mathrm{LQ}^\mathrm{u}_{22}$ signal grid for SR-HM2.

Acceptance across the $\mathrm{LQ}^\mathrm{u}_{22}$ signal grid for SR-HM3.

Acceptance across the $\mathrm{LQ}^\mathrm{u}_{22}$ signal grid for SR-HM-Disc.

Acceptance for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM1.

Acceptance for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM2.

Acceptance for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM3.

Acceptance for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM-Disc.

Acceptance for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM1.

Acceptance for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM2.

Acceptance for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM3.

Acceptance for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM-Disc.

Efficiency for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM1. Due to the limited statistical power of the signal samples, efficiencies can fluctuate and can drop to zero.

Efficiency for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM2. Due to the limited statistical power of the signal samples, efficiencies can fluctuate and can drop to zero.

Efficiency for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM3. Due to the limited statistical power of the signal samples, efficiencies can fluctuate and can drop to zero.

Efficiency for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM-Disc. Due to the limited statistical power of the signal samples, efficiencies can fluctuate and can drop to zero.

Efficiency for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM1. Due to the limited statistical power of the signal samples, efficiencies can fluctuate and can drop to zero.

Efficiency for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM2. Due to the limited statistical power of the signal samples, efficiencies can fluctuate and can drop to zero.

Efficiency for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM3. Due to the limited statistical power of the signal samples, efficiencies can fluctuate and can drop to zero.

Efficiency for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM-Disc. Due to the limited statistical power of the signal samples, efficiencies can fluctuate and can drop to zero.

Acceptance times efficiency across the SUSY signal grid for SR-Comp-1c.

Acceptance times efficiency across the SUSY signal grid for SR-HM1.

Acceptance times efficiency across the SUSY signal grid for SR-HM2.

Acceptance times efficiency across the SUSY signal grid for SR-HM3.

Acceptance times efficiency across the SUSY signal grid for SR-HM-Disc.

Acceptance times efficiency across the SUSY signal grid for SR-Comp1.

Acceptance times efficiency across the SUSY signal grid for SR-Comp2.

Acceptance times efficiency across the SUSY signal grid for SR-Comp3.

Acceptance times efficiency across the $\mathrm{LQ}^\mathrm{u}_{21}$ signal grid for SR-HM1.

Acceptance times efficiency across the $\mathrm{LQ}^\mathrm{u}_{21}$ signal grid for SR-HM2.

Acceptance times efficiency across the $\mathrm{LQ}^\mathrm{u}_{21}$ signal grid for SR-HM3.

Acceptance times efficiency across the $\mathrm{LQ}^\mathrm{u}_{21}$ signal grid for SR-HM-Disc.

Acceptance times efficiency across the $\mathrm{LQ}^\mathrm{u}_{22}$ signal grid for SR-HM1.

Acceptance times efficiency across the $\mathrm{LQ}^\mathrm{u}_{22}$ signal grid for SR-HM2.

Acceptance times efficiency across the $\mathrm{LQ}^\mathrm{u}_{22}$ signal grid for SR-HM3.

Acceptance times efficiency across the $\mathrm{LQ}^\mathrm{u}_{22}$ signal grid for SR-HM-Disc.

Acceptance times efficiency for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM1.

Acceptance times efficiency for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM2.

Acceptance times efficiency for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM3.

Acceptance times efficiency for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM-Disc.

Acceptance times efficiency for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM1.

Acceptance times efficiency for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM2.

Acceptance times efficiency for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM3.

Acceptance times efficiency for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM-Disc.

Cutflow table for SR-Comp-1c using the $m(\tilde{t}_{1},\tilde{\chi}_{1}^{0}) = (450,430)$ GeV signal point as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-Comp1 using the $m(\tilde{t}_{1},\tilde{\chi}_{1}^{0}) = (600,550)$ GeV signal point as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-Comp2 using the $m(\tilde{t}_{1},\tilde{\chi}_{1}^{0}) = (550,470)$ GeV signal point as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-Comp3 using the $m(\tilde{t}_{1},\tilde{\chi}_{1}^{0}) = (550,375)$ GeV signal point as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-HM1 using the $m(\tilde{t}_{1},\tilde{\chi}_{1}^{0}) = (1000,1)$ GeV signal point as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>40$ GeV, however, the $E_\mathrm{T}^\mathrm{miss}$ is calculated with jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-HM2 using the $m(\tilde{t}_{1},\tilde{\chi}_{1}^{0}) = (750,450)$ GeV signal point as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>40$ GeV, however, the $E_\mathrm{T}^\mathrm{miss}$ is calculated with jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-HM3 using the $m(\tilde{t}_{1},\tilde{\chi}_{1}^{0}) = (750,450)$ GeV signal point as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>40$ GeV, however, the $E_\mathrm{T}^\mathrm{miss}$ is calculated with jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-HM-Disc using the $m(\tilde{t}_{1},\tilde{\chi}_{1}^{0}) = (750,450)$ GeV signal point as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>40$ GeV, however, the $E_\mathrm{T}^\mathrm{miss}$ is calculated with jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-HM1 using the $m(\mathrm{LQ}_{21}^\mathrm{u}) = 600$ GeV signal point, with a 50% branching fraction to $c\nu$, as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>40$ GeV, however, the $E_\mathrm{T}^\mathrm{miss}$ is calculated with jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-HM2 using the $m(\mathrm{LQ}_{21}^\mathrm{u}) = 600$ GeV signal point, with a 50% branching fraction to $c\nu$, as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>40$ GeV, however, the $E_\mathrm{T}^\mathrm{miss}$ is calculated with jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-HM3 using the $m(\mathrm{LQ}_{21}^\mathrm{u}) = 600$ GeV signal point, with a 50% branching fraction to $c\nu$, as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>40$ GeV, however, the $E_\mathrm{T}^\mathrm{miss}$ is calculated with jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-HM-Disc using the $m(\mathrm{LQ}_{21}^\mathrm{u}) = 600$ GeV signal point, with a 50% branching fraction to $c\nu$, as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>40$ GeV, however, the $E_\mathrm{T}^\mathrm{miss}$ is calculated with jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the compressed SUSY simplified model with a bino-like LSP and wino-like NLSPs being considered. These are shown with $\pm1\sigma_\text{exp}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm1\sigma^{\text{SUSY}}_{\text{theory}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the ATLAS searches using the soft lepton signature is illustrated by the blue region while the limit imposed by the LEP experiments is shown in grey.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the compressed SUSY simplified model with a bino-like LSP and wino-like NLSPs being considered. These are shown with $\pm1\sigma_ ext{exp}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm1\sigma^{ ext{SUSY}}_{ ext{theory}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the ATLAS searches using the soft lepton signature is illustrated by the blue region while the limit imposed by the LEP experiments is shown in grey.

The expected upper limits on the cross-sections at 95% CL for the simplified SUSY model featuring the production of mass-degenerate pairs of wino-like $\widetilde{\chi}_{2}^{0}$ and $\widetilde{\chi}_{1}^{\pm}$ particles in association with at least two jets. The production modes considered are $\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{1}^{\pm}$, $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\mp}$, $\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{2}^{0}$, and $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\pm}$.

The observed upper limits on the cross-sections at 95% CL for the simplified SUSY model featuring the production of mass-degenerate pairs of wino-like $\widetilde{\chi}_{2}^{0}$ and $\widetilde{\chi}_{1}^{\pm}$ particles in association with at least two jets. The production modes considered are $\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{1}^{\pm}$, $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\mp}$, $\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{2}^{0}$, and $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\pm}$.

Truth-level signal acceptances in $\text{SR}_\text{2j}$ with BDT score $\in [0.88, 1.0]$. All of the considered production modes ($\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{1}^{\pm}$, $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\mp}$, $\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{2}^{0}$, and $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\pm}$) are included. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation. The generator-level selections are the same as the ones at reconstruction level.

Truth-level signal acceptances in $\text{SR}_{\geq\text{3j}}$ with BDT score $\in [0.88, 1.0]$. All of the considered production modes ($\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{1}^{\pm}$, $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\mp}$, $\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{2}^{0}$, and $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\pm}$) are included. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation. The generator-level selections are the same as the ones at reconstruction level.

Signal efficiencies in $\text{SR}_\text{2j}$ with BDT score $\in [0.88, 1.0]$. All of the considered production modes ($\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{1}^{\pm}$, $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\mp}$, $\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{2}^{0}$, and $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\pm}$) are included. The efficiency is defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level. The generator-level selections are the same as the ones at reconstruction level.

Signal efficiencies in $\text{SR}_{\geq\text{3j}}$ with BDT score $\in [0.88, 1.0]$. All of the considered production modes ($\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{1}^{\pm}$, $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\mp}$, $\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{2}^{0}$, and $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\pm}$) are included. The efficiency is defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level. The generator-level selections are the same as the ones at reconstruction level.

Event selection cutflows for signal samples with $m(\tilde\chi^0_2) =100~\text{GeV}$ and $\Delta m(\tilde\chi^0_2 / \tilde\chi^\pm_1,\tilde\chi^0_1) = 0.2~\text{GeV}, 1.0~\text{GeV},$ and $5.0~\text{GeV}$. The first number in the column for each mass point corresponds to the signal event yield after applying the associated cut while the second number in parentheses represents the efficiency with respect to the previous cut. All of the production modes are included. The total cross-section used to obtain the initial number of events ($\sigma$) includes the cuts applied at parton level as described in the main text.

Self-normalised ZN energy as a function of the self-normalised charged-particle-multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the high-ZN-energy classification (V0M+SPDClusters event classes).

Self-normalised ZN energy as a function of the self-normalised charged-particle-multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the low-ZN-energy classification (V0M+SPDClusters event classes).

$K^{0}_{S}$ transverse momentum spectrum in pp collisions at $\sqrt{s}$ = 13 TeV in the high-ZN-energy classification (V0M+SPDClusters event classes).

$\Lambda+\overline{\Lambda}$ transverse momentum spectrum in pp collisions at $\sqrt{s}$ = 13 TeV in the high-ZN-energy classification (V0M+SPDClusters event classes).

$\Xi^{-}+\overline{\Xi}^{+}$ transverse momentum spectrum in pp collisions at $\sqrt{s}$ = 13 TeV in the high-ZN-energy classification (V0M+SPDClusters event classes).

$K^{0}_{S}$ transverse momentum yield ratios to class III in the high-ZN-energy classification (V0M+SPDClusters event classes) in pp collisions at $\sqrt{s}$ = 13 TeV.

$\Lambda+\overline{\Lambda}$ transverse momentum yield ratios to class III in the high-ZN-energy classification (V0M+SPDClusters event classes) in pp collisions at $\sqrt{s}$ = 13 TeV.

$\Xi^{-}+\overline{\Xi}^{+}$ transverse momentum yield ratios to class III in the high-ZN-energy classification (V0M+SPDClusters event classes) in pp collisions at $\sqrt{s}$ = 13 TeV.

$K^{0}_{S}$ transverse momentum spectrum in pp collisions at $\sqrt{s}$ = 13 TeV in the high-local-multiplicity classification (V0M+SPDClusters event classes).

$\Lambda+\overline{\Lambda}$ transverse momentum spectrum in pp collisions at $\sqrt{s}$ = 13 TeV in the high-local-multiplicity classification (V0M+SPDClusters event classes).

$\Xi^{-}+\overline{\Xi}^{+}$ transverse momentum spectrum in pp collisions at $\sqrt{s}$ = 13 TeV in the high-local-multiplicity classification (V0M+SPDClusters event classes).

$K^{0}_{S}$ transverse momentum yield ratios to class IV in the high-local-multiplicity classification (V0M+SPDClusters event classes) in pp collisions at $\sqrt{s}$ = 13 TeV.

$\Lambda+\overline{\Lambda}$ transverse momentum yield ratios to class IV in the high-local-multiplicity classification (V0M+SPDClusters event classes) in pp collisions at $\sqrt{s}$ = 13 TeV.

$\Xi^{-}+\overline{\Xi}^{+}$ transverse momentum yield ratios to class IV in the high-local-multiplicity classification (V0M+SPDClusters event classes) in pp collisions at $\sqrt{s}$ = 13 TeV.

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised charged-particle multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the standalone classification (V0M event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised $ZN$ signal in pp collisions at $\sqrt{s}$ = 13 TeV in the standalone classification (V0M event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised charged-particle multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the high-local-multiplicity classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised charged-particle multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the low-local-multiplicity classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised charged-particle multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the high-ZN-energy classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised charged-particle multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the low-ZN-energy classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised $ZN$ signal in pp collisions at $\sqrt{s}$ = 13 TeV in the high-local-multiplicity classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised $ZN$ signal in pp collisions at $\sqrt{s}$ = 13 TeV in the low-local-multiplicity classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised $ZN$ signal in pp collisions at $\sqrt{s}$ = 13 TeV in the high-ZN-energy classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised $ZN$ signal in pp collisions at $\sqrt{s}$ = 13 TeV in the low-ZN-energy classification (V0M+SPDClusters event classes).

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{loose}}^{2\ell3j}$ signal region, which uses the invariant mass of the c-tagged jet and the geometrically closest b-tagged jet, $m_{\mathit{cb}}^{\mathrm{min}\Delta R}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{loose}}^{1\ell6ji}$ signal region, which uses the jet multiplicity, $N_{\mathrm{jets}}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{tight}}^{1\ell6ji}$ signal region, which uses the jet multiplicity, $N_{\mathrm{jets}}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{loose}}^{2\ell4ji}$ signal region, which uses the jet multiplicity, $N_{\mathrm{jets}}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{tight}}^{2\ell4ji}$ signal region, which uses the jet multiplicity, $N_{\mathrm{jets}}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Jet multiplicity and tagging criteria for the SRs and CRs in the single-lepton and dilepton channels. All single-lepton regions are defined in the 5-jet-exclusive and 6-jet-inclusive jet selections, all dilepton regions in the 3-jet-exclusive and 4-jet-inclusive jet selections. By construction, $\mathrm{SR}^{2\ell}_{\mathrm{tight}}$ only exists for the 4-jet-inclusive selection. Identical requirements on inclusive and exclusive working points, e.g., one c@22% and one c@11% in the $\mathrm{CR}_1^{1\ell}$, indicate a veto of additional tags at the inclusive working point.

Breakdown of the fiducial $t\bar{t}+{\geq}2c$ and $t\bar{t}+1c$ cross-section uncertainties. Systematic uncertainties are grouped into signal and background modeling, instrumental, and MC statistics categories. The table lists the fractional uncertainty in the measured cross-sections in percent and is estimated from the covariance matrix of the profile likelihood fit, following EPJC 84 (2024) 593. Individual groups of uncertainties can be added in quadrature to obtain the total uncertainty. JES and JER denote the jet energy scale and resolution uncertainties, respectively. The fractional uncertainties in the fitted $t\bar{t}+{\geq}1c$ and $t\bar{t}+\text{light}$ normalization factors are included in the data statistics category. For presentation purposes, all uncertainties are symmetrized.

Measured fiducial cross-section values in comparison with various NLO+PS predictions. The uncertainties in the predictions include independent and simultaneous variations of the $\mu_R$ and $\mu_F$ scales and uncertainties in the choice of the PDF set, but no uncertainties in the parton shower or hadronization. The uncertainties in the measured values include all statistical and systematic uncertainties. For presentation purposes, the quoted measurement and prediction uncertainties are symmetrized.

Measured and predicted values for the $t\bar{t}+{\geq}1b$, $t\bar{t}+{\geq}2c$, and $t\bar{t}+1c$ cross-sections and for the cross-section ratios of these processes to total $t\bar{t}+\text{jets}$ production. The quoted uncertainties in the measurements include statistical and systematic uncertainties. The predictions are taken from the nominal setup using inclusive $t\bar{t}$ Powheg+Pythia8 predictions for the $t\bar{t}+{\geq}2c$, $t\bar{t}+1c$, and $t\bar{t}+\text{light}$ processes, and $t\bar{t}+b\bar{b}$ Powheg+Pythia8 predictions for the $t\bar{t}+{\geq}1b$ process. The first use the 5FS, the latter the 4FS, where the $b$-quarks are treated as massive particles and are included in the matrix element calculation. The uncertainties in the predictions include independent and simultaneous variations of the $\mu_R$ and $\mu_F$ scales and uncertainties in the choice of the PDF set.

Nuisance parameters with a large impact on the fiducial measurement: ranking of the ten most impactful nuisance parameters for the $t\bar{t}+{\geq}2c$ cross-section measurement. The colored boxes indicate the impact of the nuisance parameters, $\Delta\mu$ (top axis), on the corresponding parameter of interest, $\mu$. They are calculated from the covariance matrix using the post-fit up (down) uncertainties of the parameter of interest and the nuisance parameter, $\sigma_{\mathrm{POI}}^{\mathrm{up}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{up}}$ ($\sigma_{\mathrm{POI}}^{\mathrm{down}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{down}}$), and their linear correlation coefficient, $\rho$, following EPJC 84 (2024) 593. The data points and error bars indicate the post-fit values and uncertainties of the nuisance parameters relative to their pre-fit values, $\theta_0$, in units of their pre-fit uncertainties, $\Delta \theta$ (bottom axis). A hollow circle indicates that the nuisance parameter has a Poisson constraint.

Nuisance parameters with a large impact on the fiducial measurement: ranking of the ten most impactful nuisance parameters for the $t\bar{t}+1c$ cross-section measurement. The colored boxes indicate the impact of the nuisance parameters, $\Delta\mu$ (top axis), on the corresponding parameter of interest, $\mu$. They are calculated from the covariance matrix using the post-fit up (down) uncertainties of the parameter of interest and the nuisance parameter, $\sigma_{\mathrm{POI}}^{\mathrm{up}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{up}}$ ($\sigma_{\mathrm{POI}}^{\mathrm{down}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{down}}$), and their linear correlation coefficient, $\rho$, following EPJC 84 (2024) 593. The data points and error bars indicate the post-fit values and uncertainties of the nuisance parameters relative to their pre-fit values, $\theta_0$, in units of their pre-fit uncertainties, $\Delta \theta$ (bottom axis). A hollow circle indicates that the nuisance parameter has a Poisson constraint.

Nuisance parameters with a large impact on the fiducial measurement: ranking of the ten most impactful nuisance parameters for the $t\bar{t}+{\geq}2c$ ratio to total $t\bar{t}+\text{jets}$ production. The colored boxes indicate the impact of the nuisance parameters, $\Delta\mu$ (top axis), on the corresponding parameter of interest, $\mu$. They are calculated from the covariance matrix using the post-fit up (down) uncertainties of the parameter of interest and the nuisance parameter, $\sigma_{\mathrm{POI}}^{\mathrm{up}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{up}}$ ($\sigma_{\mathrm{POI}}^{\mathrm{down}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{down}}$), and their linear correlation coefficient, $\rho$, following EPJC 84 (2024) 593. The data points and error bars indicate the post-fit values and uncertainties of the nuisance parameters relative to their pre-fit values, $\theta_0$, in units of their pre-fit uncertainties, $\Delta \theta$ (bottom axis). A hollow circle indicates that the nuisance parameter has a Poisson constraint.

Nuisance parameters with a large impact on the fiducial measurement: ranking of the ten most impactful nuisance parameters for the $t\bar{t}+1c$ ratio to total $t\bar{t}+\text{jets}$ production. The colored boxes indicate the impact of the nuisance parameters, $\Delta\mu$ (top axis), on the corresponding parameter of interest, $\mu$. They are calculated from the covariance matrix using the post-fit up (down) uncertainties of the parameter of interest and the nuisance parameter, $\sigma_{\mathrm{POI}}^{\mathrm{up}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{up}}$ ($\sigma_{\mathrm{POI}}^{\mathrm{down}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{down}}$), and their linear correlation coefficient, $\rho$, following EPJC 84 (2024) 593. The data points and error bars indicate the post-fit values and uncertainties of the nuisance parameters relative to their pre-fit values, $\theta_0$, in units of their pre-fit uncertainties, $\Delta \theta$ (bottom axis). A hollow circle indicates that the nuisance parameter has a Poisson constraint.

Efficiencies and rejection rates for the $b$- and $c$-jet tagging working points of the $b/c$-tagger, as well as the selection cuts on the discriminants, $\mathcal{D}_b'$ and $\mathcal{D}_c'$ that define the tagging bins. The primes indicate that a standard logistic function is applied to the discriminants. The b@70% and c@22% working points are inclusive of their respective tighter working points. The values were estimated from single-lepton and dilepton $t\bar{t}$ events simulated with Powheg + Pythia 8.

Efficiencies and rejection rates for the $b$- and $c$-jet tagging working points of the $b/c$-tagger, as well as the selection cuts on the discriminants, $\mathcal{D}_b'$ and $\mathcal{D}_c'$ that define the tagging bins. The primes indicate that a standard logistic function is applied to the discriminants. The b@70% and c@22% working points are inclusive of their respective tighter working points. The values were estimated from single-lepton and dilepton $t\bar{t}$ events simulated with Powheg + Pythia 8.

Pre-fit event yields in the twelve single-lepton and dilepton control regions (CRs) of the analysis. The quoted uncertainties include both statistical and systematic components, but all normalization factors are fixed to their pre-fit values with no uncertainties. All uncertainties are assumed to be uncorrelated. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes. The "Fakes" category includes events from misidentified leptons and non-prompt leptons, estimated through the data-driven matrix method in the single-lepton channel and through MC predictions in the dilepton channel.

Pre-fit event yields in the seven signal regions (SRs) of the analysis. The quoted uncertainties include both statistical and systematic components, but all normalization factors are fixed to their pre-fit values with no uncertainties. All uncertainties are assumed to be uncorrelated. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes. The "Fakes" category includes events from misidentified leptons and non-prompt leptons, estimated through the data-driven matrix method in the single-lepton channel and through MC predictions in the dilepton channel.

Post-fit event yields in the twelve single-lepton and dilepton control regions (CRs) of the analysis. The quoted uncertainties include both statistical and systematic components, and post-fit correlations between uncertainties are taken into account. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes. The "Fakes" category includes events from misidentified leptons and non-prompt leptons, estimated through the data-driven matrix method in the single-lepton channel and through MC predictions in the dilepton channel.

Post-fit event yields in the seven signal regions (SRs) of the analysis. The quoted uncertainties include both statistical and systematic components, and post-fit correlations between uncertainties are taken into account. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes. The "Fakes" category includes events from misidentified leptons and non-prompt leptons, estimated through the data-driven matrix method in the single-lepton channel and through MC predictions in the dilepton channel.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of the number of $b$-jets compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of the number of $l/c$-jets compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $H_{\text{T}}^{\text{had}}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $\Delta R_{\text{avg}}^{bb}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{1})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{1}^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{2}^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{3})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{1}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $m(b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $m(bb^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(bb^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $\Delta R(e\mu bb^{\text{top}}, b_{1}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets and at least one $l/c$-jet as a function of $\Delta R(e\mu bb^{\text{top}}, l/c\text{-jet}_{1})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets and at least one $l/c$-jet as a function of $p_{\text{T}}(l/c\text{-jet}_{1})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets and at least one $l/c$-jet as a function of $p_{\text{T}}(l/c\text{-jet}_{1}) - p_{\text{T}}(b_{1}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $m(bb^{\text{min}\Delta R})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(bb^{\text{min}\Delta R})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $m(bb^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(bb^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(b_{3})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(b_{1}^{\text{add}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $\Delta R(b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $m(e\mu bb^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(l/c\text{-jet}_{1})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $\Delta\eta_{\text{max}}^{jj}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $H_{\text{T}}^{\text{all}}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $m(e\mu b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(b_{1})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(b_{2})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(b_{1}^{\text{top}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(b_{2}^{\text{top}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{1})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{1}^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{2}^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{3})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{4})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{1}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{2}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $m(b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $m(bb^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(bb^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $H_{\text{T}}^{\text{all}}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $m(e\mu b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $m(e\mu bb^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $H_{\text{T}}^{\text{had}}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $\text{min}\Delta R(bb)$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $\Delta R(e\mu bb^{\text{top}}, b_{1}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $\Delta R_{\text{avg}}^{bb}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $\Delta\eta_{\text{max}}^{jj}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of the number of $l/c$-jets compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets and at least one $l/c$-jet as a function of $p_{\text{T}}(l/c\text{-jet}_{1})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets and at least one $l/c$-jet as a function of $|\eta(l/c\text{-jet}_{1})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets and at least one $l/c$-jet as a function of $\Delta R(e\mu bb^{\text{top}}, l/c\text{-jet}_{1})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets and at least one $l/c$-jet as a function of $p_{\text{T}}(l/c\text{-jet}_{1}) - p_{\text{T}}(b_{1}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{1})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{2})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{1}^{\text{top}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at leastfour $b$-jets as a function of $|\eta(b_{2}^{\text{top}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{3})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{4})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{1}^{\text{add}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{2}^{\text{add}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

The measured normalised differential cross-section as a function of $N_{b-\text{jets}}$ in the $e\mu+\geq2b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $H_{\text{T}}^{\text{had}}$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $H_{\text{T}}^{\text{all}}$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R_{\text{avg}}^{bb}$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta\eta_{\text{max}}^{jj}$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1}^{\text{top}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{2})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{2}^{\text{top}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{3})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1}^{\text{add}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{1})|$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{1}^{\text{top}})|$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{2})|$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{2}^{\text{top}})|$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{3})|$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{1}^{\text{add}})|$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(b_{1}b_{2})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1}b_{2})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(bb^{\text{top}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(bb^{\text{top}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(e\mu b_{1}b_{2})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(e\mu bb^{\text{top}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(b_{1}b_{2})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $N_{l/c-\text{jets}}$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(e\mu b_{1}b_{2},b_{3})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(e\mu bb^{\text{top}}, b_{1}^{\text{add}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(e\mu bb^{\text{top}},l/c-\text{jet})$ in the $e\mu+\geq3b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(l/c\text{-jet}_{1}) - p_{\text{T}}(b_{1}^{\text{add}})$ in the $e\mu+\geq3b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(l/c\text{-jet}_{1})|$ in the $e\mu+\geq3b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(l/c\text{-jet}_{1})$ in the $e\mu+\geq3b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $H_{\text{T}}^{\text{had}}$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $H_{\text{T}}^{\text{all}}$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R_{\text{avg}}^{bb}$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta\eta_{\text{max}}^{jj}$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1}^{\text{top}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{2})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{2}^{\text{top}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{3})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1}^{\text{add}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{4})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{2}^{\text{add}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{1})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{1}^{\text{top}})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{2})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{2}^{\text{top}})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{3})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{1}^{\text{add}})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{4})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{2}^{\text{add}})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(b_{1}b_{2})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1}b_{2})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(bb^{\text{top}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(bb^{\text{top}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(e\mu b_{1}b_{2})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(e\mu bb^{\text{top}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(bb^{\text{min}\Delta R})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(bb^{\text{min}\Delta R})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(bb^{\text{add}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(bb^{\text{add}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\text{min}\Delta R(bb)$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(b_{1}b_{2})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $N_{l/c-\text{jets}}$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(e\mu b_{1}b_{2},b_{3})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(e\mu bb^{\text{top}}, b_{1}^{\text{add}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(e\mu bb^{\text{top}}, l/c\text{-jet}_{1})$ in the $e\mu+\geq4b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(l/c\text{-jet}_{1}) - p_{\text{T}}(b_{1}^{\text{add}})$ in the $e\mu+\geq4b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(l/c\text{-jet}_{1})|$ in the $e\mu+\geq4b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(l/c\text{-jet}_{1})$ in the $e\mu+\geq4b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $N_{b-\text{jets}}$ in the phase space with at least two b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $N_{b-\text{jets}}$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $H_{\text{T}}^{\text{had}}$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $H_{\text{T}}^{\text{all}}$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R_{\text{avg}}^{bb}$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta\eta_{\text{max}}^{jj}$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1}^{\text{top}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{2})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{2}^{\text{top}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{3})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1}^{\text{add}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{1})|$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{1}^{\text{top}})|$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{2}^{\text{top}})|$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{2}^{\text{top}})|$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{3})|$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{1}^{\text{add}})|$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(b_{1}b_{2})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1}b_{2})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(bb^{\text{top}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(bb^{\text{top}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(e\mu b_{1}b_{2})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(e\mu bb^{\text{top}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(b_{1}b_{2})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $N_{l/c-\text{jets}}$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(e\mu b_{1}b_{2},b_{3})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(e\mu bb^{\text{top}}, b_{1}^{\text{add}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(e\mu bb^{\text{top}},l/c-\text{jet})$ in the phase space with at least three b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(l/c\text{-jet}_{1}) - p_{\text{T}}(b_{1}^{\text{add}})$ in the phase space with at least three b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(l/c\text{-jet}_{1})|$ in the phase space with at least three b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(l/c\text{-jet}_{1})$ in the phase space with at least three b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $H_{\text{T}}^{\text{had}}$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $H_{\text{T}}^{\text{all}}$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R_{\text{avg}}^{bb}$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta\eta_{\text{max}}^{jj}$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1}^{\text{top}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{2})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{2}^{\text{top}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{3})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1}^{\text{add}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{4})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{2}^{\text{add}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{1})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{1}^{\text{top}})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{2})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{2}^{\text{top}})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{3})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{1}^{\text{add}})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{4})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{2}^{\text{add}})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(b_{1}b_{2})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1}b_{2})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(bb^{\text{top}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(bb^{\text{top}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(e\mu b_{1}b_{2})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(e\mu bb^{\text{top}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(bb^{\text{min}\Delta R})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(bb^{\text{min}\Delta R})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(bb^{\text{add}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(bb^{\text{add}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\text{min}\Delta R(bb)$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(b_{1}b_{2})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $N_{l/c-\text{jets}}$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(e\mu b_{1}b_{2},b_{3})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(e\mu bb^{\text{top}}, b_{1}^{\text{add}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(e\mu bb^{\text{top}}, l/c\text{-jet}_{1})$ in the phase space with at least four b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(l/c\text{-jet}_{1}) - p_{\text{T}}(b_{1}^{\text{add}})$ in the phase space with at least four b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(l/c\text{-jet}_{1})|$ in the phase space with at least four b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(l/c\text{-jet}_{1})$ in the phase space with at least four b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

Correlations of the ttH signal-strength modifiers in the different Higgs pT bins of the STXS measurement.

Observed and expected upper 95% CL limit on the tH signal strength modifier $\mu_{tH}$ with $\mu_{ttH}=1$ for different channels and years, where the uncertainties are uncorrelated between the channels and years, and in their combination. The one and two standard deviation confidence intervals on the expected limit are also listed (sigma).

Likelihood-ratio test statistic as a function of the ttH and tH signal strength modifiers $\mu_{ttH}$ and $\mu_{tH}$. The observed best fit point is $(\mu_{tH}, \mu_{ttH}) = (-3.83, +0.35)$.

Observed and expected likelihood ratio test statistic as a function of $\kappa_{t}$ and $\kappa_{V}$. The observed best fit point is $(\kappa_{t}, \kappa_{V}) = (+0.59, +1.40)$.

Observed and expected likelihood ratio test statistic as a function of $\kappa_{t}$ for $\kappa_{V}=1$. The observed best fit point is $\kappa_{t} = +0.54$.

Observed and expected likelihood ratio test statistic as a function of $\kappa_{t}$ and $\tilde{\kappa}_{t}$ with $\kappa_{V}=1$. The observed best fit point is $(\kappa_{t}, \tilde{\kappa}_{t}) = (+0.53, 0.00)$.

Observed and expected likelihood ratio test statistic as a function of $f_{CP}$. The overall modifier of the top-Higgs coupling strength is profiled such that the likelihood attains its minimum at each point in the plane, and $\kappa_{V}$ is set to 1. The observed best fit point is $f_{CP} = 0.00$.

Observed and expected likelihood ratio test statistic as a function of $\cos(\alpha)$. The overall modifier of the top-Higgs coupling strength is profiled such that the likelihood attains its minimum at each point in the plane, and $\kappa_{V}$ is set to 1. The observed best fit point is $\cos(\alpha) = +1.00$.

Observed and expected likelihood ratio test statistic as a function of $\kappa_{t}$ and $\tilde{\kappa}_{t}$ with $\kappa_{V}=1$, in the $H \rightarrow bb$, $H\rightarrow\gamma\gamma$, $H\rightarrow ZZ$, $H\rightarrow WW$, and $H\rightarrow\tau\tau$ decay channels and in their combination.

Observed and expected likelihood ratio test statistic as a function of $f_{CP}$, in the $H \rightarrow bb$, $H\rightarrow\gamma\gamma$, $H\rightarrow ZZ$, $H\rightarrow WW$, and $H\rightarrow\tau\tau$ decay channel and in their combination. The overall modifier of the top-Higgs coupling strength is profiled such that the likelihood attains its minimum at each point in the plane, and $\kappa_{V}$ is set to 1.

Observed and expected likelihood ratio test statistic as a function of $\cos(\alpha)$, in the $H \rightarrow bb$, $H\rightarrow\gamma\gamma$, $H\rightarrow ZZ$, $H\rightarrow WW$, and $H\rightarrow\tau\tau$ decay channels and in their combination. The overall modifier of the top-Higgs coupling strength is profiled such that the likelihood attains its minimum at each point in the plane, and $\kappa_{V}$ is set to 1.