Negative log-likelihood contours corresponding to 68% and 95% CL in the ($\kappa_{V}$, $\kappa_{F}$) plane obtained from a combined fit, assuming no contributions from invisible or undetected non-SM Higgs boson decays.

Reduced coupling strength modifiers $\kappa_{F}\cdot m_{F}/\text{vev}$ for fermions ($F=t,\,b,\,\tau,\,\mu$) and $\sqrt{\kappa_{V}}\cdot m_{V}/\text{vev}$ for vector bosons as a function of their masses $m_{F}$ and $m_{V}$ with vev = 246 GeV. Fit scenario with $\kappa_{c}=\kappa_{t}$ (coloured circle markers) is shown. Loop-induced processes are assumed to have the SM structure, and Higgs boson decays to non-SM particles are not allowed. The $p$-value for compatibility of the combined measurement and the SM prediction is 56%. The lower panel shows the values of the coupling strength modifiers.

Reduced coupling strength modifiers $\kappa_{F}\cdot m_{F}/\text{vev}$ for fermions ($F=t,\,b,\,\tau,\,\mu$) and $\sqrt{\kappa_{V}}\cdot m_{V}/\text{vev}$ for vector bosons as a function of their masses $m_{F}$ and $m_{V}$ with vev = 246 GeV. Fit scenario with $\kappa_{c}$ left free-floating (grey cross markers) is shown. Loop-induced processes are assumed to have the SM structure, and Higgs boson decays to non-SM particles are not allowed. The $p$-value for compatibility of the combined measurement and the SM prediction is 65%. The lower panel shows the values of the coupling strength modifiers. The grey arrow points in the direction of the best-fit value and the grey uncertainty bar extends beyond the lower panel range.

Reduced coupling strength modifiers and their uncertainties per particle type with effective photon, $Z\gamma$ and gluon couplings. The scenario where $B_{inv.}=B_{u.}=0$ is assumed is shown as solid lines. The $p$-value for compatibility with the SM prediction is 61% in this case. The scenario where $B_{inv.}$ and $B_{u.}$ are allowed to contribute to the total Higgs boson decay width while assuming that $\kappa_{V}\le1$ and $B_{u.}\ge0$ is shown as dashed lines.

Reduced coupling strength modifiers and their uncertainties per particle type with effective photon, $Z\gamma$ and gluon couplings. The scenario where $B_{inv.}$ and $B_{u.}$ are allowed to contribute to the total Higgs boson decay width while assuming that $\kappa_{V}\le1$ and $B_{u.}\ge0$ is shown as dashed lines. The lower panel shows the 95% CL upper limits on $B_{inv.}$ and $B_{u.}$.

Observed and predicted Higgs boson production cross sections in different kinematic regions. The $p$-value for compatibility of the combined measurement and the SM prediction is 94%. Kinematic regions are defined separately for each production process, based on the jet multiplicity, the transverse momentum of the Higgs ($p_{T}^{H}$) and vector bosons ($p_{T}^{W}$ and $p_{T}^{Z}$) and the two-jet invariant mass ($m_{jj}$). The VH-enriched and VBF-enriched regions with the respective requirements of $m_{jj}\in[60, 120)$ GeV and $m_{jj}\notin[60, 120)$ GeV are enhanced in signal events from $VH$ and VBF productions, respectively.

Observed variations of −2 ln $\Lambda(\mu)$ as a function of $\mu$ with all systematic uncertainties included.

Observed variations of −2 ln $\Lambda(\mu)$ as a function of $\mu$ with with parameters describing theory uncertainties in background processes fixed to their best-fit values.

Observed variations of −2 ln $\Lambda(\mu)$ as a function of $\mu$ with parameters describing theory uncertainties in background and signal processes fixed to their best-fit values.

Observed variations of −2 ln $\Lambda(\mu)$ as a function of $\mu$ with all systematic uncertainties fixed to their best-fit values.

Cross sections for ggF, VBF, $WH$, $ZH$, $t\bar{t}H$ and $tH$ production modes. The cross sections are normalized to their SM predictions, measured assuming SM values for the decay branching fractions. The black error bars, blue boxes and yellow boxes show the total, systematic, and statistical uncertainties in the measurements, respectively. The gray bands indicate the theory uncertainties on the SM cross section predictions. The level of compatibility between the measurement and the SM prediction corresponds to a $p$-value of $p_{SM}=65\%$.

Correlation matrix from the measurement of ggF, VBF, $WH$, $ZH$, $ttH$ and $tH$ production cross sections without theory uncertainties.

Correlation matrix from the measurement of ggF, VBF, $WH$, $ZH$, $ttH$ and $tH$ production cross sections with the theory uncertainties.

Cross sections for ggF, VBF, $WH$, $ZH$, $t\bar{t}H+tH$ production modes, obtained without including the theory uncertainties on the predicted SM cross sections in the fit. The cross sections are normalized to their SM predictions, measured assuming SM values for the decay branching fractions. The black error bars, blue boxes and yellow boxes show the total, systematic, and statistical uncertainties in the measurements, respectively. The gray bands indicate the theory uncertainties on the SM cross section predictions. The level of compatibility between the measurement and the SM prediction corresponds to a $p$-value of $p_{SM}=89\%$.

Correlation matrix from the measurement of ggF, VBF, $WH$, $ZH$, $ttH+tH$ production cross sections without theory uncertainties.

Correlation matrix from the measurement of the branching fractions, assuming SM values for the decay production cross sections and no contributions from non-SM decays to the total Higgs boson decay width.

Correlation matrix from the measured values of the production cross sections times branching fractions of the Higgs boson, for the combinations in which sufficient sensitivity is provided by the input analyses.

Observed correlation matrix from the fit of $\kappa_{V}$ and $\kappa_{F}$ coupling modifiers. No contributions from non-SM invisible and undetected Higgs boson decays are allowed, i.e. $B_{inv.}=B_{u.}=0$.

Negative log-likelihood contour corresponding to 68% CL in the ($\kappa_{V}$, $\kappa_{F}$) plane, corresponding to $H\to\gamma\gamma$ decay, obtained from a combined fit, assuming no contributions from invisible or undetected non-SM Higgs boson decays.

Negative log-likelihood contour corresponding to 68% CL in the ($\kappa_{V}$, $\kappa_{F}$) plane, corresponding to $H\to ZZ$ decay, obtained from a combined fit, assuming no contributions from invisible or undetected non-SM Higgs boson decays.

Negative log-likelihood contour corresponding to 68% CL in the ($\kappa_{V}$, $\kappa_{F}$) plane, corresponding to $H\to\tau\tau$ decay, obtained from a combined fit, assuming no contributions from invisible or undetected non-SM Higgs boson decays.

Negative log-likelihood contour corresponding to 68% CL in the ($\kappa_{V}$, $\kappa_{F}$) plane, corresponding to $H\to WW$ decay, obtained from a combined fit, assuming no contributions from invisible or undetected non-SM Higgs boson decays.

Negative log-likelihood contour corresponding to 68% CL in the ($\kappa_{V}$, $\kappa_{F}$) plane, corresponding to $H\to bb$ decay, obtained from a combined fit, assuming no contributions from invisible or undetected non-SM Higgs boson decays.

Observed correlation matrix from the fit of $\kappa_{Z}$, $\kappa_{W}$, $\kappa_{b}$, $\kappa_{t}$, $\kappa_{\tau}$ and $\kappa_{\mu}$ coupling modifiers with $\kappa_{c}$ = $\kappa_{t}$ in the fit. All fitted parameters are assumed to be positive. No contributions from non-SM invisible and undetected Higgs boson decays are allowed, i.e. $B_{inv.}=B_{u.}=0$.

Observed correlation matrix from the fit of the $\kappa_{Z}$, $\kappa_{W}$, $\kappa_{b}$, $\kappa_{t}$, $\kappa_{\tau}$, $\kappa_{\mu}$ and effective photon, $Z \gamma$ and gluon coupling modifiers. No contributions from non-SM invisible and undetected Higgs boson decays are allowed, i.e. $B_{inv.}=B_{u.}=0$.

Coupling-strength modifiers and their uncertainties for the effective photon, $Z\gamma$ and gluon couplings as free parameters and all other coupling modifiers set to unity. The scenario where $B_{inv.}=B_{u.}=0$ is assumed is shown as solid lines. The $p$-value of the compatibility with the SM prediction is 63% in this case. The scenario where $B_{inv.}$ and $B_{u.}$ are allowed to contribute to the total Higgs width is shown as dashed lines, with the corresponding $p$-value of 49%. Both $B_{inv.}$ and $B_{u.}$ are allowed to take negative values in the latter case.

Coupling-strength modifiers and their uncertainties for the effective photon, $Z\gamma$ and gluon couplings as free parameters and all other coupling modifiers set to unity. The scenario where $B_{inv.}$ and $B_{u.}$ are allowed to contribute to the total Higgs width is shown as dashed lines, with the corresponding $p$-value of 49%. Both $B_{inv.}$ and $B_{u.}$ are allowed to take negative values. The lower pad shows the 95% CL upper limits on $B_{inv.}$ and $B_{u.}$.

Negative log-likelihood contours corresponding to 68% and 95% CL in the ($\kappa_{g}$, $\kappa_{\gamma}$) plane obtained from a combined fit of ($\kappa_{g}$, $\kappa_{\gamma}$ and $\kappa_{Z\gamma}$) .

Measured ratios of coupling modifiers. The red line indicates the SM value of unity for each parameter. The level of compatibility between the combined measurement and the SM prediction corresponds to a $p$-value of 71% with ten degrees of freedom.

Best-fit values and uncertainties for the cross sections in each measurement region, normalized to the SM predictions for the various parameters. The measurements assume SM branching fractions for all measured decays. The black error bars, blue boxes and yellow boxes show the total, systematic, and statistical uncertainties in the measurements, respectively. The gray bands show the theory uncertainties on the predictions. The level of compatibility between the combined measurement and the SM prediction corresponds to a $p$-value of 94%.

Correlation matrix for the measured values of the simplified template cross sections.

Expected negative log-likelihood scans as a function of $\kappa_{Z}$.

Observed negative log-likelihood scans as a function of $\kappa_{Z}$.

Expected negative log-likelihood scans as a function of $\kappa_{W}$.

Observed negative log-likelihood scans as a function of $\kappa_{W}$.

Expected negative log-likelihood scans as a function of $\kappa_{t}$.

Observed negative log-likelihood scans as a function of $\kappa_{t}$.

Expected negative log-likelihood scans as a function of $\kappa_{b}$.

Observed negative log-likelihood scans as a function of $\kappa_{b}$.

Expected negative log-likelihood scans as a function of $\kappa_{\tau}$.

Observed negative log-likelihood scans as a function of $\kappa_{\tau}$.

Expected negative log-likelihood scans as a function of $\kappa_{\mu}$.

Observed negative log-likelihood scans as a function of $\kappa_{\mu}$.

Expected negative log-likelihood scans as a function of $\kappa_{g}$.

Observed negative log-likelihood scans as a function of $\kappa_{g}$.

Expected negative log-likelihood scans as a function of $\kappa_{\gamma}$.

Observed negative log-likelihood scans as a function of $\kappa_{\gamma}$.

Expected negative log-likelihood scans as a function of $\kappa_{Z\gamma}$.

Observed negative log-likelihood scans as a function of $\kappa_{Z\gamma}$.

Expected negative log-likelihood scans as a function of $B_{inv.}$.

Observed negative log-likelihood scans as a function of $B_{inv.}$.

Expected negative log-likelihood scans as a function of $B_{u.}$.

Observed negative log-likelihood scans as a function of $B_{u.}$.

The acceptances of STXS stage 1.2 kinematic regions in the Higgs boson production processes.

Simultaneous coupling measurement $\kappa_V/\kappa_f$

Simultaneous coupling measurement $\kappa_V/\kappa_f$

Simultaneous coupling measurement $\kappa_V/\kappa_f$

Couplings resolved loops. Numbers are in $\kappa_i$.

Couplings effective loops.

Upper limits at 95% CL on the SM strength modifier for the HH production, per final state and combined.

Upper limits at 95% CL on the HH production cross section for different values of $\kappa_\lambda$.

Upper limits at 95% CL on the HH production cross section for different values of $\kappa_{2V}$.

Signal strength modifiers per production channel tims decay mode $\mu_i^f$.

Couplings effective loops including the invisible and undetermined branching fractions.

$\kappa_{\lambda}$ 68 and 95% confidence intervals, from single and double Higgs productions.

Simultaneous coupling measurement $\kappa_\lambda/\kappa_{2V}$

Simultaneous coupling measurement $\kappa_\lambda/\kappa_{2V}$

Simultaneous coupling measurement $\kappa_\lambda/\kappa_{2V}$

Observed correlations between the signal strength modifiers per production mode, $\mu_i$

Observed correlations between the signal strength modifiers per decay channel, $\mu^f$

Observed correlations between the signal strength modifiers per production mode times decay channel $\mu_i^f$.

Observed correlations between the coupling modifiers in the resolved loop model

Observed correlations between the coupling modifiers in the effective loop model

The $1\sigma_{\text{exp}}$ variation of the expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=1.0$.

The observed exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.1$.

The expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.1$.

The observed exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.15$.

The expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.15$.

The observed exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.2$.

The expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.2$.

The observed exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.4$.

The expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.4$.

The observed exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.6$.

The expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0p6$.

The observed exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=1.5$.

The expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=1.5$.

The observed exclusion contour at 95% CL as a function of the leptoquark mass and coupling strength.

The expected exclusion contour at 95% CL as a function of the leptoquark mass and coupling strength.

The minus $1\sigma_{\text{theory}}$ variation of the observed exclusion contour at 95% CL as a function of the leptoquark mass and coupling strength.

The plus $1\sigma_{\text{theory}}$ variation of the observed exclusion contour at 95% CL as a function of the leptoquark mass and coupling strength.

The $1\sigma_{\text{exp}}$ variation of the expected exclusion contour at 95% CL as a function of the leptoquark mass and coupling strength.

Observed yields, and fake lepton background yields in the $e^{+}\mu^{-}$ and $e^{-}\mu^{+}$ channels of SR-MET, along with the results of the $e^{+}\mu^{-}/e^{-}\mu^{+}$ ratio measurement and 1-sided p-value in SR-MET, binned in $M_{T2}$.

Observed yields, and fake lepton background yields in the $e^{+}\mu^{-}$ and $e^{-}\mu^{+}$ channels of SR-JET, along with the results of the $e^{+}\mu^{-}/e^{-}\mu^{+}$ ratio measurement and 1-sided p-value in SR-JET, binned in $H_{\text{P}}$.

Observed and expected 95% CL upper limits on the total number of signal events entering the $e^{+}\mu^{-}$ and $e^{-}\mu^{+}$ channels of each bin of SR-MET. The regions are binned in the same way as the ratio $\rho$ measurement. The limits are shown for a selection of 'z' values, where 'z' is the fraction of the total signal events entering the $e^{+}\mu^{-}$ channel.

Observed and expected 95% CL upper limits on the total number of signal events entering the $e^{+}\mu^{-}$ and $e^{-}\mu^{+}$ channels of each bin of SR-JET. The regions are binned in the same way as the ratio $\rho$ measurement. The limits are shown for a selection of 'z' values, where 'z' is the fraction of the total signal events entering the $e^{+}\mu^{-}$ channel.

Signal yields following each cut in the analysis, for representative $R$-parity-violating supersymmetry and leptoquark signals. All yields are MC generator-weighted and normalised to $139~\text{fb}^{-1}$. The cut labelled `Preselection' includes trigger requirements, and requires exactly one Baseline electron and one Baseline muon. At this point, the muon charge-bias correction weights are also applied. The $R$-parity-violating supersymmetry models were generated by specifying a top-quark in the final state and applying a two-lepton filter, hence the first row also includes events where the top quark decays to an electron.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 5-jet bHbH channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 5-jet bHbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 5-jet bZbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 5-jet bHtW channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 5-jet bZtW channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 6-jet bHbH channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 6-jet bHbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 6-jet bZbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 6-jet bHtW channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 6-jet bZtW channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the semileptonic 3-jet bHbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the semileptonic 3-jet bZbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the semileptonic 4-jet bHbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the semileptonic 4-jet bZbZ channel.

The limit at 95% CL on the cross section for VLQ pair production for the branching fraction hypothesis 0% $\mathcal{B}(B \to bH)$, 100% $\mathcal{B}(B \to bH)$, and 0% $\mathcal{B}(B \to bH)$

The limit at 95% CL on the cross section for VLQ pair production for the branching fraction hypothesis 25% $\mathcal{B}(B \to bH)$, 25% $\mathcal{B}(B \to bH)$, and 50% $\mathcal{B}(B \to bH)$

The limit at 95% CL on the cross section for VLQ pair production for the branching fraction hypothesis 50% $\mathcal{B}(B \to bH)$, 50% $\mathcal{B}(B \to bH)$, and 0% $\mathcal{B}(B \to bH)$

The limit at 95% CL on the cross section for VLQ pair production for the branching fraction hypothesis 100% $\mathcal{B}(B \to bH)$, 0% $\mathcal{B}(B \to bH)$, and 0% $\mathcal{B}(B \to bH)$

Median expected exclusion limits on the VLQ mass at 95% CL as a function of the branching fractions $\mathcal{B}(B \to bH)$ and $\mathcal{B}(B \to tW)$, with $\mathcal{B}(B \to tW) = 1 - \mathcal{B}(B \to bH) - \mathcal{B}(B \to bZ)$. The grey area corresponds to the region where the exclusion limit is less than 1000 GeV.

Median observed exclusion limits on the VLQ mass at 95% CL as a function of the branching fractions $\mathcal{B}(B \to bH)$ and $\mathcal{B}(B \to tW)$, with $\mathcal{B}(B \to tW) = 1 - \mathcal{B}(B \to bH) - \mathcal{B}(B \to bZ)$. The grey area corresponds to the region where the exclusion limit is less than 1000 GeV.

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Distribution of $\chi^{2}$/ndf for the best jet combination for simulated 1200\GeV VLQ events (red histogram) and data (black points), for 4-jet events (left), 5-jet events (center), and 6-jet events (right). The simulated signal events and data events are normalized to the same value within the displayed $\chi^{2}$/ndf range.

Distribution of $\chi^{2}$/ndf for the best jet combination for simulated 1200\GeV VLQ events (red histogram) and data (black points), for 4-jet events (left), 5-jet events (center), and 6-jet events (right). The simulated signal events and data events are normalized to the same value within the displayed $\chi^{2}$/ndf range.

Distribution of $\chi^{2}$/ndf for the best jet combination for simulated 1200\GeV VLQ events (red histogram) and data (black points), for 4-jet events (left), 5-jet events (center), and 6-jet events (right). The simulated signal events and data events are normalized to the same value within the displayed $\chi^{2}$/ndf range.

Distributions of the average reconstructed mass of VLQ candidates for the jet combination with the least $\chi^{2}$ in 4-jet (left), 5-jet (center), and 6-jet (right) multiplicity events. The red lines show the exponential fit in the range 1000-2000 GeV. The lower panels show the fractional difference, (data-fit)/fit

Distributions of the average reconstructed mass of VLQ candidates for the jet combination with the least $\chi^{2}$ in 4-jet (left), 5-jet (center), and 6-jet (right) multiplicity events. The red lines show the exponential fit in the range 1000-2000 GeV. The lower panels show the fractional difference, (data-fit)/fit

Distributions of the average reconstructed mass of VLQ candidates for the jet combination with the least $\chi^{2}$ in 4-jet (left), 5-jet (center), and 6-jet (right) multiplicity events. The red lines show the exponential fit in the range 1000-2000 GeV. The lower panels show the fractional difference, (data-fit)/fit

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Expected limits on the VLQ mass at 95% CL as a function of the branching fractions $\mathcal{B}$( B -> bH ) and $\mathcal{B}$( B -> bZ)

Observed limits on the VLQ mass at 95% CL as a function of the branching fractions $\mathcal{B}$( B -> bH ) and $\mathcal{B}$( B -> bZ)

The 95% confidence limit on the cross section for VLQ pair production as a function of VLQ mass for three branching fraction hypotheses: B( B -> bH ) = 100% (upper left), B( B -> bZ ) = 100% (upper right), and and B( B -> bH ) = B( B -> bZ ) = 50% (lower). The solid black line indicates the observed limit and the dashed line indicates the expected limit with 1 sigma (green band) and 2 sigma (yellow band) uncertainties. The theoretical cross section and its uncertainty are shown as the red line and pale red band.

The 95% confidence limit on the cross section for VLQ pair production as a function of VLQ mass for three branching fraction hypotheses: B( B -> bH ) = 100% (upper left), B( B -> bZ ) = 100% (upper right), and and B( B -> bH ) = B( B -> bZ ) = 50% (lower). The solid black line indicates the observed limit and the dashed line indicates the expected limit with 1 sigma (green band) and 2 sigma (yellow band) uncertainties. The theoretical cross section and its uncertainty are shown as the red line and pale red band.

The 95% confidence limit on the cross section for VLQ pair production as a function of VLQ mass for three branching fraction hypotheses: B( B -> bH ) = 100% (upper left), B( B -> bZ ) = 100% (upper right), and and B( B -> bH ) = B( B -> bZ ) = 50% (lower). The solid black line indicates the observed limit and the dashed line indicates the expected limit with 1 sigma (green band) and 2 sigma (yellow band) uncertainties. The theoretical cross section and its uncertainty are shown as the red line and pale red band.

A summary of benchmark simplified models, the most sensitive $n_{\mathrm{jet}}$ categories, and representative values for the corresponding experimental acceptance times efficiency ($\mathcal{A}\varepsilon$), the dominant systematic uncertainties, the theoretical production cross section ($\sigma_{\mathrm{theory}}$), and the expected and observed upper limits on the production cross section, expressed in terms of the signal strength parameter ($\mu$).

Summary of the mass limits obtained for the fourteen classes of simplified models. The limits indicate the strongest observed and expected (in parentheses) mass exclusions in $\tilde{g}$, $\tilde{q}$, $\tilde{b}$, $\tilde{t}$, and $\tilde{\chi}_1^0$. The quoted values have uncertainties of $\pm 25$ and $\pm 10$ GeV for models involving the pair production of, respectively, gluinos and squarks.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1qqqq model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1qqqq model. This line is the observed exclusion.

Covariance matrix for the SM background estimates obtained using the simplified binning scheme, determined from a simultaneous fit to data in the control regions only (CR-only fit). The uncertainties in the background estimates are correlated in such a way that the covariance is typically positive. Small positive values, as well as the few negative values, are not shown.

Correlation matrix for the SM background estimates obtained using the simplified binning scheme, determined from a simultaneous fit to data in the control regions only (CR-only fit).

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{q} } }, m_{\tilde{\chi}^0_1 })$ plane for the T2qq model. This line is the expected exclusion.

Observed data counts and "post-fit" background expectations based on the result of a combined fit to the signal region and multiple control regions under the SM-only hypothesis for the "monojet" event category. The rows labelled SM "pre-fit' show the background expectations when excluding the signal region from the fit. The uncertainties include statistical as well as systematic contributions.

Observed data counts and "post-fit" background expectations based on the result of a combined fit to the signal region and multiple control regions under the SM-only hypothesis for the "asymmetric" event categories. The rows labelled SM "pre-fit" show the background expectations when excluding the signal region from the fit. The uncertainties include statistical as well as systematic contributions.

Observed data counts and "post-fit" background expectations based on the result of a combined fit to the signal region and multiple control regions under the SM-only hypothesis for the "symmetric" event categories. The rows labelled SM "pre-fit" show the background expectations when excluding the signal region from the fit. The uncertainties include statistical as well as systematic contributions.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{q} } }, m_{\tilde{\chi}^0_1 })$ plane for the T2qq model. This line is the observed exclusion.

Summary of a simplified binning scheme comprising $H_{\mathrm{T}}^{\mathrm{miss}}$ shapes for eight representative event topologies, defined in terms of requirements on $n_{\mathrm{jet}}$ both symmetric and asymmetric categories) and $n_{\mathrm{b}}$. For each topology, event yields are integrated over the full $H_{\mathrm{T}}$ range, $H_{\mathrm{T}} > 200$ GeV, and are categorised according to eight bins in $H_{\mathrm{T}}^{\mathrm{miss}}$ -- $130 < H_{\mathrm{T}}^{\mathrm{miss}} < 200$ GeV, six bins 100 GeV wide in the region $200 < H_{\mathrm{T}}^{\mathrm{miss}} < 800$ GeV, and an open final bin, $H_{\mathrm{T}}^{\mathrm{miss}} > 800$ GeV. This categorisation scheme leads to 64 bins that are exclusive, contiguous, and provide complete coverage of the signal region. The corresponding SM background estimates for these 64 bins are obtained using the nominal likelihood model. The $H_{\mathrm{T}}^{\mathrm{miss}}$ template for each topology is determined from simulation and the normalisation for each template is given by an aggregate of estimates from several ($n_{\mathrm{jet}}$, $n_{\mathrm{b}}$, $H_{\mathrm{T}}$) bins, defined according to the nominal binning scheme.

Observed data yields and "CR-only fit" SM background expectations using the simplified binning scheme, as a function of topology and $H_{\mathrm{T}}^{\mathrm{miss}}$. The SM backgrounds are determined from a simultaneous fit to data in the control regions. The data yields in the signal region are masked and excluded from the fit. The uncertainties include statistical as well as systematic contributions.

Expected ($\mu_{\mathrm{exp}}$) and observed ($\mu_{\mathrm{obs}}$) upper limits on the production cross section, expressed in terms of the signal strength parameter, obtained using both the nominal and simplified binning schema. Comparable values are obtained from the two schema for these benchmark SUSY models because the signal events typically populate bins at high values of $H_{\mathrm{T}}$ and $H_{\mathrm{T}}^{\mathrm{miss}}$ and so are at or near the limit of the search sensitivity. Larger differences in $\mu$ can be observed for models that populate the core of these distributions.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1bbbb model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1bbbb model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1tttt model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1tttt model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1ttbb model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1qqqq model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1ttbb model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1qqqq model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1qqqq model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5tttt_DM175 model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1qqqq model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{q} } }, m_{\tilde{\chi}^0_1 })$ plane for the T2qq model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{q} } }, m_{\tilde{\chi}^0_1 })$ plane for the T2qq model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5tttt_DM175 model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{q} } }, m_{\tilde{\chi}^0_1 })$ plane for the T2qq model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{q} } }, m_{\tilde{\chi}^0_1 })$ plane for the T2qq model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1bbbb model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5ttcc model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1bbbb model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1bbbb model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1bbbb model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5ttcc model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1tttt model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1tttt model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1tttt model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ b } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bb model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1tttt model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1ttbb model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1ttbb model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ b } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bb model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1ttbb model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1ttbb model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5tttt_DM175 model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tb model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5tttt_DM175 model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5tttt_DM175 model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5tttt_DM175 model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tb model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5ttcc model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5ttcc model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5ttcc model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5ttcc model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ b } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bb model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ b } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bb model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ b } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bb model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ b } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bb model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tb model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2cc model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tb model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tb model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tb model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2cc model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt_degen model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2cc model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2cc model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt_degen model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2cc model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2cc model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt_degen model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2_mixed model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt_degen model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt_degen model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt_degen model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2_mixed model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2_mixed model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2_mixed model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2_mixed model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bW model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2_mixed model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bW model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bW model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bW model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bW model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bW model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Comparison of the $m_{miss}$ shapes for the simulated signal events within the fiducial region and those outside it, after including the effect of PU protons as describe in the text, for a generated $m_{X}$ mass of 1000 GeV. The distributions are shown for single(+z)-single(−z) proton reconstruction categories.

Product of the acceptance times the combined reconstruction and identification efficiency, as a function of $m_{X}$ , for events generated inside the fiducial volume defined in Table 1. The curves shown panel display the different final states. The definition of the fiducial region and signal model used to estimate the acceptance is provided in the text.

Product of the acceptance times the combined reconstruction and identification efficiency, as a function of $m_{X}$ , for different proton reconstruction categories in the $Z \rightarrow \mu\mu $ analysis, and for events generated inside the fiducial volume defined in Table 1. The curves shown panel display the different final states. The definition of the fiducial region and signal model used to estimate the acceptance is provided in the text.

Distributions of the reconstructed proton $\xi$ in the negative arm from the proton mixing procedure with simulated MC events, are compared to data. The ee final states are shown.The ee events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the positive arm from the proton mixing procedure with simulated MC events, are compared to data. The ee final states are shown.The ee events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed di-proton rapidity from the proton mixing procedure with simulated MC events, are compared to data. The ee final states are shown.The ee events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the negative arm from the proton mixing procedure with simulated MC events, are compared to data. The $\mu\mu$ final states are shown.The $\mu\mu$ events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the positive arm from the proton mixing procedure with simulated MC events, are compared to data. The $\mu\mu$ final states are shown.The $\mu\mu$ events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed di-proton rapidity from the proton mixing procedure with simulated MC events, are compared to data. The $\mu\mu$ final states are shown.The $\mu\mu$ events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the negative arm from the proton mixing procedure with simulated MC events, are compared to data. The $\gamma$ final states are shown.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the positive arm from the proton mixing procedure with simulated MC events, are compared to data. The $\gamma$ final states are shown.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed di-proton rapidity from the proton mixing procedure with simulated MC events, are compared to data. The $\gamma$ final states are shown.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The proton $\xi$ distribution for the positive CT-PPS arm is shown.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The proton $\xi$ distribution for the negative CT-PPS arm is shown.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The di-proton invariant mass is shown.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The $m_{miss}$ is shown.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the multi(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the multi(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the single(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the single(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the multi(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the multi(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the single(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the single(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the multi(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the multi(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the single(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the single(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $Z\rightarrow ee$ analysis.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $Z\rightarrow \mu\mu$ analysis.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $Z$ analysis.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $\gamma$ analysis.

Results of applying the BH algorithm to the mass spectra in the leading photon category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading electron category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading muon category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading small-R jet category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading bjet category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading large-R jet category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading photon category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading electron category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading muon category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the inclusive category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the inclusive category, using the fitted background estimations from the initial step

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading small-R jet category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading bjet category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading large-R jet category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading photon category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading electron category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading muon category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading small-R jet category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading bjet category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading large-R jet category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading photon category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading electron category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading muon category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the inclusive category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the inclusive category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading small-R jet category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading bjet category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading large-R jet category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading photon category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading electron category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading muon category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading small-R jet category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading bjet category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading large-R jet category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading photon category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading electron category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading muon category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the inclusive category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the inclusive category

Acceptance times efficiency in an HVT model as a function of mass in the leading large-R jet category

Upper limits at 95% CL on the cross section times the branching fraction($W^{\prime} \to ZW$) for the HVT signal as functions of mass($m_{ZX}$) in the leading large-R jet category.The full(dashed) red curves correspond to the theoretical predictions of HVT model A(B)

Comparison of the identification efficiencies using standard and merged-ee reconstruction as a function of true PT(Z)

Background rejection factor as a function of signal efficiency. The red curve shows the BDT performance whereas the blue curve corresponds to that of a cut-based analysis relying only on the jet energyfraction deposited in the EM calorimeter

BH p-values of the 100 pseudo-experiments as a function of in the leading small-R jet category for an injected Gaussian-shaped signal with a relative width value of 3%. The fractions of the pseudo-experiments that have the correctly identified interval and BH p-values below the threshold of 0.01 are indicated. The background is derived by the background-only fit in the full fitting range.

BH p-values of the 100 pseudo-experiments as a function of in the leading small-R jet category for an injected Gaussian-shaped signal with a relative width value of 3%. The fractions of the pseudo-experiments that have the correctly identified interval and BH p-values below the threshold of 0.01 are indicated. The background is derived by the background-only fit in the range excluding the BH interval.

Fractions of pseudo-experiments in which the detected BH interval agrees with the injected mass point and the BH p-value is below 0.01 as a function of mass in the leading small-R jet category for Gaussian-shaped signal with a relative width of 3%.

Fractions of pseudo-experiments in which the detected BH interval agrees with the injected mass point and the BH p-value is below 0.01 as a function of mass in the leading small-R jet category for Gaussian-shaped signal with a relative width of 3%.

Distribution of exclusion upper limits on signal event yields at 95% CL from 1000 pseudo-experiments for Gaussian-shaped signals with relative width values of 3% at the low boundary of the limit-sensitive mass range for the ZX spectrum of the leading small-R jet category. The vertical line corresponds to the expected nominal exclusion limit at 95% CL.

Distribution of exclusion upper limits on signal event yields at 95% CL from 1000 pseudo-experiments for Gaussian-shaped signals with relative width values of 3% at the high boundary of the limit-sensitive mass range for the ZX spectrum of the leading small-R jet category. The vertical line corresponds to the expected nominal exclusion limit at 95% CL.

Distribution of exclusion upper limits on signal event yields at 95% CL from 1000 pseudo-experiments for Gaussian-shaped signals with relative width values of 5% at the low boundary of the limit-sensitive mass range for the ZX spectrum of the leading small-R jet category. The vertical line corresponds to the expected nominal exclusion limit at 95% CL.

Distribution of exclusion upper limits on signal event yields at 95% CL from 1000 pseudo-experiments for Gaussian-shaped signals with relative width values of 5% at the high boundary of the limit-sensitive mass range for the ZX spectrum of the leading small-R jet category. The vertical line corresponds to the expected nominal exclusion limit at 95% CL.

Distribution of exclusion upper limits on signal event yields at 95% CL from 1000 pseudo-experiments for Gaussian-shaped signals with relative width values of 10% at the low boundary of the limit-sensitive mass range for the ZX spectrum of the leading small-R jet category. The vertical line corresponds to the expected nominal exclusion limit at 95% CL.

Distribution of exclusion upper limits on signal event yields at 95% CL from 1000 pseudo-experiments for Gaussian-shaped signals with relative width values of 10% at the high boundary of the limit-sensitive mass range for the ZX spectrum of the leading small-R jet category. The vertical line corresponds to the expected nominal exclusion limit at 95% CL.

Data yields of the six event categories in the $Z\to e^+e^-$ and $\mu^+\mu^-$ decay channels. The merged-$e^+e^-$ events are included in the $e^+e^-$ channel, increasing the event yield, mainly in the leading large-$R$-jet category, by 0.6 %.

A list of mass spectra, event categories and their corresponding fit ranges, functional forms, numbers of free parameters and global $\chi^2$ $p$-values from background-only fits. The fit range values are rounded to the nearest 5 GeV. Here $f_1(x)=p_0\left(\mathrm{e}^{-p_1x}+p_2\mathrm{e}^{-(p_1+p_3)x}+p_4\mathrm{e}^{-(p_1+p_3+p_5)x}+\cdots\right)$ and $f_2(x)=p_0(1-x)^{p_1}x^{p_2+p_3\ln\!x+p_4\ln^2\!x+\cdots}$

A list of mass spectra, event categories and their corresponding signal-sensitive mass ranges. The initial BH $p$-value is obtained by using the background derived from the background-only fit in the full fit range, whereas the new BH $p$-value uses the background derived from the background-only fit in the range excluding the initial BH interval.

A list of mass spectra, event categories, relative width values of Gaussian-shaped signals and limit-sensitive mass ranges and fractions, corresponding to the mass values in the previous column, of pseudo-experiments having exclusion upper limits higher than the nominal exclusion limit at 95% CL

Cutflow of HVT model $A$ signals ($W^\prime \to ZW \to \ell\ell qq$) with $m_{W^\prime} = 1$ TeV and $m_{W^\prime} = 4$ TeV based on MC simulations.

Acceptance times efficiency ($\mathcal{A} \times \epsilon$) values in % in the dominant event category for $p^Z_{\mathrm{T}} > 100$ GeV in the $Z \to \ell^{+}\ell^{-}$ decay channel.

NMSSM 95% CL upper limit on cross section times branching ratio

NMSSM 95% CL upper limit on cross section times branching ratio

NMSSM 95% CL upper limit on cross section times branching ratio

NMSSM 95% CL upper limit on cross section times branching ratio

90% CL upper limit on MSSMD kinetic mixing parameter epsilon

90% CL upper limit on MSSMD kinetic mixing parameter epsilon

90% CL upper limit on MSSMD kinetic mixing parameter epsilon

90% CL upper limit on MSSMD kinetic mixing parameter epsilon