The distributions of the variables used in the simultaneous fit for the SR. The black points with error bars represent the data and their statistical uncertainties, whereas the shaded band represents the predicted uncertainties. The bottom panel in each figure shows the ratio of the number of events observed in data to that of the total SM prediction. The last bin of each plot has been extended to include the overflow contribution.

Expected and observed 95% upper limits on the product of the cross section and branching fraction$\sigma(\mathrm{pp} \ ightarrow W\mathrm{a})\mathcal{B}(W ightarrow \ell^{+} v_{\ell})\mathcal{B}(\mathrm{a} \rightarrow Z\gamma)\mathcal{B}(Z \rightarrow \ell^{+} \ell^{-})$ s a function of the ALP mass from 110 to 400 GeV

Expected and observed 95% upper limits on the photophobic ALP model parameter $1/f_{\mathrm{a}}$ as a function of ALP mass reinterpreted from $1/f_{\mathrm{a}}=2~\mathrm{TeV}^{-1}$

Comparison of the shapes extrapolated with a linear or quadratic extrapolation of the QCD contribution in the W+ signal region.

Comparison of the shapes extrapolated with a linear or quadratic extrapolation of the QCD contribution in the W- signal region.

Post-fit distributions of the transverse mass in the W+ signal region.

Post-fit distributions of the transverse mass in the W- signal region.

Post-fit distributions of the dimuon mass in the Z boson signal region.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W+ signal region in the transverse mass bin 12 < mT < 18 GeV. The values corresponding to the smallest relative muon isolation are extrapolated with a linear (smaller value) and quadratic (larger value) function.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W- signal region in the transverse mass bin 12 < mT < 18 GeV.The values corresponding to the smallest relative muon isolation are extrapolated with a linear (smaller value) and quadratic (larger value) function.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W+ signal region in the transverse mass bin 102 < mT < 108 GeV.The values corresponding to the smallest relative muon isolation are extrapolated with a linear (larger value) and quadratic (smaller value) function.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W- signal region in the transverse mass bin 102 < mT < 108 GeV.The values corresponding to the smallest relative muon isolation are extrapolated with a linear (smaller value) and quadratic (larger value) function.

The table compares theoretical and measured values of the product of the fiducial cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value. The relative uncertainties on the measured values are given separately with and without the contribution from the luminosity uncertainty.

The table compares theoretical and measured values of the product of the total cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value. The relative uncertainties on the measured values are given separately with and without the contribution from the luminosity uncertainty.

The table compares theoretical and measured ratios of the product of the fiducial cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value.

The table compares theoretical and measured ratios of the product of the total cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value.

Comparison of the measured and theoretical cross sections for Z, W+, W-, and W production at different center-of-mass energies. The theoretical predictions are generated with DYTurbo in combination with the NNPDF 3.1 PDF set.

The relative systematic uncertainties of the fitted number of neutrino interactions.

The relative systematic uncertainties of the fitted number of anti-neutrino interactions.

The relative systematic uncertainties of the fitted number of neutrino and anti-neutrino interactions.

Number of observed CC $\nu_{\mu}$ neutrino candidate events in data and simulation.

Number of observed CC $\nu_{\mu}$ neutrino candidate events per inv. bin width.

Number of CC $\nu_{\mu$} neutrino interactions in the fiducial volume.

Number of CC $\nu_{\mu$} neutrino interactions in the fiducial volume per bin width.

Number of CC $\bar{nu}_{\mu$} neutrino interactions in the fiducial volume per bin width.

Number of CC $\nu_{\mu}$ and $\bar{nu}_{\mu$} neutrino interactions in the fiducial volume per bin width.

Neutrino cross section

Anti-neutrino cross section

Neutrino and anti-neutrino cross section

Neutrino cross section divided by energy.

Anti-neutrino cross section divided by energy.

Neutrino and anti-neutrino cross section divided by energy.

Neutrino flux divided by the cross sectional area of the fiducial volume.

Anti-Neutrino flux divided by the cross sectional area of the fiducial volume.

Neutrino and anti-Neutrino flux divided by the cross sectional area of the fiducial volume.

Neutrino flux per bin width divided by the cross sectional area of the fiducial volume.

Anti-neutrino flux per bin width divided by the cross sectional area of the fiducial volume.

Neutrino and anti-neutrino flux per bin width divided by the cross sectional area of the fiducial volume.

Number of neutrino interactions from pion and kaon decays.

Covariance matrix of the fitted number of neutrinos from pion and kaon decays.

Response matrix to unfold the reconstructed, calibrated muon momentum $q/p^{'}_{\mu}$ to the true neutrino energy $-L/E_{\nu}$.

Distribution of the number of clusters in the IFT for neutrino-like data, muon-like data, and simulation.

Distribution of the muon angle for neutrino-like data, muon-like data, and simulation.

Distribution of the muon charge over momentum for neutrino-like data, muon-like data, and simulation.

Efficiency in bins of the neutrino energy $-L/E_{\nu}$.

Covariance matrix from the likelihood of the number of neutrino interactions in bins of the neutrino energy $-L/E_{\nu}$.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Prediction from the HVT model.

Prediction from the HVT model.

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

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

Distributions for data and expected signal and background contributions of the MVA score for the µ + jets channel in the 3j1b category, before the maximum likelihood fit. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty in the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-prediction ratio. The first and last bins in each distribution include underflow and overflow events, respectively.

Observed and predicted number of events in each of the eight categories of the signal region, before the maximum likelihood fit. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty in the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-prediction ratio.

Distributions for the e + jets final state before the maximum likelihood fit, MVA score bins for the 3j1b category and ∆R med ( j, j 0 ) bins for the other categories. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty of the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-simulation ratio, where the simulations are adjusted according to their normalization before the fit. The first and last bins in each distribution include underflow and overflow events, respectively.

Distributions for the e + jets final state after the maximum likelihood fit, MVA score bins for the 3j1b category and ∆R med ( j, j 0 ) bins for the other categories. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty of the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-simulation ratio, where the simulations are adjusted according to their normalization after the fit. The first and last bins in each distribution include underflow and overflow events, respectively.

Distributions for the µ + jets final state before the maximum likelihood fit, MVA score bins for the 3j1b category and ∆R med ( j, j 0 ) bins for the other categories. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty of the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-simulation ratio, where the simulations are adjusted according to their normalization before the fit. The first and last bins in each distribution include underflow and overflow events, respectively.

Distributions for the µ + jets final state after the maximum likelihood fit, MVA score bins for the 3j1b category and ∆R med ( j, j 0 ) bins for the other categories. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty of the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-simulation ratio, where the simulations are adjusted according to their normalization after the fit. The first and last bins in each distribution include underflow and overflow events, respectively.

Covariance matrix for the lepton+jets ML fit.

Covariance matrix for the combination ML fit.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Acceptance across the SUSY signal grid for SR-HM1.

Acceptance across the SUSY signal grid for SR-HM2.

Acceptance across the SUSY signal grid for SR-HM3.

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

Acceptance across the SUSY signal grid for SR-Comp1.

Acceptance across the SUSY signal grid for SR-Comp2.

Acceptance across the SUSY signal grid for SR-Comp3.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Selection efficieny for GMSB supersymmetric $\tau$.

Selection efficiency for gluino R-hardon samples.

Selection efficiency for supersymmetric top R-hadrons.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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