The expected exclusion limit contours at 95% CL for the $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ production in the $m(\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2})$-$m(\tilde{\chi}^{0}_{1})$ plane.

The +1 $\sigma$ expected exclusion limit contours at 95% CL for the $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ production in the $m(\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2})$-$m(\tilde{\chi}^{0}_{1})$ plane.

The -1 $\sigma$ expected exclusion limit contours at 95% CL for the $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ production in the $m(\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2})$-$m(\tilde{\chi}^{0}_{1})$ plane.

The +1 $\sigma$ observed exclusion limit contours at 95% CL for the $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ production in the $m(\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2})$-$m(\tilde{\chi}^{0}_{1})$ plane.

The -1 $\sigma$ observed exclusion limit contours at 95% CL for the $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ production in the $m(\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2})$-$m(\tilde{\chi}^{0}_{1})$ plane.

The observed and expected exclusion limit at 95% CL for the $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ production in the $m(\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2})$-$m(\tilde{\chi}^{0}_{1})$ plane.

Expected and observed 95% CL exclusion upper limits on the higgsino production ($\tilde{\chi}\tilde{\chi} \equiv \tilde{\chi}^{0}_{1}\tilde{\chi}^{0}_{2}, \tilde{\chi}^{0}_{1}\tilde{\chi}^{\pm}_{1}, \tilde{\chi}^{0}_{2}\tilde{\chi}^{\pm}_{1}, \tilde{\chi}^{\pm}_{1}\tilde{\chi}^{\mp}_{1}$) cross-section in the signal of $h \tilde{G} h \tilde{G} $ as a function of the higgsino mass.

The distribution of $S_{E_\mathrm{T}^{\mathrm{miss}}}$ after the selection of diphoton candidates with $120~GeV < m_{\gamma\gamma} < 130~GeV$. Expected distributions are shown for the $\tilde{\chi}^{\pm}_{1} \tilde{\chi}^{0}_{2} \to W^{\pm} \tilde{\chi}^{0}_{1} h \tilde{\chi}^{0}_{1} $ signal with $m(\tilde{\chi}^{\pm}_{1} \tilde{\chi}^{0}_{2})=200~GeV$ and $m(\tilde{\chi}^{0}_{1})=0.5~GeV$, and the $h\tilde{G} h\tilde{G}$ signal with $m(\tilde{\chi}^{0}_{1})=150~GeV$ and $m(\tilde{G})=1~MeV$. These overlaid signal points are representative of the model kinematics. The sum in quadrature of the MC statistical and experimental systematic uncertainties in the total background is shown as the hatched bands, while the theoretical uncertainties in the background normalisation are not included. The $t\bar{t}\gamma$ and $t\bar{t}\gamma\gamma$ processes have a negligible contribution and are not represented. Overflow events are included in the rightmost bin. The lower panel shows the ratio of the data to the background.

Acceptance and efficiency for the $W^{\pm} \tilde{\chi}^{0}_{1} h \tilde{\chi}^{0}_{1}$ signal in each category.

Acceptance and efficiency for the $h \tilde{G} h \tilde{G}$ signal in each category.

Cut flow for the $W^{\pm} \tilde{\chi}^{0}_{1} h \tilde{\chi}^{0}_{1}$ signal sample at $m(\tilde{\chi}^{\pm}_{1} \tilde{\chi}^{0}_{2})=200~GeV$, $m(\tilde{\chi}^{0}_{1})=0.5~GeV$, with 260000 entries.

Cut flow for the $h \tilde{G} h \tilde{G}$ signal sample at $m(\tilde{\chi}^{0}_{1})=150~GeV$, $m(\tilde{G})=1~MeV$, with 50000 entries.

Event yields in the range 120 <$m\gamma\gamma$< 130 GeV for data, signal models, the SM Higgs boson background and non-resonant background in each analysis category after the fit.

Event yields in the range 120 <$m\gamma\gamma$< 130 GeV for data, signal models, the SM Higgs boson background and non-resonant background in each analysis category before the fit.

Acceptance for the $W^{\pm} \tilde{\chi}^{0}_{1} h \tilde{\chi}^{0}_{1}$ signals in the leptonic channel.

Efficiency for the $W^{\pm} \tilde{\chi}^{0}_{1} h \tilde{\chi}^{0}_{1}$ signals in the leptonic channel.

Acceptance for the $W^{\pm} \tilde{\chi}^{0}_{1} h \tilde{\chi}^{0}_{1}$ signals in the hadronic channel.

Efficiency for the $W^{\pm} \tilde{\chi}^{0}_{1} h \tilde{\chi}^{0}_{1}$ signals in the hadronic channel.

Acceptance for the $W^{\pm} \tilde{\chi}^{0}_{1} h \tilde{\chi}^{0}_{1}$ signals in the rest channel.

Efficiency for the $W^{\pm} \tilde{\chi}^{0}_{1} h \tilde{\chi}^{0}_{1}$ signals in the rest channel.

Acceptance for the $h \tilde{G} h \tilde{G}$ signals in the each channel.

Efficiency for the $h \tilde{G} h \tilde{G}$ signals in the each channel.

Data and SM predictions in SRs for the $1lb\bar{b}$ analysis for $E_{\mathrm{T}}^{\mathrm{miss}}$ in SR1Lbb-Medium. All SRs selections but the one on the quantity shown are applied. All uncertainties are included in the uncertainty band. Example SUSY models are superimposed for illustrative purposes.

Distributions of $m_{\gamma\gamma}$ before the final requirement on $m_{\gamma\gamma}$ in SR1Lyy-a. The expected contributions from both the peaking and non-peaking backgrounds are shown as stacked colored histograms. Two example SUSY models are superimposed for illustrative purposes.

Distributions of $m_{\gamma\gamma}$ before the final requirement on $m_{\gamma\gamma}$ in SR1Lyy-b. The expected contributions from both the peaking and non-peaking backgrounds are shown as stacked colored histograms. Two example SUSY models are superimposed for illustrative purposes.

Observed and predicted distributions for $m_{jj}$ in SRSS-j1. All SRs selections but the one on the quantity shown are applied. All uncertainties are included in the uncertainty band. An example SUSY model is superimposed for illustrative purposes.

Observed and predicted distributions for $m_{T2}$ in SRSS-j23. All SRs selections but the one on the quantity shown are applied. All uncertainties are included in the uncertainty band. An example SUSY model is superimposed for illustrative purposes.

The observed exclusion for the $0lb\bar{b}$ channel. Experimental and theoretical systematic uncertainties are applied to background and signal samples and illustrated by the yellow band and the red dotted contour lines, respectively. The red dotted lines indicate the $\pm$ 1 standard-deviation variation on the observed exclusion limit due to theoretical uncertainties in the signal cross-section.

The expected exclusion for the $0lb\bar{b}$ channel. Experimental and theoretical systematic uncertainties are applied to background and signal samples and illustrated by the yellow band and the red dotted contour lines, respectively. The red dotted lines indicate the $\pm$ 1 standard-deviation variation on the observed exclusion limit due to theoretical uncertainties in the signal cross-section.

The observed exclusion for the $1lb\bar{b}$ channel. Experimental and theoretical systematic uncertainties are applied to background and signal samples and illustrated by the yellow band and the red dotted contour lines, respectively. The red dotted lines indicate the $\pm$ 1 standard-deviation variation on the observed exclusion limit due to theoretical uncertainties in the signal cross-section.

The expected exclusion for the $1lb\bar{b}$ channel. Experimental and theoretical systematic uncertainties are applied to background and signal samples and illustrated by the yellow band and the red dotted contour lines, respectively. The red dotted lines indicate the $\pm$ 1 standard-deviation variation on the observed exclusion limit due to theoretical uncertainties in the signal cross-section.

The expected exclusion for the $1l\gamma\gamma$ channel. Experimental and theoretical systematic uncertainties are applied to background and signal samples and illustrated by the yellow band and the red dotted contour lines, respectively. The red dotted lines indicate the $\pm$ 1 standard-deviation variation on the observed exclusion limit due to theoretical uncertainties in the signal cross-section.

The observed exclusion for the $l^{\pm}l^{\pm}$ channel. Experimental and theoretical systematic uncertainties are applied to background and signal samples and illustrated by the yellow band and the red dotted contour lines, respectively. The red dotted lines indicate the $\pm$ 1 standard-deviation variation on the observed exclusion limit due to theoretical uncertainties in the signal cross-section.

The expected exclusion for the $l^{\pm}l^{\pm}$ channel. Experimental and theoretical systematic uncertainties are applied to background and signal samples and illustrated by the yellow band and the red dotted contour lines, respectively. The red dotted lines indicate the $\pm$ 1 standard-deviation variation on the observed exclusion limit due to theoretical uncertainties in the signal cross-section.

The observed cross-section exclusion as a function of the $\tilde{\chi}_{1}^{\pm}$/$\tilde{\chi}_{2}^{0}$ masses for the $3l$ analysis for signal models with $m(\tilde{\chi}_{2}^{0}) - m(\tilde{\chi}_{1}^{0}) = 130$ GeV.

The expected cross-section exclusion as a function of the $\tilde{\chi}_{1}^{\pm}$/$\tilde{\chi}_{2}^{0}$ masses for the $3l$ analysis for signal models with $m(\tilde{\chi}_{2}^{0}) - m(\tilde{\chi}_{1}^{0}) = 130$ GeV.

The observed cross-section exclusion as a function of the $\tilde{\chi}_{1}^{\pm}$/$\tilde{\chi}_{2}^{0}$ masses for the $1lb\bar{b}$ analysis for signal models with $m(\tilde{\chi}_{2}^{0}) - m(\tilde{\chi}_{1}^{0}) = 130$ GeV.

The expected cross-section exclusion as a function of the $\tilde{\chi}_{1}^{\pm}$/$\tilde{\chi}_{2}^{0}$ masses for the $1lb\bar{b}$ analysis for signal models with $m(\tilde{\chi}_{2}^{0}) - m(\tilde{\chi}_{1}^{0}) = 130$ GeV.

The observed cross-section exclusion as a function of the $\tilde{\chi}_{1}^{\pm}$/$\tilde{\chi}_{2}^{0}$ masses for the $1l\gamma\gamma$ analysis for signal models with $m(\tilde{\chi}_{2}^{0}) - m(\tilde{\chi}_{1}^{0}) = 130$ GeV.

The expected cross-section exclusion as a function of the $\tilde{\chi}_{1}^{\pm}$/$\tilde{\chi}_{2}^{0}$ masses for the $1l\gamma\gamma$ analysis for signal models with $m(\tilde{\chi}_{2}^{0}) - m(\tilde{\chi}_{1}^{0}) = 130$ GeV.

The observed cross-section exclusion as a function of the $\tilde{\chi}_{1}^{\pm}$/$\tilde{\chi}_{2}^{0}$ masses for the $l^{\pm}l^{\pm}$ analysis for signal models with $m(\tilde{\chi}_{2}^{0}) - m(\tilde{\chi}_{1}^{0}) = 130$ GeV.

The expected cross-section exclusion as a function of the $\tilde{\chi}_{1}^{\pm}$/$\tilde{\chi}_{2}^{0}$ masses for the $l^{\pm}l^{\pm}$ analysis for signal models with $m(\tilde{\chi}_{2}^{0}) - m(\tilde{\chi}_{1}^{0}) = 130$ GeV.

The observed cross-section exclusion as a function of the $\tilde{\chi}_{1}^{\pm}$/$\tilde{\chi}_{2}^{0}$ masses for the $0lb\bar{b}$ analysis for signal models assuming $m(\tilde{\chi}_{1}^{0}) = 0$ GeV.

The expected cross-section exclusion as a function of the $\tilde{\chi}_{1}^{\pm}$/$\tilde{\chi}_{2}^{0}$ masses for the $0lb\bar{b}$ analysis for signal models assuming $m(\tilde{\chi}_{1}^{0}) = 0$ GeV.

The observed cross-section exclusion as a function of the $\tilde{\chi}_{1}^{\pm}$/$\tilde{\chi}_{2}^{0}$ masses for the $1lb\bar{b}$ analysis for signal models assuming $m(\tilde{\chi}_{1}^{0}) = 0$ GeV.

The expected cross-section exclusion as a function of the $\tilde{\chi}_{1}^{\pm}$/$\tilde{\chi}_{2}^{0}$ masses for the $1lb\bar{b}$ analysis for signal models assuming $m(\tilde{\chi}_{1}^{0}) = 0$ GeV.

The observed cross-section exclusion as a function of the $\tilde{\chi}_{1}^{\pm}$/$\tilde{\chi}_{2}^{0}$ masses for the $1l\gamma\gamma$ analysis for signal models assuming $m(\tilde{\chi}_{1}^{0}) = 0$ GeV.

The expected cross-section exclusion as a function of the $\tilde{\chi}_{1}^{\pm}$/$\tilde{\chi}_{2}^{0}$ masses for the $1l\gamma\gamma$ analysis for signal models assuming $m(\tilde{\chi}_{1}^{0}) = 0$ GeV.

The observed cross-section exclusion as a function of the $\tilde{\chi}_{1}^{\pm}$/$\tilde{\chi}_{2}^{0}$ masses for the $l^{\pm}l^{\pm}$ analysis for signal models assuming $m(\tilde{\chi}_{1}^{0}) = 0$ GeV.

The expected cross-section exclusion as a function of the $\tilde{\chi}_{1}^{\pm}$/$\tilde{\chi}_{2}^{0}$ masses for the $l^{\pm}l^{\pm}$ analysis for signal models assuming $m(\tilde{\chi}_{1}^{0}) = 0$ GeV.

The observed cross-section exclusion as a function of the $\tilde{\chi}_{1}^{\pm}$/$\tilde{\chi}_{2}^{0}$ masses for the $3l$ analysis for signal models assuming $m(\tilde{\chi}_{1}^{0}) = 0$ GeV.

The expected cross-section exclusion as a function of the $\tilde{\chi}_{1}^{\pm}$/$\tilde{\chi}_{2}^{0}$ masses for the $3l$ analysis for signal models assuming $m(\tilde{\chi}_{1}^{0}) = 0$ GeV.

Acceptance for $0lb\bar{b}$ SRHad-Low signal region. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Efficiencies for $0lb\bar{b}$ SRHad-Low signal region.

Acceptance for $0lb\bar{b}$ SRHad-High signal region. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Efficiencies for $0lb\bar{b}$ SRHad-High signal region.

Upper cross section limits for the $0lb\bar{b}$ analysis.

Acceptance for $1lb\bar{b}$ SR1Lbb-Low signal region. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Efficiency for $1lb\bar{b}$ SR1Lbb-Low signal region. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Acceptance for $1lb\bar{b}$ SR1Lbb-Medium signal region. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Efficiency for $1lb\bar{b}$ SR1Lbb-Medium signal region. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Acceptance for $1lb\bar{b}$ SR1Lbb-High signal region. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Efficiency for $1lb\bar{b}$ SR1Lbb-High signal region. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Upper cross section limits for the $1lb\bar{b}$ analysis.

Acceptance for $1l\gamma\gamma$ signal region SR1Lyy-a. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Efficiencies for $1l\gamma\gamma$ signal region SR1Lyy-a.

Acceptance for $1l\gamma\gamma$ signal region SR1Lyy-b. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Efficiencies for $1l\gamma\gamma$ signal region SR1Lyy-b.

Upper cross section limits for the $1l\gamma\gamma$ analysis.

Acceptance for $l^{\pm}l^{\pm}$ signal region SRSS-j1. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Efficiencies for $l^{\pm}l^{\pm}$ signal region SRSS-j1. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Acceptance for $l^{\pm}l^{\pm}$ signal region SRSS-j23. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Efficiencies for $l^{\pm}l^{\pm}$ signal region SRSS-j23. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Acceptance for $3L$ SFOS signal region SR3L-SFOS-1J. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Efficiencies for $3L$ SFOS signal region SR3L-SFOS-1J. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Acceptance for $3L$ SFOS signal region SR3L-SFOS-0Ja. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Efficiencies for $3L$ SFOS signal region SR3L-SFOS-0Ja. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Acceptance for $3L$ SFOS signal region SR3L-SFOS-0Jb. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Efficiencies for $3L$ SFOS signal region SR3L-SFOS-0Jb. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Acceptance for $3L$ DFOS signal region SR3L-DFOS-0J. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Efficiencies for $3L$ DFOS signal region SR3L-DFOS-0J. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Acceptance for $3L$ DFOS signal region SR3L-DFOS-1Ja. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Efficiencies for $3L$ DFOS signal region SR3L-DFOS-1Ja. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Acceptance for $3L$ DFOS signal region SR3L-DFOS-1Jb. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Efficiencies for $3L$ DFOS signal region SR3L-DFOS-1Jb. Note that the acceptance is relative to the total chargino-neutralino production cross section.

Upper cross section limits for the $l^{\pm}l^{\pm}$ analysis.

Upper cross section limits for the $3l$ analysis.

Event selection cutflow for SM background (pre-fit) and representative signal samples for the $0lb\bar{b}$ SR. The masses of next-lightest-neutralinos and LSPs are reported. Only statistical uncertainties are shown. The MC statistics for the signal samples is 50K and 44K events, respectively. Samples are produced with generator filters which selects $h \rightarrow b\bar{b}$ and $W \rightarrow qq'$ decays. ``All events" for SUSY scenarios represent the total number of signal events for the models considered taking into account the BR of $W \rightarrow qq'$ (0.674) and $h \rightarrow b\bar{b}$ (0.582).

Event selection cutflow for SM background (pre-fit) and representative signal samples for the $1lb\bar{b}$ channel. The masses of next-lightest-neutralinos and LSPs are reported. Only statistical uncertainties are shown. The MC statistics for the signal samples is 28K, 29K and 29K events, respectively. Samples are produced with generator filters which selects $h \rightarrow b\bar{b}$ and $W \rightarrow l \nu$ decays. ``All events" for SUSY scenarios represent the total number of signal events for the models considered taking into account the BR of $W \rightarrow l \nu$ (0.324) and $h \rightarrow b\bar{b}$ (0.582).

Event selection cutflow is presented for SM background (pre-fit) and three signal scenarios for the $1l\gamma\gamma$ channel SRs. The masses of next-lightest-neutralinos and LSPs are reported. The uncertainties only include statistical uncertainties. The MC statistics for the signal samples is 10K events in all cases. Samples are produced with generator filters which selects $h \rightarrow \gamma\gamma$ and $W \rightarrow l \nu$ decays. The data-driven-estimated non-peaking background is not evaluated at each stage of the selection, hence only the peaking background cutflow is presented.

Event selection cutflow is presented for one signal scenario for SRSS-j23. The model with next-to-lightest neutralino mass of 175 GeV and massless $\tilde{\chi}^{0}_{1})$ with 100% BR is considered. The MC statistics for the signal samples is 75K events. Samples are produced with a generator filter which selects events with at least two leptons with $p_{\mathrm{T}} > 7$ GeV. ``2SS leptons" indicates the selection of two same-sign electrons or muons as described in the text (Section 6.4).

Event selection cutflow is presented for one signal scenario for SRSS-j1. The model with next-to-lightest neutralino mass of 175 GeV and massless $\tilde{\chi}^{0}_{1})$ with 100% BR is considered. The MC statistics for the signal samples is 75K events. Samples are produced with a generator filter which selects events with at least two leptons with $p_{\mathrm{T}} > 7$ GeV. ``2SS leptons" indicates the selection of two same-sign electrons or muons as described in the text (Section 6.4).

Event selection cut-flow is presented for one signal scenario for each signal regions. The model with next-to-lightest neutralino mass of 150 GeV and massless LSP with 100 \% BR is considered. The MC statistics for the signal samples is 123K events. Samples are produced with a generator filter which selects events with at least two leptons with pT $>$ 7 GeV. ``3L+Trigger" indicates three-lepton and trigger requirements. ``Preselection" includes the veto on bjets, ETmiss above 20 GeV and invariant mass of the three leptons above 20 GeV.