$a_2(1320)$ differential cross section in the m=-2 projection split into different reflectivity components, $\frac{d\sigma^+_{-2}}{dt}$ and $\frac{d\sigma^-_{-2}}{dt}$, in bins of four-momentum transfer. The first uncertainty is statistical, the second systematic.

$a_2(1320)$ differential cross section in the m=0 projection split into different reflectivity components, $\frac{d\sigma^+_{0}}{dt}$ and $\frac{d\sigma^-_{0}}{dt}$, in bins of four-momentum transfer. The first uncertainty is statistical, the second systematic.

$a_2(1320)$ differential cross section in the m=+1 projection split into different reflectivity components, $\frac{d\sigma^+_{+1}}{dt}$ and $\frac{d\sigma^-_{+1}}{dt}$, in bins of four-momentum transfer. The first uncertainty is statistical, the second systematic.

$a_2(1320)$ differential cross section in the m=+2 projection split into different reflectivity components, $\frac{d\sigma^+_{+2}}{dt}$ and $\frac{d\sigma^-_{+2}}{dt}$, in bins of four-momentum transfer. The first uncertainty is statistical, the second systematic.

Total invariant mass distribution $m_{4\ell}$ in Signal Region 2.

Selected events in the ($m_{4\ell}$, $\left\langle m_{\ell\ell}\right\rangle$) plane, SR1

Selected events in the ($m_{4\ell}$, $\left\langle m_{\ell\ell}\right\rangle$) plane, SR2

Expected background events in the $(m_{4\ell}, \langle m_{\ell\ell}\rangle)$ plane, SR1

Expected signal shape for $(m_{S}, m_{Z_d}) = (110~\text{GeV},30~\text{GeV}$) in the $(m_{4\ell}, \langle m_{\ell\ell}\rangle)$ plane, SR1. The signal distributions are normalized to unit volume.

Expected background events in the $(m_{4\ell}, \langle m_{\ell\ell}\rangle)$ plane, SR2. Events that satisfy $|m_{ab,cd} - m_Z| < 8~\text{GeV}$ are removed by the Z-veto requirement, however not all events lying behind the mask in the plot, corresponding to $\langle m_{\ell\ell}\rangle$ values between 83 and 99 GeV, are deleted.

Expected signal shape for $(m_{S}, m_{Z_{d}}) = (750~\text{GeV},150~\text{GeV}$) in the $(m_{4\ell}, \langle m_{\ell\ell}\rangle)$ plane, SR2. The signal distributions are normalized to unit volume. Events that satisfy $|m_{ab,cd} - m_Z| < 8~\text{GeV}$ are removed by the Z-veto requirement, however not all events lying behind the mask in the plot, corresponding to $\langle m_{\ell\ell}\rangle$ values between 83 and 99 GeV, are deleted.

Expected 95% CL limit on $\sigma \mathcal{B}(S\rightarrow Z_{d}Z_{d} \rightarrow 4\ell)$, SR1

Observed 95% CL limit on $\sigma \mathcal{B}(S\rightarrow Z_{d}Z_{d} \rightarrow 4\ell)$, SR1

Observed / Expected 95% CL limit ratio in SR1

Expected 95% CL limit on $\mathcal{B}(S\rightarrow Z_{d}Z_{d} \rightarrow 4\ell)$ with Z-veto, SR2

Observed 95% CL limit on $\sigma \mathcal{B}(S\rightarrow Z_{d}Z_{d} \rightarrow 4\ell)$ with Z-veto, SR2

Observed / Expected 95% CL limit ratio with Z-veto, SR2

Local $p_0$-values in the ($m_{S}, m_{Z_{d}}$) plane, SR1

Local $p_0$-values in the ($m_{S}, m_{Z_{d}}$) plane, SR2 with Z-veto applied

Observed data in Signal Region 1 in the $(m_{4\ell}, \langle m_{\ell\ell}\rangle)$ plane.

Observed data in Signal Region 2 in the $(m_{4\ell}, \langle m_{\ell\ell}\rangle)$ plane. Events that satisfy $|m_{ab,cd} - m_Z| < 8~\text{GeV}$ are removed by the Z-veto requirement, however not all events lying behind the mask in the plot, corresponding to <m_ll> values between 83 and 99 GeV, are deleted.

Distributions of the average dilepton mass $\left\langle m_{\ell\ell}\right\rangle = \frac{1}{2}\left(m_{ab} + m_{cd}\right)$ for Signal Region 1 in the $4e$ final state.

Distributions of the average dilepton mass $\left\langle m_{\ell\ell}\right\rangle = \frac{1}{2}\left(m_{ab} + m_{cd}\right)$ for Signal Region 1 in the $2e2\mu$ final state.

Distributions of the average dilepton mass $\left\langle m_{\ell\ell}\right\rangle = \frac{1}{2}\left(m_{ab} + m_{cd}\right)$ for Signal Region 1 in the $4\mu$ final state.

Distributions of the average dilepton mass $\left\langle m_{\ell\ell}\right\rangle = \frac{1}{2}\left(m_{ab} + m_{cd}\right)$ for Signal Region 2 in the $4e$ final state.

Distributions of the average dilepton mass $\left\langle m_{\ell\ell}\right\rangle = \frac{1}{2}\left(m_{ab} + m_{cd}\right)$ for Signal Region 2 in the $2e2\mu$ final state.

Distributions of the average dilepton mass $\left\langle m_{\ell\ell}\right\rangle = \frac{1}{2}\left(m_{ab} + m_{cd}\right)$ for Signal Region 2 in the $4\mu$ final state.

Distributions of the total invariant mass $m_{4\ell}$ for Signal Region 1 in the $4e$ final state. The SR1 requirement $m_{4\ell} < 115~\text{GeV}$ is not applied here.

Distributions of the total invariant mass $m_{4\ell}$ for Signal Region 1 in the $2e2\mu$ final state. The SR1 requirement $m_{4\ell} < 115~\text{GeV}$ is not applied here.

Distributions of the total invariant mass $m_{4\ell}$ for Signal Region 1 in the $4\mu$ final state. The SR1 requirement $m_{4\ell} < 115~\text{GeV}$ is not applied here.

Distributions of the total invariant mass $m_{4\ell}$ for Signal Region 2 in the $4e$ final state.

Distributions of the total invariant mass $m_{4\ell}$ for Signal Region 2 in the $2e2\mu$ final state.

Distributions of the total invariant mass $m_{4\ell}$ for Signal Region 2 in the $4\mu$ final state.

Data selection efficiency as a function of nuclear recoil energy

Data points used in analysis in log_10(S2)-S1 space

Isoscalar interaction coupling limit for Lagrangian 7

Isoscalar interaction coupling limit for Lagrangian 10

Isovector interaction coupling limit for Lagrangian 14

Isoscalar interaction coupling limit for Lagrangian 6

Isovector interaction coupling limit for Lagrangian 7

Isovector interaction coupling limit for Lagrangian 9

Isoscalar interaction coupling limit for Lagrangian 18

Isoscalar interaction coupling limit for Lagrangian 4

Isoscalar interaction coupling limit for Lagrangian 2

Isovector interaction coupling limit for Lagrangian 18

Isoscalar interaction coupling limit for Lagrangian 9

Isovector interaction coupling limit for Lagrangian 1

Isovector interaction coupling limit for Lagrangian 8

Isovector interaction coupling limit for Lagrangian 5

Isovector interaction coupling limit for Lagrangian 17

Isoscalar interaction coupling limit for Lagrangian 5

Isoscalar interaction coupling limit for Lagrangian 8

Isovector interaction coupling limit for Lagrangian 12

Isoscalar interaction coupling limit for Lagrangian 11

Isovector interaction coupling limit for Lagrangian 15

Isoscalar interaction coupling limit for Lagrangian 3

Isoscalar interaction coupling limit for Lagrangian 13

Isoscalar interaction coupling limit for Lagrangian 16

Isovector interaction coupling limit for Lagrangian 3

Isovector interaction coupling limit for Lagrangian 20

Isovector interaction coupling limit for Lagrangian 4

Isoscalar interaction coupling limit for Lagrangian 14

Isoscalar interaction coupling limit for Lagrangian 12

Isovector interaction coupling limit for Lagrangian 13

Isovector interaction coupling limit for Lagrangian 11

Isoscalar interaction coupling limit for Lagrangian 20

Isoscalar interaction coupling limit for Lagrangian 17

Isoscalar interaction coupling limit for Lagrangian 15

Isovector interaction coupling limit for Lagrangian 6

Isovector interaction coupling limit for Lagrangian 2

Isovector interaction coupling limit for Lagrangian 10

Isovector interaction coupling limit for Lagrangian 16

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

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the combined fit, when all other coupling modifiers are fixed to their SM predictions.

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

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

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

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{gghh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross. The best-fit value corresponds to the very first entry in the table.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{gghh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{tthh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross. The best-fit value corresponds to the very first entry in the table.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{tthh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{H}$, $c_{H◻}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross. The best-fit value corresponds to the very first entry in the table.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{H}$, $c_{H◻}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross.

The acceptance times efficiency [in %] for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ SLT channel. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ LTT channel. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ SLT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ LTT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ SLT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ LTT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

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

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

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

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

Distribution of the diphoton invariant mass in validation region VR2. The solid histograms are stacked to show the SM expectations after the 2&times;2D background estimation technique is applied. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. Statistical and systematic uncertainties are indicated by the shaded area. The lower panel of each plot shows the ratio of the data to the SM prediction for the respective bin. The first and last bins include the underflows and overflows respectively.

Distribution of the missing transverse momentum in validation region VR2. The solid histograms are stacked to show the SM expectations after the 2&times;2D background estimation technique is applied. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. Statistical and systematic uncertainties are indicated by the shaded area. The lower panel of each plot shows the ratio of the data to the SM prediction for the respective bin. The first and last bins include the underflows and overflows respectively.

Distribution of the diphoton invariant mass with all selections of the signal regions applied, except on m_{&gamma;&gamma;} itself, for signal region SR1h. The background estimation techniques described in the text are applied. The different backgrounds are stacked to add up to the total SM prediction in each bin. The predicted yields for signal benchmark models of varying &chi;&#771;⁰₁ mass are also overlaid (not stacked) assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; ) to equal 100%. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The sizes of the statistical and systematic uncertainties are indicated by the shaded areas. The lower panels show the ratio of the data to the SM prediction. Arrows indicate the borders of the signal region (|m_{&gamma;&gamma;}-125 GeV|&lt;5 GeV). The first and last bins include the underflows and overflows respectively.

Distribution of the diphoton invariant mass with all selections of the signal regions applied, except on m_{&gamma;&gamma;} itself, for signal region SR1Z. The background estimation techniques described in the text are applied. The different backgrounds are stacked to add up to the total SM prediction in each bin. The predicted yields for signal benchmark models of varying &chi;&#771;⁰₁ mass are also overlaid (not stacked) assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; ) to equal 50%. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The sizes of the statistical and systematic uncertainties are indicated by the shaded areas. The lower panels show the ratio of the data to the SM prediction. Arrows indicate the borders of the signal region (|m_{&gamma;&gamma;}-125 GeV|&lt;5 GeV). The first and last bins include the underflows and overflows respectively.

Distribution of the diphoton invariant mass with all selections of the signal regions applied, except on m_{&gamma;&gamma;} itself, for signal region SR2. The background estimation techniques described in the text are applied. The different backgrounds are stacked to add up to the total SM prediction in each bin. The predicted yields for signal benchmark models of varying &chi;&#771;⁰₁ mass are also overlaid (not stacked) assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; ) to equal 100%. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The sizes of the statistical and systematic uncertainties are indicated by the shaded areas. The lower panels show the ratio of the data to the SM prediction. Arrows indicate the borders of the signal region (|m_{&gamma;&gamma;}-125 GeV|&lt;5 GeV). The first and last bins include the underflows and overflows respectively.

Observed and expected limits on the pure higgsino cross-section at 95% CL assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; )=100% for different &chi;&#771;⁰₁ masses, obtained by a statistical combination of the three signal regions SR1h, SR1Z and SR2. The inner and outer bands indicate the 1&sigma; and 2&sigma; variation on the expected limit respectively.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(&chi;&#771;⁰₁ &rarr; hG&#771; ) as a function of the higgsino mass m(&chi;&#771;⁰₁) assuming it decays via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The &plusmn; 1&sigma; variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by &plusmn;1 &sigma;. Values of B(&chi;&#771;⁰₁ &rarr; hG&#771; ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(&chi;&#771;⁰₁ &rarr; hG&#771; ) as a function of the higgsino mass m(&chi;&#771;⁰₁) assuming it decays via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The &plusmn; 1&sigma; variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by &plusmn;1 &sigma;. Values of B(&chi;&#771;⁰₁ &rarr; hG&#771; ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(&chi;&#771;⁰₁ &rarr; hG&#771; ) as a function of the higgsino mass m(&chi;&#771;⁰₁) assuming it decays via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The &plusmn; 1&sigma; variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by &plusmn;1 &sigma;. Values of B(&chi;&#771;⁰₁ &rarr; hG&#771; ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(&chi;&#771;⁰₁ &rarr; hG&#771; ) as a function of the higgsino mass m(&chi;&#771;⁰₁) assuming it decays via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The &plusmn; 1&sigma; variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by &plusmn;1 &sigma;. Values of B(&chi;&#771;⁰₁ &rarr; hG&#771; ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(&chi;&#771;⁰₁ &rarr; hG&#771; ) as a function of the higgsino mass m(&chi;&#771;⁰₁) assuming it decays via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The &plusmn; 1&sigma; variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by &plusmn;1 &sigma;. Values of B(&chi;&#771;⁰₁ &rarr; hG&#771; ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(&chi;&#771;⁰₁ &rarr; hG&#771; ) as a function of the higgsino mass m(&chi;&#771;⁰₁) assuming it decays via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The &plusmn; 1&sigma; variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by &plusmn;1 &sigma;. Values of B(&chi;&#771;⁰₁ &rarr; hG&#771; ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Numbers of signal and background events in the signal regions. The respective background estimation techniques are applied. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The different backgrounds are stacked to add up to the total Standard Model prediction in each bin. The predicted yields for signal benchmark models of varying &chi;&#771;⁰₁ mass are also plotted (not stacked), assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; )=100% and a &chi;&#771;⁰₁ mass of 130 or 200 GeV. The statistical and systematic uncertainties are indicated by the shaded areas in the top plot. The bottom panel shows the statistical significance&nbsp;<a href="https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2020-025/">[Ref]</a> of the difference between the SM prediction and the observed data in each region.

Observed 95% CL limits in pb on the pure higgsino cross-section, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Limits are obtained by a statistical combination of the three signal regions SR1h, SR1Z and SR2, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Acceptances for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Acceptances are provided separately for signal region SR1h, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Acceptances for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Acceptances are provided separately for signal region SR1Z, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Acceptances for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Acceptances are provided separately for signal region SR2, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Efficiencies for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Efficiencies are provided separately for signal region SR1h, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Efficiencies for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Efficiencies are provided separately for signal region SR1Z, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Efficiencies for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Efficiencies are provided separately for signal region SR2, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Observed and expected numbers of events in the three signal regions. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The table also includes model-independent 95% CL upper limits on the visible number of BSM events (S⁹⁵_obs), the number of BSM events given the expected number of background events (S⁹⁵_exp) and the visible BSM cross-section (&lang;&epsilon; &sigma;&rang;_obs⁹⁵), all calculated from pseudo-experiments. The discovery p-value (p₀) is also shown and its value is capped at 0.5 if the observed number of events is below the expected number of events.

Cutflows of two benchmark signal points assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; )=100% for all three discovery signal regions. The initial selection includes the leptons veto. Only statistical uncertainties are included. Expected yields are normalised to a luminosity of 139 fb^-1.