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 (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 (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 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.

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.

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-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.

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)$ = 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.

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.

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.

Distributions of representative kinematic variables in the misidentified background-enriched same-sign region: (a) the Higgs boson transverse momentum ($p_\text{T}^H$) in the $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ category, (b) the missing transverse momentum (${E_{\mathrm{T}}^{\mathrm{miss}}}$) in the $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ category, (c) the radial distance (dR$(\ell,\ell)$) between the two light leptons associated to the $Z\to\ell{}\ell$ decay process in the $ZH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ category, and (d) the invariant mass ($m_{\ell\ell}$) of the two light leptons associated to the $Z\to\ell{}\ell$ decay in the $ZH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ category. The hatched band represents the pre-fit statistical, experimental and theoretical uncertainties. The signal contributions are considered as part of the predictions and are normalized as predicted by the Standard Model.

Post-fit distributions for NN-based analysis of the NN-scores in the $ZH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ (a), $ZH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ (b), $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ (c) and $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ (d) categories. The hatched band indicates the total post-fit uncertainty of the total predicted yields. The post-fit signal contributions are considered as part of the predictions.

Post-fit distributions for NN-based analysis of the NN-scores in the $ZH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ (a), $ZH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ (b), $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ (c) and $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ (d) categories. The hatched band indicates the total post-fit uncertainty of the total predicted yields. The post-fit signal contributions are considered as part of the predictions.

Post-fit distributions for NN-based analysis of the NN-scores in the $ZH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ (a), $ZH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ (b), $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ (c) and $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ (d) categories. The hatched band indicates the total post-fit uncertainty of the total predicted yields. The post-fit signal contributions are considered as part of the predictions.

Post-fit distributions for NN-based analysis of the NN-scores in the $ZH(\\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ (a), $ZH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ (b), $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ (c) and $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ (d) categories. The hatched band indicates the total post-fit uncertainty of the total predicted yields. The post-fit signal contributions are considered as part of the predictions.

Post-fit yields from the NN-based fit performed with $\mu_{VH}^{\tau\tau}=\mu_{ZH}^{\tau\tau}=\mu_{WH}^{\tau\tau}$. The symbol $-$ is used when no events or $<10^{-2}$ events are present.

Event yields as a function of $\log_{10}(S/B)$ for data, background, and a Higgs boson signal with $m_{H} = 125$ GeV. Final-discriminant bins in all regions are combined into bins of $\log_{10}(S/B)$, where $S$ is the fitted signal and $B$ is the fitted background from the NN-based fit. The Higgs boson signal contribution is shown after re-scaling the SM cross-section according to the value of the signal strength extracted from data ($\mu = 1.28$). In the lower panel, the pull of the data relative to the background-only expectation, evaluated as the difference between the data and the background-only expectation divided by the square-root sum of the data and background statistical uncertainty, is shown. The solid line shows the expected pull in each bin for the best-fit signal value.

The fitted values of the Higgs boson signal strength $\mu_{VH}^{\tau\tau}$ for $m_{H} = 125$ GeV for the $WH$ and $ZH$ processes and their combination from the NN-based fit. The individual $\mu_{VH}^{\tau\tau}$ values for the $(W/Z)H$ processes are obtained from a simultaneous fit with the signal strength for each of the $WH$ and $ZH$ processes floating independently. The probability of compatibility of the individual signal strengths is 56\%.

Post-fit distributions for mass-based analysis for $m_{\mathrm{MMC}}$ in the (a) $ZH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$, (b) $ZH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$, and for $m_{\mathrm{2T}}$ in the (c) $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ and (d) $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ categories. The hatched band indicates the total post-fit uncertainty of the total predicted yields. The post-fit signal contributions are considered as part of the predictions.

Post-fit distributions for mass-based analysis for $m_{\mathrm{MMC}}$ in the (a) $ZH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$, (b) $ZH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$, and for $m_{\mathrm{2T}}$ in the (c) $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ and (d) $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ categories. The hatched band indicates the total post-fit uncertainty of the total predicted yields. The post-fit signal contributions are considered as part of the predictions.

Post-fit distributions for mass-based analysis for $m_{\mathrm{MMC}}$ in the (a) $ZH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$, (b) $ZH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$, and for $m_{\mathrm{2T}}$ in the (c) $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ and (d) $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ categories. The hatched band indicates the total post-fit uncertainty of the total predicted yields. The post-fit signal contributions are considered as part of the predictions.

Post-fit distributions for mass-based analysis for $m_{\mathrm{MMC}}$ in the (a) $ZH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$, (b) $ZH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$, and for $m_{\mathrm{2T}}$ in the (c) $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})# and (d) $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ categories. The hatched band indicates the total post-fit uncertainty of the total predicted yields. The post-fit signal contributions are considered as part of the predictions.

The fitted values from the NN-based fit of the Higgs boson signal strength $\mu_{VH}^{\tau\tau}$ for $m_{H} = 125$ GeV are shown separately for the $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$, $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$, $ZH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$, $ZH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ categories, the $WH$ and $ZH$ processes, and their full combination. The individual $\mu_{VH}^{\tau\tau}$ values for the $(W/Z)H$ processes are obtained from a simultaneous fit with the signal strength for each of the $WH$ and $ZH$ processes floating independently.

Summary table with the expected (exp) and observed (obs) values of the $\mu^{\tau\tau}_{\text{VH}}$ and the significance resulting from the mass-based fit.

Event yields as a function of $\log_{10}(S/B)$ for data, background, and a Higgs boson signal with $m_{H} = 125$ GeV. Final-discriminant bins in all regions are combined into bins of $\log_{10}(S/B)$, where $S$ is the fitted signal and $B$ is the fitted background from the NN-based fit. The Higgs boson signal contribution is shown after re-scaling the SM cross-section according to the value of the signal strength extracted from data ($\mu = 1.28$). In the lower panel, the pull of the data relative to the background-only expectation, evaluated as the difference between the data and the background-only expectation divided by the square-root sum of the data and background statistical uncertainty, is shown. The solid line shows the expected pull in each bin for the best-fit signal value.

Summary table with the expected (exp) and observed (obs) values of the $\mu^{\tau\tau}_{\text{VH}}$ and the significance resulting from the mass-based fit.

The fitted values from the mass-based fit of the Higgs boson signal strength $\mu_{VH}^{\tau\tau}$ for $m_{H} = 125$ GeV are shown for the $WH$ and $ZH$ processes and their combination. The individual $\mu_{VH}^{\tau\tau}$ values for the $(W/Z)H$ processes are obtained from a simultaneous fit with the signal strength for each of the $WH$ and $ZH$ processes floating independently. The probability of compatibility of the individual signal strengths is 85\%.

The fitted values from the mass-based fit of the Higgs boson signal strength $\mu_{VH}^{\tau\tau}$ for $m_{H} = 125$ GeV are shown separately for the $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$, $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$, $ZH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$, $ZH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ categories, the $WH$ and $ZH$ processes, and their full combination. The individual $\mu_{VH}^{\tau\tau}$ values for the $(W/Z)H$ processes are obtained from a simultaneous fit with the signal strength for each of the $WH$ and $ZH$ processes floating independently.

Distributions of $d\phi({E_{\mathrm{T}}^{\mathrm{miss}}},\tau_{1})$ (left) and $\eta$ of the leading $p_\text{T}$ $\tau_{\mathrm{had}-\mathrm{vis}}$ (right) in the $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ $Z\to\tau\tau$ validation region. The hatched band represents the pre-fit statistical, experimental and theoretical uncertainties. The signal contributions are considered as part of the predictions.

Distributions of $d\phi({E_{\mathrm{T}}^{\mathrm{miss}}},\tau_{1})$ (left) and $\eta$ of the leading $p_\text{T}$ $\tau_{\mathrm{had}-\mathrm{vis}}$ (right) in the $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ $Z\to\tau\tau$ validation region. The hatched band represents the pre-fit statistical, experimental and theoretical uncertainties. The signal contributions are considered as part of the predictions.

Distributions of the radial distance $dR(\tau_{1},\tau_{2})$ between the two $\tau_{\mathrm{had}-\mathrm{vis}}$ (left) and of the ${E_{\mathrm{T}}^{\mathrm{miss}}}$ (right) in the $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ $W\mathrm{ +jets}$ validation region. The hatched band represents the pre-fit statistical, experimental and theoretical uncertainties. The signal contributions are considered as part of the predictions.

Distributions of the radial distance $dR(\tau_{1},\tau_{2})$ between the two $\tau_{\mathrm{had}-\mathrm{vis}}$ (left) and of the ${E_{\mathrm{T}}^{\mathrm{miss}}}$ (right) in the $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ $W\mathrm{ +jets}$ validation region. The hatched band represents the pre-fit statistical, experimental and theoretical uncertainties. The signal contributions are considered as part of the predictions.

Distributions of the ${m_{\mathrm{MMC}}}$ (left) and of the $p_\text{T}$ of the leading $p_\text{T}$ light lepton (right) in the $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ $t\bar{t}$ validation region. The hatched band represents the pre-fit statistical, experimental and theoretical uncertainties. The signal contributions are considered as part of the predictions.

Distributions of the ${m_{\mathrm{MMC}}}$ (left) and of the $p_\text{T}$ of the leading $p_\text{T}$ light lepton (right) in the $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ $t\bar{t}$ validation region. The hatched band represents the pre-fit statistical, experimental and theoretical uncertainties. The signal contributions are considered as part of the predictions.

Distributions of the $m_{\mathrm{2T}}$ (left) and of the $\eta$ of the leading $p_\text{T}$ $\tau_{\mathrm{had}-\mathrm{vis}}$ in the $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ $Z\to\tau\tau$ validation region. The hatched band represents the pre-fit statistical, experimental and theoretical uncertainties. The signal contributions are considered as part of the predictions.

Distributions of the $m_{\mathrm{2T}}$ (left) and of the $\eta$ of the leading $p_\text{T}$ $\tau_{\mathrm{had}-\mathrm{vis}}$ in the $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ $Z\to\tau\tau$ validation region. The hatched band represents the pre-fit statistical, experimental and theoretical uncertainties. The signal contributions are considered as part of the predictions.

Distributions of the $dR({\tau_{\mathrm{lep}}},{\tau_{\mathrm{had}}})$ between the two objects associated to the Higgs boson decay and of the invariant mass $m_{\ell,\ell}$ of the two light leptons in the $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ all-same-sign validation region. The hatched band represents the pre-fit statistical, experimental and theoretical uncertainties. The signal contributions are considered as part of the predictions.

Distributions of the $dR({\tau_{\mathrm{lep}}},{\tau_{\mathrm{had}}})$ between the two objects associated to the Higgs boson decay and of the invariant mass $m_{\ell,\ell}$ of the two light leptons in the $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ all-same-sign validation region. The hatched band represents the pre-fit statistical, experimental and theoretical uncertainties. The signal contributions are considered as part of the predictions.

Distributions of the NN-score in the signal regions sidebands of the $ZH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ (a), $ZH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ (b), $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ (c), $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ (d) categories.The hatched band represents the pre-fit statistical, experimental and theoretical uncertainties. The signal contributions are considered as part of the predictions and are normalized to the fit result performed with one parameter of interest.

Distributions of the NN-score in the signal regions sidebands of the $ZH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ (a), $ZH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ (b), $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ (c), $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ (d) categories.The hatched band represents the pre-fit statistical, experimental and theoretical uncertainties. The signal contributions are considered as part of the predictions and are normalized to the fit result performed with one parameter of interest.

Distributions of the NN-score in the signal regions sidebands of the $ZH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ (a), $ZH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ (b), $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ (c), $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ (d) categories.The hatched band represents the pre-fit statistical, experimental and theoretical uncertainties. The signal contributions are considered as part of the predictions and are normalized to the fit result performed with one parameter of interest.

Distributions of the NN-score in the signal regions sidebands of the $ZH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ (a), $ZH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ (b), $WH(\tau_{\mathrm{lep}}\tau_{\mathrm{had}})$ (c), $WH(\tau_{\mathrm{had}}\tau_{\mathrm{had}})$ (d) categories.The hatched band represents the pre-fit statistical, experimental and theoretical uncertainties. The signal contributions are considered as part of the predictions and are normalized to the fit result performed with one parameter of interest.

Expected CLs values in the High mA SR, mA vs ma signal grid.

Observed CLs values in the High mA SR, mA vs ma signal grid.

Expected CLs values in the Low mA SR, mA vs ma signal grid.

Observed CLs values in the High mA SR, mA vs ma signal grid.

CLs+b values in the Low mA SR, mA vs tanB signal grid.

CLs+b values in the High mA SR, mA vs ma signal grid.

CLs+b values in the Low mA SR, mA vs ma signal grid.

Cut flow of the 2HDM+a signal points, gluon–gluon fusion production, Low mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

Cut flow of the 2HDM+a signal points, gluon–gluon fusion production, High mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

Cut flow of the 2HDM+a signal points, bb annihilation production, Low mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

Cut flow of the 2HDM+a signal points, bb annihilation production, High mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

A comparison of the observed and expected yields in the four bins of the Low mA SR.

A comparison of the observed and expected yields in the two bins of the High mA SR.

Expected exclusion contours at 95% CL as a function of mA and ma.

Observed exclusion contours at 95% CL as a function of mA and ma.

Expected +- 1sigma exclusion contours at 95% CL as a function of mA and ma.

Expected exclusion contours at 95% CL as a function of mA and ma.

Observed exclusion contours at 95% CL as a function of mA and ma.

Expected +- 1sigma exclusion contours at 95% CL as a function of mA and ma.

Acceptance, High mA SR, mA vs tanB grid, 400-750 GeV, bb prod

Acceptance, High mA SR, mA vs tanB grid, >750 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, 100-250 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, 250-400 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, 400-550 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, >550 GeV, bb prod

Acceptance, High mA SR, mA vs ma grid, 400-750 GeV, bb prod

Acceptance, High mA SR, mA vs ma grid, >750 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, 100-250 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, 250-400 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, 400-550 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, >550 GeV, bb prod

Acceptance, High mA SR, mA vs tanB grid, 400-750 GeV, ggF prod

Acceptance, High mA SR, mA vs tanB grid, >750 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, 100-250 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, 250-400 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, 400-550 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, >550 GeV, ggF prod

Acceptance, High mA SR, mA vs ma grid, 400-750 GeV, ggF prod

Acceptance, High mA SR, mA vs ma grid, >750 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, 100-250 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, 250-400 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, 400-550 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, >550 GeV, ggF prod

Efficiency, High mA SR, mA vs tanB grid, 400-750 GeV, bb prod

Efficiency, High mA SR, mA vs tanB grid, >750 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, 100-250 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, 250-400 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, 400-550 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, >550 GeV, bb prod

Efficiency, High mA SR, mA vs ma grid, 400-750 GeV, bb prod

Efficiency, High mA SR, mA vs ma grid, >750 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, 100-250 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, 250-400 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, 400-550 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, >550 GeV, bb prod

Efficiency, High mA SR, mA vs tanB grid, 400-750 GeV, ggF prod

Efficiency, High mA SR, mA vs tanB grid, >750 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, 100-250 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, 250-400 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, 400-550 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, >550 GeV, ggF prod

Efficiency, High mA SR, mA vs ma grid, 400-750 GeV, ggF prod

Efficiency, High mA SR, mA vs ma grid, >750 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, 100-250 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, 250-400 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, 400-550 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, >550 GeV, ggF prod