Selection acceptance times efficiency for VBF H -> ZZ -> vvqq as a function of the Higgs boson mass, combining the HP and LP signal regions.

Selection acceptance times efficiency for DY W' -> ZW -> llqq as a function of the W' mass, combining the HP and LP signal regions of the ZV -> llJ selection and the b-tagged and untagged regions of the ZV -> lljj selection.

Selection acceptance times efficiency for VBF W' -> ZW -> llqq as a function of the W' mass, combining the HP and LP signal regions of the ZV -> llJ selection and the b-tagged and untagged regions of the ZV -> lljj selection.

Selection acceptance times efficiency for DY W' -> ZW -> vvqq as a function of the W' mass, combining the HP and LP signal regions.

Selection acceptance times efficiency for VBF W' -> ZW -> vvqq as a function of the W' mass, combining the HP and LP signal regions.

Selection acceptance times efficiency for G* -> ZZ -> llqq as a function of the W' mass, combining the HP and LP signal regions of the ZV -> llJ selection and the b-tagged and untagged regions of the ZV -> lljj selection.

Selection acceptance times efficiency for G* -> ZZ -> vvqq as a function of the W' mass, combining the HP and LP signal regions.

Observed and expected 95% CL upper limits on $\sigma\times$BR($H\to ZZ$) at $\sqrt{s} = 13$ TeV for ggF production of a heavy Higgs boson as a function of its mass, combining $\ell\ell qq$ and $\nu\nu qq$ searches.

Observed and expected 95% CL upper limits on $\sigma\times$BR($H\to ZZ$) at $\sqrt{s} = 13$ TeV for VBF production of a heavy Higgs boson as a function of its mass, combining $\ell\ell qq$ and $\nu\nu qq$ searches.

Observed and expected 95% CL upper limits on $\sigma\times$BR($W'\to ZW$) at $\sqrt{s} = 13$ TeV for Drell-Yan production of a $W'$ boson in the HVT model as a function of its mass, combining $\ell\ell qq$ and $\nu\nu qq$ searches.

Observed and expected 95% CL upper limits on $\sigma\times$BR($W'\to ZW$) at $\sqrt{s} = 13$ TeV for VBF production of a $W'$ boson in the HVT model as a function of its mass, combining $\ell\ell qq$ and $\nu\nu qq$ searches.

Observed and expected 95% CL upper limits on $\sigma\times$BR($G_{\mathrm{KK}}\to ZZ$) at $\sqrt{s} = 13$ TeV for the production of a bulk RS graviton with couplings of $k/\overline{M}_{\mathrm{Pl}}=1$.

Observed and expected 95% CL upper limits on $\sigma\times$BR($G_{\mathrm{KK}}\to ZZ$) at $\sqrt{s} = 13$ TeV for the production of a bulk RS graviton with couplings of $k/\overline{M}_{\mathrm{Pl}}=0.5$.

Fit results of the simultaneous measurements of the $H\to e\tau$ and $H\to \mu\tau$ signals (2POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to e\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the simultaneous measurements of the $H\to e\tau$ and $H\to \mu\tau$ signals (2POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to \mu\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the simultaneous measurements of the $H\to e\tau$ and $H\to \mu\tau$ signals (2POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to \mu\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the simultaneous measurements of the $H\to e\tau$ and $H\to \mu\tau$ signals (2POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to \mu\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the simultaneous measurements of the $H\to e\tau$ and $H\to \mu\tau$ signals (2POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to \mu\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to e\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to e\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to e\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to e\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to \mu\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to \mu\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to \mu\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to \mu\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $e\tau_{\mu}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $e\tau_{\mu}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $\mu\tau_{e}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $\mu\tau_{e}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $e\tau_{\mu}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $e\tau_{\mu}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $\mu\tau_{e}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $\mu\tau_{e}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $e\tau_{had}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $e\tau_{had}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $\mu\tau_{had}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $\mu\tau_{had}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $e\tau_{had}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $e\tau_{had}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $\mu\tau_{had}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $\mu\tau_{had}$ final state in the VBF category.

NN score distribution, after an independent fit of the $H\to \mu\tau$ signal, obtained by fitting the data of the Symmetry $\ell\tau_{\ell'}$ channel, for the $\mu\tau_{e}$ final state in the non-VBF category.

NN score distribution, after an independent fit of the $H\to \mu\tau$ signal, obtained by fitting the data of the Symmetry $\ell\tau_{\ell'}$ channel, for the $\mu\tau_{e}$ final state in the non-VBF category.

NN score distribution, after an independent fit of the $H\to \mu\tau$ signal, obtained by fitting the data of the Symmetry $\ell\tau_{\ell'}$ channel, for the $\mu\tau_{e}$ final state in the VBF category.

NN score distribution, after an independent fit of the $H\to \mu\tau$ signal, obtained by fitting the data of the Symmetry $\ell\tau_{\ell'}$ channel, for the $\mu\tau_{e}$ final state in the VBF category.

Best-fit $\mathcal{B}(H\to \mu\tau)- \mathcal{B}(H\to e\tau)$ values, given in %, obtained from the standalone fits of the $\ell\tau_\ell'$ final state with either background estimation method. For the MC-template method, the correlated uncertainties are fixed to their best-fit values and the uncertainties are re-derived considering only the uncorrelated sources. For the Symmetry method, the full uncertainty is considered.

Best-fit $\mathcal{B}(H\to \mu\tau)- \mathcal{B}(H\to e\tau)$ values, given in %, obtained from the standalone fits of the $\ell\tau_\ell'$ final state with either background estimation method. For the MC-template method, the correlated uncertainties are fixed to their best-fit values and the uncertainties are re-derived considering only the uncorrelated sources. For the Symmetry method, the full uncertainty is considered.

Best-fit value of the branching ratios $\mathcal{B}(H\to e\tau)$ and $\mathcal{B}(H\to \mu\tau)$, given in %, and likelihood contours at 68% and 95% C.L. Obtained from the simultaneous fit of $H\to e\tau$ and $H\to \mu\tau$ signals based on the MC-template method, compared to the SM expectation.

Best-fit value of the branching ratios $\mathcal{B}(H\to e\tau)$ and $\mathcal{B}(H\to \mu\tau)$, given in %, and likelihood contours at 68% and 95% C.L. Obtained from the simultaneous fit of $H\to e\tau$ and $H\to \mu\tau$ signals based on the MC-template method, compared to the SM expectation.

Breakdown of expected and observed yields in the electroweak search Int, Low, and OffShell signal regions after a simultaneous fit to the signal regions and control regions. All statistical and systematic uncertainties are included.

Breakdown of expected and observed yields in the four edge signal regions, integrated over the $m_{\ell\ell}$ distribution after a separate simultaneous fit to each signal region and control region pair. The uncertainties include both the statistical and systematic sources.

Breakdown of expected and observed yields in the three on-$Z$ signal regions after a separate simultaneous fit to each signal region and control region pair. The uncertainties include both the statistical and systematic sources.

Distributions of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in VR-High-Sideband-EWK (top-left), VR-High-R-EWK (top-right), VR-1J-High-EWK (bottom-left), and VR-$\ell\ell bb$-EWK (bottom-right) from the EWK search after a simultaneous fit of the control regions. The hatched band includes both the systematic and statistical uncertainties. The last bin includes the overflow.

Distributions of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in VR-High-Sideband-EWK (top-left), VR-High-R-EWK (top-right), VR-1J-High-EWK (bottom-left), and VR-$\ell\ell bb$-EWK (bottom-right) from the EWK search after a simultaneous fit of the control regions. The hatched band includes both the systematic and statistical uncertainties. The last bin includes the overflow.

Distributions of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in VR-High-Sideband-EWK (top-left), VR-High-R-EWK (top-right), VR-1J-High-EWK (bottom-left), and VR-$\ell\ell bb$-EWK (bottom-right) from the EWK search after a simultaneous fit of the control regions. The hatched band includes both the systematic and statistical uncertainties. The last bin includes the overflow.

Distributions of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in VR-High-Sideband-EWK (top-left), VR-High-R-EWK (top-right), VR-1J-High-EWK (bottom-left), and VR-$\ell\ell bb$-EWK (bottom-right) from the EWK search after a simultaneous fit of the control regions. The hatched band includes both the systematic and statistical uncertainties. The last bin includes the overflow.

Distributions of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in VR-Int-EWK (top-left), VR-Low-EWK (top-right), VR-Low-2-EWK (bottom-left), and VR-OffShell-EWK (bottom-right) from the EWK search after a simultaneous fit of the control regions. The hatched band includes both the systematic and statistical uncertainties. The last bin includes the overflow.

Distributions of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in VR-Int-EWK (top-left), VR-Low-EWK (top-right), VR-Low-2-EWK (bottom-left), and VR-OffShell-EWK (bottom-right) from the EWK search after a simultaneous fit of the control regions. The hatched band includes both the systematic and statistical uncertainties. The last bin includes the overflow.

Distributions of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in VR-Int-EWK (top-left), VR-Low-EWK (top-right), VR-Low-2-EWK (bottom-left), and VR-OffShell-EWK (bottom-right) from the EWK search after a simultaneous fit of the control regions. The hatched band includes both the systematic and statistical uncertainties. The last bin includes the overflow.

Distributions of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in VR-Int-EWK (top-left), VR-Low-EWK (top-right), VR-Low-2-EWK (bottom-left), and VR-OffShell-EWK (bottom-right) from the EWK search after a simultaneous fit of the control regions. The hatched band includes both the systematic and statistical uncertainties. The last bin includes the overflow.

Observed and expected dilepton mass distributions in VRC-STR (top-left), VRLow-STR (top-right), VRMed-STR (bottom-left), and VRHigh-STR (bottom-right). Each validation region is fit separately with the corresponding control region. All statistical and systematic uncertainties are included in the hatched band. The entries are normalized to the bin width, and the last bin is the overflow.

Observed and expected dilepton mass distributions in VRC-STR (top-left), VRLow-STR (top-right), VRMed-STR (bottom-left), and VRHigh-STR (bottom-right). Each validation region is fit separately with the corresponding control region. All statistical and systematic uncertainties are included in the hatched band. The entries are normalized to the bin width, and the last bin is the overflow.

Observed and expected dilepton mass distributions in VRC-STR (top-left), VRLow-STR (top-right), VRMed-STR (bottom-left), and VRHigh-STR (bottom-right). Each validation region is fit separately with the corresponding control region. All statistical and systematic uncertainties are included in the hatched band. The entries are normalized to the bin width, and the last bin is the overflow.

Observed and expected dilepton mass distributions in VRC-STR (top-left), VRLow-STR (top-right), VRMed-STR (bottom-left), and VRHigh-STR (bottom-right). Each validation region is fit separately with the corresponding control region. All statistical and systematic uncertainties are included in the hatched band. The entries are normalized to the bin width, and the last bin is the overflow.

Observed and expected jet multiplicity in VRLow-STR (top-left), VRMed-STR (top-right), and VRHigh-STR (bottom) after a fit performed on the $m_{\ell\ell}$ distribution and corresponding control region. All statistical and systematic uncertainties are included in the hatched band. The last bin contains the overflow.

Observed and expected jet multiplicity in VRLow-STR (top-left), VRMed-STR (top-right), and VRHigh-STR (bottom) after a fit performed on the $m_{\ell\ell}$ distribution and corresponding control region. All statistical and systematic uncertainties are included in the hatched band. The last bin contains the overflow.

Observed and expected jet multiplicity in VRLow-STR (top-left), VRMed-STR (top-right), and VRHigh-STR (bottom) after a fit performed on the $m_{\ell\ell}$ distribution and corresponding control region. All statistical and systematic uncertainties are included in the hatched band. The last bin contains the overflow.

Observed and expected dilepton mass distributions in VR3L-STR without a fit to the data. The 'Other' category includes the negligible contributions from $t\bar{t}$ and $Z/\gamma^*$+jets processes. The hatched band contains the statistical uncertainty and the theoretical systematic uncertainties of the $WZ$/$ZZ$ prediction, which are the dominant sources of uncertainty. No fit is performed. The last bin contains the overflow.

Observed and expected distributions in five EWK search regions after a simultaneous fit to the SR and CR. In the top row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-High_8-EWK and $m_{bb}$ in SR-$\ell\ell bb$-EWK. In the middle row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Int-EWK and $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Low-EWK. In the bottom row is $m_{\ell\ell}$ in SR-OffShell-EWK. Overlaid are example C1N2 and GMSB signal models, where the numbers in the brackets indicate the masses, in $\mathrm{GeV}$, of the $\tilde{\chi}_1^\pm$ and $\tilde{\chi}_2^0$ or the mass of the $\tilde{\chi}_1^0$ and branching ratio to the Higgs boson respectively. All statistical and systematic uncertainties are included in the hatched bands. The last bin includes the overflow.

Observed and expected distributions in five EWK search regions after a simultaneous fit to the SR and CR. In the top row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-High_8-EWK and $m_{bb}$ in SR-$\ell\ell bb$-EWK. In the middle row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Int-EWK and $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Low-EWK. In the bottom row is $m_{\ell\ell}$ in SR-OffShell-EWK. Overlaid are example C1N2 and GMSB signal models, where the numbers in the brackets indicate the masses, in $\mathrm{GeV}$, of the $\tilde{\chi}_1^\pm$ and $\tilde{\chi}_2^0$ or the mass of the $\tilde{\chi}_1^0$ and branching ratio to the Higgs boson respectively. All statistical and systematic uncertainties are included in the hatched bands. The last bin includes the overflow.

Observed and expected distributions in five EWK search regions after a simultaneous fit to the SR and CR. In the top row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-High_8-EWK and $m_{bb}$ in SR-$\ell\ell bb$-EWK. In the middle row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Int-EWK and $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Low-EWK. In the bottom row is $m_{\ell\ell}$ in SR-OffShell-EWK. Overlaid are example C1N2 and GMSB signal models, where the numbers in the brackets indicate the masses, in $\mathrm{GeV}$, of the $\tilde{\chi}_1^\pm$ and $\tilde{\chi}_2^0$ or the mass of the $\tilde{\chi}_1^0$ and branching ratio to the Higgs boson respectively. All statistical and systematic uncertainties are included in the hatched bands. The last bin includes the overflow.

Observed and expected distributions in five EWK search regions after a simultaneous fit to the SR and CR. In the top row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-High_8-EWK and $m_{bb}$ in SR-$\ell\ell bb$-EWK. In the middle row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Int-EWK and $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Low-EWK. In the bottom row is $m_{\ell\ell}$ in SR-OffShell-EWK. Overlaid are example C1N2 and GMSB signal models, where the numbers in the brackets indicate the masses, in $\mathrm{GeV}$, of the $\tilde{\chi}_1^\pm$ and $\tilde{\chi}_2^0$ or the mass of the $\tilde{\chi}_1^0$ and branching ratio to the Higgs boson respectively. All statistical and systematic uncertainties are included in the hatched bands. The last bin includes the overflow.

Observed and expected distributions in five EWK search regions after a simultaneous fit to the SR and CR. In the top row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-High_8-EWK and $m_{bb}$ in SR-$\ell\ell bb$-EWK. In the middle row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Int-EWK and $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Low-EWK. In the bottom row is $m_{\ell\ell}$ in SR-OffShell-EWK. Overlaid are example C1N2 and GMSB signal models, where the numbers in the brackets indicate the masses, in $\mathrm{GeV}$, of the $\tilde{\chi}_1^\pm$ and $\tilde{\chi}_2^0$ or the mass of the $\tilde{\chi}_1^0$ and branching ratio to the Higgs boson respectively. All statistical and systematic uncertainties are included in the hatched bands. The last bin includes the overflow.

Observed and expected dilepton mass distributions in SRC-STR (top-left), SRLow-STR (top-right), SRMed-STR (bottom-left), and SRHigh-STR (bottom-right), with the binning used for interpretations after a separate simultaneous fit to each signal region and control region pair. The red dashed lines are example signal models overlaid on the figure. All statistical and systematic uncertainties are included in the hatched bands. The last bins are the overflow.

Observed and expected dilepton mass distributions in SRC-STR (top-left), SRLow-STR (top-right), SRMed-STR (bottom-left), and SRHigh-STR (bottom-right), with the binning used for interpretations after a separate simultaneous fit to each signal region and control region pair. The red dashed lines are example signal models overlaid on the figure. All statistical and systematic uncertainties are included in the hatched bands. The last bins are the overflow.

Observed and expected dilepton mass distributions in SRC-STR (top-left), SRLow-STR (top-right), SRMed-STR (bottom-left), and SRHigh-STR (bottom-right), with the binning used for interpretations after a separate simultaneous fit to each signal region and control region pair. The red dashed lines are example signal models overlaid on the figure. All statistical and systematic uncertainties are included in the hatched bands. The last bins are the overflow.

Observed and expected dilepton mass distributions in SRC-STR (top-left), SRLow-STR (top-right), SRMed-STR (bottom-left), and SRHigh-STR (bottom-right), with the binning used for interpretations after a separate simultaneous fit to each signal region and control region pair. The red dashed lines are example signal models overlaid on the figure. All statistical and systematic uncertainties are included in the hatched bands. The last bins are the overflow.

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294]. The grey numbers indicate the observed 95\% CLs upper limit on the cross section.

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294]. The grey numbers indicate the observed 95$\%$ CLs upper limit on the cross section.

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$ ilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$ ilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$ ilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$ ilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].The grey numbers indicated the observed 95\% CL upper limit on the cross section.

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].The grey numbers indicated the observed 95\% CL upper limit on the cross section.

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].The grey numbers indicated the observed 95\% CL upper limit on the cross section.

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

The combined $E_{\mathrm{T}}^{\mathrm{miss}}$ distribution of VRLow-STR and SRLow-STR (left), and the same region with the $\Delta\phi(\boldsymbol{j}_{1,2},\boldsymbol{\mathit{p}}_{ ext{T}}^{ ext{miss}})<0.4$ requirement, used as a control region to normalize the $Z/\gamma^*+\mathrm{jets}$ process (right). Separate fits for the SR and VR are performed, as for the results in the paper, and the resulting distributions are merged. All statistical and systematic uncertainties are included in the hatched bands. The last bins contain the overflow.

Cutflow of expected events in the four Strong search edge signal regions. `Leptons' refers to electrons and muons only. The gluino-$Z^{(*)}$ model with $m_{ ilde{g}}=800~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for SRC-STR with 60000 Monte Carlo (MC) events generated. The slepton-$Z^{(*)}$ model with $m_{ ilde{\ell}}=1200~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for SRMed-STR with 30000 (MC) events generated. The gluino-slepton model with $m_{ ilde{g}}=2~TeV$ and $m_{ ilde{\ell}}=1.3~TeV$ is used for SRLow-STR and SRHigh-STR with 30000 MC events generated. The Generator Filter requires two 5~GeV leptons and 100~GeV of \met. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~GeV$ or at least one lepton with $p_{\mathrm{T}}>25~GeV$ and a photon with $p_{\mathrm{T}}>40~GeV$, with all objects within $|\eta|=2.6$.

Cutflow of expected events in the four Strong search edge signal regions. `Leptons' refers to electrons and muons only. The gluino-$Z^{(*)}$ model with $m_{ ilde{g}}=800~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for SRC-STR with 60000 Monte Carlo (MC) events generated. The slepton-$Z^{(*)}$ model with $m_{ ilde{\ell}}=1200~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for SRMed-STR with 30000 (MC) events generated. The gluino-slepton model with $m_{ ilde{g}}=2~TeV$ and $m_{ ilde{\ell}}=1.3~TeV$ is used for SRLow-STR and SRHigh-STR with 30000 MC events generated. The Generator Filter requires two 5~GeV leptons and 100~GeV of \met. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~GeV$ or at least one lepton with $p_{\mathrm{T}}>25~GeV$ and a photon with $p_{\mathrm{T}}>40~GeV$, with all objects within $|\eta|=2.6$.

Cutflow of expected events in the four Strong search edge signal regions. `Leptons' refers to electrons and muons only. The gluino-$Z^{(*)}$ model with $m_{ ilde{g}}=800~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for SRC-STR with 60000 Monte Carlo (MC) events generated. The slepton-$Z^{(*)}$ model with $m_{ ilde{\ell}}=1200~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for SRMed-STR with 30000 (MC) events generated. The gluino-slepton model with $m_{ ilde{g}}=2~TeV$ and $m_{ ilde{\ell}}=1.3~TeV$ is used for SRLow-STR and SRHigh-STR with 30000 MC events generated. The Generator Filter requires two 5~GeV leptons and 100~GeV of \met. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~GeV$ or at least one lepton with $p_{\mathrm{T}}>25~GeV$ and a photon with $p_{\mathrm{T}}>40~GeV$, with all objects within $|\eta|=2.6$.

Cutflow of expected events in the four Strong search edge signal regions. `Leptons' refers to electrons and muons only. The gluino-$Z^{(*)}$ model with $m_{ ilde{g}}=800~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for SRC-STR with 60000 Monte Carlo (MC) events generated. The slepton-$Z^{(*)}$ model with $m_{ ilde{\ell}}=1200~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for SRMed-STR with 30000 (MC) events generated. The gluino-slepton model with $m_{ ilde{g}}=2~TeV$ and $m_{ ilde{\ell}}=1.3~TeV$ is used for SRLow-STR and SRHigh-STR with 30000 MC events generated. The Generator Filter requires two 5~GeV leptons and 100~GeV of \met. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~GeV$ or at least one lepton with $p_{\mathrm{T}}>25~GeV$ and a photon with $p_{\mathrm{T}}>40~GeV$, with all objects within $|\eta|=2.6$.

Cutflow of expected events in the three Strong search on-$Z$ signal regions. The cutflow up to the signal region specific requirements is the same as in the Strong search edge cutflow. The slepton-$Z^{(*)}$ model with $m_{ ilde{\ell}}=1200~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for all of the on-$Z$ signal regions with 30000 (MC) events generated.

Cutflow of expected events in the three Strong search on-$Z$ signal regions. The cutflow up to the signal region specific requirements is the same as in the Strong search edge cutflow. The slepton-$Z^{(*)}$ model with $m_{ ilde{\ell}}=1200~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for all of the on-$Z$ signal regions with 30000 (MC) events generated.

Cutflow of expected events in the three Strong search on-$Z$ signal regions. The cutflow up to the signal region specific requirements is the same as in the Strong search edge cutflow. The slepton-$Z^{(*)}$ model with $m_{ ilde{\ell}}=1200~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for all of the on-$Z$ signal regions with 30000 (MC) events generated.

Table 36: Cutflow of expected events in the region SR-OffShell_a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 36: Cutflow of expected events in the region SR-OffShell_a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 37: Cutflow of expected events in the region SR-OffShell_b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 37: Cutflow of expected events in the region SR-OffShell_b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 38: Cutflow of expected events in the region SR-Low_a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 38: Cutflow of expected events in the region SR-Low_a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 39: Cutflow of expected events in the region SR-Low_b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 39: Cutflow of expected events in the region SR-Low_b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 40: Cutflow of expected events in the region SR-Low-2-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 40: Cutflow of expected events in the region SR-Low-2-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 41: Cutflow of expected events in the region SR-Int_a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 41: Cutflow of expected events in the region SR-Int_a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 42: Cutflow of expected events in the region SR-Int_b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 42: Cutflow of expected events in the region SR-Int_b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 43: Cutflow of expected events in the region SR-High_16a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 43: Cutflow of expected events in the region SR-High_16a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 44: Cutflow of expected events in the region SR-High_16b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 44: Cutflow of expected events in the region SR-High_16b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 45: Cutflow of expected events in the region SR-High_8a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 45: Cutflow of expected events in the region SR-High_8a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 46: Cutflow of expected events in the region SR-High_8b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 46: Cutflow of expected events in the region SR-High_8b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 47: Cutflow of expected events in the region SR-1J-High-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 47: Cutflow of expected events in the region SR-1J-High-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 48: Cutflow of expected events in the region SR-llbb-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 48: Cutflow of expected events in the region SR-llbb-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

The combined $m_{jj}$ distribution of CR-Z-EWK and SR-Low-EWK (left), and the $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ distribution in CR-Z-met-EWK (right), which removes the upper limit of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}}) < 9$ from the definition of CR-Z-EWK. This $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ tail overlaps other control and validation regions, but not signal regions. The arrows indicate the signal region SR-Low-EWK (left), and the $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ phase space which is not included in CR-Z-EWK (right). All EWK search control and signal regions are included in the fit. All statistical and systematic uncertainties are included in the hatched bands. The theoretical uncertainties from CR-Z-EWK are applied to CR-Z-met-EWK. The last bins contain the overflow.

The combined $m_{jj}$ distribution of CR-Z-EWK and SR-Low-EWK (left), and the $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ distribution in CR-Z-met-EWK (right), which removes the upper limit of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}}) < 9$ from the definition of CR-Z-EWK. This $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ tail overlaps other control and validation regions, but not signal regions. The arrows indicate the signal region SR-Low-EWK (left), and the $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ phase space which is not included in CR-Z-EWK (right). All EWK search control and signal regions are included in the fit. All statistical and systematic uncertainties are included in the hatched bands. The theoretical uncertainties from CR-Z-EWK are applied to CR-Z-met-EWK. The last bins contain the overflow.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-1J-High-EWK and SR-$\ell\ell bb$-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-1J-High-EWK and SR-$\ell\ell bb$-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-1J-High-EWK and SR-$\ell\ell bb$-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-1J-High-EWK and SR-$\ell\ell bb$-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-1J-High-EWK and SR-$\ell\ell bb$-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-1J-High-EWK and SR-$\ell\ell bb$-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-1J-High-EWK and SR-$\ell\ell bb$-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-1J-High-EWK and SR-$\ell\ell bb$-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_SLN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_SLN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_SLN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_SLN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_SLN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_SLN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_SLN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_SLN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Differential cross-section and uncertainty from each source, as a function of leading $p_{T}^j$ for the dominant process in the $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of leading $p_{T}^j$ for the dominant process in the extreme $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of leading $p_{T}^j$ for the dominant process in the extreme $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of leading $p_{T}^j$ for the dominant process in the extreme $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of subleading $p_{T}^j$ for the dominant process in the $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of subleading $p_{T}^j$ for the dominant process in the $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of subleading $p_{T}^j$ for the dominant process in the $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of subleading $p_{T}^j$ for the dominant process in the extreme $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of subleading $p_{T}^j$ for the dominant process in the extreme $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of subleading $p_{T}^j$ for the dominant process in the extreme $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $min\Delta\phi(j_0,l)$ for the dominant process in the $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $min\Delta\phi(j_0,l)$ for the dominant process in the $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $min\Delta\phi(j_0,l)$ for the dominant process in the $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $min\Delta\phi(j_0,l)$ for the dominant process in the extreme $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $min\Delta\phi(j_0,l)$ for the dominant process in the extreme $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $min\Delta\phi(j_0,l)$ for the dominant process in the extreme $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $min\Delta\phi(j_1,l)$ for the dominant process in the $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $min\Delta\phi(j_1,l)$ for the dominant process in the $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $min\Delta\phi(j_1,l)$ for the dominant process in the $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $min\Delta\phi(j_1,l)$ for the dominant process in the extreme $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $min\Delta\phi(j_1,l)$ for the dominant process in the extreme $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $min\Delta\phi(j_1,l)$ for the dominant process in the extreme $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $\Delta\eta_{jj}$ for the dominant process in the $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $\Delta\eta_{jj}$ for the dominant process in the $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $\Delta\eta_{jj}$ for the dominant process in the $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $\Delta\eta_{jj}$ for the dominant process in the extreme $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $\Delta\eta_{jj}$ for the dominant process in the extreme $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $\Delta\eta_{jj}$ for the dominant process in the extreme $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $\Delta\phi_{jj}$ for the dominant process in the $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $\Delta\phi_{jj}$ for the dominant process in the $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $\Delta\phi_{jj}$ for the dominant process in the $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $\Delta\phi_{jj}$ for the dominant process in the extreme $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $\Delta\phi_{jj}$ for the dominant process in the extreme $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $\Delta\phi_{jj}$ for the dominant process in the extreme $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $\Delta\phi_{ll}$ for the dominant process in the $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $\Delta\phi_{ll}$ for the dominant process in the $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $\Delta\phi_{ll}$ for the dominant process in the $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $\Delta\phi_{ll}$ for the dominant process in the extreme $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $\Delta\phi_{ll}$ for the dominant process in the extreme $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $\Delta\phi_{ll}$ for the dominant process in the extreme $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $m_{jj}$ for the dominant process in the $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $m_{jj}$ for the dominant process in the $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $m_{jj}$ for the dominant process in the $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $m_{jj}$ for the dominant process in the extreme $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $m_{jj}$ for the dominant process in the extreme $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $m_{jj}$ for the dominant process in the extreme $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $p_{T}^{ee}$ for the dominant process in the $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $p_{T}^{\mu\mu}$ for the dominant process in the $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $p_{T}^{e\mu}$ for the dominant process in the $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $p_{T}^{ee}$ for the dominant process in the extreme $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $p_{T}^{\mu\mu}$ for the dominant process in the extreme $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $p_{T}^{e\mu}$ for the dominant process in the extreme $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $H_{T}$ for the dominant process in the $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $H_{T}$ for the dominant process in the $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $H_{T}$ for the dominant process in the $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $H_{T}$ for the dominant process in the extreme $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $H_{T}$ for the dominant process in the extreme $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $H_{T}$ for the dominant process in the extreme $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $S_{T}$ for the dominant process in the $ee jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $S_{T}$ for the dominant process in the $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $S_{T}$ for the dominant process in the $e\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $S_{T}$ for the dominant process in the extreme $eejj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $S_{T}$ for the dominant process in the extreme $\mu\mu jj$ measurement region.

Differential cross-section and uncertainty from each source, as a function of $S_{T}$ for the dominant process in the extreme $e\mu jj$ measurement region.

Expected and observed 95% CL lower limits on first- and second-generation leptoquark masses for different values of $\beta$.

Event yields in the dimuon channel control regions with total uncertainties. The observed number of events is given in the first row. The background event numbers as obtained from the fit are shown together with the total uncertainties. The second row shows the total background expectation, the further rows show the breakdown into different background components.

Event yields in the dielectron channel control regions with total uncertainties. The observed number of events is given in the first row. The background event numbers as obtained from the fit are shown together with the total uncertainties. The second row shows the total background expectation, the further rows show the breakdown into different background components.

Distribution of $m_{LQ}^{min}$ in the training region for the BDT for the $ee jj$ and $\mu\mu jj$ channels. Data are shown together with predicted total background expectation.

Distribution of $m_{LQ}^{T}$ in the training region for the BDT for the $e\nu jj$ and $\mu\nu jj$ channels. Data are shown together with predicted total background expectation.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$, when all other coupling modifiers are fixed to their SM predictions.

Likelihood observed contours at 68% CL (solid line) and 95% CL (dashed line) in the $\kappa_{\lambda}$, $\kappa_{2V}$ parameter space, when all other coupling modifiers are fixed to their SM predictions. The best-fit value, denoted by a cross in the plot, corresponds to the very first entry in the table.

Likelihood expected contours at 68% CL (teal shaded region) and 95% CL (yellow shaded region) in the $\kappa_{\lambda}$, $\kappa_{2V}$ parameter space, when all other coupling modifiers are fixed to their SM predictions.In the plot, the SM prediction is indicated by the star.

Likelihood observed contours at 68% CL (solid line) and 95% CL (dashed line) in the $c_{gghh}$ versus $c_{hhh}$ HEFT parameter space, with the remaining coefficient fixed to its SM value. The best-fit value, denoted by a cross in the plot, corresponds to the very first entry in the table.

Likelihood expected contours at 68% CL (teal shaded region) 95% CL (yellow shaded region) in the $c_{gghh}$ versus $c_{hhh}$ HEFT parameter space, with the remaining coefficient fixed to its SM value. In the plot, the SM prediction is indicated by the star.

Likelihood observed contours at 68% CL (solid line) and 95% CL (dashed line) in the $c_{tthh}$ versus $c_{hhh}$ HEFT parameter space, with the remaining coefficient fixed to its SM value. The best-fit value, denoted by a cross in the plot, corresponds to the very first entry in the table.

Likelihood exptected contours at 68% CL (teal shaded region) and 95% CL (yellow shaded region) in the $c_{tthh}$ versus $c_{hhh}$ HEFT parameter space, with the remaining coefficient fixed to its SM value. In the plot, the SM prediction is indicated by the star.

The observed (filled circles) and expected (hollow circles) 95% CL upper limits on the $HH$ ggF production cross-section in the Standard Model and for seven HEFT benchmark points. The teal and yellow colored bands indicate the $\pm 1\sigma$ and $\pm 2\sigma$ variations on the expected limit due to statistical and systematic uncertainties. The predicted cross-sections of each of the models under consideration are shown by the red crosses. The contribution from VBF production to the total $HH$ production cross-section is neglected.

Likelihood observed contours at 68% CL (solid line) and 95% CL (dashed line) in the $c_{{H}\boxed{}}$ versus $c_{H}$ SMEFT parameter space, with the remaining coefficient fixed to its SM value. The best-fit value, denoted by a cross in the plot, corresponds to the very first entry in the table.

Likelihood expected contours at 68% CL (teal shaded region) 95% CL (yellow shaded region) in the $c_{{H}\boxed{}}$ versus $c_{H}$ SMEFT parameter space, with the remaining coefficient fixed to its SM value. In the plot, the SM prediction is indicated by the star.

The acceptance times efficiency for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions. The bands indicate the simulation statistical uncertainty.

The acceptance times efficiency for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category. When varying $\kappa_{\lambda}$, the other coupling modifiers affecting the VBF $HH$ production mode are set to their SM predictions. The bands indicate the simulation statistical uncertainty.

The acceptance times efficiency for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category. When varying $\kappa_{2V}$, the other coupling modifiers affecting the VBF $HH$ production mode are set to their SM predictions. The bands indicate the simulation statistical uncertainty.

Expected yields of the signal $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ after applying the analysis preselection and in each analysis category after the final selection. The $HH$ yield is obtained considering both the $ggF$ and $VBF$ contributions. The other coupling modifiers affecting $HH$ production are fixed to their SM predictions. The bands indicate the simulation statistical uncertainty.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{hhh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{tthh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{gghh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{H}$ SMEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{{H}\boxed{}}$ SMEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{hhh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{tthh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{gghh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{H}$ SMEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{{H}\boxed{}}$ SMEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the low mass 1 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the low mass 2 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the low mass 3 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the low mass 4 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the high mass 1 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the high mass 2 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the high mass 3 region.

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{hhh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{hhh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{gghh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{gghh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{tthh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{tthh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{H}$ SMEFT coefficients when all other coupling modifiers are fixed to the SM value.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{H}$ SMEFT coefficients when all other coupling modifiers are fixed to the SM value.

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{{H}\boxed{}}$ SMEFT coefficients when all other coupling modifiers are fixed to the SM value.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{{H}\boxed{}}$ SMEFT coefficients when all other coupling modifiers are fixed to the SM value.

Summary of the expected (blue) and observed (orange) one-dimensional constraints at 68% (solid error bars) and 95% (dashed error bars) CL on the coefficients of selected HEFT (top) and SMEFT (bottom) operators describing Higgs boson interactions affecting $HH$ production through gluon--gluon fusion. The results are obtained from one dimensional scans of the profile log-likelihood assuming that all other Wilson coefficients are fixed to their SM values. The contribution from VBF $HH$ production is neglected.

Summary of the total efficiencies. Quoted uncertainties correspond to statistical uncertainties of the simulated samples used. For the $B_c^+\to J/\psi \pi^+$ channel, the efficiency $\epsilon_{B_c^+\to J/\psi \pi^+}$ entering the equations for $R_{D_s^{(*)+}/\pi^+}$ is shown, while the efficiency for the full dataset is not defined.

Post-fit BDT distributions after the VBF event selection for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. The VBF signal is shown for $\mu = 0.73$ and $\tilde d = -0.01$. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given.

The observed $\Delta\mathrm{NLL}$ curve as a function of $\tilde d$ values. For comparison, expected $\Delta\mathrm{NLL}$ curves are also shown. The constraints on the nuisance parameters and normalization factors are first determined in a CR-only fit, and then a fit to pseudo-data corresponding to these nuisance parameters, normalization factors, and to $\tilde d=0, \mu = 1$ or $\tilde d =0, \mu = 0.73$ is performed to obtain these $\Delta\mathrm{NLL}$ curves. Moreover, a pre-fit expected $\Delta\mathrm{NLL}$ is shown, using pseudo-data corresponding to $\tilde d =0$ and $\mu = 1$ in the signal and control regions.

The expected $\Delta\mathrm{NLL}$ curves comparing the sensitivity of the fit with and without systematic uncertainties. For comparison, other curves are shown which remove the effect of jet-based systematic uncertainties, $\tau$-based systematic uncertainties, and MC statistical uncertainties.

The observed $\Delta\mathrm{NLL}$ curves for each analysis channel as a function of $\tilde d$, compared to the combined result. For the individual analysis channel $\Delta\mathrm{NLL}$ curves, only event yield information in the other SRs is used, ensuring that the Optimal Observable distributions in the other SRs do not influence the preferred value of $\tilde d$. The signal strength is constrained to be positive in these individual channel $\Delta\mathrm{NLL}$ curves. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit BDT distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit BDT distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit BDT distributions in the $Z\to \ell\ell$ CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit Optimal Observable distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit Optimal Observable distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit Optimal Observable distributions in the $Z\to \ell\ell$ CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit distribution of weighted event yields as a function of the Optimal Observable for all four SRs combined. The contributions of the different SRs are weighted by a factor of ln(1 + S/B), where S and B are the post-fit expected numbers of signal and background events in that region, respectively. The size of the combined statistical, experimental, and theoretical uncertainties is given.

The number of data and fitted simulated events in each bin of the $|d_{0}^{\mu}|$ distribution in the $10<p_{\textrm{T}}^{\mu}<20$ GeV selection in the $\mu-\mu$ channel.

The number of data and fitted simulated events in each bin of the $|d_{0}^{\mu}|$ distribution in the $20<p_{\textrm{T}}^{\mu}<250$ GeV selection in the $e-\mu$ channel.

The number of data and fitted simulated events in each bin of the $|d_{0}^{\mu}|$ distribution in the $20<p_{\textrm{T}}^{\mu}<250$ GeV selection in the $\mu-\mu$ channel.

The measurement of the ratio of the rate of decay of W bosons to τ-leptons and muons, $R(\tau/\mu)=B(W\rightarrow\tau\nu_\tau)/B(W\rightarrow \mu\nu_\mu)$.