Observed 95% CL upper limits on the production cross section times branching ratio of the A → ZH → Ztt process in the (mA , mH) plane.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tan β, mA) parameter space at mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tan β, mA) parameter space at mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tanβ, cos(β−α)) parameter space at mA = 600 GeV and mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tanβ, cos(β−α)) parameter space at mA = 600 GeV and mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hq}^{(3)}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hu}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hu}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hd}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hd}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hu}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hu}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hd}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hd}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $c_{Hd}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $c_{Hd}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $g^{(2)}_{ZZ}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $g^{(2)}_{ZZ}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Results for the $P_3$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the $P_4^\prime$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the $P_5^\prime$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the $P_6^\prime$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the $P_8^\prime$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the CP averaged observables $F_L$ and $P_1$–$P_8^\prime$. The first uncertainties are statistical and the second systematic.

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 1.1 < $q^2$ < 2 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 2 < $q^2$ < 4.3 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 4.3 < $q^2$ < 6 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 6 < $q^2$ < 8.68 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 10.09 < $q^2$ < 12.86 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 14.18 < $q^2$ < 16 GeV$^2$

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.

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 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})$.

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.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH3 at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 and ISR starting scale SpaceShower:pT0Ref = 1 GeV at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH3 and ISR starting scale SudakovCommon:pTmin=0.7 GeV at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-proton center-of-mass energy 38.8 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-proton center-of-mass energy 8000 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-nucleon center-of-mass energy 8160 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-nucleon center-of-mass energy 13000 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-proton center-of-mass energy 38.8 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-proton center-of-mass energy 8000 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-nucleon center-of-mass energy 8160 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-proton center-of-mass energy 13000 GeV versus the dilepton mass range of the process.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 38.8 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 62 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 200 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-anti-proton center-of-mass energy 1800 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-anti-proton center-of-mass energy 1960 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 2760 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 8000 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-nucleon center-of-mass energy 8160 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 13000 GeV when fitting to the CMS measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 13000 GeV when fitting to the LHCb measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 38.8 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 62 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 200 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-anti-proton center-of-mass energy 1800 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-anti-proton center-of-mass energy 1960 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 2760 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 8000 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-nucleon center-of-mass energy 8160 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 13000 GeV when fitting to the CMS measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 13000 GeV when fitting to the LHCb measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 38.8 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 62 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 200 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-anti-proton center-of-mass energy 1800 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-anti-proton center-of-mass energy 1960 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 2760 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 8000 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-nucleon center-of-mass energy 8160 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 13000 GeV when fitting to the CMS measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 13000 GeV when fitting to the LHCb measurement.

Distributions of the $\mathrm{NN^{SvsB}}$ output for data and the expected background after the likelihood fit in the $SR_{cQu --}$ signal region. The post-fit background expectations are shown as filled histograms, the combined pre-fit background expectations are shown as dashed lines. The signal distribution using the Wilson coefficient values $c_{tu}^{(1)}=0.04$, $c_{Qu}^{(1)}=0.1$, $c_{Qu}^{(8)}=0.1$ is shown with a dotted line, normalized to the same number of events as the background.

Summary of the observed and predicted number of events in the five 2$\ell$ control regions. The background prediction is shown after the combined likelihood fit to data under the signal-plus-background hypothesis across all control regions and signal regions. The uncertainties in the total yields are smaller than the sum in quadrature of the uncertainties in the individual contributions due to anti-correlations resulting from the likelihood fit. Int Conv refers to events with photon conversions in the inner detector, while Mat Conv refers to events with photon conversions in the detector material. QMisID refers to events with charge-misidentified lepton.

Summary of the observed and predicted number of events in the four 3$\ell$ control regions. The background prediction is shown after the combined likelihood fit to data under the signal-plus-background hypothesis across all control regions and signal regions. The uncertainties in the total yields are smaller than the sum in quadrature of the uncertainties in the individual contributions due to anti-correlations resulting from the likelihood fit. Int Conv refers to events with photon conversions in the inner detector, while Mat Conv refers to events with photon conversions in the detector material. QMisID refers to events with charge-misidentified lepton.

Comparison of the observed limits on the three WCs to the limit obtained by the ATLAS Run~1 analysis of events with $b$-tagged jets and a same-sign lepton pair [JHEP10(2015)150]. The cross-section limits from that analysis are converted to the limits on the WCs by using the fitted EFT parametrization. The bounds on the WCs are shown at the 68$\%$ CL (solid) and/or 95$\%$ CL (dashed) levels. For the ATLAS Run~1 analysis only the bounds at 95$\%$ CL are available. The vertical bar represents the SM prediction. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$.

Observed lower limits at 95$\%$ confidence level on the scale of new physics $\Lambda$ for Wilson coefficient values of 0.01, 1 and $4\pi^2$. The limits are compared with the limits obtained by the ATLAS Run~1 analysis of events with $b$-tagged jets and a same-sign lepton pair.

Observed and expected limits at 95$\%$ confidence level on the cross-section of $\sigma(pp \rightarrow tt)$.

2D likelihood scans of the Wilson coefficients $c_{tu}^{(1)}$ versus $c_{Qu}^{(1)}$. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$. In the table the $2\Delta$NLL values are shown which correspond to the minimized negative log likelihood for a given Wilson coefficient pair.

2D likelihood scans of the Wilson coefficients $c_{tu}^{(1)}$ versus $c_{Qu}^{(8)}$. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$. In the table the $2\Delta$NLL values are shown which correspond to the minimized negative log likelihood for a given Wilson coefficient pair.

2D likelihood scans of the Wilson coefficients $c_{Qu}^{(1)}$ versus $c_{Qu}^{(8)}$. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$. In the table the $2\Delta$NLL values are shown which correspond to the minimized negative log likelihood for a given Wilson coefficient pair.

1D-likelihood scan for $c_{tu}^{(1)}$. The observed and expected limits at 68\% and 95\% confidence level (CL) can be read off by finding the intersection of the likelihood scan curve with the 1$\sigma$ and 2$\sigma$ lines. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$. In the table the $2\Delta$NLL values are shown which correspond to the minimized negative log likelihood for a given Wilson coefficient value.

1D-likelihood scan for $c_{Qu}^{(1)}$. The observed and expected limits at 68\% and 95\% confidence level (CL) can be read off by finding the intersection of the likelihood scan curve with the 1$\sigma$ and 2$\sigma$ lines. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$. In the table the $2\Delta$NLL values are shown which correspond to the minimized negative log likelihood for a given Wilson coefficient value.

1D-likelihood scan for $c_{Qu}^{(8)}$. The observed and expected limits at 68\% and 95\% confidence level (CL) can be read off by finding the intersection of the likelihood scan curve with the 1$\sigma$ and 2$\sigma$ lines. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$. In the table the $2\Delta$NLL values are shown which correspond to the minimized negative log likelihood for a given Wilson coefficient value.