The normalized distributions of the BDT input variable $D_{p_{\text{T}}}$ for the two highest $p_{\text{T}}$ jets from the simulated SM backgrounds and several signal models. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized BDT discriminator distribution for the two highest $p_{\text{T}}$ jets from the simulated SM backgrounds and several signal models.

The BDT ROC curves for the two highest $p_{\text{T}}$ jets, comparing the simulated SM backgrounds with one signal model with $m_{\text{Z}^{\prime}} = 3.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, $r_{\text{inv}} = 0.3$, $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$.

The BDT ROC curves for the two highest $p_{\text{T}}$ jets, comparing the simulated SM backgrounds with one signal model with $m_{\text{Z}^{\prime}} = 3.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, $r_{\text{inv}} = 0.3$, $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$.

The BDT ROC curves for the two highest $p_{\text{T}}$ jets, comparing the simulated SM backgrounds with one signal model with $m_{\text{Z}^{\prime}} = 3.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, $r_{\text{inv}} = 0.3$, $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$.

The BDT ROC curves for the two highest $p_{\text{T}}$ jets, comparing the simulated SM backgrounds with one signal model with $m_{\text{Z}^{\prime}} = 3.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, $r_{\text{inv}} = 0.3$, $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$.

The $m_{\text{T}}$ distribution for the high-$R_{\text{T}}$ signal region, comparing the observed data to the background prediction from the analytic fit ($g_{3}(x) = \exp(p_{1}x)x^{p_{2}(1+p_{3}\ln(x))}$, $x = m_{\text{T}}/\sqrt{s}$). The distributions from several example signal models, with cross sections corresponding to the observed limits, are superimposed.

The $m_{\text{T}}$ distribution for the low-$R_{\text{T}}$ signal region, comparing the observed data to the background prediction from the analytic fit ($g_{2}(x) = \exp(p_{1}x)x^{p_{2}}$, $x = m_{\text{T}}/\sqrt{s}$). The distributions from several example signal models, with cross sections corresponding to the observed limits, are superimposed.

The $m_{\text{T}}$ distribution for the high-SVJ2 signal region, comparing the observed data to the background prediction from the analytic fit ($g_{2}(x) = \exp(p_{1}x)x^{p_{2}}$, $x = m_{\text{T}}/\sqrt{s}$). The distributions from several example signal models, with cross sections corresponding to the observed limits, are superimposed.

The $m_{\text{T}}$ distribution for the low-SVJ2 signal region, comparing the observed data to the background prediction from the analytic fit ($g_{2}(x) = \exp(p_{1}x)x^{p_{2}}$, $x = m_{\text{T}}/\sqrt{s}$). The distributions from several example signal models, with cross sections corresponding to the observed limits, are superimposed.

The 95% CL observed upper limits on the product of the cross section and branching fraction from the inclusive search for variations of the mediator mass and the dark hadron mass.

The expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The observed exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The observed exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The 95% CL observed upper limits on the product of the cross section and branching fraction from the inclusive search for variations of the mediator mass and the invisible fraction.

The expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The observed exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The 95% CL observed upper limits on the product of the cross section and branching fraction from the BDT-based search for variations of the mediator mass and the dark hadron mass.

The expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The observed exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The observed exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The observed exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The observed exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The 95% CL observed upper limits on the product of the cross section and branching fraction from the BDT-based search for variations of the mediator mass and the invisible fraction.

The expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The observed exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The observed exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The 95% CL upper limits on the product of the cross section and branching fraction from the inclusive search for the $\alpha_{\text{dark}}$ variations.

The 95% CL upper limits on the product of the cross section and branching fraction from the BDT-based search for the $\alpha_{\text{dark}}$ variations.

The three two-dimensional signal model parameter scans.

Metrics representing the performance of the BDT for the benchmark signal model ($m_{\text{Z}^{\prime}} = 3.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, $r_{\text{inv}} = 0.3$, $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$), compared to each of the major SM background processes.

The range of effects on the signal yield for each systematic uncertainty and the total. Values less than 0.01% are rounded to 0.0%.

The normalized distribution of the variable $m_{\text{T}}$ for the simulated SM backgrounds and several signal models. The requirement on this variable is omitted, but all other preselection requirements are applied. The last bin of each histogram includes the overflow events.

The normalized distribution of the variable $\Delta\eta(\text{J}_{1},\text{J}_{2})$ for the simulated SM backgrounds and several signal models. The requirement on this variable is omitted, but all other preselection requirements are applied. The last bin of each histogram includes the overflow events.

The normalized distribution of the variable $p_{\text{T}}^{\text{miss}}$ for the simulated SM backgrounds and several signal models. The $R_{\text{T}}$ requirement is omitted, but all other preselection requirements are applied. The last bin of each histogram includes the overflow events.

The normalized distribution of the variable $N_{\text{e}}$ for the simulated SM backgrounds and several signal models. The requirement on this variable is omitted, but all other preselection requirements are applied. The last bin of each histogram includes the overflow events.

The normalized distribution of the variable $N_{\mu}$ for the simulated SM backgrounds and several signal models. The requirement on this variable is omitted, but all other preselection requirements are applied. The last bin of each histogram includes the overflow events.

The normalized distribution of $\Delta\eta(\text{J}_{1},\text{J}_{2})$ vs. $R_{\text{T}}$ for the simulated QCD background. The preselection requirements on both variables are omitted, but all other preselection requirements are applied.

The normalized distribution of $p_{\text{T}}^{\text{miss}}$ vs. $m_{\text{T}}$ for the simulated QCD background. All selection requirements are omitted, except for the requirement of two high-$p_{\text{T}}$ wide jets.

The normalized distribution of $R_{\text{T}}$ vs. $m_{\text{T}}$ for the simulated QCD background. All selection requirements are omitted, except for the requirement of two high-$p_{\text{T}}$ wide jets.

The normalized distributions of the BDT input variable $\tau_{21}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $m_{\text{dark}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $\tau_{32}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $m_{\text{dark}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $N_{2}^{(1)}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $m_{\text{dark}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $N_{3}^{(1)}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $m_{\text{dark}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $g_{\text{jet}}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $r_{\text{inv}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $\sigma_{\text{major}}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $r_{\text{inv}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $\sigma_{\text{minor}}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $r_{\text{inv}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $\Delta\phi(\vec{J},\vec{p}_{\text{T}}^{\text{miss}})$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $r_{\text{inv}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $f_{\text{h}^{\pm}}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $m_{\text{Z}^{\prime}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $f_{\text{e}}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $m_{\text{Z}^{\prime}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $f_{\mu}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $m_{\text{Z}^{\prime}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $f_{\text{h}^{0}}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $m_{\text{Z}^{\prime}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $f_{\gamma}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $m_{\text{Z}^{\prime}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The product of signal acceptance and efficiency in the high-$R_{\text{T}}$ signal region, for variations of the mediator mass and the dark hadron mass.

The product of signal acceptance and efficiency in the high-$R_{\text{T}}$ signal region, for variations of the mediator mass and the invisible fraction.

The product of signal acceptance and efficiency in the high-$R_{\text{T}}$ signal region, for variations of the mediator mass and the dark coupling strength.

The product of signal acceptance and efficiency in the low-$R_{\text{T}}$ signal region, for variations of the mediator mass and the dark hadron mass.

The product of signal acceptance and efficiency in the low-$R_{\text{T}}$ signal region, for variations of the mediator mass and the invisible fraction.

The product of signal acceptance and efficiency in the low-$R_{\text{T}}$ signal region, for variations of the mediator mass and the dark coupling strength.

The product of signal acceptance and efficiency in the high-SVJ2 signal region, for variations of the mediator mass and the dark hadron mass.

The product of signal acceptance and efficiency in the high-SVJ2 signal region, for variations of the mediator mass and the invisible fraction.

The product of signal acceptance and efficiency in the high-SVJ2 signal region, for variations of the mediator mass and the dark coupling strength.

The product of signal acceptance and efficiency in the low-SVJ2 signal region, for variations of the mediator mass and the dark hadron mass.

The product of signal acceptance and efficiency in the low-SVJ2 signal region, for variations of the mediator mass and the invisible fraction.

The product of signal acceptance and efficiency in the low-SVJ2 signal region, for variations of the mediator mass and the dark coupling strength.

Comparison of different the dijet mass $m_{\text{J}\text{J}}$, the transverse mass $m_{\text{T}}$, and the Monte Carlo (MC) mass $m_{\text{MC}}$ for a signal model with $m_{\text{Z}^{\prime}} = 2.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, $r_{\text{inv}} = 0.3$, $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$. No selection is applied, except that there must be at least two jets. $m_{\text{MC}}$ is computed by adding the generator-level four-vectors for invisible particles to the dijet system, to represent the achievable resolution if the invisible component were fully measured. The last bin of each histogram includes the overflow events.

$m_{\text{T}}$ distributions for signal models with different $m_{\text{dark}}$ values for the high-$R_{\text{T}}$ inclusive signal region.

$m_{\text{T}}$ distributions for signal models with different $m_{\text{dark}}$ values for the low-$R_{\text{T}}$ inclusive signal region.

$m_{\text{T}}$ distributions for signal models with different $r_{\text{inv}}$ values for the high-$R_{\text{T}}$ inclusive signal region.

$m_{\text{T}}$ distributions for signal models with different $r_{\text{inv}}$ values for the low-$R_{\text{T}}$ inclusive signal region.

$m_{\text{T}}$ distributions for signal models with different $\alpha_{\text{dark}}$ values for the high-$R_{\text{T}}$ inclusive signal region.

$m_{\text{T}}$ distributions for signal models with different $\alpha_{\text{dark}}$ values for the low-$R_{\text{T}}$ inclusive signal region.

The proportions of each SM background process in the high-$R_{\text{T}}$ signal region.

The proportions of each SM background process in the low-$R_{\text{T}}$ signal region.

The proportions of each SM background process in the high-SVJ2 signal region.

The proportions of each SM background process in the low-SVJ2 signal region.

The 95% CL expected upper limits on the product of the cross section and branching fraction from the inclusive search for variations of the mediator mass and the dark hadron mass.

The 95% CL expected upper limits on the product of the cross section and branching fraction from the inclusive search for variations of the mediator mass and the invisible fraction.

The 95% CL expected upper limits on the product of the cross section and branching fraction from the BDT-based search for variations of the mediator mass and the dark hadron mass.

The 95% CL expected upper limits on the product of the cross section and branching fraction from the BDT-based search for variations of the mediator mass and the invisible fraction.

Relative efficiencies in % for each step of the event selection process for the major background processes. Statistical uncertainties, at most 1.8%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Relative efficiencies in % for each step of the event selection process for signals with $m_{\text{Z}^{\prime}} = 2.1\,\text{TeV}$, varying $m_{\text{dark}}$ values, $r_{\text{inv}} = 0.3$, and $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$. Statistical uncertainties, at most 0.5%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Relative efficiencies in % for each step of the event selection process for signals with $m_{\text{Z}^{\prime}} = 3.1\,\text{TeV}$, varying $m_{\text{dark}}$ values, $r_{\text{inv}} = 0.3$, and $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$. Statistical uncertainties, at most 0.4%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Relative efficiencies in % for each step of the event selection process for signals with $m_{\text{Z}^{\prime}} = 4.1\,\text{TeV}$, varying $m_{\text{dark}}$ values, $r_{\text{inv}} = 0.3$, and $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$. Statistical uncertainties, at most 0.4%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Relative efficiencies in % for each step of the event selection process for signals with $m_{\text{Z}^{\prime}} = 2.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, varying $r_{\text{inv}}$ values, and $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$. Statistical uncertainties, at most 2.6%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Relative efficiencies in % for each step of the event selection process for signals with $m_{\text{Z}^{\prime}} = 3.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, varying $r_{\text{inv}}$ values, and $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$. Statistical uncertainties, at most 1.2%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Relative efficiencies in % for each step of the event selection process for signals with $m_{\text{Z}^{\prime}} = 4.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, varying $r_{\text{inv}}$ values, and $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$. Statistical uncertainties, at most 0.9%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Relative efficiencies in % for each step of the event selection process for signals with $m_{\text{Z}^{\prime}} = 2.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, $r_{\text{inv}} = 0.3$, and varying $\alpha_{\text{dark}}$ values. Statistical uncertainties, at most 0.4%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Relative efficiencies in % for each step of the event selection process for signals with $m_{\text{Z}^{\prime}} = 3.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, $r_{\text{inv}} = 0.3$, and varying $\alpha_{\text{dark}}$ values. Statistical uncertainties, at most 0.4%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Relative efficiencies in % for each step of the event selection process for signals with $m_{\text{Z}^{\prime}} = 4.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, $r_{\text{inv}} = 0.3$, and varying $\alpha_{\text{dark}}$ values. Statistical uncertainties, at most 0.4%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Distribution of the three leading leptons flavour in the SR-WZ with uncertainties evaluated after the inclusive cross section fit

Efficiency, acceptance, and proportion of events with leptonic tau decays in WZ production

WZ fiducial cross section in the four flavour exclusive and the flavour inclusive channels

WZ total cross section extrapolated from the four flavour exclusive and the flavour inclusive channels

Distribution of the total lepton charge in the SR-WZ with uncertainties evaluated after the inclusive cross section fit

W$^{+}$Z fiducial cross section in the four flavour exclusive and the flavour inclusive channels

W$^{-}$Z fiducial cross section in the four flavour exclusive and the flavour inclusive channels

WZ charge asymmetry ratio measured on each of the four flavour exclusive and the flavour inclusive channels

Distribution of the cosine of the W polarization angle times total lepton charge in the SR-WZ with uncertainties evaluated after the W polarization fit

Distribution of the cosine of the Z polarization angle in the SR-WZ with uncertainties evaluated after the Z polarization fit

Best fits to the W and Z polarization fractions

2D confidence regions at the 68, 95, and 99% CL in the $f_O^W$-$f_{L}^W-f_R^W$ plane

2D confidence regions at the 68, 95, and 99% CL in the $f_O^Z$-$f_{L}^Z-f_R^Z$ plane

Distribution of the invariant mass of the WZ system in the SR-WZ with uncertainties evaluated after the inclusive cross section fit

Best fit values and one dimensional confidence regions in several EFT coefficients obtained from the EFT fit considering both the SM interferences and purely BSM (order $\Lambda^{-2}$ and $\Lambda^{-4}$) terms

2D confidence regions at the 68, 95, and 99% CL in the $c_{www}$-$c_{w}$ plane

2D confidence regions at the 68, 95, and 99% CL in the $c_{w}$-$c_{b}$ plane

2D confidence regions at the 68, 95, and 99% CL in the $c_{www}$-$c_{w}$ plane

Best fit values and one dimensional confidence regions in several EFT coefficients obtained from the EFT fit considering only the SM-EFT interference (order $\Lambda^{-2}$) terms

Evolution of the best fit and expected and observed 95% CI for the $c_{w}$ parameter as a function of the cutoff scale

Evolution of the best fit and expected and observed 95% CI for the $c_{b}$ parameter as a function of the cutoff scale

Evolution of the best fit and expected and observed 95% CI for the $c_{www}$ parameter as a function of the cutoff scale

Evolution of the best fit and expected and observed 95% CI for the $\tilde{c}_{www}$ parameter as a function of the cutoff scale

Evolution of the best fit and expected and observed 95% CI for the $\tilde{c}_{w}$ parameter as a function of the cutoff scale

Differential cross section with respect to the transverse momentum of the Z boson

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the $p_{T}$ of the Z boson

Response matrix for the $p_{T}$ of the Z boson obtained with POWHEG

Differential cross section with respect to the transverse momentum of the leading jet

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the $p_{T}$ of the leading jet

Response matrix for the $p_{T}$ of the leading jet obtained with POWHEG

Differential cross section with respect to the jet multiplicity

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the jet multiplicity

Response matrix for the jet multiplicity obtained with POWHEG

Differential cross section with respect to the invariant mass of the WZ system

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the invariant mass of the WZ system

Response matrix for the invariant mass of the WZ system obtained with POWHEG

Differential cross section with respect to the transverse momentum of the lepton associated to the W boson

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the $p_{T}$ of the lepton associated to the W boson

Response matrix for the $p_{T}$ of the lepton associated to the W boson obtained with POWHEG

Differential cross section with respect to the transverse momentum of the lepton associated to the W boson, W$^{+}$Z only

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the $p_{T}$ of the lepton associated to the W boson, W$^{+}$Z only

Differential cross section with respect to the transverse momentum of the lepton associated to the W boson, W$^{-}$Z only

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the $p_{T}$ of the lepton associated to the W boson, W$^{-}$Z only

Differential cross section with respect to the cosine of the W polarization angle times total lepton charge

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the cosine of the W polarization angle times total lepton charge

Response matrix for the cosine of the W polarization angle times total lepton charge obtained with POWHEG

Differential cross section with respect to the cosine of the W polarization angle times total lepton charge, W$^{+}$Z only

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the cosine of the W polarization angle times total lepton charge, W$^{+}$Z only

Differential cross section with respect to the cosine of the W polarization angle times total lepton charge, W$^{-}$Z only

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the cosine of the W polarization angle times total lepton charge, W$^{-}$Z only

Differential cross section with respect to the cosine of the Z polarization angle

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the cosine of the Z polarization angle

Response matrix for the cosine of the Z polarization angle obtained with POWHEG

Differential cross section with respect to the cosine of the Z polarization angle, W$^{+}$Z only

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the cosine of the Z polarization angle, W$^{+}$Z only

Differential cross section with respect to the cosine of the Z polarization angle, W$^{-}$Z only

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the cosine of the Z polarization angle, W$^{-}$Z only

Differential signal yields for various signal hypotheses.

Differential signal yields for various signal hypotheses.

Differential signal yields for various signal hypotheses.

Differential signal yields for various signal hypotheses.

Differential signal yields for various signal hypotheses.

Background and data yields in the control and signal region bins. The prediction before ("prefit") and after the background only fit ("b-only") are given separately.

Background and data yields in the control and signal region bins. The prediction before ("prefit") and after the background only fit ("b-only") are given separately.

Background and data yields in the control and signal region bins. The prediction before ("prefit") and after the background only fit ("b-only") are given separately.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit to the control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit to the corresponding control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit to the control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit to the corresponding control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit to the control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit to the corresponding control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Median Expected exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with axial couplings.

Observed exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with axial couplings.

Expected plus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with axial couplings.

Expected minus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with axial couplings.

Median Expected exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with vector couplings.

Observed exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with vector couplings.

Expected plus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with vector couplings.

Expected minus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with vector couplings.

Upper limits on the coupling $g_{\chi}$ in the simplified model with a axial mediator.

Upper limits on the coupling $g_{q}$ in the simplified model with a axial mediator.

Upper limits on the coupling $g_{\chi}$ in the simplified model with a vector mediator.

Upper limits on the coupling $g_{q}$ in the simplified model with a vector mediator.

Exclusion limits on the signal strength in the simplified model with scalar couplings.

Exclusion limits on the signal strength in the simplified model with pseudoscalar couplings.

Exclusion limits on the fundamental Planck scale $M_{D}$ as a function of the number of extra dimensions $d$.

Median Expected exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the fermion portal model.

Observed exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the fermion portal model.

Expected plus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the fermion portal model.

Expected minus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the fermion portal model.

Tagging efficiency for AK8 jets. The efficiency includes the effect of the machine-learning based DeepAK8 tagger, as well as the application of the mass window requirement on the jet. The efficiency is split depending on the matching generator-level object. For reinterpretation purposes, this efficiency can directly be applied to any AK8 jet that is matched to a given type of generator-level object, with no prior selection on the jet mass. Other acceptance requirements on jet $\p_{T}$ and $\eta$ should still be applied.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 7 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The cluster efficiency in bins of hadronic and EM energy in region A. Region A is defined as 391 cm $< r <$ 695.5 cm and 400 cm $< |z| <$ 671 cm. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The cluster efficiency in bins of hadronic and EM energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The efficiency of $N_{station} > 1$ requirement in bins of hadronic energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bin corresponds to LLPs that decayed leptonically with 0 hadronic energy. The efficiency is calculated with respect to clusters that pass all cluster-level cuts described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 10%.

The geometric signal acceptance as a function of $c\tau$. The acceptance is shown for four different LLP mass hypotheses: 7, 15, 40, and 55 GeV. The acceptance is defined by requiring at least 1 LLP to decay in the region defined as 400 cm $< |z| <$ 1100 cm, $r < $ 695.5 cm, and $|\eta| <$ 2.4.

The product of signal acceptance and efficiency in the 2l categories for the signal produced via vector boson fusion.

$m_{X}^{T}$ distribution in data in the 0l 2b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}^{T}$ distribution in data in the 0l $\leq$1b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}$ distribution in data in the 2e 2b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2e $\leq$1b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ 2b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ $\leq$1b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}^{T}$ distribution in data in the 0l 2b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}^{T}$ distribution in data in the 0l $\leq$1b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}$ distribution in data in the 2e 2b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2e $\leq$1b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ 2b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ $\leq$1b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

Observed and expected 95% CL upper limit on $\sigma \mathcal{B}$(Z'-> ZH) with all categories combined for the non-VBF signal, including all statistical and systematic uncertainties. The inner green band and the outer yellow band indicate the regions containing 68 and 95%, respectively, of the distribution of expected limits under the background-only hypothesis. The CMS search for a heavy resonance using 2016 data and the same final state [JHEP 11 (2018) 172] is shown as a comparison.

Observed and expected 95% CL upper limit on $\sigma \mathcal{B}$(Z'-> ZH) with all categories combined for the VBF signal, including all statistical and systematic uncertainties. The inner green band and the outer yellow band indicate the regions containing 68 and 95%, respectively, of the distribution of expected limits under the background-only hypothesis.

Observed exclusion limit in the space of the HVT model parameters [$g_{V}c_{H}$, $g^{2}c_{F}/g_{V}$] for mass hypotheses of 2 TeV for the non-VBF signal.

Observed exclusion limit in the space of the HVT model parameters [$g_{V}c_{H}$, $g^{2}c_{F}/g_{V}$] for mass hypotheses of 3 TeV for the non-VBF signal.

Observed exclusion limit in the space of the HVT model parameters [$g_{V}c_{H}$, $g^{2}c_{F}/g_{V}$] for mass hypotheses of 4 TeV for the non-VBF signal.

The $L_{T}$ distribution in the MisID-enriched region. The last bin contains the overflow events.

The $L_{T}+p^{miss}_{T}$ distribution in the 3L below-Z signal region. The last bin contains the overflow events.

The $M_{T}$ distribution in the 3L on-Z signal region. The last bin contains the overflow events.

The $L_{T}+p^{miss}_{T}$ distribution in the 3L above-Z signal region. The last bin contains the overflow events.

The $L_{T}+p^{miss}_{T}$ distribution in the 3L OSSF0 signal region. The last bin contains the overflow events.

The $L_{T}+p^{miss}_{T}$ distribution in the 4L OSSF0 signal region. The last bin contains the overflow events.

The $L_{T}+p^{miss}_{T}$ distribution in the 4L OSSF1 signal region. The last bin contains the overflow events.

The $L_{T}+p^{miss}_{T}$ distribution in the 4L OSSF2 signal region. The last bin contains the overflow events.

The dielectron $M_{OSSF}^{20}$ distribution in the 3L(ee) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 3L(ee) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 3L(ee) 0B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 3L(ee) 0B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 3L(ee) 0B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 3L(ee) 0B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 3L(ee) 1B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 3L(ee) 1B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 3L(ee) 1B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 3L(ee) 1B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 3L(ee) 1B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 3L(ee) 1B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 4L(ee) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 4L(ee) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 4L(ee) 0B, $S_{T}$>400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 4L(ee) 0B, $S_{T}$>400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 4L(ee) 1B signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 4L(ee) 1B signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 3L($\mu\mu$) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 3L($\mu\mu$) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 3L($\mu\mu$) 0B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 3L($\mu\mu$) 0B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 3L($\mu\mu$) 0B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 3L($\mu\mu$) 0B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 3L($\mu\mu$) 1B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 3L($\mu\mu$) 1B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 3L($\mu\mu$) 1B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 3L($\mu\mu$) 1B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 3L($\mu\mu$) 1B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 3L($\mu\mu$) 1B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 4L($\mu\mu$) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 4L($\mu\mu$) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 4L($\mu\mu$) 0B, $S_{T}$>400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 4L($\mu\mu$) 0B, $S_{T}$>400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 4L($\mu\mu$) 1B signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 4L($\mu\mu$) 1B signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The 95% confidence level exclusion limits for the flavor-democratic scenario on the total production cross section of heavy fermion pairs.

The 95% confidence level exclusion limits on the product of the total production cross section of tt$\phi$(S) and $\mathcal{B}(\phi - {ee})$.

The 95% confidence level exclusion limits on the product of the total production cross section of tt$\phi$(PS) and $\mathcal{B}(\phi - {ee})$.

The 95% confidence level exclusion limits on the product of the total production cross section of tt$\phi$(S) and $\mathcal{B}(\phi - {\mu\mu})$.

The 95% confidence level exclusion limits on the product of the total production cross section of tt$\phi$(PS) and $\mathcal{B}(\phi - {\mu\mu})$.

The 95% confidence level exclusion limits on $g_{t}^2\mathcal{B}(\phi - {ee})$ for tt$\phi$(S).

The 95% confidence level exclusion limits on $g_{t}^2\mathcal{B}(\phi - {ee})$ for tt$\phi$(PS).

The 95% confidence level exclusion limits on $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$ for tt$\phi$(S).

The 95% confidence level exclusion limits on $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$ for tt$\phi$(PS).

Product of the fiducial acceptance and the event selection efficiency for the type-III Seesaw signal model at various signal mass hypotheses calculated after all analysis selection requirements.

Product of the fiducial acceptance and the event selection efficiency for the tt$\phi$ signal models at various signal mass hypotheses calculated after all analysis selection requirements.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for bbA prodcution and the branching fraction $\mathrm{A} \rightarrow \tau\tau$, obtained from the combination of the $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels in the 1 b tag category, as a function of the A boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for bbA prodcution and the branching fraction $\mathrm{A} \rightarrow \tau\tau$, obtained in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1 b tag category, as a function of the A boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for bbA prodcution and the branching fraction $\mathrm{A} \rightarrow \tau\tau$, obtained in the $\mu\tau_\mathrm{h}$ channel of the 1 b tag category, as a function of the A boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

The product of acceptance, efficiency, and branching fraction of the qbX signal with $\mathrm{X} \rightarrow \tau\tau$ and $m_{\mathrm{B}}=170~\mathrm{GeV}$ in $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels of the 1b1f and 1b1c event categories, as a function of the X boson mass. The selections are as described in the paper. The uncertainty refers to the statistical uncertainty only.

Observed $m_{\tau\tau}$ distribution in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1b1f category, compared to the expected SM background contributions. The distributions of the qbX signal with $m_{\mathrm{B}}=170~\mathrm{GeV}$, and X boson mass 40 and 60 GeV are overlaid to illustrate the sensitivity. They are normalized to the cross section times branching fraction of 20 pb. The uncertainty band represents the sum in quadrature of statistical and systematic uncertainties obtained from the fit. The lower panel shows the ratio between the observed and expected events in each bin.

Observed $m_{\tau\tau}$ distribution in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1b1c category, compared to the expected SM background contributions. The distributions of the qbX signal with $m_{\mathrm{B}}=170~\mathrm{GeV}$, and X boson mass 40 and 60 GeV are overlaid to illustrate the sensitivity. They are normalized to the cross section times branching fraction of 20 pb. The uncertainty band represents the sum in quadrature of statistical and systematic uncertainties obtained from the fit. The lower panel shows the ratio between the observed and expected events in each bin.

Observed $m_{\tau\tau}$ distribution in the $\mu\tau_\mathrm{h}$ channel of the 1b1f category, compared to the expected SM background contributions. The distributions of the qbX signal with $m_{\mathrm{B}}=170~\mathrm{GeV}$, and X boson mass 40 and 60 GeV are overlaid to illustrate the sensitivity. They are normalized to the cross section times branching fraction of 20 pb. The uncertainty band represents the sum in quadrature of statistical and systematic uncertainties obtained from the fit. The lower panel shows the ratio between the observed and expected events in each bin.

Observed $m_{\tau\tau}$ distribution in the $\mu\tau_\mathrm{h}$ channel of the 1b1c category, compared to the expected SM background contributions. The distributions of the qbX signal with $m_{\mathrm{B}}=170~\mathrm{GeV}$, and X boson mass 40 and 60 GeV are overlaid to illustrate the sensitivity. They are normalized to the cross section times branching fraction of 20 pb. The uncertainty band represents the sum in quadrature of statistical and systematic uncertainties obtained from the fit. The lower panel shows the ratio between the observed and expected events in each bin.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained from the combination of the $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels, and the 1b1f and 1b1c categories, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 300 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained from the combination of the $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels, and the 1b1f and 1b1c categories, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 450 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained from the combination of the $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels, and the 1b1f and 1b1c categories, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1b1f category, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1b1c category, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained in the $\mu\tau_\mathrm{h}$ channel of the 1b1f category, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained in the $\mu\tau_\mathrm{h}$ channel of the 1b1c category, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Missing transverse momentum distribution after SR3b selection, beside the $E_\mathrm{T}^\mathrm{miss}$ requirement. The results in the signal region correspond to the last inclusive bin. The systematic uncertainties include theory uncertainties for the backgrounds with prompt SS/3L and the full systematic uncertainties for data-driven backgrounds. For illustration the distribution for a benchmark SUSY scenario ($pp\to \tilde g\tilde g$, $\tilde g\to t\bar t\tilde\chi_1^0$, $m_{\tilde g}=1.2$ TeV, $m_{\tilde\chi_1^0}=0.7$ TeV) is also shown.

Observed exclusion limits on the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to qq(\tilde\ell\ell/\tilde\nu\nu)$ decays. All limits are computed at 95% CL.

Expected exclusion limits on the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to qq(\tilde\ell\ell/\tilde\nu\nu)$ decays. All limits are computed at 95% CL.

Upper limits on signal cross-sections as function of the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to qq(\tilde\ell\ell/\tilde\nu\nu)$ decays, obtained using the signal efficiency and acceptance specific to each model. All limits are computed at 95% CL.

Observed exclusion limits on the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to qqWZ\tilde\chi_1^0$ decays. All limits are computed at 95% CL.

Expected exclusion limits on the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to qqWZ\tilde\chi_1^0$ decays. All limits are computed at 95% CL.

Upper limits on signal cross-sections as function of the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to qqWZ\tilde\chi_1^0$ decays, obtained using the signal efficiency and acceptance specific to each model. All limits are computed at 95% CL.

Observed exclusion limits on the $\tilde b_1$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde b_1\tilde b_1^*$ pair production with exclusive $\tilde b_1\to t\tilde\chi_1^-$ decays. All limits are computed at 95% CL.

Expected exclusion limits on the $\tilde b_1$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde b_1\tilde b_1^*$ pair production with exclusive $\tilde b_1\to t\tilde\chi_1^-$ decays. All limits are computed at 95% CL.

Upper limits on signal cross-sections as function of the $\tilde b_1$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde b_1\tilde b_1^*$ pair production with exclusive $\tilde b_1\to t\tilde\chi_1^-$ decays, obtained using the signal efficiency and acceptance specific to each model. All limits are computed at 95% CL.

Observed exclusion limits on the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to t\bar t\tilde\chi_1^0$ decays. All limits are computed at 95% CL.

Expected exclusion limits on the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to t\bar t\tilde\chi_1^0$ decays. All limits are computed at 95% CL.

Upper limits on signal cross-sections as function of the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to t\bar t\tilde\chi_1^0$ decays, obtained using the signal efficiency and acceptance specific to each model. All limits are computed at 95% CL.

SUSY scenario with $\tilde g\tilde g$ production and $\tilde g\to q\bar q(\tilde\ell\ell/\tilde\nu\nu)$ decay: signal acceptance (in %) in the signal region SR0b3j. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde g\tilde g$ production and $\tilde g\to q\bar q(\tilde\ell\ell/\tilde\nu\nu)$ decay: reconstruction efficiency (in %) in the signal region SR0b3j. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde g\tilde g$ production and $\tilde g\to q\bar qWZ\tilde\chi_1^0$ decay: signal acceptance (in %) in the signal region SR0b5j. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde g\tilde g$ production and $\tilde g\to q\bar qWZ\tilde\chi_1^0$ decay: reconstruction efficiency (in %) in the signal region SR0b5j. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde b_1\tilde b_1^*$ production and $\tilde b_1\to tW\tilde\chi_1^0$ decay: signal acceptance (in %) in the signal region SR1b. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde b_1\tilde b_1^*$ production and $\tilde b_1\to tW\tilde\chi_1^0$ decay: reconstruction efficiency (in %) in the signal region SR1b. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde g\tilde g$ production and $\tilde g\to t\bar t\tilde\chi_1^0$ decay: signal acceptance (in %) in the signal region SR3b. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde g\tilde g$ production and $\tilde g\to t\bar t\tilde\chi_1^0$ decay: reconstruction efficiency (in %) in the signal region SR3b. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

Observed 95% CL limit for the pMSSM grid.

Observed 95% CL limit for the pMSSM grid when the signal cross section is increased by one standard deviation.

Observed 95% CL limit for the pMSSM grid when the signal cross section is decreased by one standard deviation.

Expected 95% CL limit for the pMSSM grid.

+1 sigma excursion of the expected 95% CL limit for the pMSSM grid.

-1 sigma excursion of the expected 95% CL limit for the pMSSM grid.

Observed 95% CL limit for the 2Step grid.

Observed 95% CL limit for the 2Step grid when the signal cross section is increased by one standard deviation.

Observed 95% CL limit for the 2Step grid when the signal cross section is decreased by one standard deviation.

Expected 95% CL limit for the 2Step grid.

+1 sigma excursion of the expected 95% CL limit for the 2Step grid.

-1 sigma excursion of the expected 95% CL limit for the 2Step grid.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j50 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j50 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 9j50 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 9j50 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 9j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 7j80 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 7j80 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 7j80 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

Degree of multijet closure for signal and vaidation regions with at no b-jet requirement. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The bins labelled in bold are signal regions, while the others are validation regions. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Degree of multijet closure for signal and vaidation regions with at least 1 b-jet. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The bins labelled in bold are signal regions, while the others are validation regions. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Degree of multijet closure for signal and vaidation regions with at least 2 b-jets. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The bins labelled in bold are signal regions, while the others are validation regions. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Summary of all 15 signal regions (post-fit).

Signal region yielding the best-expected CLs value, the best expected CLs value, and the corresponding observed CLs value for the 2Step grid.

Signal region yielding the best-expected CLs value, the best expected CLs value, and the corresponding observed CLs value for the pMSSM grid.

95% CLs observed upper limit on model cross-section for 2-step signal points for the best-expected signal region.

Performance of the 8j50-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.