Overall efficiencies in the STA-STA (green) and TMS-TMS (red) dimuon categories, as well as their combination (black) as a function of $c\tau$ for the HAHM signal events with $m(Z_D) = 20\ GeV$. The solid curves show efficiencies achieved with the 2022 Run 3 triggers, whereas dashed curves show efficiencies for the subset of events selected by the triggers used in the 2018 Run 2 analysis. The efficiency is defined as the fraction of signal events that satisfy the criteria of the indicated trigger as well as the full set of offline selection criteria. The lower panel shows the relative improvement of the overall signal efficiency brought in by improvements in the trigger.

Comparison of the observed (black points) and expected (histograms) numbers of events in nonoverlapping $m_{\mu \mu}$ intervals in the STA-STA dimuon category, in the signal region optimized for the HAHM model. Yellow and green stacked histograms represent mean expected background contributions from QCD and DY, respectively, while statistical uncertainties in the total expected background are shown as hatched histograms. Signal contributions expected from simulated signals indicated in the legends are shown in red and blue. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow. All uncertainties shown are statistical only.

Comparison of the observed (black points) and expected (histograms) numbers of events in nonoverlapping $m_{\mu \mu}$ intervals in the STA-STA dimuon category, in the signal region optimized for the HAHM model. Yellow and green stacked filled histograms represent mean expected background contributions from QCD and DY, respectively, while statistical uncertainties in the total expected background are shown as hatched histograms. Signal contributions expected from simulated signals indicated in the legends are shown in red and blue. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Comparison of the observed (black points) and expected (histograms) numbers of events in nonoverlapping $m^{corr}_{\mu\mu}$ intervals in the STA-STA dimuon category, in the signal region optimized for the RPV SUSY model. Yellow and green stacked histograms represent mean expected background contributions from QCD and DY, respectively, while statistical uncertainties in the total expected background are shown as hatched histograms. Signal contributions expected from simulated signals indicated in the legends are shown in red and blue. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow. All uncertainties shown are statistical only.

Comparison of the observed (black points) and expected (histograms) numbers of events in nonoverlapping $m^{corr}_{\mu\mu}$ intervals in the STA-STA dimuon category, in the signal region optimized for the RPV SUSY model. Yellow and green stacked filled histograms represent mean expected background contributions from QCD and DY, respectively, while statistical uncertainties in the total expected background are shown as hatched histograms. Signal contributions expected from simulated signals indicated in the legends are shown in red and blue. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Distribution of min($d_0 / \sigma_{d_0}$) for TMS-TMS dimuons with $|\Delta\Phi| < \pi/30$, for events in all mass intervals combined. Events are required to satisfy all nominal selection criteria with the exception of the $d_0 / \sigma_{d_0}$ requirement. Notations are as in the Fig. 10 caption.

Distribution of min($d_0 / \sigma_{d_0}$) for TMS-TMS dimuons with $|\Delta\Phi| < \pi/30$, for events in all mass intervals combined, for both the validation (min($d_0 / \sigma_{d_0}$) < 6) and signal (min($d_0 / \sigma_{d_0}$) > 6) regions. The number of observed events (black circles) is overlaid with the stacked histograms showing the expected numbers of QCD (yellow) and DY (green) background events. Statistical uncertainties in the total expected background are shown as hatched histograms. Signal contributions expected from simulated signals indicated in the legends are shown in red and blue. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. Events are required to satisfy all nominal selection criteria with the exception of the $d_0 / \sigma_{d_0}$ requirement. The last bin includes events in the histogram overflow.

Distribution of min($d_0 / \sigma_{d_0}$) for TMS-TMS dimuons with $|\Delta\Phi| < \pi/4$, for events in all mass intervals combined. Events are required to satisfy all nominal selection criteria with the exception of the $d_0 / \sigma_{d_0}$ requirement. Notations are as in the Fig. 10 caption.

Distribution of min($d_0 / \sigma_{d_0}$) for TMS-TMS dimuons with $|\Delta\Phi| < \pi/4$, for events in all mass intervals combined, for both the validation (min($d_0 / \sigma_{d_0}$) < 6) and signal (min($d_0 / \sigma_{d_0}$) > 6) regions. The number of observed events (black circles) is overlaid with the stacked histograms showing the expected numbers of QCD (yellow) and DY (green) background events. Statistical uncertainties in the total expected background are shown as hatched histograms. Signal contributions expected from simulated signals indicated in the legends are shown in red and blue. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. Events are required to satisfy all nominal selection criteria with the exception of the $d_0 / \sigma_{d_0}$ requirement. The last bin includes events in the histogram overflow.

Comparison of observed and expected numbers of events in the TMS-TMS dimuon category, in the RPV SUSY study that requires $|\Delta\Phi| < \pi/4$, in bins of $m^{corr}_{\mu\mu}$. The number of observed events (black circles) is overlaid with the stacked histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m^{corr}_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: 6-10. Contributions expected from signal events predicted by the RPV SUSY model with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow. All uncertainties shown are statistical only.

Comparison of observed and expected numbers of events in bins of $m^{corr}_{\mu\mu}$ in the TMS-TMS dimuon category, in the signal regions optimized for the RPV SUSY model. The number of observed events (black circles) is overlaid with the stacked filled histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m^{corr}_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: 6-10. Hatched histograms show statistical uncertainties in the total expected background. Contributions expected from signal events predicted by the RPV SUSY model with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Comparison of observed and expected numbers of events in the TMS-TMS dimuon category, in the RPV SUSY study that requires $|\Delta\Phi| < \pi/4$, in bins of $m^{corr}_{\mu\mu}$. The number of observed events (black circles) is overlaid with the stacked histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m^{corr}_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: 10-20. Contributions expected from signal events predicted by the RPV SUSY model with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow. All uncertainties shown are statistical only.

Comparison of observed and expected numbers of events in bins of $m^{corr}_{\mu\mu}$ in the TMS-TMS dimuon category, in the signal regions optimized for the RPV SUSY model. The number of observed events (black circles) is overlaid with the stacked filled histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m^{corr}_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: 10-20. Hatched histograms show statistical uncertainties in the total expected background. Contributions expected from signal events predicted by the RPV SUSY model with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Comparison of observed and expected numbers of events in the TMS-TMS dimuon category, in the RPV SUSY study that requires $|\Delta\Phi| < \pi/4$, in bins of $m^{corr}_{\mu\mu}$. The number of observed events (black circles) is overlaid with the stacked histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m^{corr}_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: > 20. Contributions expected from signal events predicted by the RPV SUSY model with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow. All uncertainties shown are statistical only.

Comparison of observed and expected numbers of events in bins of $m^{corr}_{\mu\mu}$ in the TMS-TMS dimuon category, in the signal regions optimized for the RPV SUSY model. The number of observed events (black circles) is overlaid with the stacked filled histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m^{corr}_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: >20. Hatched histograms show statistical uncertainties in the total expected background. Contributions expected from signal events predicted by the RPV SUSY model with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Comparison of observed and expected numbers of events in the TMS-TMS dimuon category, in the HAHM study that requires $|\Delta\Phi| < \pi/30$, in bins of $m_{\mu\mu}$. The number of observed events (black circles) is overlaid with the stacked histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: 6-10. Signal contributions expected from simulated $H \rightarrow Z_DZ_D$ events with the parameters indicated in the legends are shown as red and blue histograms. Other notations are as in the Fig. 12 caption.

Comparison of observed and expected numbers of events in bins of $m_{\mu\mu}$ in the TMS-TMS dimuon category, in the signal regions optimized for the HAHM. The number of observed events (black circles) is overlaid with the stacked filled histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: 6-10. Hatched histograms show statistical uncertainties in the total expected background. Signal contributions expected from simulated $H \rightarrow Z_DZ_D$ events with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Comparison of observed and expected numbers of events in the TMS-TMS dimuon category, in the HAHM study that requires $|\Delta\Phi| < \pi/30$, in bins of $m_{\mu\mu}$. The number of observed events (black circles) is overlaid with the stacked histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: 10-20. Signal contributions expected from simulated $H \rightarrow Z_DZ_D$ events with the parameters indicated in the legends are shown as red and blue histograms. Other notations are as in the Fig. 12 caption.

Comparison of observed and expected numbers of events in bins of $m_{\mu\mu}$ in the TMS-TMS dimuon category, in the signal regions optimized for the HAHM. The number of observed events (black circles) is overlaid with the stacked filled histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: 10-20. Hatched histograms show statistical uncertainties in the total expected background. Signal contributions expected from simulated $H \rightarrow Z_DZ_D$ events with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Comparison of observed and expected numbers of events in the TMS-TMS dimuon category, in the HAHM study that requires $|\Delta\Phi| < \pi/30$, in bins of $m_{\mu\mu}$. The number of observed events (black circles) is overlaid with the stacked histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: > 20. Signal contributions expected from simulated $H \rightarrow Z_DZ_D$ events with the parameters indicated in the legends are shown as red and blue histograms. Other notations are as in the Fig. 12 caption.

Comparison of observed and expected numbers of events in bins of $m_{\mu\mu}$ in the TMS-TMS dimuon category, in the signal regions optimized for the HAHM. The number of observed events (black circles) is overlaid with the stacked filled histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: >20. Hatched histograms show statistical uncertainties in the total expected background. Signal contributions expected from simulated $H \rightarrow Z_DZ_D$ events with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 10\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 10\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 20\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 20\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 30\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 30\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 40\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 40\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 50\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 50\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 60\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 60\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 10\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 10\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 20\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 20\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 30\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 30\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 40\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 40\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 50\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 50\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 60\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 60\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 125\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The predicted cross section for $m(\tilde{q}) = 125\ GeV$ is 7200 pb, and falls outside the y-axis range.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 125\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The predicted cross section for $m(\tilde{q}) = 125\ GeV$ is 7200 pb, and falls outside the y-axis range.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 200\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The predicted cross section for $m(\tilde{q}) = 200 GeV$ is 840 pb, and falls outside the y-axis range.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 200\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The predicted cross section for $m(\tilde{q}) = 200 GeV$ is 840 pb, and falls outside the y-axis range.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 350\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The predicted cross section for $m(\tilde{q}) = 350\ GeV$ is 50 pb, and falls outside the y-axis range.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 350\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The predicted cross section for $m(\tilde{q}) = 350\ GeV$ is 50 pb, and falls outside the y-axis range.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 700\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The gray horizontal line indicates the theoretical value of the squark-antisquark production cross section with the uncertainties shown as the gray shaded band.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 700\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The gray horizontal line indicates the theoretical value of the squark-antisquark production cross section with the uncertainties shown as the gray shaded band.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 1150\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The gray horizontal line indicates the theoretical value of the squark-antisquark production cross section with the uncertainties shown as the gray shaded band.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 1150\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The gray horizontal line indicates the theoretical value of the squark-antisquark production cross section with the uncertainties shown as the gray shaded band.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 1600\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The gray horizontal line indicates the theoretical value of the squark-antisquark production cross section with the uncertainties shown as the gray shaded band.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 1600\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The gray horizontal line indicates the theoretical value of the squark-antisquark production cross section with the uncertainties shown as the gray shaded band.

Fractions of signal events with zero (green), one (blue), and two (red) STA muons matched to TMS muons by the STA to TMS association procedure, as a function of generated $L_{xy}$, in all HAHM signal samples combined.

Efficiencies of the Run 2 and Run 3 displaced dimuon triggers as a function of $c\tau$ for the HAHM signal events with $m(Z_D) = 50\ GeV$. The efficiency is defined as the fraction of simulated events that satisfy the requirements of the following sets of trigger paths: the Run 2 (2018) triggers (dashed black); the Run 3 (2022, L3) triggers (blue); the Run 3 (2022, L2) triggers (red); and the OR of all these triggers (Run 3 (2022), black). The lower panel shows the ratio of the overall Run 3 (2022) efficiency to the Run 2 (2018) efficiency.

Efficiencies of the Run 2 (2018) (red) and Run 3 (2022) (black) sets of displaced dimuon triggers as a function of $m(Z_D)$ for the HAHM signal events with $c\tau = 1\ cm$. The efficiency is defined as the fraction of simulated events that satisfy the detector acceptance and the requirements of the indicated set of trigger paths. The lower panel shows the ratio of the Run 3 (2022) efficiency to the Run 2 (2018) efficiency.

Efficiencies of the Run 2 (2018) (red) and Run 3 (2022) (black) sets of displaced dimuon triggers as a function of $m(Z_D)$ for the HAHM signal events with $c\tau = 10\ m$. The efficiency is defined as the fraction of simulated events that satisfy the detector acceptance and the requirements of the indicated set of trigger paths. The lower panel shows the ratio of the Run 3 (2022) efficiency to the Run 2 (2018) efficiency.

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM signal with $m(Z_D) = 20\ GeV$ in different years of data taking. Efficiencies are computed as the ratios of the number of simulated signal events in which at least one dimuon candidate passes all 2016 (dashed green), 2018 (dashed red), and 2022 (solid black) trigger and offline selection criteria to the total number of simulated signal events. The lower panel shows the ratio of the 2022 efficiency to the 2018 efficiency (dashed red) and to the 2016 efficiency (dashed green).

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM signal with $m(Z_D) = 50\ GeV$ in different years of data taking. Efficiencies are computed as the ratios of the number of simulated signal events in which at least one dimuon candidate passes all 2016 (dashed green), 2018 (dashed red), and 2022 (solid black) trigger and offline selection criteria to the total number of simulated signal events. The lower panel shows the ratio of the 2022 efficiency to the 2018 efficiency (dashed red) and to the 2016 efficiency (dashed green).

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM model with $m(Z_D) = 10\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM model with $m(Z_D) = 20\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM model with $m(Z_D) = 30\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM model with $m(Z_D) = 40\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM model with $m(Z_D) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM model with $m(Z_D) = 60\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 125\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 200\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 350\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 700\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 1150\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 1600\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 700\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 500\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 1150\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 500\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 1600\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 500\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the TMS-TMS dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ smaller than 20 cm in the HAHM signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the STA-STA dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ smaller than 20 cm in the HAHM signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the TMS-TMS dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ 20-70 cm in the HAHM signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the STA-STA dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ 20-70 cm in the HAHM signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the STA-STA dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ 70-500 cm in the HAHM signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the TMS-TMS dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ smaller than 20 cm in the RPV SUSY signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the STA-STA dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ smaller than 20 cm in the RPV SUSY signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the TMS-TMS dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ 20-70 cm in the RPV SUSY signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the STA-STA dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ 20-70 cm in the RPV SUSY signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the STA-STA dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ 70-500 cm in the RPV SUSY signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Observed upper limits at 95% CL on B($\text{H} \rightarrow \text{a}_{1}\text{a}_{1}$) in percent, obtained from the combination of the $\mu\mu$bb and $\tau\tau$bb channels. The results are obtained as functions of $m_{\text{a}_{1}}$ for 2HDM+S Type I (independent of tan$\beta$), Type II (tan$\beta$=2.0), Type III (tan$\beta$=2.0), and Type IV (tan$\beta$=0.6), respectively.

Observed exclusion limits at the 95% CL in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$ and a $-1 \sigma$ deviation of the NNLO+NNLL theoretical cross-section of a $\tilde{t}_1$ pair-production.

Expected exclusion limits at the 95% CL in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

The $1\sigma$ variation of expected 95% CL exclusion contour obtained by varying MC statistical uncertainties, detector-related systematic uncertainties, and theoretical uncertainties in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

The $2\sigma$ variation of expected 95%CL exclusion contour obtained by varying MC statistical uncertainties, detector-related systematic uncertainties, and theoretical uncertainties in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Observed exclusion limits at the 95% CL in the $m_{\tilde{t}_1} - $BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1)$ plane, assuming $m_{\tilde\chi^0_1} = 1 \mathrm{GeV}$.

Observed exclusion limits at the 95% CL in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming $m_{\tilde\chi^0_1} = 1 \mathrm{ GeV}$ and a $+1 \sigma$ deviation of the NNLO+NNLL theoretical cross-section of a $\tilde{t}_1$ pair-production.

Observed exclusion limits at the 95% CL in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming $m_{\tilde\chi^0_1} = 1 \mathrm{ GeV}$ and a $-1 \sigma$ deviation of the NNLO+NNLL theoretical cross-section of a $\tilde{t}_1$ pair-production.

Expected exclusion limits at the 95% CL in the $m_{\tilde{t}_1} - $BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1)$ plane, assuming $m_{\tilde\chi^0_1} = 1 \mathrm{GeV}$.

The $1\sigma$ variation of expected 95% CL exclusion contour obtained by varying MC statistical uncertainties, detector-related systematic uncertainties, and theoretical uncertainties in the $m_{\tilde{t}_1} - $BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1)$ plane, assuming $m_{\tilde\chi^0_1} = 1 \mathrm{GeV}$.

The $2\sigma$ variation of expected 95%CL exclusion contour obtained by varying MC statistical uncertainties, detector-related systematic uncertainties, and theoretical uncertainties in the $m_{\tilde{t}_1} - $BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1)$ plane, assuming $m_{\tilde\chi^0_1} = 1 \mathrm{GeV}$.

Observed upper limits on the top-spartner pair production cross-section in fb at the 95% CL in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$

Expected upper limits on the top-spartner pair production cross-section in fb at the 95% CL in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$

Observed upper limits on the top-spartner pair production cross-section in fb at the 95% CL in the $m_{\tilde{t}_1} - $BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1)$ plane, assuming $m_{\tilde\chi^0_1} = 1 \mathrm{GeV}$.

Expected upper limits on the top-spartner pair production cross-section in fb at the 95% CL in the $m_{\tilde{t}_1} - $BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1)$ plane, assuming $m_{\tilde\chi^0_1} = 1 \mathrm{GeV}$.

Post-fit distribution of $N_{\mathrm{tops}}^{\mathrm{DNN}} (\mathrm{R=1.0})$ in the SRA signal region presented without the associated SRA applied to the variable. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Post-fit distribution of $m_{\mathrm{T2}}(j^{b}_{R=1.0}, c)$ in the SRA signal region presented without the associated SRA applied to the variable. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Post-fit distribution of $p_{\mathrm{T}}(c_{1})$ in the SRB signal region presented without the associated SRB applied to the variable. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Post-fit distribution of $m_{\mathrm{T}}(j, \mathrm{E}_{\mathrm{T}}^{\mathrm{miss}})_{\mathrm{close}}$ in the SRB signal region presented without the associated SRB applied to the variable. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Post-fit distribution of $\mathrm{E}_{\mathrm{T}}^{\mathrm{miss}} \textrm{Significance}$ in the SRC signal region presented without the associated SRC applied to the variable. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Post-fit distribution of $m_{\mathrm{T}}(j, \mathrm{E}_{\mathrm{T}}^{\mathrm{miss}})_{\mathrm{close}}$ in the SRC signal region presented without the associated SRC applied to the variable. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Post-fit distribution of NN signal score in the SRD signal region presented without the associated SRD applied to the variable. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Post-fit distribution of $m_{\mathrm{eff}}$ in the SRD signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Pull plots showing the SRA, SRB and SRC post-fit data and SM agreement using the background-only fit configuration

Pull plots showing the SRD post-fit data and SM agreement using the background-only fit configuration

Cutflow for the reference point $(m_{\tilde{t}_1}, m_{\tilde\chi^0_1})=$ (900,1) , (700,300), (550,375) in the Signal region A. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(m_{\tilde{t}_1}, m_{\tilde\chi^0_1})=$ (900,1) , (700,300), (550,375) in the Signal region B. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(m_{\tilde{t}_1}, m_{\tilde\chi^0_1})=$ (900,1) , (700,300), (550,375) in the Signal region C. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(m_{\tilde{t}_1}, m_{\tilde\chi^0_1})=$ (900,1) , (700,300), (550,375) in the Signal region D 750. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(m_{\tilde{t}_1}, m_{\tilde\chi^0_1})=$ (900,1) , (700,300), (550,375) in the Signal region D 1000. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(m_{\tilde{t}_1}, m_{\tilde\chi^0_1})=$ (900,1) , (700,300), (550,375) in the Signal region D 1250. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(m_{\tilde{t}_1}, m_{\tilde\chi^0_1})=$ (900,1) , (700,300), (550,375) in the Signal region D 1500. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(m_{\tilde{t}_1}, m_{\tilde\chi^0_1})=$ (900,1) , (700,300), (550,375) in the Signal region D 1750. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(m_{\tilde{t}_1}, m_{\tilde\chi^0_1})=$ (900,1) , (700,300), (550,375) in the Signal region D 2000. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Signal acceptance in the SRA$^{[450,575]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRA$^{[450,575]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRA$^{\geq 575}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRA$^{\geq 575}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRB$^{[100,150]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRB$^{[100,150]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRB$^{[150,400]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRB$^{[150,400]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRB$^{\geq 400}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRB$^{\geq 400}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRC$^{[100,150]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRC$^{[100,150]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRC$^{[150,300]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRC$^{[150,300]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRC$^{[300,500]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRC$^{[300,500]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRC$^{\geq 500}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRC$^{\geq 500}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRD750 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRD750 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRD1000 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRD1000 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRD1250 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRD1250 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRD1500 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRD1500 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRD1750 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRD1750 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRD2000 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRD2000 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Observed 95% CL exclusion limits on the simplified models of chargino-pair production decaying via W bosons.

$+1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino-pair production decaying via W bosons.

$-1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino-pair production decaying via W bosons.

Expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$+1\sigma$ expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$-1\sigma$ expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

Observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$+1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$-1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

Expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and h bosons.

$+1\sigma$ expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and h bosons.

$-1\sigma$ expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and h bosons.

Observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and h bosons.

$+1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and h bosons.

$-1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and h bosons.

Expected 95% CL exclusion limits on the simplified models of higgsino GGM scenarios.

$+1\sigma$ expected 95% CL exclusion limits on the simplified models of higgsino GGM scenarios.

$-1\sigma$ expected 95% CL exclusion limits on the simplified models of higgsino GGM scenarios.

Observed 95% CL exclusion limits on the simplified models of higgsino GGM scenarios.

$+1\sigma$ observed 95% CL exclusion limits on the simplified models of higgsino GGM scenarios.

$-1\sigma$ observed 95% CL exclusion limits on the simplified models of higgsino GGM scenarios.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\chi}_1^{+}\tilde{\chi}_{1}^{-}$.

The analyses used in combination for each scenario to set limits in models of the production of $\tilde{\chi}_1^{+}\tilde{\chi}_{1}^{-}$.

Observed upper limit on the signal cross section in fb for chargino--neutralino production decaying via W and Z bosons.

The analyses used in combination for each scenario to set limits in models of chargino--neutralino production decaying via W and Z bosons.

Expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$+1\sigma$ expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$-1\sigma$ expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

Observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$+1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$-1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

Observed upper limit on the signal cross section in fb for chargino--neutralino production decaying via W and h bosons.

The analyses used in combination for each scenario to set limits in models of chargino--neutralino production decaying via W and h bosons.

Observed upper limit on the signal cross section in fb for higgsino GGM scenarios.

The analyses used in combination for each scenario to set limits in higgsino GGM scenarios.

Upper limit on the WIMP-nucleus scattering cross section from the single-scatter analysis.

Total acceptance evaluated on a simulated dataset of MIMP events. The mean and upper and lower bounds of the 1 sigma uncertainty band are provided.

The charged multiplicity normalzied excess yield of 54.4 GeV in 0-80% centrality.

T vs. uB from 27 GeV and 54.4 GeV

Data compared to the SM prediction in each region of the positive \lam{} analysis, after the fit to data. The shaded bands represent the total post-fit uncertainty on the prediction.

Fit results in the ($\kappa_{Z}$, $\kappa_{W}$) plane, using the negative-$\lambda_{WZ}$ analysis. The results are overlaid with the confidence regions (shown in blue) from a separate fit using a combination of the Higgs boson measurements collected in Nature 607 (2022) 52. This fit assumes that all Higgs boson couplings besides the ones represented are positive, and that only SM particles contribute to loop processes. Confidence regions are constructed from the log-likelihood ratio $\Lambda_{LR}=-2\ln(L/L_{max})$, where $L_{\max}$ is the likelihood for the best-fit point, which is shown as a black dot for the VBF $WH$ analysis or a blue dot for the Higgs combination. The $1\sigma$, $2\sigma$, and $5\sigma$ regions are defined by $\Lambda_{LR}$ values smaller than 2.30, 6.18, and 28.7, respectively. The SM value is marked with a star, while green triangles mark the points with $\kappa_{Z} = \pm1$, $\kappa_{W} = \mp1$.

The post-fit BDT score distribution for the direct stau channel, showing the scores for BDT4, before the selections on the BDT score is made. The black arrow depicts the BDT score selection for the SR-BDT. A few example SUSY scenarios targeted by each BDT are overlaid for illustration.

The post-fit kinematic distribution for the Intermediate stau channel, showing the $m_{\mathrm{T}2}$ distribution in SR-C1C1-LM, before the selection on $m_{\mathrm{T}2}$ is made. The black arrow depicts the selection for the signal region. An example SUSY scenario is overlaid for illustration.

The post-fit kinematic distribution for the Intermediate stau channel, showing the $m_{\mathrm{T}2}$ distribution in SR-C1N2OS-LM, before the selection on $m_{\mathrm{T}2}$ is made. The black arrow depicts the selection for the signal region. An example SUSY scenario is overlaid for illustration.

The post-fit kinematic distribution for the Intermediate stau channel, showing the $m_{\mathrm{T}2}$ distribution in SR-C1N2SS-LM, before the selection on $m_{\mathrm{T}2}$ is made. The black arrow depicts the selection for the signal region. An example SUSY scenario is overlaid for illustration.

The post-fit kinematic distribution for the Intermediate $Wh$ channel, showing the $m_{\mathrm{T}2}$ distribution in SR-Wh-LM, before the selection on $m_{\mathrm{T}2}$ is made. The black arrow depicts the selection for the signal region. An example SUSY scenario is overlaid for illustration.

The post-fit kinematic distribution for the Intermediate $Wh$ channel, showing the $m_{\mathrm{Tsum}}$ distribution in SR-Wh-HM, before the selection on $m_{\mathrm{Tsum}}$ is made. The black arrow depicts the selection for the signal region. An example SUSY scenario is overlaid for illustration.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$ production.

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$ production.

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$ production.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$ production.

$+1\sigma$ observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$ production.

$-1\sigma$ observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$ production.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L}\tilde{\tau}_{L}$ production.

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L}\tilde{\tau}_{L}$ production.

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L}\tilde{\tau}_{L}$ production.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L}\tilde{\tau}_{L}$ production.

$+1\sigma$ observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L}\tilde{\tau}_{L}$ production.

$-1\sigma$ observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L}\tilde{\tau}_{L}$ production.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{R}\tilde{\tau}_{R}$ production.

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{R}\tilde{\tau}_{R}$ production.

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{R}\tilde{\tau}_{R}$ production.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{R}\tilde{\tau}_{R}$ production.

$+1\sigma$ observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{R}\tilde{\tau}_{R}$ production.

$-1\sigma$ observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{R}\tilde{\tau}_{R}$ production.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ production decaying via staus.

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ production decaying via staus.

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ production decaying via staus.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ production decaying via staus.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$).

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$).

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$).

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$).

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$) using the OS signal regions.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$) using the OS signal regions.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$) using the SS signal regions.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$) using the SS signal regions.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$).

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$).

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$).

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$).

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$) using the OS signal regions.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$) using the OS signal regions.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$) using the SS signal regions.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$) using the SS signal regions.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$.

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$.

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\tau}_{L}\tilde{\tau}_{L}$.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\tau}_{R}\tilde{\tau}_{R}$.

The best expected signal region used to set limits in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The best expected signal region used to set limits in models of the production of $\tilde{\tau}_{L}\tilde{\tau}_{L}$.

The best expected signal region used to set limits in models of the production of $\tilde{\tau}_{R}\tilde{\tau}_{R}$.

The acceptance $A$ of the direct stau production signal region SR-BDT1 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The efficiency $\epsilon$ of the direct stau production signal region SR-BDT1 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The acceptance $A$ of the direct stau production signal region SR-BDT2 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The efficiency $\epsilon$ of the direct stau production signal region SR-BDT2 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The acceptance $A$ of the direct stau production signal region SR-BDT3 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The efficiency $\epsilon$ of the direct stau production signal region SR-BDT3 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The acceptance $A$ of the direct stau production signal region SR-BDT4 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The efficiency $\epsilon$ of the direct stau production signal region SR-BDT4 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

Cutflow event yields of the direct stau production SRs in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$. All yields correspond to weighted events, so that effects from reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalised to the integrated luminosity of the data sample, $\int L dt = 139$ fb$^{-1}$. The preliminary event reduction is a centralised stage requiring the event to: have at least two hadronically-decaying taus; be selected and matched by the asymmetric ditau trigger (with plateau cuts); have at least two medium taus of opposite sign and have $m_{\mathrm{T2}} > 30$ GeV. Cuts except the BDT score cuts are applied consecutively, while the BDT score cuts are applied one at a time.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$.

The acceptance $A$ of the Intermediate stau signal region SR-C1C1-LM in models of the production of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ decaying via staus.

The efficiency $\epsilon$ of the Intermediate stau signal region SR-C1C1-LM in models of the production of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ decaying via staus.

The acceptance $A$ of the Intermediate stau signal region SR-C1C1-HM in models of the production of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ decaying via staus.

The efficiency $\epsilon$ of the Intermediate stau signal region SR-C1C1-HM in models of the production of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ decaying via staus.

The acceptance $A$ of the Intermediate stau signal region SR-C1N2OS-LM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The efficiency $\epsilon$ of the Intermediate stau signal region SR-C1N2OS-LM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The acceptance $A$ of the Intermediate stau signal region SR-C1N2SS-LM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The efficiency $\epsilon$ of the Intermediate stau signal region SR-C1N2SS-LM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The acceptance $A$ of the Intermediate stau signal region SR-C1N2OS-HM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The efficiency $\epsilon$ of the Intermediate stau signal region SR-C1N2OS-HM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The acceptance $A$ of the Intermediate stau signal region SR-C1N2SS-HM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The efficiency $\epsilon$ of the Intermediate stau signal region SR-C1N2SS-HM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

Cutflow event yields of the Intermediate stau SRs in models of the production of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$. All yields correspond to weighted events, so that effects from reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalised to the integrated luminosity of the data sample, $\int L dt = 139$ fb$^{-1}$. The preliminary event reduction is a centralised stage requiring the event to: have at least two hadronically-decaying medium taus; be selected and matched by the asymmetric ditau trigger or ditau+MET trigger with offline cuts using HLT threshold values. The different yields of preliminary event reduction are caused by different trigger scale factor.

Cutflow event yields of the Intermediate stau SRs in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$. All yields correspond to weighted events, so that effects from reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalised to the integrated luminosity of the data sample, $\int L dt = 139$ fb$^{-1}$. The preliminary event reduction is a centralised stage requiring the event to: have at least two hadronically-decaying medium taus; be selected and matched by the asymmetric ditau trigger or ditau+MET trigger with offline cuts using HLT threshold values. The different yields of preliminary event reduction are caused by different trigger scale factor.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$.

The best expected signal region used to set limits in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$.

The acceptance $A$ of the Intermediate $Wh$ signal region SR-Wh-LM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$.

The efficiency $\epsilon$ of the Intermediate $Wh$ signal region SR-Wh-LM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$.

The acceptance $A$ of the Intermediate $Wh$ signal region SR-Wh-HM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$.

The efficiency $\epsilon$ of the Intermediate $Wh$ signal region SR-Wh-HM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$.

Cutflow event yields of the Intermediate $Wh$ SRs in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$. All yields correspond to weighted events, so that effects from reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalised to the integrated luminosity of the data sample, $\int L dt = 139$ fb$^{-1}$. The preliminary event reduction is a centralised stage requiring the event to: have one single lepton and two hadronically decaying medium taus; be selected and matched by the single lepton trigger with plateau cuts.