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.

Jet pairing efficiencies over the parameter space for the wide-width heavy resonance signals. The pairing efficiency is evaluated in the 6$b$ region when a correct pairing is possible — that is, the six leading jets are geometrically matched to truth-level b-quarks.

Distributions of scores for nonresDNN in 6b data and the background prediction after background-only fits to the observed data. The Low- and High-Score regions are included. Benchmark signal models (normalized to the background) are overlaid: SM HHH signal and non-resonant TRSM signal with $(m_X, m_S)=(400,200)$ GeV.

Distributions of scores for resDNN in 6b data and the background prediction after background-only fits to the observed data. The Low- and High-Score regions are included. A benchmark signal model (normalized to the background) is overlaid: resonant TRSM signal with $(m_X, m_S)=(500,350)$ GeV.

Distributions of scores for heavyresDNN in 6b data and the background prediction after background-only fits to the observed data. The Low- and High-Score regions are included. Benchmark signal models (normalized to the background) are overlaid: heavy-resonant signals with (mX, mS)=(900,325) GeV with a 1% or 20% width.

Expected 95% CL $HHH$ cross section upper limits for the phase space within the perturbative unitarity bounds of the TRSM signals. While the points were generated under Benchmark point 3 of the TRSM, they can also be interpreted in the DM-CPV model.

Observed 95% CL $HHH$ cross section upper limits for the phase space within the perturbative unitarity bounds of the TRSM signals. While the points were generated under Benchmark point 3 of the TRSM, they can also be interpreted in the DM-CPV model.

Expected 95% CL $HHH$ cross section upper limits for the narrow-width (1%) heavy resonance signals.

Expected 95% CL $HHH$ cross section upper limits for the wide-width (20%) heavy resonance signals.

Observed 95% CL $HHH$ cross section upper limits for the narrow-width (1%) heavy resonance signals.

Observed 95% CL $HHH$ cross section upper limits for the wide-width (20%) heavy resonance signals.

Expected 95% CL constraints on $\kappa_3$ and $\kappa_4$, the ratios of the Higgs tri-linear and quartic self-couplings to their predicted SM values.

Expected 68% CL constraints on $\kappa_3$ and $\kappa_4$, the ratios of the Higgs tri-linear and quartic self-couplings to their predicted SM values.

Observed 95% CL constraints on $\kappa_3$ and $\kappa_4$, the ratios of the Higgs tri-linear and quartic self-couplings to their predicted SM values.

Observed 68% CL constraints on $\kappa_3$ and $\kappa_4$, the ratios of the Higgs tri-linear and quartic self-couplings to their predicted SM values.

Unitarity limits, calculated in Eur. Phys. J. C 84 (2024) 366, are overlaid in the region bounded by the gray dashed line.

Overall acceptance x efficiency for events entering the 6$b$ region over the parameter space for the TRSM signals. Generator weights, and weights associated to pileup modeling corrections are applied. Jet trigger scale factors, NNJVT scale factors and scale factors associated to the $b$-tagging calibration are applied.

Overall acceptance x efficiency for events entering the 6$b$ region over the parameter space for the narrow-width heavy resonance signals. Generator weights, and weights associated to pileup modeling corrections are applied. Jet trigger scale factors, NNJVT scale factors and scale factors associated to the $b$-tagging calibration are applied.

Overall acceptance x efficiency for events entering the 6$b$ region over the parameter space for the wide-width heavy resonance signals. Generator weights, and weights associated to pileup modeling corrections are applied. Jet trigger scale factors, NNJVT scale factors and scale factors associated to the $b$-tagging calibration are applied.

Expected one-dimensional likelihood scan over the Higgs self-coupling modifier $\kappa_3$. Both of the modifiers are profiled in the profile-likelihood fit, and the likelihood scan keeps $\kappa_4$ fixed at one.

Observed one-dimensional likelihood scan over the Higgs self-coupling modifier $\kappa_3$. Both of the modifiers are profiled in the profile-likelihood fit, and the likelihood scan keeps $\kappa_4$ fixed at one.

Expected one-dimensional likelihood scan over the Higgs self-coupling modifier $\kappa_4$. Both of the modifiers are profiled in the profile-likelihood fit, and the likelihood scan keeps $\kappa_3$ fixed at one.

Observed one-dimensional likelihood scan over the Higgs self-coupling modifier $\kappa_4$. Both of the modifiers are profiled in the profile-likelihood fit, and the likelihood scan keeps $\kappa_3$ fixed at one.

Summary of the impact on the analysis of different sources of systematic uncertainty. The impact is estimated as the absolute value of the difference between the expected upper limit with all sources of uncertainty included and the upper limit with a subset of systematic uncertainties excluded. In the left-hand column, more heavily indented rows are a sub-category of the less indented row above.

Cutflow for selected signal points. Generator weights, and weights associated to pileup modeling corrections are applied from the second row onwards. These weight the SM events assuming the SM cross section, whilst the TRSM and Generic resonance events are weighted to 50 fb. Jet trigger scale factors are applied from the third row onwards. NNJVT scale factors are applied from the fourth row onwards. Scale factors associated to the b-tagging calibration are applied in the last three rows.

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.

Total invariant mass distribution $m_{4\ell}$ in Signal Region 2.

Selected events in the ($m_{4\ell}$, $\left\langle m_{\ell\ell}\right\rangle$) plane, SR1

Selected events in the ($m_{4\ell}$, $\left\langle m_{\ell\ell}\right\rangle$) plane, SR2

Expected background events in the $(m_{4\ell}, \langle m_{\ell\ell}\rangle)$ plane, SR1

Expected signal shape for $(m_{S}, m_{Z_d}) = (110~\text{GeV},30~\text{GeV}$) in the $(m_{4\ell}, \langle m_{\ell\ell}\rangle)$ plane, SR1. The signal distributions are normalized to unit volume.

Expected background events in the $(m_{4\ell}, \langle m_{\ell\ell}\rangle)$ plane, SR2. Events that satisfy $|m_{ab,cd} - m_Z| < 8~\text{GeV}$ are removed by the Z-veto requirement, however not all events lying behind the mask in the plot, corresponding to $\langle m_{\ell\ell}\rangle$ values between 83 and 99 GeV, are deleted.

Expected signal shape for $(m_{S}, m_{Z_{d}}) = (750~\text{GeV},150~\text{GeV}$) in the $(m_{4\ell}, \langle m_{\ell\ell}\rangle)$ plane, SR2. The signal distributions are normalized to unit volume. Events that satisfy $|m_{ab,cd} - m_Z| < 8~\text{GeV}$ are removed by the Z-veto requirement, however not all events lying behind the mask in the plot, corresponding to $\langle m_{\ell\ell}\rangle$ values between 83 and 99 GeV, are deleted.

Expected 95% CL limit on $\sigma \mathcal{B}(S\rightarrow Z_{d}Z_{d} \rightarrow 4\ell)$, SR1

Observed 95% CL limit on $\sigma \mathcal{B}(S\rightarrow Z_{d}Z_{d} \rightarrow 4\ell)$, SR1

Observed / Expected 95% CL limit ratio in SR1

Expected 95% CL limit on $\mathcal{B}(S\rightarrow Z_{d}Z_{d} \rightarrow 4\ell)$ with Z-veto, SR2

Observed 95% CL limit on $\sigma \mathcal{B}(S\rightarrow Z_{d}Z_{d} \rightarrow 4\ell)$ with Z-veto, SR2

Observed / Expected 95% CL limit ratio with Z-veto, SR2

Local $p_0$-values in the ($m_{S}, m_{Z_{d}}$) plane, SR1

Local $p_0$-values in the ($m_{S}, m_{Z_{d}}$) plane, SR2 with Z-veto applied

Observed data in Signal Region 1 in the $(m_{4\ell}, \langle m_{\ell\ell}\rangle)$ plane.

Observed data in Signal Region 2 in the $(m_{4\ell}, \langle m_{\ell\ell}\rangle)$ plane. Events that satisfy $|m_{ab,cd} - m_Z| < 8~\text{GeV}$ are removed by the Z-veto requirement, however not all events lying behind the mask in the plot, corresponding to <m_ll> values between 83 and 99 GeV, are deleted.

Distributions of the average dilepton mass $\left\langle m_{\ell\ell}\right\rangle = \frac{1}{2}\left(m_{ab} + m_{cd}\right)$ for Signal Region 1 in the $4e$ final state.

Distributions of the average dilepton mass $\left\langle m_{\ell\ell}\right\rangle = \frac{1}{2}\left(m_{ab} + m_{cd}\right)$ for Signal Region 1 in the $2e2\mu$ final state.

Distributions of the average dilepton mass $\left\langle m_{\ell\ell}\right\rangle = \frac{1}{2}\left(m_{ab} + m_{cd}\right)$ for Signal Region 1 in the $4\mu$ final state.

Distributions of the average dilepton mass $\left\langle m_{\ell\ell}\right\rangle = \frac{1}{2}\left(m_{ab} + m_{cd}\right)$ for Signal Region 2 in the $4e$ final state.

Distributions of the average dilepton mass $\left\langle m_{\ell\ell}\right\rangle = \frac{1}{2}\left(m_{ab} + m_{cd}\right)$ for Signal Region 2 in the $2e2\mu$ final state.

Distributions of the average dilepton mass $\left\langle m_{\ell\ell}\right\rangle = \frac{1}{2}\left(m_{ab} + m_{cd}\right)$ for Signal Region 2 in the $4\mu$ final state.

Distributions of the total invariant mass $m_{4\ell}$ for Signal Region 1 in the $4e$ final state. The SR1 requirement $m_{4\ell} < 115~\text{GeV}$ is not applied here.

Distributions of the total invariant mass $m_{4\ell}$ for Signal Region 1 in the $2e2\mu$ final state. The SR1 requirement $m_{4\ell} < 115~\text{GeV}$ is not applied here.

Distributions of the total invariant mass $m_{4\ell}$ for Signal Region 1 in the $4\mu$ final state. The SR1 requirement $m_{4\ell} < 115~\text{GeV}$ is not applied here.

Distributions of the total invariant mass $m_{4\ell}$ for Signal Region 2 in the $4e$ final state.

Distributions of the total invariant mass $m_{4\ell}$ for Signal Region 2 in the $2e2\mu$ final state.

Distributions of the total invariant mass $m_{4\ell}$ for Signal Region 2 in the $4\mu$ final state.

90% CL WIMP SDn cross sections, including sensitivities and nuclear structure uncertainties

90% CL WIMP SDp cross sections, including sensitivities and nuclear structure uncertainties

90% CL WIMP SDp cross sections, including sensitivities and nuclear structure uncertainties

Data selection efficiency as a function of nuclear recoil energy

Data selection efficiency as a function of nuclear recoil energy

Data points used in analysis in log_10(S2)-S1 space

Data points used in analysis in log_10(S2)-S1 space

The $G_{\text{i}}^{\text{Strips}}$ distribution in the PASS region for events passing the event selection and with $55 < p_{\mathrm{T}} < 200$ GeV.

The $G_{\text{i}}^{\text{Strips}}$ distribution in the FAIL region for events passing the event selection and with $p_{\mathrm{T}} > 200$ GeV.

The $G_{\text{i}}^{\text{Strips}}$ distribution in the PASS region for events passing the event selection and with $p_{\mathrm{T}} > 200$ GeV.

Mass spectrum predicted in the signal region defined by $G_{\text{i}}^{\text{Strips}} > 0.22$ and $p_{\mathrm{T}} > 70$ GeV.

Cross section limits for gluino and supersymmetric top R-hadrons for both background prediction methods.

Cross section limits for supersymmetric tau models for both background prediction methods.

Cross section limits for DY-produced tau prime models (single and multicharged) for both background prediction methods.

Cross section limits for Z prime to multicharged tau prime model for both background prediction methods. All Z prime models assume a narrow width, and a 100% branching fraction to 600 GeV tau primes.

2D exclusion showing the observed cross section limit as a function of the multicharged tau prime mass and Z prime mass for the ionization method.

2D exclusion showing the observed cross section limit as a function of the multicharged tau prime mass and Z prime mass for the mass method.

Cumulative selection efficiency for the data and for two signal hypotheses.

Expected and observed mass limits obtained using 2017-2018 data for various HSCP candidate models,for the two background estimate methods.

Mass windows used in the mass method as a function of the signal target mass for the signal samples assuming a charge of $1e$ and $2e$, as appropriate to the signal model.

Mass spectrum predicted in the validation region defined by $0.018 < G_{\text{i}}^{\text{Strips}} < 0.057$ and $p_{\mathrm{T}} > 70$ GeV. The thresholds used in the $G_{\text{i}}^{\text{Strips}}$ requirement represent the 50% and 90% quantile of the distribution.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $\text{abs(}\eta\text{)}<0.3$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $0.3<\text{abs(}\eta\text{)}<0.6$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $0.6<\text{abs(}\eta\text{)}<0.9$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $0.9<\text{abs(}\eta\text{)}<1.2$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $1.2<\text{abs(}\eta\text{)}<2.1$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $2.1<\text{abs(}\eta\text{)}<2.4$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $\text{abs(}\eta\text{)}<0.3$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $0.3<\text{abs(}\eta\text{)}<0.6$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $0.6<\text{abs(}\eta\text{)}<0.9$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $0.9<\text{abs(}\eta\text{)}<1.2$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $1.2<\text{abs(}\eta\text{)}<2.1$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $2.1<\text{abs(}\eta\text{)}<2.4$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency as a function of $\beta$, for the $\tilde{g}$ signals.

Trigger efficiency as a function of $\beta$, for the $\tilde{\tau}$ signals.

Cross section limits for $\tilde{g}$ R-hadrons obtained with the ionization method.

Cross section limits for $\tilde{g}$ R-hadrons obtained with the mass method.

Cross section limits for $\tilde{t}$ R-hadrons obtained with the ionization method.

Cross section limits for $\tilde{t}$ R-hadrons obtained with the mass method.

Cross section limits for $\tilde{\tau}$ obtained with the ionization method.

Cross section limits for $\tilde{\tau}$ obtained with the mass method.

Cross section limits for $\tilde{\tau}$ production within the GMSB SPS7 modelobtained with the ionization method.

Cross section limits for $\tilde{\tau}$ production within the GMSB SPS7 model obtained with the mass method.

Cross section limits for DY-produced $\tau'$ with $abs(Q) = 1e$ with the ionization method.

Cross section limits for DY-produced $\tau'$ with $abs(Q) = 1e$ with the mass method.

Cross section limits for DY-produced $\tau'$ with $abs(Q) = 2e$ with the ionization method.

Cross section limits for DY-produced $\tau'$ with $abs(Q) = 2e$ with the mass method.

Cross section limits for the production of $Z'$ boson decaying into a pair of $\tau'$ fermions of charge $2e$ (with a branching fraction equal to 1 and a fixed $\tau'$ mass of 600 GeV), obtained with the ionization method.

Cross section limits for the production of $Z'$ boson decaying into a pair of $\tau'$ fermions of charge $2e$ (with a branching fraction equal to 1 and a fixed $\tau'$ mass of 600 GeV), obtained with the mass method.

Two-dimensional exclusion showing the observed cross section limit as a function of the masses of the $\tau'$ (on the $x$ axis) and of the $Z'$ boson (on the $y$ axis), for the ionization method.

Two-dimensional exclusion showing the observed cross section limit as a function of the masses of the $\tau'$ (on the $x$ axis) and of the $Z'$ boson (on the $y$ axis), for the mass method.

Expected background yield, expected signal yield, and observed data for the mass method for 2017.

Expected background yield, expected signal yield, and observed data for the mass method for 2018.

Selection efficieny for GMSB supersymmetric $\tau$.

Selection efficiency for gluino R-hardon samples.

Selection efficiency for supersymmetric top R-hadrons.

Selection efficiency for pair-produced supersymmetric $\tau$.

Selection efficiency for $\tau'$ with $|Q| = 1e$.

Selection efficiecny for $\tau'$ with $|Q| = 2e$.

Selection efficiency for $Z'$ (mass = 3 TeV) going to $\tau'$ with $|Q| = 2e$.

Selection efficiency for $Z'$ (mass 5 TeV) goes to $\tau'$ with $|Q| = 2e$.

Selection efficiency for $Z'$ going to $\tau'$ (mass 600 GeV) with $|Q| = 2e$.

Selection efficiency for $Z'$ going to $\tau'$ (mass 1 TeV) with charge $|Q| = 2e$.

Post-fit $m_\mathrm{MMC}$ distribution in the high-mass SR for the $m_A = 90\,\mathrm{GeV}$ signal mass hypothesis. $m_\mathrm{MMC}$ is the mass reconstructed by the Missing Mass Calculator. Processes contributing to the background Others are $Z/\gamma^* \rightarrow ee/\mu\mu$ and SM Higgs. The subscript on the $A\to\tau\tau$ process indicates the mass of the $A$ boson. otal includes all backgrounds and the signal process. The high-mass Signal Region is defined as: - 1 electron and 1 muon with opposite charge - $p_\mathrm{T}$ requirements of the leptons are a combination of the following: - $p_\mathrm{T}^e > 18\,\mathrm{GeV}$ and $p_\mathrm{T}^\mu > 15\,\mathrm{GeV}$ or - $p_\mathrm{T}^e > 10\,\mathrm{GeV}$ and $p_\mathrm{T}^\mu > 25\,\mathrm{GeV}$ or - $p_\mathrm{T}^e > 27\,\mathrm{GeV}$ and $p_\mathrm{T}^\mu > 10\,\mathrm{GeV}$ - $\vert \eta_e \vert < 2.47$, excluding $1.37 < \vert \eta_e \vert < 1.52$ - $\vert \eta_\mu \vert < 2.7$ - no jets with $b$-quarks - $\Delta R_{\ell\ell} < 1.0$ - $E_\mathrm{T}^\mathrm{miss} > 30\,\mathrm{GeV}$ - $m_\mathrm{T}^\mathrm{tot} = \sqrt{\left(p_\mathrm{T}^e+p_\mathrm{T}^\mu+E_\mathrm{T}^\mathrm{miss}\right)^2-\left(\vec{p}_\mathrm{T}^{\,e}+\vec{p}_\mathrm{T}^{\,\mu}+\vec{E}_\mathrm{T}^{\,\mathrm{miss}}\right)^2} < 65\,\mathrm{GeV}$ - $35\,\mathrm{GeV} < m_\mathrm{MMC} < 130\,\mathrm{GeV}$

Expected and observed $95\%$ CL limits on the production cross-section for $gg\rightarrow A$ times the branching ratio for $A$ decaying into two $\tau$-leptons for $A$ boson masses ranging from $20$ to $90\,\mathrm{GeV}$.

Expected and observed $95\%$ CL limits on the production cross-section for $gg\rightarrow A$ times the branching ratio for $A$ decaying into two $\tau$-leptons for $A$ boson masses ranging from $20$ to $90\,\mathrm{GeV}$.

Expected and observed $95\%$ CL limits on the absolute value of the up-type quark coupling parameter $|\zeta_\mathrm{u}|$ for masses of the $A$ boson ranging from $20$ to $90\,\mathrm{GeV}$.

Expected and observed $95\%$ CL limits on the absolute value of the up-type quark coupling parameter $|\zeta_\mathrm{u}|$ for masses of the $A$ boson ranging from $20$ to $90\,\mathrm{GeV}$.

Observed $p_0$ values for the fits on the production cross-section for $gg\rightarrow A$ times the branching ratio for $A$ decaying into two $\tau$-leptons for $A$ boson masses from $20$ to $90\,\mathrm{GeV}$. For the fit on the up-type quark coupling parameter the $p_0$ values agree up to a precision of $O(10^{-5})$.

Observed $p_0$ values for the fits on the production cross-section for $gg\rightarrow A$ times the branching ratio for $A$ decaying into two $\tau$-leptons for $A$ boson masses from $20$ to $90\,\mathrm{GeV}$. For the fit on the up-type quark coupling parameter the $p_0$ values agree up to a precision of $O(10^{-5})$.