Postfit pulls, constraints, and impacts (both nominal and 'global') for all nuisance parameters in the W-like Z boson mass fit (flip odd/even events to choose muon charge), sorted by the absolute value of the nominal impact.

Postfit pulls, constraints, and impacts (both nominal and 'global') for all nuisance parameters in the W-like Z boson mass fit (charge difference), sorted by the absolute value of the nominal impact.

Postfit pulls, constraints, and impacts (only nominal) for all nuisance parameters in the Z boson mass fit (using dimuon mass), sorted by the absolute value of the nominal impact.

Postfit pulls, constraints, and impacts (only nominal) for all nuisance parameters in the W boson mass fit (helicity cross section fit), sorted by the absolute value of the nominal impact.

$d\sigma^{W}/d\mathit{p}_{T}\ (pb\,/\,GeV)$: Generator-level $\mathit{p}_{T}^{W}$ distribution and uncertainty. The last column reports the result of a simultaneous fit to the single muon $(\mathit{p}_{T}^{μ}, \mathit{\eta}^{μ}, \mathit{q}^{μ})$ distribution in $W$ boson decays and the dimuon $(\mathit{p}_{T}^{μμ},\mathit{y}^{μμ})$ distribution in $Z$ boson events. The measured cross section in each bin is reported after dividing the value by the corresponding bin width.

$d\sigma^{Z}/d|\mathit{y}|\ (pb)$: Unfolded measured $\mathit{y}^{Z}$ distribution compared with the generator-level predictions before the likelihood fit or after the W-like $\mathit{m}_{Z}$ fit or from the direct fit to the $(\mathit{p}_{T}^{μμ},\mathit{y}^{μμ})$ distribution. The measured cross section in each bin is reported after dividing the value by the corresponding bin width.

$d\sigma^{Z}/d\mathit{p}_{T}\ (pb\,/\,GeV)$: Unfolded measured $\mathit{p}_{T}^{Z}$ distribution compared with the generator-level predictions before the likelihood fit or after the W-like $\mathit{m}_{Z}$ fit or from the direct fit to the $(\mathit{p}_{T}^{μμ},\mathit{y}^{μμ})$ distribution. The measured cross section in each bin is reported after dividing the value by the corresponding bin width.

$d\sigma^{W}/d\mathit{p}_{T}\ (pb\,/\,GeV)$: Generator-level $\mathit{p}_{T}^{W}$ distribution and uncertainty, before and after running the helicity fit. The measured cross section in each bin is reported after dividing the value by the corresponding bin width.

$d\sigma^{W}/d|\mathit{y}|\ (pb)$: Generator-level $\mathit{y}^{Z}$ distribution and uncertainty, before and after running the helicity fit. The measured cross section in each bin is reported after dividing the value by the corresponding bin width.

Comparison of the nominal $\mathit{m}_{W}$ measurement and its uncertainty (total or only from $\mathit{p}_{T}^{W}$ modeling), with the alternative measurements using different approaches to the $\mathit{p}_{T}^{W}$ modeling and its uncertainty. The first column reports the measured mass for each configuration, while the second one shows the difference between the measurement and the main result.

Comparison of the nominal $\mathit{m}_{W}$ measurement and its uncertainty (total or only from the CT18Z PDF set), with the alternative measurements using different PDFs and their uncertainty before the scaling procedure described in the paper text. The measured mass values are those reported in Table A.7 of the paper. The first column reports the measured mass for each configuration, while the second one shows the difference between the measurement and the main result.

Comparison of the nominal $\mathit{m}_{W}$ measurement and its uncertainty (total or only from the CT18Z PDF set), with the alternative measurements using different PDFs and their uncertainty after the scaling procedure described in the paper text. The measured mass values are those reported in Table A.7 of the paper. The scale factors for the PDF uncertainties from each set are reported in Table A.3 of the paper. The first column reports the measured mass for each configuration, while the second one shows the difference between the measurement and the main result.

Values of the test statistic $t_{\kappa_{g,\mathrm{off-shell}}}$ as a function of $\kappa_{g,\mathrm{off-shell}}$ while profiling $\kappa_{V,\mathrm{off-shell}}$ obtained with an Asimov dataset and with data. The values with all nuisance parameters fixed at their best-fit values (stat-only) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

Values of the test statistic $t_{\kappa_{V,\mathrm{off-shell}}}$ as a function of $\kappa_{V,\mathrm{off-shell}}$ while profiling $\kappa_{g,\mathrm{off-shell}}$ obtained with an Asimov dataset and with data. The values with all nuisance parameters fixed at their best-fit values (stat-only) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

Values of the test statistic $t_{\kappa_{H}}$ as a function of $\kappa_{H}=\Gamma_H/\Gamma_H^{\mathrm{SM}}$ obtained with an Asimov dataset and with data. The values with all nuisance parameters fixed at their best-fit values (stat-only) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

Values of the test statistic $t_{R_{gg}}$ as a function of $R_{gg}=\kappa_{g,\mathrm{on-shell}}^2/\kappa_{g,\mathrm{off-shell}}^2$ while profiling $R_{VV}=\kappa_{V,\mathrm{on-shell}}^2/\kappa_{V,\mathrm{off-shell}}^2$ and fixing $\kappa_H=1$ obtained with an Asimov dataset and with data. The values with all nuisance parameters fixed at their best-fit values (stat-only) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

Values of the test statistic $t_{R_{VV}}$ as a function of $R_{VV}=\kappa_{V,\mathrm{on-shell}}^2/\kappa_{V,\mathrm{off-shell}}^2$ while profiling $R_{gg}=\kappa_{g,\mathrm{on-shell}}^2/\kappa_{g,\mathrm{off-shell}}^2$ and fixing $\kappa_H=1$ obtained with an Asimov dataset and with data. The values with all nuisance parameters fixed at their best-fit values (stat-only) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

Values of the test statistic $t_{\mu_{\mathrm{off-shell}}}$ assuming a single parameter of interest $\mu_{\mathrm{off-shell}}$ obtained with an Asimov dataset and with data when combining the $H^*\rightarrow ZZ\rightarrow 4\ell$ and $H^*\rightarrow ZZ\rightarrow 2\ell 2\nu$ decay channels. The values from the histogram-based analysis (Phys. Lett. B 846 (2023) 138223) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

Values of the test statistic $t_{\kappa_{H}}$ as a function of $\kappa_{H}=\Gamma_H/\Gamma_H^{\mathrm{SM}}$ obtained with an Asimov dataset and with data in the $H^*\rightarrow ZZ\rightarrow 4\ell$ decay channel. The values from the histogram-based analysis (Phys. Lett. B 846 (2023) 138223) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

Values of the test statistic $t_{\kappa_{H}}$ as a function of $\kappa_{H}=\Gamma_H/\Gamma_H^{\mathrm{SM}}$ obtained with an Asimov dataset and with data when combining the $H^*\rightarrow ZZ\rightarrow 4\ell$ and $H^*\rightarrow ZZ\rightarrow 2\ell 2\nu$ decay channels. The values from the histogram-based analysis (Phys. Lett. B 846 (2023) 138223) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

Values of the test statistic $t_{R_{gg}}$ as a function of $R_{gg}$ while profiling $R_{VV}$ and fixing $\kappa_H=1$ obtained with an Asimov dataset and with data. The values from the histogram-based analysis (Phys. Lett. B 846 (2023) 138223) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

Values of the test statistic $t_{R_{VV}}$ as a function of $R_{VV}$ while profiling $R_{gg}$ and fixing $\kappa_H=1$ obtained with an Asimov dataset and with data. The values from the histogram-based analysis (Phys. Lett. B 846 (2023) 138223) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

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.

Invariant mass of the lepton pair plus jet system, for data, background (pre- and post-reweighting) and three $H\rightarrow Za$ signal hypotheses (for $a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection but not the classification NN requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Regression NN output variable, for data, reweighted background, and three $H\rightarrow Za$ signal hypotheses ($a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection, including a $120<m_{\ell\ell j}<140$ GeV requirement, but not the classification NN output variable requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Classification NN output variable, for data, reweighted background, and three $H\rightarrow Za$ signal hypotheses ($a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection, including a $120<m_{\ell\ell j}<140$ GeV requirement, but not the classification NN output variable requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Invariant mass of the lepton pair plus jet system for data, reweighted background, and various signal hypotheses. Events are required to pass the complete event selection (preselection and classification NN requirement). The pre-fit background distribution is shown along with three $H\rightarrow Za$ signal hypotheses coming from the pure MC sample ($a\rightarrow q\bar{q}/gg$ inclusively). The background normalization is set equal to that of the data for events passing the preselection and being in the region 100-180 GeV. The error bars represent the data sample's statistical uncertainty and the grey bands represent the full background model uncertainty before the fit (assuming all the uncertainties are uncorrelated). The lower panel shows the ratio of the data to the pre-fit background prediction.

Invariant mass of the lepton pair plus jet system for data, reweighted background, and 2.0 GeV signal hypotheses. Events are required to pass the complete event selection (preselection and classification NN requirement). The background shape and normalization come from the fit to the data. The signal includes only $gg$ decays, its shape comes from the fit, and it is normalized to $\mathcal{B}(H\rightarrow Za)=100\%$. The error bars represent the data sample's statistical uncertainty and the grey bands represent the full background model uncertainty after the fit. The lower panel shows the per-bin signal significance, $(N_\text{data}-N_\text{MC})\,\text{\Large{/}}\sqrt{\sigma_\text{data}^2+\sigma_\text{MC}^2}$. The post-fit background uncertainty is calculated from a fit under the background-only hypothesis. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Observed 95$\%$ CL upper limits on $\mathcal{B}(H \to Za)=100\%$ for $\mathcal{B}(a \to gg)=100\%$, and the expectation under the background-only hypothesis together with its $\pm1\sigma$ and $\pm2\sigma$ intervals. The weaker limits from the previous version of the analysis (PRL 125(2020) 221802) are also shown. Linear interpolation is used to set limits between the mass points for which MC signal samples were generated.

Observed 95$\%$ confidence level upper limits on $\mathcal{B}(H \to Za)=100\%$ for $\mathcal{B}(a\rightarrow q\bar{q})=100\%$, and the expectation under the background-only hypothesis together with its $\pm1\sigma$ and $\pm2\sigma$ intervals. The weaker limits from the previous version of the analysis (PRL 125(2020) 221802) are also shown. Linear interpolation is used to set limits between the mass points for which MC signal samples were generated.

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