Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected and observed CLs values per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Expected and observed CLs values per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Expected and observed cross-section upper-limit per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Expected and observed cross-section upper-limit per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Event selection cutflows for signal samples with $m(\tilde{\chi}_{1}^0)$ = 150 GeV and $\Delta m(\tilde{\chi}_{1}^\pm, \tilde{\chi}_{1}^0)$ = 1.5, 1.0, and 0.75 GeV, including all six production processes ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$). The cross-section used to obtain the initial number of events ($\sigma(\mathrm{n}_{\mathrm{jets}}) \geq 1$) refers to an emission of at least one gluon or quark with $p_{\mathrm{T}} > 50$ GeV at the parton level.

Event selection cutflows for signal samples with $m(\tilde{\chi}_{1}^0)$ = 150 GeV and $\Delta m(\tilde{\chi}_{1}^\pm, \tilde{\chi}_{1}^0)$ = 0.5, 0.35, and 0.25 GeV, including all six production processes ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$). The cross-section used to obtain the initial number of events ($\sigma(\mathrm{n}_{\mathrm{jets}}) \geq 1$) refers to an emission of at least one gluon or quark with $p_{\mathrm{T}} > 50$ GeV at the parton level.

The 95% CL limits on the effective coupling constants. The scale of new physics is set to $\Lambda = 1$ TeV.

The 95% CL limits on the branching fractions. The scale of new physics is set to $\Lambda = 1$ TeV.

The pre-fit distribution of the NN discriminant for the FCNC $tu\gamma$ LH signal scaled to 10 times the observed limit and all SM backgrounds as well as measured data in the SR. The last bin of the distribution contains the overflow and the first bin also contains the underflow. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction.

The pre-fit distribution of the photon transverse momentum for SM backgrounds as well as measured data in the $t\bar{t}\gamma$ CR. The last bin of the distribution contains the overflow. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction.

The pre-fit distribution of the photon transverse momentum for the all SM backgrounds as well as measured data in the $W\gamma$+jets CR. The last bin of the distribution contains the overflow. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction.

Observed and expected distributions of the invariant mass of the four-lepton system in the $A\to ZH\to 4\ell+X$ search for SR1 under a background-only fit to data. The total background includes the $q\overline{q}\to ZZ$, $gg\to ZZ$, $q\overline{q}\to ZZ$ (EW), $VVV$, $t\overline{t}V$, $t\overline{t}$, $Z$+jets and $WZ$ processes. The distribution of the $(m_{A}, m_{H}) = (510, 380)$ GeV signal is normalised to the observed upper limit on the cross-section (29.9 fb).

Observed and expected distributions of the invariant mass of the four-lepton system in the $A\to ZH\to 4\ell+X$ search for SR2 under a background-only fit to data. The total background includes the $q\overline{q}\to ZZ$, $gg\to ZZ$, $q\overline{q}\to ZZ$ (EW), $VVV$, $t\overline{t}V$, $t\overline{t}$, $Z$+jets and $WZ$ processes. The distribution of the $(m_{A}, m_{H}) = (510, 380)$ GeV signal is normalised to the observed upper limit on the cross-section (29.9 fb).

Observed and expected distributions of the invariant mass of the four-lepton system in the $A\to ZH\to 4\ell+X$ search for SR3 under a background-only fit to data. The total background includes the $q\overline{q}\to ZZ$, $gg\to ZZ$, $q\overline{q}\to ZZ$ (EW), $VVV$, $t\overline{t}V$, $t\overline{t}$, $Z$+jets and $WZ$ processes. The distribution of the $(m_{A}, m_{H}) = (510, 380)$ GeV signal is normalised to the observed upper limit on the cross-section (29.9 fb).

Observed and expected distributions of the invariant mass of the four-lepton system in the $A\to ZH\to 4\ell+X$ search for SR4 under a background-only fit to data. The total background includes the $q\overline{q}\to ZZ$, $gg\to ZZ$, $q\overline{q}\to ZZ$ (EW), $VVV$, $t\overline{t}V$, $t\overline{t}$, $Z$+jets and $WZ$ processes. The distribution of the $(m_{A}, m_{H}) = (510, 380)$ GeV signal is normalised to the observed upper limit on the cross-section (29.9 fb).

Observed and expected distributions of the invariant mass of the four-lepton system in the $A\to ZH\to 4\ell+X$ search for SR5 under a background-only fit to data. The total background includes the $q\overline{q}\to ZZ$, $gg\to ZZ$, $q\overline{q}\to ZZ$ (EW), $VVV$, $t\overline{t}V$, $t\overline{t}$, $Z$+jets and $WZ$ processes. The distribution of the $(m_{A}, m_{H}) = (510, 380)$ GeV signal is normalised to the observed upper limit on the cross-section (29.9 fb).

Observed and expected distributions of the invariant mass of the four-lepton system in the $A\to ZH\to 4\ell+X$ search for SR6 under a background-only fit to data. The total background includes the $q\overline{q}\to ZZ$, $gg\to ZZ$, $q\overline{q}\to ZZ$ (EW), $VVV$, $t\overline{t}V$, $t\overline{t}$, $Z$+jets and $WZ$ processes. The distribution of the $(m_{A}, m_{H}) = (510, 380)$ GeV signal is normalised to the observed upper limit on the cross-section (29.9 fb).

Observed and expected distributions of the invariant mass of the four-lepton system in the $A\to ZH\to 4\ell+X$ search for SR7 under a background-only fit to data. The total background includes the $q\overline{q}\to ZZ$, $gg\to ZZ$, $q\overline{q}\to ZZ$ (EW), $VVV$, $t\overline{t}V$, $t\overline{t}$, $Z$+jets and $WZ$ processes. The distribution of the $(m_{A}, m_{H}) = (510, 380)$ GeV signal is normalised to the observed upper limit on the cross-section (29.9 fb).

Local p0-values in the $(m_{H}, m_{R})$ plane for the $R\to SH\to 4\ell+E^{\textrm{miss}}_{\textrm{T}}$ search with $m_{S} = 160$ GeV.

Local p0-values in the $(m_{H}, m_{A})$ plane for the $A\to ZH\to 4\ell+X$ search.

The observed upper limits at 95% confidence level on $\sigma(gg\to R)\times \mathcal{B}(R\to SH)\times (H\to ZZ)$ across the $(m_{H}, m_{R})$ plane with $m_{S} = 160$ GeV for the $R\to SH\to 4\ell+E^{\textrm{miss}}_{\textrm{T}}$ search.

The expected upper limits at 95% confidence level on $\sigma(gg\to R)\times \mathcal{B}(R\to SH)\times (H\to ZZ)$ across the $(m_{H}, m_{R})$ plane with $m_{S} = 160$ GeV for the $R\to SH\to 4\ell+E^{\textrm{miss}}_{\textrm{T}}$ search.

The observed upper limits at 95% confidence level on $\sigma(gg\to A)\times \mathcal{B}(A\to ZH)\times (H\to ZZ)$ across the $(m_{H}, m_{A})$ plane for the $A\to ZH\to 4\ell+X$ search.

The expected upper limits at 95% confidence level on $\sigma(gg\to A)\times \mathcal{B}(A\to ZH)\times (H\to ZZ)$ across the $(m_{H}, m_{A})$ plane for the $A\to ZH\to 4\ell+X$ search.

Cut-flow of the raw events at each selection stage for the $R\rightarrow SH\rightarrow 4\ell + E^{\textrm{miss}}_{\textrm{T}}$ signal with a mass point of $(m_R, m_H ) = (390, 220)$ GeV and $m_S = 160$ GeV. The events are shown for the individual and combined four-lepton channels. The SFOS denotes same flavour and opposite sign lepton pairs selection and the final selection is shown for each signal region. The final selection criteria for SR1 to SR3 are defined in Table 2 of the paper.

Cut-flow of the raw events at each selection stage for the $A\rightarrow Z(\rightarrow jj/\ell^+\ell^-/\textrm{invisible})H(\rightarrow 4\ell)$ signal with a mass point of $(m_A, m_H ) = (330, 220)$ GeV. The events are shown for the individual and combined four-lepton channels. The SFOS denotes same flavour and opposite sign lepton pairs selection and the final selection is shown for each signal region. The final selection criteria for SR1 to SR7 are defined in Table 2 of the paper.

Cut-flow of the raw events at each selection stage for the $A\rightarrow Z(\rightarrow 2\ell)H(\rightarrow 2\ell+jj/\textrm{invisible})$ signal with a mass point of $(m_A, m_H ) = (330, 220)$ GeV. The events are shown for the individual and combined four-lepton channels. The SFOS denotes same flavour and opposite sign lepton pairs selection and the final selection is shown for each signal region. The final selection criteria for SR1 to SR7 are defined in Table 2 of the paper.

Summary of observed and expected background yields in the SR after the likelihood fit under the background-only hypothesis.

The expected minus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case with mass thresholds.

The expected plus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case with mass thresholds.

The expected minus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case with mass thresholds.

The observed exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The median expected exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The expected plus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The expected minus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The expected plus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The expected minus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The observed exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The median expected exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The expected plus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The expected minus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The expected plus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The expected minus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The observed exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The median expected exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The expected plus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The expected minus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The expected plus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The expected minus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

Distributions of kinematic variables in the signal regions for the $2\ell$ channels after applying all selection requirements. The histograms show the post-fit background predictions. The last bin includes the overflow. The distribution for $R_{\textrm{ISR}}$ in SR$2\ell$_ISR is plotted. The expected distribution for a benchmark signal model, normalized to the NLO+NLL cross-section times integrated luminosity, is also shown for comparison.

Distributions of kinematic variables in the signal regions for the $3\ell$ channels after applying all selection requirements. The histograms show the post-fit background predictions. The last bin includes the overflow. The distribution for $H_{3,1}^{\textrm{PP}}$ in SR$3\ell$_Low is plotted. The expected distribution for a benchmark signal model, normalized to the NLO+NLL cross-section times integrated luminosity, is also shown for comparison.

Distributions of kinematic variables in the signal regions for the $3\ell$ channels after applying all selection requirements. The histograms show the post-fit background predictions. The last bin includes the overflow. The distribution of $p_{\textrm{T}}^{\ell_{1}}$ in SR$3\ell$_Low is plotted. The expected distribution for a benchmark signal model, normalized to the NLO+NLL cross-section times integrated luminosity, is also shown for comparison.

Distributions of kinematic variables in the signal regions for the $3\ell$ channels after applying all selection requirements. The histograms show the post-fit background predictions. The last bin includes the overflow. The distribution of $p_{\mathrm{T\ ISR}}^{~\textrm{CM}}$ in SR$3\ell$_ISR is plotted. The expected distribution for a benchmark signal model, normalized to the NLO+NLL cross-section times integrated luminosity, is also shown for comparison.

Distributions of kinematic variables in the signal regions for the $3\ell$ channels after applying all selection requirements. The histograms show the post-fit background predictions. The last bin includes the overflow. The distribution of $R_{\textrm{ISR}}$ in SR$3\ell$_ISR is plotted. The expected distribution for a benchmark signal model, normalized to the NLO+NLL cross-section times integrated luminosity, is also shown for comparison.

Expected exclusion limits at 95% CL on the masses of C1/N2 and N1 from the analysis of 36.1fb$^{-1}$ of 13 TeV $pp$ collision data obtained from the 2$\ell$ search, assuming 100\% branching ratio of the sparticles to decay to SM $W/Z$ bosons and LSP.

Observed exclusion limits at 95% CL on the masses of C1/N2 and N1 from the analysis of 36.1fb$^{-1}$ of 13 TeV $pp$ collision data obtained from the 2$\ell$ search, assuming 100% branching ratio of the sparticles to decay to SM $W/Z$ bosons and LSP.

Expected exclusion limits at 95% CL on the masses of C1/N2 and N1 from the analysis of 36.1fb$^{-1}$ of 13 TeV $pp$ collision data obtained from the the $3\ell$ search, assuming 100% branching ratio of the sparticles to decay to SM $W/Z$ bosons and LSP.

Observed exclusion limits at 95% CL on the masses of C1/N2 and N1 from the analysis of 36.1fb$^{-1}$ of 13 TeV $pp$ collision data obtained from the the $3\ell$ search, assuming 100% branching ratio of the sparticles to decay to SM $W/Z$ bosons and LSP.

Expected exclusion limits at 95% CL on the masses of C1/N2 and N1 from the analysis of 36.1fb$^{-1}$ of 13 TeV $pp$ collision data obtained from the statistical combination of the $2\ell$ and 3$\ell$ search channels, assuming 100\% branching ratio of the sparticles to decay to SM $W/Z$ bosons and LSP.

Observed exclusion limits at 95% CL on the masses of C1/N2 and N1 from the analysis of 36.1fb$^{-1}$ of 13 TeV $pp$ collision data obtained from the statistical combination of the $2\ell$ and 3$\ell$ search channels, assuming 100% branching ratio of the sparticles to decay to SM $W/Z$ bosons and LSP.

Signal region acceptance for chargino-neutralino production in SR2L_High. The acceptance is with respect to all possible $W/Z$ decays.

Signal region efficiency for chargino-neutralino production in SR2L_High.

Signal region acceptanceXefficiency for chargino-neutralino production in SR2L_High.

Signal region acceptance for chargino-neutralino production in SR2L_Int. The acceptance is with respect to all possible $W/Z$ decays.

Signal region efficiency for chargino-neutralino production in SR2L_Int.

Signal region acceptanceXefficiency for chargino-neutralino production in SR2L_Int.

Signal region acceptance for chargino-neutralino production in SR2L_Low. The acceptance is with respect to all possible $W/Z$ decays.

Signal region efficiency for chargino-neutralino production in SR2L_Low.

Signal region acceptanceXefficiency for chargino-neutralino production in SR2L_Low.

Signal region acceptance for chargino-neutralino production in SR2L_ISR. The acceptance is with respect to all possible $W/Z$ decays.

Signal region efficiency for chargino-neutralino production in SR2L_ISR.

Signal region acceptanceXefficiency for chargino-neutralino production in SR2L_ISR.

Signal region acceptance for chargino-neutralino production in SR3L_High. The acceptance is with respect to all possible $W/Z$ decays.

Signal region efficiency for chargino-neutralino production in SR3L_High.

Signal region acceptanceXefficiency for chargino-neutralino production in SR3L_High.

Signal region acceptance for chargino-neutralino production in SR3L_Int. The acceptance is with respect to all possible $W/Z$ decays.

Signal region efficiency for chargino-neutralino production in SR3L_Int.

Signal region acceptanceXefficiency for chargino-neutralino production in SR3L_Int.

Signal region acceptance for chargino-neutralino production in SR3L_Low. The acceptance is with respect to all possible $W/Z$ decays.

Signal region efficiency for chargino-neutralino production in SR3L_Low.

Signal region acceptanceXefficiency for chargino-neutralino production in SR3L_Low.

Signal region acceptance for chargino-neutralino production in SR3L_ISR. The acceptance is with respect to all possible $W/Z$ decays.

Signal region efficiency for chargino-neutralino production in SR3L_ISR.

Signal region acceptanceXefficiency for chargino-neutralino production in SR3L_ISR.

Observed 95% CL upper limit on the production cross-section of C1/N2 from the analysis of 36.1fb$^{-1}$ of 13 TeV $pp$ collision data obtained from the 2$\ell$ search, assuming 100% branching ratio of the sparticles to decay to SM $W/Z$ bosons and N1.

Observed 95% CL upper limit on the production cross-section of C1/N2 from the analysis of 36.1fb$^{-1}$ of 13 TeV $pp$ collision data obtained from the 3$\ell$ search, assuming 100% branching ratio of the sparticles to decay to SM $W/Z$ bosons and N1.

Observed 95% CL upper limit on the production cross-section of C1/N2 from the analysis of 36.1fb$^{-1}$ of 13 TeV $pp$ collision data obtained from the statistical combination of 2 and 3$\ell$ searches, assuming 100% branching ratio of the sparticles to decay to SM $W/Z$ bosons and N1.

Signal cutflow for SR2L_High and m[C1/N2,N1] = [500,0] GeV. 5000 events were generated for this point. The unweighted number of events correspond to the number of selected MC events in the produced sample, while the weighted yield is normalized to the luminosity of the data.

Signal cutflow for SR2L_Int and m[C1/N2,N1] = [400,200] GeV. 10000 events were generated for this point. The unweighted number of events correspond to the number of selected MC events in the produced sample, while the weighted yield is normalized to the luminosity of the data.

Signal cutflow for SR2L_Low and m[C1/N2,N1] = [200,100] GeV. 20000 events were generated for this point. The unweighted number of events correspond to the number of selected MC events in the produced sample, while the weighted yield is normalized to the luminosity of the data.

Signal cutflow for SR2L_ISR and m[C1/N2,N1] = [200,100] GeV. 20000 events were generated for this point. The unweighted number of events correspond to the number of selected MC events in the produced sample, while the weighted yield is normalized to the luminosity of the data.

Signal cutflow for SR3L_High and m[C1/N2,N1] = [500,0] GeV. 5000 events were generated for this point. The unweighted number of events correspond to the number of selected MC events in the produced sample, while the weighted yield is normalized to the luminosity of the data.

Signal cutflow for SR3L_Int and m[C1/N2,N1] = [400,200] GeV. 10000 events were generated for this point. The unweighted number of events correspond to the number of selected MC events in the produced sample, while the weighted yield is normalized to the luminosity of the data.

Signal cutflow for SR3L_Low and m[C1/N2,N1] = [200,100] GeV. 20000 events were generated for this point. The unweighted number of events correspond to the number of selected MC events in the produced sample, while the weighted yield is normalized to the luminosity of the data.

Signal cutflow for SR3L_ISR and m[C1/N2,N1] = [200,100] GeV. 20000 events were generated for this point. The unweighted number of events correspond to the number of selected MC events in the produced sample, while the weighted yield is normalized to the luminosity of the data.

Observed and expected 95% CL upper limits on the cross-section of gluon-fusion initialted resonant production of the mass of the resonance times the branch ratio of X to HH and the BR of HH to gammagamma WW, without assuming the BR of H to gammagamma and H to WW.

Differential fiducial cross section of the Z boson $|y|$ in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the leading b-jet $|y|$ in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta \phi$ between Z boson and leading $b$-jet in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta y$ between Z boson and leading $b$-jet in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta R$ between Z boson and leading $b$-jet in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta \phi$ between the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta y$ between the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta R$ between the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the invariant mass of the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the Z boson $p_{\text{T}}$ in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $p_{\text{T}}$ of the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the ratio between the $p_{\text{T}}$ and the invariant mass of the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.