Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2016 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 2.3%.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2017 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 3.7%.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2018 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 1.8%.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the high-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the high-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Results of the background-only fit in the low-mass channel discovery region SR_LM_150. Both pre-fit and post-fit values are shown.

Results of the background-only fit in the low-mass channel discovery region SR_LM_300. Both pre-fit and post-fit values are shown.

The experimental efficiency of the low-mass channel for the exclusion and discovery signal regions as a function of higgsino mass. The experimental efficiency is defined as the number of events passing the detector-level event selections divided by the number of events passing the event selections for a perfect detector. The denominator is obtained by implementing particle-level event selections that emulate the detector-level selections. This treats the lack of availability of $b$-jet triggers as an inefficiency.

The particle-level acceptance for the low-mass exclusion and discovery signal regions, shown as a function of higgsino mass. The acceptance is defined as the fraction of signal events passing the particle-level event selection that emulates the detector-level selection. The acceptance calculation considers only those signal events where both higgsinos decay to Higgs bosons.

The experimental efficiency of the high-mass channel discovery regions as a function of higgsino mass. For each higgsino mass, the efficiency is shown for the SR-1 region corresponding to the mass. For masses above 1100 GeV, SR-1-1100 is used. The experimental efficiency is defined as the number of events passing the detector-level event selections divided by the number of events passing the event selections for a perfect detector. The denominator is obtained by implementing particle-level event selections that emulate the detector-level selections. The efficiency calculation considers only those signal events where both higgsinos decay to Higgs bosons.

The particle-level acceptance for the high-mass signal regions, shown as a function of higgsino mass. For each higgsino mass, the acceptance is shown for the SR-1 region corresponding to the mass. For masses above 1100 GeV, SR-1-1100 is used. The acceptance is defined as the fraction of signal events passing the particle-level event selection that emulates the detector-level selection. The acceptance calculation considers only those signal events where both higgsinos decay to Higgs bosons.

Cutflow for the low-mass channel for a representative 130 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 150 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 250 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 600 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 700 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 800 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 900 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 1000 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 1100 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 250 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 600 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 700 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 800 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 900 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1000 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1100 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Distribution of the transverse mass of the $W$-candidate $m_{\rm T}^{W}$ for the one lepton and two $b$-jets channel in SRlbb-1 and SRlbb-2, $m_{bb}$ sidebands. The background histograms are obtained from the background-only fit, and their uncertainty represents the total background uncertainty after the fit. The last bin includes overflow.

Distribution of the number of $b$-jets for the one lepton and two $b$-jets channel in SRlbb-1 and SRlbb-2 without the b-jet multiplicity requirement, central $m_{bb}$ bin. The background histograms are obtained from the background-only fit, and their uncertainty represents the total background uncertainty after the fit.

Distribution of the invariant mass of the $b$-jets $m_{bb}$ for the one lepton and two $b$-jets channel in the SRlbb-1 and SRlbb-2. The background histograms are obtained from the background-only fit, and their uncertainty represents the total background uncertainty after the fit.

Distribution of missing transverse momentum $E_{\rm T}^{\rm miss}$ in the one lepton and two photons signal regions SR$\ell\gamma\gamma$-1 and SR$\ell\gamma\gamma$-2 for the Higgs-mass window ($120\lt m_{\gamma\gamma} \lt 130$ GeV). The final column (background) is a simulation-based cross check. The contributions from non-Higgs backgrounds are scaled by 10 GeV / 50 GeV = 0.2 from the $m_{\gamma\gamma}$ sideband ($100 \lt m_{\gamma\gamma} \lt 120$ GeV and $130 \lt m_{\gamma\gamma} \lt 160$ GeV) into the Higgs-mass window. The last bin includes overflow. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

Distribution of the azimuth difference between the $W$ and Higgs boson candidates $\Delta\phi(W,h)$ in the one lepton and two photons signal regions SR$\ell\gamma\gamma$-1 and SR$\ell\gamma\gamma$-2 for the Higgs-mass window ($120\lt m_{\gamma\gamma} \lt 130$ GeV). The final column (background) is a simulation-based cross check. The contributions from non-Higgs backgrounds are scaled by 10 GeV / 50 GeV = 0.2 from the $m_{\gamma\gamma}$ sideband ($100 \lt m_{\gamma\gamma} \lt 120$ GeV and $130 \lt m_{\gamma\gamma} \lt 160$ GeV) into the Higgs-mass window. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

Distribution of the transverse mass of the $W$ and photon system $m_{\rm{T}}^{W\gamma_1}$ in the one lepton and two photons signal regions SR$\ell\gamma\gamma$-1 and SR$\ell\gamma\gamma$-2 for the Higgs-mass window ($120\lt m_{\gamma\gamma} \lt 130$ GeV). The final column (background) is a simulation-based cross check. The contributions from non-Higgs backgrounds are scaled by 10 GeV / 50 GeV = 0.2 from the $m_{\gamma\gamma}$ sideband ($100 \lt m_{\gamma\gamma} \lt 120$ GeV and $130 \lt m_{\gamma\gamma} \lt 160$ GeV) into the Higgs-mass window. The last bin includes overflow. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

Distribution of the transverse mass of the $W$ and photon system $m_{\rm{T}}^{W\gamma_2}$ in the one lepton and two photons signal regions SR$\ell\gamma\gamma$-1 and SR$\ell\gamma\gamma$-2 for the Higgs-mass window ($120\lt m_{\gamma\gamma} \lt 130$ GeV). The final column (background) is a simulation-based cross check. The contributions from non-Higgs backgrounds are scaled by 10 GeV / 50 GeV = 0.2 from the $m_{\gamma\gamma}$ sideband ($100 \lt m_{\gamma\gamma} \lt 120$ GeV and $130 \lt m_{\gamma\gamma} \lt 160$ GeV) into the Higgs-mass window. The last bin includes overflow. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

Results of the background-only fit to the diphoton invariant mass, $m_{\gamma\gamma}$, distribution in the one lepton and two photons signal region SR$l\gamma\gamma$-1. The contributions from SM Higgs boson production are constrained to the MC prediction and associated systematic uncertainties. The fit is performed on events with 100 GeV $ \lt m_{\gamma\gamma} \lt $ 160 GeV, with events in SR$l\gamma\gamma$-1 or SR$l\gamma\gamma$-2 in the Higgs-mass window (120 GeV $\le m_{\gamma\gamma} \le$ 130 GeV) excluded from the fit.

Results of the background-only fit to the diphoton invariant mass, $m_{\gamma\gamma}$, distribution in the one lepton and two photons signal region SR$l\gamma\gamma$-2. The contributions from SM Higgs boson production are constrained to the MC prediction and associated systematic uncertainties. The fit is performed on events with 100 GeV $ \lt m_{\gamma\gamma} \lt $ 160 GeV, with events in SR$l\gamma\gamma$-1 or SR$l\gamma\gamma$-2 in the Higgs-mass window (120 GeV $\le m_{\gamma\gamma} \le$ 130 GeV) excluded from the fit.

Results of the background-only fit to the diphoton invariant mass, $m_{\gamma\gamma}$, distribution in the one lepton and two photons validation region VR$l\gamma\gamma$-1. The contributions from SM Higgs boson production are constrained to the MC prediction and associated systematic uncertainties. The fit is performed on events with 100 GeV $ \lt m_{\gamma\gamma} \lt $ 160 GeV, with events in SR$l\gamma\gamma$-1 or SR$l\gamma\gamma$-2 in the Higgs-mass window (120 GeV $\le m_{\gamma\gamma} \le$ 130 GeV) excluded from the fit.

Results of the background-only fit to the diphoton invariant mass, $m_{\gamma\gamma}$, distribution in the one lepton and two photons validation region VR$l\gamma\gamma$-2. The contributions from SM Higgs boson production are constrained to the MC prediction and associated systematic uncertainties. The fit is performed on events with 100 GeV $ \lt m_{\gamma\gamma} \lt $ 160 GeV, with events in SR$l\gamma\gamma$-1 or SR$l\gamma\gamma$-2 in the Higgs-mass window (120 GeV $\le m_{\gamma\gamma} \le$ 130 GeV) excluded from the fit.

Distribution of effective mass $m_{\rm eff}$ in the validation region of the same-sign $e\mu$ channel. This validation region is defined by requiring one, two, or three jets, and reversing the $m_{lj}$, $m_{ljj}$ criteria. The distribution of a signal hypothesis is also shown.

Distribution of effective mass $m_{\rm eff}$ for the same-sign dilepton channel in the signal regions with one jet SR$ll$-1. SR$ll$-1 is the sum of SR$ee$-1, SR$e\mu$-1, and SR$\mu\mu$-1. All selection criteria are applied, except for the one on $m_{\rm eff}$. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(130,0)$ GeV. The last bin includes overflow.

Distribution of effective mass $m_{\rm eff}$ for the same-sign dilepton channel in the signal regions with two or three jets SR$ll$-2. SR$ll$-2 is the sum of SR$ee$-2, SR$e\mu$-2, and SR$\mu\mu$-2. All selection criteria are applied, except for the one on $m_{\rm eff}$. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(130,0)$ GeV. The last bin includes overflow.

Distribution of largest transverse mass $m_{\rm T}^{\rm max}$ for the same-sign dilepton channel in the signal regions with one jet SR$ll$-1. SR$ll$-1 is the sum of SR$ee$-1, SR$e\mu$-1, and SR$\mu\mu$-1. All selection criteria are applied, except for the one on $m_{\rm T}^{\rm max}$. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(130,0)$ GeV. The last bin includes overflow.

Distribution of largest transverse mass $m_{\rm T}^{\rm max}$ for the same-sign dilepton channel in the signal regions with two or three jets SR$ll$-2. SR$ll$-2 is the sum of SR$ee$-2, SR$e\mu$-2, and SR$\mu\mu$-2. All selection criteria are applied, except for the one on $m_{\rm T}^{\rm max}$. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(130,0)$ GeV. The last bin includes overflow.

Distribution of invariant mass of lepton and jet $m_{lj}$ for the same-sign dilepton channel in the signal regions with one jet SR$ll$-1. SR$ll$-1 is the sum of SR$ee$-1, SR$e\mu$-1, and SR$\mu\mu$-1. All selection criteria are applied, except for the one on $m_{lj}$. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(130,0)$ GeV. The last bin includes overflow.

Distribution of invariant mass of lepton and di-jet $m_{ljj}$ for the same-sign dilepton channel in the signal regions with two or three jets SR$ll$-2. SR$ll$-2 is the sum of SR$ee$-2, SR$e\mu$-2, and SR$\mu\mu$-2. All selection criteria are applied, except for the one on $m_{ljj}$. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(130,0)$ GeV. The last bin includes overflow.

One lepton and two $b$-jets channel: observed and expected 95% CL upper limits on the cross section normalised by the simplified model prediction as a function of the common mass $m_{\tilde{\chi}_1^\pm \tilde{\chi}^0_2}$ for $m_{\tilde{\chi}^0_1}=0$.

One lepton and two photons channel: observed and expected 95% CL upper limits on the cross section normalised by the simplified model prediction as a function of the common mass $m_{\tilde{\chi}_1^\pm \tilde{\chi}^0_2}$ for $m_{\tilde{\chi}^0_1}=0$.

Same-sign dilepton channel: observed and expected 95% CL upper limits on the cross section normalised by the simplified model prediction as a function of the common mass $m_{\tilde{\chi}_1^\pm \tilde{\chi}^0_2}$ for $m_{\tilde{\chi}^0_1}=0$.

Combination: observed and expected 95% CL upper limits on the cross section normalised by the simplified model prediction as a function of the common mass $m_{\tilde{\chi}_1^\pm \tilde{\chi}^0_2}$ for $m_{\tilde{\chi}^0_1}=0$. This combination is obtained using the result from the ATLAS three-lepton search, J. High Energy Phys. 04 (2014) 169, in addition to the three channels reported in this paper.

One lepton and two b-jets channel: Expected 95% CL exclusion contour for chargino neutralino production via Wh.

One lepton and two b-jets channel: Observed 95% CL exclusion contour for chargino neutralino production via Wh.

One lepton and two photons channel: Expected 95% CL exclusion contour for chargino neutralino production via Wh.

One lepton and two photons channel: Observed 95% CL exclusion contour for chargino neutralino production via Wh.

Same-sign dilepton channel: Expected 95% CL exclusion contour for chargino neutralino production via Wh.

Same-sign dilepton channel: Observed 95% CL exclusion contour for chargino neutralino production via Wh.

Combination: Expected 95% CL exclusion contour for chargino neutralino production via Wh.

Combination: Observed 95% CL exclusion contour for chargino neutralino production via Wh.

Combination: excluded model cross-section at 95% CL IN PB.

Acceptance SR$\ell\gamma\gamma$-1.

Acceptance SR$\ell\gamma\gamma$-2.

Efficiency SR$\ell\gamma\gamma$-1.

Efficiency SR$\ell\gamma\gamma$-2.

$m_{\gamma\gamma}$ distribution in the SR$l\gamma\gamma$-1 region for the full $m_{\gamma\gamma}$ window. The last column (background) is from a simulation-based cross check. $Z\gamma$ events, with $Z\rightarrow{ee}$ and one electron mis-identified as a photon, are included in the $Z\gamma(\gamma)$ background. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

$m_{\gamma\gamma}$ distribution in the SR$l\gamma\gamma$-2 region for the full $m_{\gamma\gamma}$ window. The last column (background) is from a simulation-based cross check. $Z\gamma$ events, with $Z\rightarrow{ee}$ and one electron mis-identified as a photon, are included in the $Z\gamma(\gamma)$ background. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

$m_{\gamma\gamma}$ distribution in the VR$l\gamma\gamma$-1 region for the full $m_{\gamma\gamma}$ window. The last column (background) is from a simulation-based cross check. $Z\gamma$ events, with $Z\rightarrow{ee}$ and one electron mis-identified as a photon, are included in the $Z\gamma(\gamma)$ background. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

$m_{\gamma\gamma}$ distribution in the VR$l\gamma\gamma$-2 region for the full $m_{\gamma\gamma}$ window. The last column (background) is from a simulation-based cross check. $Z\gamma$ events, with $Z\rightarrow{ee}$ and one electron mis-identified as a photon, are included in the $Z\gamma(\gamma)$ background. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

Acceptance for SRlbb-1, central $m_{bb}$ bin.

Acceptance for SRlbb-1, $m_{bb}$ sideband.

Acceptance for SRlbb-2, central $m_{bb}$ bin.

Acceptance for SRlbb-2, $m_{bb}$ sideband.

Efficiency for SRlbb-1, central $m_{bb}$ bin.

Efficiency for SRlbb-1, $m_{bb}$ sideband.

Efficiency for SRlbb-2, central $m_{bb}$ bin.

Efficiency for SRlbb-2, $m_{bb}$ sideband.

Acceptance for the same-sign $ee$ channel with one jet.

Acceptance for the same-sign $ee$ channel with two or three jets.

Acceptance for the same-sign $e\mu$ channel with one jet.

Acceptance for the same-sign $e\mu$ channel with two or three jets.

Acceptance for the same-sign $\mu\mu$ channel with one jet.

Acceptance for the same-sign $\mu\mu$ channel with two or three jets.

Efficiency for the same-sign $ee$ channel with one jet.

Efficiency for the same-sign $ee$ channel with two or three jets.

Efficiency for the same-sign $e\mu$ channel with one jet.

Efficiency for the same-sign $e\mu$ channel with two or three jets.

Efficiency for the same-sign $\mu\mu$ channel with one jet.

Efficiency for the same-sign $\mu\mu$ channel with two or three jets.

The observed 95\% C.L. limit for $pp\rightarrow G^{*}_{KK}\rightarrow hh\rightarrow b\bar{b}b\bar{b}$ in the bulk RS model with $k/\bar{M}_{Pl} = 2$, as a function of resonance mass.

The observed 95\% C.L. limit for $pp\rightarrow H\rightarrow hh\rightarrow b\bar{b}b\bar{b}$ with fixed $\Gamma_{H} = 1$ GeV, as a function of resonance mass.

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

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

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

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

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

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

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

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

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

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

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

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

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

Observed and predicted mTtot distribution in the b-tag category of the 2tau_h channel. Despite listing this as an exclusive final state (as there must be at least one b-jets), there is no explicit selection on the presence of additional light-flavour jets. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 300, 500 and 800 GeV and $\tan\beta$ = 10 in the hMSSM scenario are also provided.

Observed and predicted mTtot distribution for the b-inclusive selection in the 1l1tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. In the paper, the first bin is cut off at 60 GeV for aesthetics but contains underflows down to 50 GeV as in the HepData table. The last bin includes overflows. The prediction for a SSM Zprime with masses of 1500, 2000 and 2500 GeV are also provided.

Observed and predicted mTtot distribution for the b-inclusive selection in the 2tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The prediction for a SSM Zprime with masses of 1500, 2000 and 2500 GeV are also provided.

Observed and expected 95% CL upper limits on the b-associated Higgs boson production cross section times ditau branching fraction as a function of the boson mass.

Observed and expected 95% CL upper limits on the Drell Yan production cross section times ditau branching fraction as a function of the Zprime boson mass.

Observed and expected 95% CL upper limits on the Higgs boson production cross section times ditau branching fraction as a function of the boson mass and the relative strength of the b-associated production.

Ratio of the 95% CL upper limits on the production cross section times branching fraction for alternate Zprime models with respect to the SSM, both observed and expected are shown.

Acceptance, acceptance times efficiency and b-tag category fraction for a scalar boson produced by gluon-gluon fusion as a function of the scalar boson mass.

Acceptance, acceptance times efficiency and b-tag category fraction for a scalar boson produced by b-associated production as a function of the scalar boson mass.

Acceptance and acceptance times efficiency for a heavy gauge boson produced by Drell Yan as a function of the gauge boson mass.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Observed and expected 95% CL upper limits on the gluon-gluon fusion Higgs boson production cross section times ditau branching fraction as a function of the boson mass.

Observed correlations between the measured simplified template cross sections, including both the statistical and systematic uncertainties.

DM exclusion limit in the two-dimensional phase space of WIMP mass m_&chi; vs mediator mass m_med determined using the combined ee+&mu;&mu; channel. Both the observed and expected limits are presented, and the 1&sigma; uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_&chi;'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs, and the other grey line gives the perturbative limit (arXiv 1603.04156). The relic density line (arXiv 1603.04156) illustrates the combination of m_&chi; and m_med that would explain the observed DM relic density.

DM exclusion limit in the two-dimensional phase space of WIMP mass m_&chi; vs mediator mass m_med determined using the combined ee+&mu;&mu; channel. Both the observed and expected limits are presented, and the 1&sigma; uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_&chi;'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs, and the other grey line gives the perturbative limit (arXiv 1603.04156). The relic density line (arXiv 1603.04156) illustrates the combination of m_&chi; and m_med that would explain the observed DM relic density.

DM exclusion limit in the two-dimensional phase space of WIMP mass m_&chi; vs mediator mass m_med determined using the combined ee+&mu;&mu; channel. Both the observed and expected limits are presented, and the 1&sigma; uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_&chi;'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs, and the other grey line gives the perturbative limit (arXiv 1603.04156). The relic density line (arXiv 1603.04156) illustrates the combination of m_&chi; and m_med that would explain the observed DM relic density.

The 95% CL upper limits on the production cross-section of the ZH process with the prompt Z&rarr; ee and Z&rarr; &mu;&mu; decays and the invisible Higgs boson decays as a function of m_H, obtained from the combined ee+&mu;&mu; channel. The observed and expected limits are given, as well as the &pm;1&sigma; and &pm;2&sigma; error bands on the expected limits. The signal process is only modelled for the qq initial state, and the theory uncertainties on the signal prediction are not considered.

DM exclusion limit in the two-dimensional phase space of m_&chi; vs m_med set using the combined ee+&mu;&mu; channel. The DM model assumes a vector mediator, a fermionic WIMP, and the coupling parameters g_q = 0.25 and g_&chi; = 1. Both the observed and expected limits are presented, and the 1&sigma; uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_&chi;'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs. The relic density line [arXiv 1603.04156] illustrates the combination of m_&chi; and m_med that would explain the observed DM relic density.

DM exclusion limit in the two-dimensional phase space of m_&chi; vs m_med set using the combined ee+&mu;&mu; channel. The DM model assumes a vector mediator, a fermionic WIMP, and the coupling parameters g_q = 0.25 and g_&chi; = 1. Both the observed and expected limits are presented, and the 1&sigma; uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_&chi;'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs. The relic density line [arXiv 1603.04156] illustrates the combination of m_&chi; and m_med that would explain the observed DM relic density.

DM exclusion limit in the two-dimensional phase space of m_&chi; vs m_med set using the combined ee+&mu;&mu; channel. The DM model assumes a vector mediator, a fermionic WIMP, and the coupling parameters g_q = 0.25 and g_&chi; = 1. Both the observed and expected limits are presented, and the 1&sigma; uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_&chi;'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs. The relic density line [arXiv 1603.04156] illustrates the combination of m_&chi; and m_med that would explain the observed DM relic density.

DM exclusion limit in the two-dimensional phase space of m_&chi; vs m_med set using the combined ee+&mu;&mu; channel. The DM model assumes a vector mediator, a fermionic WIMP, and the coupling parameters g_q = 0.25 and g_&chi; = 1. Both the observed and expected limits are presented, and the 1&sigma; uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_&chi;'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs. The relic density line [arXiv 1603.04156] illustrates the combination of m_&chi; and m_med that would explain the observed DM relic density.

Limits on the WIMP and proton scattering cross section as a function of WIMP mass, obtained for simplified dark matter models with an axial-vector mediator (left) and a vector mediator (right). The solid black line shows the observed limit at the 90% confidence level from this search. The left plot incorporates the latest results from these competitive direct-search experiments, LUX [arXiv 1705.03380], PICO-2L [Phys. Rev. Lett. 114 (2015) 231302], and PICO-60 [arXiv 1702.07666]. The right plot incorporates the latest results from the following experiments, CRESST-II [Eur. Phys. J. C76 (2016) 25], CDMSlite [Phys. Rev. Lett. 116 (2016) 071301], PandaX-II [Phys. Rev. Lett. 117 (2016) 121303], LUX [Phys. Rev. Lett. 118 (2017) 021303], and XENON1T [arXiv 1705.06655].

Limits on the WIMP and proton scattering cross section as a function of WIMP mass, obtained for simplified dark matter models with an axial-vector mediator (left) and a vector mediator (right). The solid black line shows the observed limit at the 90% confidence level from this search. The left plot incorporates the latest results from these competitive direct-search experiments, LUX [arXiv 1705.03380], PICO-2L [Phys. Rev. Lett. 114 (2015) 231302], and PICO-60 [arXiv 1702.07666]. The right plot incorporates the latest results from the following experiments, CRESST-II [Eur. Phys. J. C76 (2016) 25], CDMSlite [Phys. Rev. Lett. 116 (2016) 071301], PandaX-II [Phys. Rev. Lett. 117 (2016) 121303], LUX [Phys. Rev. Lett. 118 (2017) 021303], and XENON1T [arXiv 1705.06655].