The observed and expected 95% CL upper limits on the minimal coupling vector LQ pair production cross-sections assuming B = 1 as a function of m$_\text{LQ}$.

The observed and expected 95% CL upper limits on the Yang-Mills coupling vector LQ pair production cross-sections assuming B = 1 as a function of m$_\text{LQ}$.

The observed and expected 95% CL upper limits on the branching ratio into charged leptons as a function of m$_\text{LQ}$ for the scalar LQ.

The observed and expected 95% CL upper limits on the branching ratio into charged leptons as a function of m$_\text{LQ}$ for the vector LQ minimal coupling scenario.

The observed and expected 95% CL upper limits on the branching ratio into charged leptons as a function of m$_\text{LQ}$ for the vector LQ Yang-Mills coupling scenario.

Summary of the observed and expected 95% CL upper limits on the cross section for $\mathrm{LQ}_{3}^{\mathrm{d}}$ pair production as a function of $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$ under the assumptions of B=1.

Summary of the observed and expected 95% CL upper limits on B as a function of $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$.

Summary of the observed and expected 95% CL upper limits on B as a function of $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$.

Cutflow of the preselection requirements (see Section 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the preselection requirements (see Section 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in 1$\ell+\geq 1\tau$ channel (see Table 3) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in 1$\ell+\geq 1\tau$ channel (see Table 3) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the 2$\ell$OS+$\geq 1\tau$ channel (see Table 4) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the 2$\ell$OS+$\geq 1\tau$ channel (see Table 4) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the $2\ell$SS/$3\ell+\geq 1\tau$ channel (see Table 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. In this channel, two signal regions (SR-L and SR-H) are defined based on $p_{\mathrm{T}, 1}^{\tau}$, with SR-L and SR-H requiring $125< p_{\mathrm{T}, 1}^{\tau} < 225$ GeV and $p_{\mathrm{T}, 1}^{\tau}>225$ GeV, respectively. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the $2\ell$SS/$3\ell+\geq 1\tau$ channel (see Table 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. In this channel, two signal regions (SR-L and SR-H) are defined based on $p_{\mathrm{T}, 1}^{\tau}$, with SR-L and SR-H requiring $125< p_{\mathrm{T}, 1}^{\tau} < 225$ GeV and $p_{\mathrm{T}, 1}^{\tau}>225$ GeV, respectively. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Expected 95% CL lower limit on the VLQ $T$ mass as a function of the decay branching ratios into $W b$ and $Ht$.

Observed 95% CL lower limit on the VLQ $T$ mass as a function of the decay branching ratios into $W b$ and $Ht$.

Comparison of the distribution of the scalar sum of small-$R$ jet and lepton transverse momenta, $S_T$, between data and the background prediction in the signal region of the pair-production (PP) $\geq 3\ell$ channel. The background prediction is shown post-fit, i.e. after the fit to the data $S_T$ distributions under the background-only hypothesis. The last bin contains the overflow. An example distribution for a $B\bar B$ signal in the singlet model with $m_B$ = 900 GeV is overlaid.

Comparison of the distribution of the invariant mass of the $Z$ boson candidate and the highest-$p_T$ top-tagged large-$R$ jet, $m_{Zt}$, between data and the background prediction in the signal region of the single-production (SP) 2$\ell$ channel. The background prediction is shown post-fit, i.e. after the fit to the data $m_{Zt}$ distributions under the background-only hypothesis. The last bin contains the overflow. An example distribution for a single-$T$-quark signal with $m_T$ = 900 GeV and $\kappa_T = 0.5$ is overlaid.

Comparison of the distribution of the scalar sum of small-$R$ jet and lepton transverse momenta, $S_T$, between data and the background prediction in the signal region of the single-production (SP) $\geq 3\ell$ channel. The background prediction is shown post-fit, i.e. after the fit to the data $S_T$ distributions under the background-only hypothesis. The last bin contains the overflow. An example distribution for a single-$T$-quark signal with $m_T$ = 900 GeV and $\kappa_T$ = 0.5 is overlaid.

Upper limits at 95$\%$ CL on the cross section of vector-like quark pair production (PP) for $T\bar T$ in the singlet model. The expected limits are shown for the combination of the channels, as are the observed limits.

Upper limits at 95$\%$ CL on the cross section of vector-like quark pair production (PP) for $B\bar B$ in the singlet model. The expected limits are shown for the combination of the channels, as are the observed limits.

Upper limits at 95$\%$ CL on the cross section of vector-like quark pair production (PP) for $T\bar T$ in the doublet model. The expected limits are shown for the combination of the channels, as are the observed limits.

Upper limits at 95$\%$ CL on the cross section of vector-like quark pair production (PP) for $B\bar B$ in the doublet model. The expected limits are shown for the combination of the channels, as are the observed limits.

Upper limits at 95$\%$ CL on the cross section of vector-like quark pair production (PP) for $T\bar T$ with a BR of 100$\%$ to $Zt$. The expected limits are shown for the combination of the channels, as are the observed limits.

Upper limits at 95$\%$ CL on the cross section of vector-like quark pair production (PP) for $B\bar B$ with a BR of 100$\%$ to $Zb$. The expected limits are shown for the combination of the channels, as are the observed limits.

Expected 95$\%$ CL lower limits from the combination of the pair-production channels on the mass of vector-like quarks for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$, adding up to unity.

Observed 95$\%$ CL lower limits from the combination of the pair-production channels on the mass of vector-like quarks for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$, adding up to unity.

Expected 95$\%$ CL lower limits from the combination of the pair-production channels on the mass of vector-like quarks for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$, adding up to unity.

Observed 95$\%$ CL lower limits from the combination of the pair-production channels on the mass of vector-like quarks for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$, adding up to unity.

Upper limits at 95$\%$ CL on the cross section times BR to $Zt$ of single production (SP) of a $T$-quark. The expected limits are shown for the combination of the channels, as are the observed limits. The expected cross section times BR to $Zt$ for single-$T$-quark production is also shown for a coupling $\kappa_T = 0.5$, which corresponds to a coupling of $c_W = \sqrt{c^2_{W,L} + c^2_{W,R}} = 0.45$ from Ref. [14]. The BR assumed here corresponds to the singlet benchmark model, i.e. $\approx 25\%$.

Expected and observed 95$\%$ CL limits from the combination of the single-production channels on the coupling of the $T$ quark to SM particles, $c_W = \sqrt c^2_L + c^2_R$, assuming the singlet-model BR of $\approx 25\%$, as a function of the mass of the $T$ quark, $m_T$. Values of $c_W$ larger than the observed limit are excluded.

Expected and observed 95$\%$ CL limits from the combination of the single-production channels on the mixing angle in the singlet model between the $T$ quark and the top quark, $|\sin\theta_L|$, as a function of the mass of the $T$ quark, $m_T$. Values of $|\sin\theta_L|$ enclosed by the observed limit are excluded, i.e. for $m_T$ larger than $\approx 1200$ GeV, no value of $|\sin\theta_L|$ is excluded.

Expected lower limit from the combination of the single-production channels on the mass of the $T$ quark as a function of the couplings of the $T$ quark to the $W$ boson, $\sqrt c^2_L + c^2_R$, and to the $Z$ boson, $\sqrt c^2_L + c^2_R$ with the assumption of equal BRs for $T\rightarrow Zt$ and $T\rightarrow Ht$ in the limit of large $T$-quark masses.

Observed lower limit from the combination of the single-production channels on the mass of the $T$ quark as a function of the couplings of the $T$ quark to the $W$ boson, $\sqrt c^2_L + c^2_R$, and to the $Z$ boson, $\sqrt c^2_L + c^2_R$ with the assumption of equal BRs for $T\rightarrow Zt$ and $T\rightarrow Ht$ in the limit of large $T$-quark masses.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 500 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 600 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 700 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 750 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 800 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 850 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 900 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 950 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 1000 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 1050 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 1100 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 1150 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 1200 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 1300 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 1400 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 500 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 600 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 700 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 750 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 800 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 850 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 900 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 950 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 1000 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 1050 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 1100 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 1150 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 1200 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 1300 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 1400 GeV.

Observed upper limits at 95$\%$ CL on the cross section times BR to $Zt$ as a function of the vector-like $T$-quark mass and coupling $\kappa_T$ for the combination of the SP $2\ell$ and SP $\geq 3\ell$ channels.

Signal cutflows for the $B\bar B$ process in the singlet model in the PP $2\ell$ $0-1$J channel in the 0-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the singlet model in the PP $2\ell$ $0-1$J channel in the 1-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the doublet model in the PP $2\ell$ $0-1$J channel in the 0-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the doublet model in the PP $2\ell$ $0-1$J channel in the 1-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the case of 100$\%$ BR for $B \rightarrow Zb$ in the PP $2\ell$ $0-1$J channel in the 0-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the case of 100$\%$ BR for $B \rightarrow Zb$ in the PP $2\ell$ $0-1$J channel in the 1-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the singlet model in the PP $2\ell$ $0-1$J channel in the 0-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the singlet model in the PP $2\ell$ $0-1$J channel in the 1-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the doublet model in the PP $2\ell$ $0-1$J channel in the 0-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the doublet model in the PP $2\ell$ $0-1$J channel in the 1-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the case of 100$\%$ BR for $T \rightarrow Zt$ in the PP $2\ell$ $0-1$J channel in the 0-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the case of 100$\%$ BR for $T \rightarrow Zt$ in the PP $2\ell$ $0-1$J channel in the 1-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the singlet model in the PP $2\ell$ $\geq 2$J channel. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the singlet model in the PP $2\ell$ $\geq 2$J channel. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the doublet model in the PP $2\ell$ $\geq 2$J channel. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the doublet model in the PP $2\ell$ $\geq 2$J channel. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the case of 100$\%$ BR for $T \rightarrow Zt$ in the PP $2\ell$ $\geq 2$J channel. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the case of 100$\%$ BR for $B \rightarrow Zb$ in the PP $2\ell$ $\geq 2$J channel. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the singlet model in the PP $\geq 3\ell$ channel. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the singlet model in the PP $\geq 3\ell$ channel. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the doublet model in the PP $\geq 3\ell$ channel. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the doublet model in the PP $\geq 3\ell$ channel. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the case of 100$\%$ BR for $T \rightarrow Zt$ in the PP $\geq 3\ell$ channel. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the case of 100$\%$ BR for $B \rightarrow Zb$ in the PP $\geq 3\ell$ channel. Only statistical uncertainties are shown.

Signal cutflows for single-$T$-quark production with $\kappa_T = 0.5$ in the SP $2\ell$ channel. Only statistical uncertainties are shown.

Signal cutflows for single-$T$-quark production with $\kappa_T = 0.5$ in the SP $\geq 3\ell$ channel. Only statistical uncertainties are shown.

Signal efficiencies in $\%$ in the PP $2\ell$ $0-1$J channel in the 0-large-$R$ jet SR. Uncertainties are statistical only.

Signal efficiencies in $\%$ in the PP $2\ell$ $0-1$J channel in the 1-large-$R$ jet SR. Uncertainties are statistical only.

Signal efficiencies in $\%$ in the PP $2\ell$ $\geq 2$J channel. Uncertainties are statistical only.

Signal efficiencies in $\%$ in the PP $\geq 3\ell$ channel. Uncertainties are statistical only.

Signal efficiencies in $\%$ in the SP $2\ell$ channel for $\kappa_T = 0.5$. Uncertainties are statistical only.

Signal efficiencies in $\%$ in the SP $\geq 3\ell$ channel for $\kappa_T = 0.5$. Uncertainties are statistical only.

Distribution of the diphoton invariant mass in validation region VR2. The solid histograms are stacked to show the SM expectations after the 2&times;2D background estimation technique is applied. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. Statistical and systematic uncertainties are indicated by the shaded area. The lower panel of each plot shows the ratio of the data to the SM prediction for the respective bin. The first and last bins include the underflows and overflows respectively.

Distribution of the missing transverse momentum in validation region VR2. The solid histograms are stacked to show the SM expectations after the 2&times;2D background estimation technique is applied. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. Statistical and systematic uncertainties are indicated by the shaded area. The lower panel of each plot shows the ratio of the data to the SM prediction for the respective bin. The first and last bins include the underflows and overflows respectively.

Distribution of the diphoton invariant mass with all selections of the signal regions applied, except on m_{&gamma;&gamma;} itself, for signal region SR1h. The background estimation techniques described in the text are applied. The different backgrounds are stacked to add up to the total SM prediction in each bin. The predicted yields for signal benchmark models of varying &chi;&#771;⁰₁ mass are also overlaid (not stacked) assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; ) to equal 100%. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The sizes of the statistical and systematic uncertainties are indicated by the shaded areas. The lower panels show the ratio of the data to the SM prediction. Arrows indicate the borders of the signal region (|m_{&gamma;&gamma;}-125 GeV|&lt;5 GeV). The first and last bins include the underflows and overflows respectively.

Distribution of the diphoton invariant mass with all selections of the signal regions applied, except on m_{&gamma;&gamma;} itself, for signal region SR1Z. The background estimation techniques described in the text are applied. The different backgrounds are stacked to add up to the total SM prediction in each bin. The predicted yields for signal benchmark models of varying &chi;&#771;⁰₁ mass are also overlaid (not stacked) assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; ) to equal 50%. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The sizes of the statistical and systematic uncertainties are indicated by the shaded areas. The lower panels show the ratio of the data to the SM prediction. Arrows indicate the borders of the signal region (|m_{&gamma;&gamma;}-125 GeV|&lt;5 GeV). The first and last bins include the underflows and overflows respectively.

Distribution of the diphoton invariant mass with all selections of the signal regions applied, except on m_{&gamma;&gamma;} itself, for signal region SR2. The background estimation techniques described in the text are applied. The different backgrounds are stacked to add up to the total SM prediction in each bin. The predicted yields for signal benchmark models of varying &chi;&#771;⁰₁ mass are also overlaid (not stacked) assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; ) to equal 100%. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The sizes of the statistical and systematic uncertainties are indicated by the shaded areas. The lower panels show the ratio of the data to the SM prediction. Arrows indicate the borders of the signal region (|m_{&gamma;&gamma;}-125 GeV|&lt;5 GeV). The first and last bins include the underflows and overflows respectively.

Observed and expected limits on the pure higgsino cross-section at 95% CL assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; )=100% for different &chi;&#771;⁰₁ masses, obtained by a statistical combination of the three signal regions SR1h, SR1Z and SR2. The inner and outer bands indicate the 1&sigma; and 2&sigma; variation on the expected limit respectively.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(&chi;&#771;⁰₁ &rarr; hG&#771; ) as a function of the higgsino mass m(&chi;&#771;⁰₁) assuming it decays via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The &plusmn; 1&sigma; variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by &plusmn;1 &sigma;. Values of B(&chi;&#771;⁰₁ &rarr; hG&#771; ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(&chi;&#771;⁰₁ &rarr; hG&#771; ) as a function of the higgsino mass m(&chi;&#771;⁰₁) assuming it decays via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The &plusmn; 1&sigma; variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by &plusmn;1 &sigma;. Values of B(&chi;&#771;⁰₁ &rarr; hG&#771; ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(&chi;&#771;⁰₁ &rarr; hG&#771; ) as a function of the higgsino mass m(&chi;&#771;⁰₁) assuming it decays via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The &plusmn; 1&sigma; variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by &plusmn;1 &sigma;. Values of B(&chi;&#771;⁰₁ &rarr; hG&#771; ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(&chi;&#771;⁰₁ &rarr; hG&#771; ) as a function of the higgsino mass m(&chi;&#771;⁰₁) assuming it decays via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The &plusmn; 1&sigma; variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by &plusmn;1 &sigma;. Values of B(&chi;&#771;⁰₁ &rarr; hG&#771; ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(&chi;&#771;⁰₁ &rarr; hG&#771; ) as a function of the higgsino mass m(&chi;&#771;⁰₁) assuming it decays via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The &plusmn; 1&sigma; variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by &plusmn;1 &sigma;. Values of B(&chi;&#771;⁰₁ &rarr; hG&#771; ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(&chi;&#771;⁰₁ &rarr; hG&#771; ) as a function of the higgsino mass m(&chi;&#771;⁰₁) assuming it decays via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The &plusmn; 1&sigma; variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by &plusmn;1 &sigma;. Values of B(&chi;&#771;⁰₁ &rarr; hG&#771; ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Numbers of signal and background events in the signal regions. The respective background estimation techniques are applied. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The different backgrounds are stacked to add up to the total Standard Model prediction in each bin. The predicted yields for signal benchmark models of varying &chi;&#771;⁰₁ mass are also plotted (not stacked), assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; )=100% and a &chi;&#771;⁰₁ mass of 130 or 200 GeV. The statistical and systematic uncertainties are indicated by the shaded areas in the top plot. The bottom panel shows the statistical significance&nbsp;<a href="https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2020-025/">[Ref]</a> of the difference between the SM prediction and the observed data in each region.

Observed 95% CL limits in pb on the pure higgsino cross-section, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Limits are obtained by a statistical combination of the three signal regions SR1h, SR1Z and SR2, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Acceptances for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Acceptances are provided separately for signal region SR1h, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Acceptances for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Acceptances are provided separately for signal region SR1Z, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Acceptances for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Acceptances are provided separately for signal region SR2, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Efficiencies for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Efficiencies are provided separately for signal region SR1h, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Efficiencies for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Efficiencies are provided separately for signal region SR1Z, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Efficiencies for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Efficiencies are provided separately for signal region SR2, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Observed and expected numbers of events in the three signal regions. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The table also includes model-independent 95% CL upper limits on the visible number of BSM events (S⁹⁵_obs), the number of BSM events given the expected number of background events (S⁹⁵_exp) and the visible BSM cross-section (&lang;&epsilon; &sigma;&rang;_obs⁹⁵), all calculated from pseudo-experiments. The discovery p-value (p₀) is also shown and its value is capped at 0.5 if the observed number of events is below the expected number of events.

Cutflows of two benchmark signal points assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; )=100% for all three discovery signal regions. The initial selection includes the leptons veto. Only statistical uncertainties are included. Expected yields are normalised to a luminosity of 139 fb^-1.

Observed 95% CL limits on the production cross-section for scalar up-type LQs decaying into a top quark and a neutrino or a bottom quark and an electron.

Expected 95% CL limits on the production cross-section for scalar up-type LQs decaying into a top quark and a neutrino or a bottom quark and an electron.

Observed 95% CL limits on the production cross-section for scalar down-type LQs decaying into a bottom quark and a neutrino or a top quark and a muon.

Expected 95% CL limits on the production cross-section for scalar down-type LQs decaying into a bottom quark and a neutrino or a top quark and a muon.

Observed 95% CL limits on the production cross-section for scalar down-type LQs decaying into a bottom quark and a neutrino or a top quark and an electron.

Expected 95% CL limits on the production cross-section for scalar down-type LQs decaying into a bottom quark and a neutrino or a top quark and an electron.

Observed 95% CL limits on the production cross-section for vector up-type LQs in the Yang-Mills coupling scenario decaying into a top quark and a neutrino or a bottom quark and a muon.

Expected 95% CL limits on the production cross-section for vector up-type LQs in the Yang-Mills coupling scenario decaying into a top quark and a neutrino or a bottom quark and a muon.

Observed 95% CL limits on the production cross-section for vector up-type LQs in the Yang-Mills coupling scenario decaying into a top quark and a neutrino or a bottom quark and an electron.

Expected 95% CL limits on the production cross-section for vector up-type LQs in the Yang-Mills coupling scenario decaying into a top quark and a neutrino or a bottom quark and an electron.

Observed 95% CL limits on the production cross-section for vector up-type LQs in the minimal coupling scenario decaying into a top quark and a neutrino or a bottom quark and a muon.

Expected 95% CL limits on the production cross-section for vector up-type LQs in the minimal coupling scenario decaying into a top quark and a neutrino or a bottom quark and a muon.

Observed 95% CL limits on the production cross-section for vector up-type LQs in the minimal coupling scenario decaying into a top quark and a neutrino or a bottom quark and an electron.

Expected 95% CL limits on the production cross-section for vector up-type LQs in the minimal coupling scenario decaying into a top quark and a neutrino or a bottom quark and an electron.

Product of signal acceptance and efficiency in the training region for scalar up-type LQs decaying into top quarks and neutrinos or bottom quarks and muons depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(LQ_{mix}^u\rightarrow b\mu)$.

Product of signal acceptance and efficiency in the training region for scalar up-type LQs decaying into top quarks and neutrinos or bottom quarks and electrons depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(LQ_{mix}^u\rightarrow be)$.

Product of signal acceptance and efficiency in the training region for scalar down-type LQs decaying into top quarks and muons or bottom quarks and neutrinos depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(LQ_{mix}^d\rightarrow t\mu)$.

Product of signal acceptance and efficiency in the training region for scalar down-type LQs decaying into top quarks and electrons or bottom quarks and neutrinos depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(LQ_{mix}^d\rightarrow te)$.

Product of signal acceptance and efficiency in the training region for vector LQs in the Yang-Mills coupling scenario decaying into top quarks and neutrinos or bottom quarks and muons depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(vLQ_{mix}^{YM}\rightarrow b\mu)$.

Product of signal acceptance and efficiency in the training region for vector LQs in the Yang-Mills coupling scenario decaying into top quarks and neutrinos or bottom quarks and electrons depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(vLQ_{mix}^{YM}\rightarrow be)$.

Product of signal acceptance and efficiency in the training region for vector LQs in the minimal coupling scenario decaying into top quarks and neutrinos or bottom quarks and muons depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(vLQ_{mix}^{min}\rightarrow b\mu)$.

Product of signal acceptance and efficiency in the training region for vector LQs in the minimal coupling scenario decaying into top quarks and neutrinos or bottom quarks and electrons depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(vLQ_{mix}^{min}\rightarrow be)$.

Cut-flow for the four models containing scalar LQs at a mass of $m_{LQ}=1300$ GeV and $\mathcal{B}(LQ_{mix}\rightarrow q_3\ell)=0.5$.

Cut-flow for the four models containing vector LQs at a mass of $m_{LQ}=1300$ GeV and $\mathcal{B}(LQ_{mix}\rightarrow q_3\ell)=0.5$.

The measured leading jet $p_T$ distribution in the $W \rightarrow \mu \nu$ control region, for the M1 selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit.

The measured $E_T^{miss}$ distribution in the $W \rightarrow e \nu$ control region, for the M1 selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit.

The measured leading jet $p_T$ distribution in the $W \rightarrow e \nu$ control region, for the M1 selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit.

The measured $E_T^{miss}$ distribution in the $Z \rightarrow \mu \mu$ control region, for the M1 selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit.

The measured leading jet $p_T$ distribution in the $Z \rightarrow \mu \mu$ control region, for the M1 selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit.

The measured $E_T^{miss}$ distribution in the $W\rightarrow \mu \nu$ control region, for the c-tagged selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit.

The measured leading jet $p_T$ distribution in the $W\rightarrow \mu \nu$ control region, for the c-tagged selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit.

The measured $E_T^{miss}$ distribution in the $W\rightarrow e \nu$ control region, for the c-tagged selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit.

The measured leading jet $p_T$ distribution in the $W\rightarrow e \nu$ control region, for the c-tagged selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit.

The measured $E_T^{miss}$ distribution in the $Z\rightarrow ll$ control region, for the c-tagged selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit.

The measured leading jet $p_T$ distribution in the $Z\rightarrow ll$ control region, for the c-tagged selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit.

The measured $E_T^{miss}$ distribution in the t-(anti)-t control region, for the c-tagged selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit.

The measured leading jet $p_T$ distribution in the t-(anti)-t control region, for the c-tagged selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit.

Measured leading jet $p_T$ distribution for the V3 selections compared to the SM predictions.

Measured $E_T^{miss}$ distribution for the V4 selections compared to the SM predictions.

Measured leading jet $p_T$ distribution for the V5 selections compared to the SM predictions.

Measured $E_T^{miss}$ distribution for the V5 selections compared to the SM predictions.

Measured $E_T^{miss}$ distribution for the M1 selection compared to the SM predictions. For illustration purposes, the distribution of two different SUSY scenarios are included.

Measured leading jet $p_T$ distribution for the M1 selection compared to the SM predictions. For illustration purposes, the distribution of two different SUSY scenarios are included.

Measured $E_T^{miss}$ distribution for the M2 selection compared to the SM predictions. For illustration purposes, the distribution of two different SUSY scenarios are included.

Measured leading jet $p_T$ distribution for the M2 selection compared to the SM predictions. For illustration purposes, the distribution of two different SUSY scenarios are included.

Measured $E_T^{miss}$ distribution for the M3 selection compared to the SM predictions. For illustration purposes, the distribution of two different SUSY scenarios are included.

Measured leading jet $p_T$ distribution for the M3 selection compared to the SM predictions. For illustration purposes, the distribution of two different SUSY scenarios are included.

Measured $E_T^{miss}$ distribution for the C1 selection before the cut in the variable shown is applied. The data are compared to the SM predictions. For illustration purposes, the distribution of two different SUSY scenarios are included.

Measured leading jet $p_T$ distribution for the C1 selection before the cut in the variable shown is applied. The data are compared to the SM predictions. For illustration purposes, the distribution of two different SUSY scenarios are included.

Measured leading jet $p_T$ for the C2 selection. The data are compared to the SM predictions. For illustration purposes, the distribution of two different SUSY scenarios are included.

Measured jet multiplicity for the C2 selection. The data are compared to the SM predictions. For illustration purposes, the distribution of two different SUSY scenarios are included.

The measured $W$ transverse mass distribution in the $W \rightarrow \mu \nu$ control region for the M1 monojet-like selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit.

The measured $W$ transverse mass distribution in the $W \rightarrow e \nu$ control region for the M1 monojet-like selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit.

The measured dimuon invariant mass distribution in the $Z \rightarrow \mu \mu$ control region for the M1 monojet-like selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit.

The measured $E_T^{miss}$ distribution in a multijet control sample for the monojet-like analysis. This region is defined with similar $E_T^{miss}$ and leading jet $p_T$ cuts as M1, but releasing the jet-veto and inverting the $\Delta \phi (jet, p_T^{miss})$ cut, now requiring it to be less than 0.4. The data are compared to predictions.

Measured jet multiplicity in the $W \rightarrow e \nu$ control region for the c-tagged analysis, compared to background predictions. The latter include the global normalization factors extracted from the fit.

Measured jet multiplicity in the $W \rightarrow \mu \nu$ control region for the c-tagged analysis, compared to background predictions. The latter include the global normalization factors extracted from the fit.

Measured jet multiplicity in the $Z \rightarrow ll$ control region for the c-tagged analysis, compared to background predictions. The latter include the global normalization factors extracted from the fit.

Measured jet multiplicity in the t-(anti)-t control region for the c-tagged analysis, compared to background predictions. The latter include the global normalization factors extracted from the fit.

Measured dilepton invariant mass in the $Z \rightarrow ll$ control region for the c-tagged analysis. The predictions include the global normalization factors extracted from the fit.

Measured jet multiplicity distribution for the M1 selection compared to the SM predictions. For illustration purposes, the impact of two different SUSY scenarios are included.

Measured leading jet $\eta$ distribution for the M1 selection compared to the SM predictions. For illustration purposes, the impact of two different SUSY scenarios are included.

Measured $\Delta \phi (jet, p_T^{miss})$ distribution for the M1 selection compared to the SM predictions. For illustration purposes, the impact of two different SUSY scenarios are included.

Measured $E_T^{miss}$ distribution after preselection cuts compared to background predictions.

Measured jet multiplicity distribution in the C1 c-tagged selection compared to background predictions.

Multiplicity of jets with loose heavy-flavor tags in the C1 selection.

Multiplicity of jets with medium heavy-flavor tags in the C1 selection.

Expected and observed CLs, number of events, acceptance and efficiency for M1 SR for $\tilde{t}_1 \rightarrow c + \tilde{\chi}^0_1$. The expected event yield in 20.3 fb$^{-1}$, $N_{\text{events}}$, is used to calculate $\sigma\times A\times\varepsilon$ and $A\times\varepsilon$. The acceptance $A$ is defined as the percentage of events passing the selection at truth level, and the efficiency $\varepsilon$ is the ratio of events passing the selection at reconstructed level over the number of events passing the selection at truth level. The modified efficiency $\varepsilon\prime$ is the ratio of events passing the selection at truth and reconstructed level over the number of events passing the selection at truth level.

Expected and observed CLs, number of events, acceptance and efficiency for M2 SR for $\tilde{t}_1 \rightarrow c + \tilde{\chi}^0_1$. The expected event yield in 20.3 fb$^{-1}$, $N_{\text{events}}$, is used to calculate $\sigma\times A\times\varepsilon$ and $A\times\varepsilon$. The acceptance $A$ is defined as the percentage of events passing the selection at truth level, and the efficiency $\varepsilon$ is the ratio of events passing the selection at reconstructed level over the number of events passing the selection at truth level. The modified efficiency $\varepsilon\prime$ is the ratio of events passing the selection at truth and reconstructed level over the number of events passing the selection at truth level.

Expected and observed CLs, number of events, acceptance and efficiency for M3 SR for $\tilde{t}_1 \rightarrow c + \tilde{\chi}^0_1$. The expected event yield in 20.3 fb$^{-1}$, $N_{\text{events}}$, is used to calculate $\sigma\times A\times\varepsilon$ and $A\times\varepsilon$. The acceptance $A$ is defined as the percentage of events passing the selection at truth level, and the efficiency $\varepsilon$ is the ratio of events passing the selection at reconstructed level over the number of events passing the selection at truth level. The modified efficiency $\varepsilon\prime$ is the ratio of events passing the selection at truth and reconstructed level over the number of events passing the selection at truth level.

Expected and observed CLs, number of events, acceptance and efficiency for C1 SR for $\tilde{t}_1 \rightarrow c + \tilde{\chi}^0_1$. The expected event yield in 20.3 fb$^{-1}$, $N_{\text{events}}$, is used to calculate $\sigma\times A\times\varepsilon$ and $A\times\varepsilon$. The acceptance $A$ is defined as the percentage of events passing the selection at truth level, and the efficiency $\varepsilon$ is the ratio of events passing the selection at reconstructed level over the number of events passing the selection at truth level.

Expected and observed CLs, number of events, acceptance and efficiency for C2 SR for $\tilde{t}_1 \rightarrow c + \tilde{\chi}^0_1$. The expected event yield in 20.3 fb$^{-1}$, $N_{\text{events}}$, is used to calculate $\sigma\times A\times\varepsilon$ and $A\times\varepsilon$. The acceptance $A$ is defined as the percentage of events passing the selection at truth level, and the efficiency $\varepsilon$ is the ratio of events passing the selection at reconstructed level over the number of events passing the selection at truth level.

Expected and observed CLs, number of events, acceptance and efficiency for M1 SR for $\tilde{b}_1 \rightarrow b + \tilde{\chi}^0_1$. The expected event yield in 20.3 fb$^{-1}$, $N_{\text{events}}$, is used to calculate $\sigma\times A\times\varepsilon$ and $A\times\varepsilon$. The acceptance $A$ is defined as the percentage of events passing the selection at truth level, and the efficiency $\varepsilon$ is the ratio of events passing the selection at reconstructed level over the number of events passing the selection at truth level. The modified efficiency $\varepsilon\prime$ is the ratio of events passing the selection at truth and reconstructed level over the number of events passing the selection at truth level.

Expected and observed CLs, number of events, acceptance and efficiency for M2 SR for $\tilde{b}_1 \rightarrow b + \tilde{\chi}^0_1$. The expected event yield in 20.3 fb$^{-1}$, $N_{\text{events}}$, is used to calculate $\sigma\times A\times\varepsilon$ and $A\times\varepsilon$. The acceptance $A$ is defined as the percentage of events passing the selection at truth level, and the efficiency $\varepsilon$ is the ratio of events passing the selection at reconstructed level over the number of events passing the selection at truth level. The modified efficiency $\varepsilon\prime$ is the ratio of events passing the selection at truth and reconstructed level over the number of events passing the selection at truth level.

Expected and observed CLs, number of events, acceptance and efficiency for M3 SR for $\tilde{b}_1 \rightarrow b + \tilde{\chi}^0_1$. The expected event yield in 20.3 fb$^{-1}$, $N_{\text{events}}$, is used to calculate $\sigma\times A\times\varepsilon$ and $A\times\varepsilon$. The acceptance $A$ is defined as the percentage of events passing the selection at truth level, and the efficiency $\varepsilon$ is the ratio of events passing the selection at reconstructed level over the number of events passing the selection at truth level. The modified efficiency $\varepsilon\prime$ is the ratio of events passing the selection at truth and reconstructed level over the number of events passing the selection at truth level.

Expected and observed CLs, number of events, acceptance and efficiency for M1 SR for $\tilde{t}_1 \rightarrow b + ff' + \tilde{\chi}^0_1$. The expected event yield in 20.3 fb$^{-1}$, $N_{\text{events}}$, is used to calculate $\sigma\times A\times\varepsilon$ and $A\times\varepsilon$. The acceptance $A$ is defined as the percentage of events passing the selection at truth level, and the efficiency $\varepsilon$ is the ratio of events passing the selection at reconstructed level over the number of events passing the selection at truth level. The modified efficiency $\varepsilon\prime$ is the ratio of events passing the selection at truth and reconstructed level over the number of events passing the selection at truth level.

Expected and observed CLs, number of events, acceptance and efficiency for M2 SR for $\tilde{t}_1 \rightarrow b + ff' + \tilde{\chi}^0_1$. The expected event yield in 20.3 fb$^{-1}$, $N_{\text{events}}$, is used to calculate $\sigma\times A\times\varepsilon$ and $A\times\varepsilon$. The acceptance $A$ is defined as the percentage of events passing the selection at truth level, and the efficiency $\varepsilon$ is the ratio of events passing the selection at reconstructed level over the number of events passing the selection at truth level. The modified efficiency $\varepsilon\prime$ is the ratio of events passing the selection at truth and reconstructed level over the number of events passing the selection at truth level.

Expected and observed CLs, number of events, acceptance and efficiency for M3 SR for $\tilde{t}_1 \rightarrow b + ff' + \tilde{\chi}^0_1$. The expected event yield in 20.3 fb$^{-1}$, $N_{\text{events}}$, is used to calculate $\sigma\times A\times\varepsilon$ and $A\times\varepsilon$. The acceptance $A$ is defined as the percentage of events passing the selection at truth level, and the efficiency $\varepsilon$ is the ratio of events passing the selection at reconstructed level over the number of events passing the selection at truth level. The modified efficiency $\varepsilon\prime$ is the ratio of events passing the selection at truth and reconstructed level over the number of events passing the selection at truth level.

Expected and observed upper limits at 95% CL on the cross section of vector-like quark pair production for $B\bar{B}$ in the singlet model.

Expected and observed upper limits at 95% CL on the cross section of vector-like quark pair production for $T\bar{T}$ in the doublet model.

Expected and observed upper limits at 95% CL on the cross section of vector-like quark pair production, considereing both $T\bar{T}$ and $B\bar{B}$, in the $(T,B)$ doublet model.

Expected 95% CL lower limits on the vector-like top quark mass as a function of the branching ratios into $Wb$ and $Ht$.

Observed 95% CL lower limits on the vector-like top quark mass as a function of the branching ratios into $Wb$ and $Ht$.

Expected 95% CL lower limits on the vector-like bottom quark mass as a function of the branching ratios into $Hb$ and $Wt$.

Observed 95% CL lower limits on the vector-like top quark mass as a function of the branching ratios into $Hb$ and $Wt$.

The expected upper limits on the production cross-section times the product of branching ratios for the benchmark signal scenario involving a scalar particle $X$ with narrow width decaying via $X\rightarrow aa\rightarrow 4\gamma$, $\sigma_X\times B(X\rightarrow aa)\times B(a\rightarrow\gamma\gamma)^2$. The limits for $m_{a}$ = 5 GeV and 10 GeV do not cover as large a range as the other mass points, since the region of interest is limited to $ m_{a} < 0.01 \times m_{X}$. Additionally, the expected limits are not provided for a small number of points, indicated with a hyphen, because of a technical failure with the computation.

The observed upper limits on the production cross-section times the product of branching ratios for the benchmark signal scenario involving a scalar particle $X$ with narrow width decaying via $X\rightarrow aa\rightarrow 6\pi^0$, $\sigma_X\times B(X\rightarrow aa)\times B(a\rightarrow 3\pi^0)^2$. The limits for $m_{a}$ = 5 GeV and 10 GeV do not cover as large a range as the other mass points, since the region of interest is limited to $ m_{a} < 0.01 \times m_{X}$.

The expected upper limits on the production cross-section times the product of branching ratios for the benchmark signal scenario involving a scalar particle $X$ with narrow width decaying via $X\rightarrow aa\rightarrow 6\pi^0$, $\sigma_X\times B(X\rightarrow aa)\times B(a\rightarrow 3\pi^0)^2$. The limits for $m_{a}$ = 5 GeV and 10 GeV do not cover as large a range as the other mass points, since the region of interest is limited to $ m_{a} < 0.01 \times m_{X}$. Additionally, the expected limits are not provided for a small number of points, indicated with a hyphen, because of a technical failure with the computation.

Observed 95% CL upper limits on the visible cross section as a function of $m_X$ and the fraction of events in the low-$\Delta E$ category.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 0.1 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 0.5 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 0.7 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 1 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 2 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 5 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 10 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 0.5 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 0.7 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 1 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 2 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 5 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 10 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 0.1 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 0.5 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 0.7 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 1 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 2 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 5 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 10 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 0.5 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 0.7 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 1 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 2 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 5 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 10 GeV.

Selection efficiency for photons originating from the BSM process $X\rightarrow\gamma\gamma$, where the $X$ particle is a high-mass narrow-width scalar particle originating from the gluon--gluon fusion process.

Fraction of photons with a value of shower shape variable $\Delta E$ lower than the threshold, for photons originating from the BSM process $X\rightarrow\gamma\gamma$, where the $X$ particle is a high-mass narrow-width scalar particle originating from the gluon--gluon fusion process.

Distribution of the resonance mass in the c-tag Signal Region of the $ ce$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the b-tag Signal Region of the $ ce$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the untagged Signal Region of the $ c\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the c-tag Signal Region of the $ c\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the b-tag Signal Region of the $ c\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 0-tag Signal Region of the $ be$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 1-tag Signal Region of the $ be$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 2-tag Signal Region of the $ be$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 0-tag Signal Region of the $ b\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 1-tag Signal Region of the $ b\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 2-tag Signal Region of the $ b\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into electrons, shown as a function of $m_{LQ}$ for the $qe$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into muons, shown as a function of $m_{LQ}$ for the $q\mu$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into electrons, shown as a function of $m_{LQ}$ for the $ce$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into muons, shown as a function of $m_{LQ}$ for the $c\mu$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into electrons, shown as a function of $m_{LQ}$ for the $be$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into muons, shown as a function of $m_{LQ}$ for the $b\mu$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $qe$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $q\mu$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $ce$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $c\mu$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $be$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $b\mu$ channel.

The signal selection efficiency x acceptance summed over all signal regions, for all masses and LQ decay channels considered.

The observed and expected limits for all masses and LQ decay channels considered.

Cutflow Table in the electron channel, considering signal samples with LQ mass of 1 TeV.

Cutflow Table in the muon channel, considering signal samples with LQ mass of 1 TeV.