Signal region detector-level distribution for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$.

Signal region detector-level distribution for the observable $m_{e\mu}$.

Signal region detector-level distribution for the observable $p_{\mathrm{T}}^{e\mu}$.

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

Observed (points with error bars) and expected background (solid histograms) distributions for $E_{T}^{miss}$ in the signal region (a) SRL, (b) SRM and (c) SRH after the background-only fit applied to the CRs. The predicted signal distributions for the two models with a gluino mass of 2000 GeV and neutralino mass of 250 GeV (SRL), 1050 GeV (SRM) or 1950 GeV (SRH) are also shown for comparison. The uncertainties in the SM background are only statistical.

Observed and expected exclusion limit in the gluino-neutralino mass plane at 95% CL combined using the signal region with the best expected sensitivity at each point, for the full Run-2 dataset corresponding to an integrated luminosity of $139~\mathrm{fb}^{-1}$, for $\gamma/Z$ (a) and $\gamma/h$ (b) signal models. The black solid line corresponds to the expected limits at 95% CL, with the light (yellow) bands indicating the 1$\sigma$ exclusions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves, the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties. For each point in the higgsino-bino parameter space, the labels indicate the best-expected signal region, where L, M and H mean SRL, SRM and SRH, respectively.

Observed and expected exclusion limit in the gluino-neutralino mass plane at 95% CL combined using the signal region with the best expected sensitivity at each point, for the full Run-2 dataset corresponding to an integrated luminosity of $139~\mathrm{fb}^{-1}$, for $\gamma/Z$ (a) and $\gamma/h$ (b) signal models. The black solid line corresponds to the expected limits at 95% CL, with the light (yellow) bands indicating the 1$\sigma$ exclusions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves, the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties. For each point in the higgsino-bino parameter space, the labels indicate the best-expected signal region, where L, M and H mean SRL, SRM and SRH, respectively.

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Cutflow for the SRL selection, for two relevant signal points for both $\gamma/Z$ and $\gamma/h$ models, where the gluinos have mass of 2000 GeV and the neutralinos have a mass of 250 GeV (10000 generated events). The numbers are normalized to a luminosity of 139 $fb^{-1}$.

Cutflow for the SRM selection, for two relevant signal points for both $\gamma/Z$ and $\gamma/h$ models, where the gluinos have mass of 2000 GeV and the neutralinos have a mass of 1050 GeV (10000 generated events). The numbers are normalized to a luminosity of 139 $fb^{-1}$.

Cutflow for the SRH selection, for two relevant signal points for both $\gamma/Z$ and $\gamma/h$ models, where the gluinos have mass of 2000 GeV and the neutralinos have a mass of 1950 GeV (10000 generated events). The numbers are normalized to a luminosity of 139 $fb^{-1}$.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

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.

Corrected $m(HH)$ distribution in the resolved $4b$ signal region (dots), after the fit under the background-only hypothesis. The error bars on the $4b$ points represent the Poisson uncertainties corresponding to their event yields. The background model (teal histogram) is shown with its total post-fit uncertainty (gray band). The final bin includes overflow. Representative spin-0 signal hypotheses (dashed, dotted, and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. The bottom panel shows the difference between the $4b$ distribution and the background model, relative to the background model. No significant excess of data relative to the SM background is observed.

Corrected $m(HH)$ distribution in the resolved $4b$ signal region (dots), after the fit under the background-only hypothesis. The error bars on the $4b$ points represent the Poisson uncertainties corresponding to their event yields. The background model (teal histogram) is shown with its total post-fit uncertainty (gray band). The final bin includes overflow. Representative spin-2 signal hypotheses (dashed, dotted, and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. The bottom panel shows the difference between the $4b$ distribution and the background model, relative to the background model. No significant excess of data relative to the SM background is observed.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The steps up to the $b$-tag categorization are shown for the spin-0.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The steps up to the $b$-tag categorization are shown for the spin-2.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The efficiencies of the three b-tag categories are shown for the spin-0 scenario; this efficiency is obtained after the other selection steps including the SR definition. The signal efficiency in the 4b region has a maximum around 1.5 TeV. Above that value the track-jets start to merge together, and for the highest resonance masses the 2b category becomes the most efficient.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The efficiencies of the three b-tag categories are shown for the spin-2 scenario; this efficiency is obtained after the other selection steps including the SR definition. The signal efficiency in the 4b region has a maximum around 1.5 TeV. Above that value the track-jets start to merge together, and for the highest resonance masses the 2b category becomes the most efficient.

Comparison of the background model (stacked histograms) with data (dots) in the $2b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

Comparison of the background model (stacked histograms) with data (dots) in the $3b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

Comparison of the background model (stacked histograms) with data (dots) in the $4b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

The $m(HH)$ distributions in the boosted $2b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $2b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $3b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $3b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $4b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $4b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

Expected (dashed black lines) and observed (solid black lines) 95% CL upper limits on the cross-section of resonant $HH$ production in the spin-0 signal models. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits (colored bands) are shown. Expected limits using each of the resolved and boosted channels individually (dashed colored lines) are shown. The nominal $H\rightarrow b\bar{b}$ branching ratio is taken as 0.582.

Expected (dashed black lines) and observed (solid black lines) 95% CL upper limits on the cross-section of resonant $HH$ production in the spin-2 signal models. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits (colored bands) are shown. Expected limits using each of the resolved and boosted channels individually (dashed colored lines) are shown. The theoretical prediction for the bulk RS model with $k/\bar{M}_{\text{Pl}} = 1$ (solid red line) is shown; the decrease below 350 GeV is due to a sharp reduction in the $G^{*}_{\text{KK}} \rightarrow HH$ branching ratio. The nominal $H\rightarrow b\bar{b}$ branching ratio is taken as 0.582.

Observed and expected 95% CL lower limits on the $T$-quark mass as a function of the $T$-quark width-to-mass ratio and the branching fraction of the $T \rightarrow Ht$ decay ($\Gamma_{T}$ is the $T$-quark width).

Cutflow table listing the number of events passing each criterion for a $T$-quark hypothesis with a mass of 1.6 TeV and $\kappa_{T} = 0.5$. The initial signal event yield is the predicted number of $T$-quark events inclusive in the Higgs-boson and top-quark decays for 139 fb$^{-1}$.

Observed 95% CL upper limits on the single $T$-quark production cross-section as a function of the $T$-quark coupling $\kappa_{T}$ and $m_{T}$.

Expected 95% CL upper limits on the single $T$-quark production cross-section as a function of the $T$-quark coupling $\kappa_{T}$ and $m_{T}$.

Observed and expected 95% CL lower limits on the $T$-quark mass as a function of the $T$-quark width-to-mass ratio and the branching fraction of the $T \rightarrow Wb$ decay ($\Gamma_{T}$ is the $T$-quark width).

The BDT distribution (with x-axis zoomed in) of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the high mass region ($m^{*}_{b\bar{b}\gamma\gamma} > 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.857 in the high mass region are discarded. The range of BDT scores is from 0.8 to 1.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 300 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.85 for $m_{X}$ = 300 GeV are discarded.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 500 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.75 for $m_{X}$ = 500 GeV are discarded.

Distributions of $m_{\gamma\gamma}$ in high mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in high mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 300 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 370 fb.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 500 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 67 fb.

The number of data events observed in the 120 GeV < $m_{\gamma\gamma}$ < 130 GeV window, the number of $HH$ signal events expected for $\kappa_{\lambda}$ = 1 and for $\kappa_{\lambda}$ = 10, and events expected for single Higgs boson production (estimated using MC simulation), as well as for continuum background. For the single Higgs boson, Rest includes VBF, $WH$, $tHqb$, and $tHW$. The values are obtained from a fit of the Asimov data set generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$ = 1. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in $HH$ signals and single Higgs boson background include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

Observed and expected limits at 95% CL on the cross section of nonresonant Higgs boson pair production as a function of the Higgs boson self-coupling modifier $\kappa_{\lambda}= \lambda_{HHH}/\lambda^{\textrm{SM}}_{HHH}$. The expected constraints on $\kappa_{\lambda}$ are obtained with a background hypothesis excluding $pp \rightarrow HH$ production. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown. The theory prediction curve represents the scenario where all parameters and couplings are set to their SM values except for $\kappa_{\lambda}$. The uncertainty band of the theory prediction curve shows the cross-section uncertainty.

Values of the negative log-profile-likelihood ratio ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties in the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence-level intervals.

The number of events observed in the 120 < $m_{\gamma\gamma}$ < 130 GeV window in data, the number of events expected for scalar resonance signals of masses $m_{X}$ = 300 GeV and $m_{X}$ = 500 GeV assuming a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ equal to the observed exclusion limits of Figure 15, and events expected for SM $HH$ and single Higgs boson production (estimated using MC simulation), as well as for continuum background. The values are obtained from a fit of the Asimov data set generated under the signal-plus-background hypothesis. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in the resonant signals and the SM $HH$ and single-Higgs-boson backgrounds include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

Observed and expected limits at 95% CL on the production cross section of a narrow-width scalar resonance $X$ as a function of the mass $m_{X}$ of the hypothetical scalar particle. The black solid line represents the observed upper limits. The dashed line represents the expected upper limits. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown.

Breakdown of the dominant systematic uncertainties. The impact of the uncertainties is defined according to the statistical analysis described in Section 7. It corresponds to the relative variation of the expected upper limit on the cross section when re-evaluating the profile likelihood ratio after fixing the nuisance parameter in question to its best-fit value, while all remaining nuisance parameters remain free to float. The impact is shown in %. Only systematic uncertainties with an impact of at least 0.2% are shown. Uncertainties of the "Norm. + Shape" type affect both the normalization and the parameters of the functional form. The rest of the uncertainties affect only the yields.

Cutflow for nonresonant di-Higgs ggF signal sample, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 300 GeV, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 500 GeV, yields are normalized to 139 $fb^{-1}$.

Comparison of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

Fit results of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

The relative amount (purity) of expected events from SM $HH$ and single Higgs boson production processes for each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b).

The expected significance in each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b). The single Higgs boson processes and the di-photon continuum spectrum are considered as background.

Spurious signal result for the exponential pdf for the various ggF nonresonant di-Higgs categories. In each category, the spurious signal value ($n_{sp}$) and its ratio to the expected statistical error ($Z_{spur}$) from data are shown.

Spurious signal result for the exponential pdf as function of the resonant di-Higgs signal mass. The spurious signal value and its ratio to the expected statistical error from data are shown.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson ggF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson VBF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Values of the negative log-profile-likelihood ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties on the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence level intervals. The best fit value corresponds to $\kappa_{\lambda}$ = $2.8^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval. The expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.5}_{-2.4}(^{+7.3}_{-4.2})$ for the 1$\sigma$(2$\sigma$) confidence interval. The dashed curves represent values of the negative log-profile-likelihood where the Higgs boson branching fractions and the cross section of the production modes are varied as a function of $\kappa_{\lambda}$. In this case,the best fit value corresponds to $\kappa_{\lambda}$ = $2.7^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ and the expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.4}_{-2.5}(^{+7.3}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval.

Minimum BDT value of the events passing the selection criteria of the resonant search. The combined BDT score is formed using as coefficients $C_{1}$ = 0.65 and $C_{2}$ = 1 − $C_{1}$. The selection efficiency for the resonant $X \rightarrow HH$ signal is also shown.

Upper limits on the production cross section times branching fraction of the b* RH hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Upper limits on the production cross section times branching fraction of the b* VL hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Upper limits on the production cross section times branching fraction of the B+b hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Upper limits on the production cross section times branching fraction of the B+t hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Post-fit $m_\text{T}(\gamma, E_\text{T}^\text{miss})$ distribution in the inclusive signal region for the dark-photon search with the 125 GeV mass $\mathcal{B}(H \to \gamma \gamma_\text{d})$ signal normalization set to zero. A $H \to \gamma \gamma_\text{d}$ decay signal is shown for two different mass hypotheses, 125 GeV and 500 GeV, and scaled to a $\mathcal{B}(H \to \gamma \gamma_\text{d})$ of 2% and 1%, respectively. Events with $m_\text{T}(\gamma, E_\text{T}^\text{miss})$ larger than the rightmost bin boundary are added to that bin.

The 95% CL upper limit on the Higgs boson production cross-section times branching ratio to $\gamma \gamma_\text{d}$ is shown for different VBF-produced scalar-mediator-mass hypotheses in the NWA. The theoretically predicted cross-section of a Higgs boson produced via VBF and with the $\mathcal{B}(H \to \gamma \gamma_\text{d}) =$ 5% is superimposed on the $\pm 1\sigma$ and $\pm 2\sigma$ NNLO QCD + NLO EW uncertainty band of the expected production cross-section limit.

Post-fit $m_\text{jj}$ distribution in the inclusive signal region. The Higgs boson invisible decay signal is scaled to a $\mathcal{B}_\text{inv}$ of 37%. Events with $m_\text{jj}$ larger than the rightmost bin boundary are added to that bin.

Post-fit $m_\text{jj}$ distribution in the one-lepton control region $W_{\ell \nu}^\gamma$ CR. Events with $m_\text{jj}$ larger than the rightmost bin boundary are added to that bin.

Post-fit $m_\text{T}$ distribution in the one lepton control region. Events with $m_\text{T}$ larger than the rightmost bin boundary are added to that bin.

Post-fit photon centrality distribution in the zero lepton signal plus control region with the $\mathcal{B}_\text{inv}$ signal normalization set to zero in the fit.

Post-fit photon $E_\text{T}$ distribution in the zero lepton signal region with the $\mathcal{B}_\text{inv}$ signal normalization set to zero in the fit.

Post-fit photon centrality distribution in the zero lepton signal plus control region resulting from the fit to the $m_\text{jj}$ distribution for EW $Z \gamma + \text{jets}$. The post-fit uncertainties include statistical, experimental, and theory contributions.

Post-fit photon $E_\text{T}$ distribution in the zero lepton signal region resulting from the fit to the $m_\text{jj}$ distribution for EW $Z \gamma + \text{jets}$. The post-fit uncertainties include statistical, experimental, and theory contributions.

Post-fit DNN output score distribution in the one lepton control region.

Yields for the EW $Z \gamma + \text{jets}$ process are shown after each selection along with relative and absolute signal acceptance efficiencies.

Yields for the 125 GeV Higgs boson with $\mathcal{B}_\text{inv.} =$ 1 signal produced by the vector boson fusion process in association with a final state photon are shown after each selection along with relative and absolute signal acceptance efficiencies.

Yields for the 125 GeV Higgs boson with $\mathcal{B}(H \to \gamma \gamma_\text{d}) =$ 1 signal produced by the vector boson fusion process are shown after each selection along with relative and absolute signal acceptance efficiencies.

Expected +- 1 sigma 95% CL exclusion limit for the Zprime-2HDM model.

Expected +- 2 sigma 95% CL exclusion limit for the Zprime-2HDM model.

Observed 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected +- 1 sigma 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected +- 2 sigma 95% CL exclusion limit for the 2HDM+a model ggF production.

Observed 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected +- 1 sigma 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected +- 2 sigma 95% CL exclusion limit for the 2HDM+a model bbA production.

Observed 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected +- 1 sigma 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected +- 2 sigma 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected and observed upper limits at 95% CL on cross-section for Zprime-2HDM model.

Expected and observed upper limits at 95% CL on cross-section for ggF producton in the 2HDM+a model.

Expected and observed upper limits at 95% CL on cross-section for bbA producton in the 2HDM+a model.

Model-independent upper limits on the visible cross-section $σ_{vis, $h(\bar{b})+DM} ≡ σ_{h+DM} \times B(h \to b\bar{b}) \times \mathcal{A} \times \epsilon$ in the different signal regions.

Theory cross-section for Zprime-2HDM model.

Theory cross-section for bbA production in the 2HDM+a model.

Theory cross-section for ggF production in the 2HDM+a model.

Distribution of Higgs boson candidate mass in 2b region with MET=150-200 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=200-350 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=350-500 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=500-750 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET > 750 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=150-200 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=200-350 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=350-500 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET > 500 GeV.

Yields in 1-lepton control region.

Yields in 2-lepton control region.

MET distribution in 2b region of the 0-lepton channel.

MET distribution in 3b region of the 0-lepton channel.

Expected signal yields after certain selection cuts in 2b region with MET=150-200 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=200-350 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=350-500 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=500-750 GeV.

Expected signal yields after certain selection cuts in 2b region with MET > 750 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=150-200 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=200-350 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=350-500 GeV.

Expected signal yields after certain selection cuts in 3b region with MET > 500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET > 750 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET>500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET > 750 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET > 500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET > 750 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET > 500 GeV.