The impact of the eight most important systematic uncertainties on \(R_t\), in %.

The post-fit yields in the two SRs. All uncertainties applied in the analysis are included.

Limits on EFT parameters and Vtb extracted from the analysis. The upper and lower limits correspond to 95% CL.

Results of the main analysis.

Pre-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR plus.

Post-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR plus.

Post-fit distribution of the \(\Delta\phi(\vec{p}_{T}^{miss},\mu)\) variable in the muon-plus CR.

Pre-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR minus.

Post-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR minus.

Post-fit distribution of the \(\Delta\phi(\vec{p}_{T}^{miss},\mu)\) variable in the muon-minus CR.

Pre-fit distribution of the invariant mass of the untagged jet and the b-tagged jet, \(m(jb)\), in SR plus.

Pre-fit distribution of the absolute value of the pseudorapidity of the untagged jet, \(|\eta(j)|\), in SR plus.

Pre-fit distribution of the absolute value of the difference in \(p_{\text{T}}\) between the reconstructed W boson and the jet pair, \(|\Delta p_{\text{T}}(W, jb)|\), in SR plus.

Pre-fit distribution of the difference in azimuth angle between the reconstructed W boson and the jet pair, \(|\Delta\phi(W, jb)|\), in SR plus.

Pre-fit distribution of the invariant mass of the reconstructed top quark, \(m(t)\), in SR plus.

Pre-fit distribution of the absolute value of the difference in pseudorapidity between the charged lepton and the untagged jet, \(|\Delta\eta(\ell ,j)|\), SR plus.

Pre-fit distribution of the angular distance of the charged lepton and the untagged jet, \(\Delta R(\ell ,j)\), in SR plus.

Pre-fit distribution of the absolute value of the difference in pseudorapidity between the b-tagged jet and the charged lepton, \(|\Delta\eta(b,\ell )|\), in SR plus.

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, 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, 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 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 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 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 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 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 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) = 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) = 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) = 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) = 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) = 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) = 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}) = 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}) = 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}) = 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.

Corrected $m(HH)$ distribution in the resolved $4b$ signal region (dots), after the fit under the background-only hypothesis. The error bars on the $4b$ points represent the Poisson uncertainties corresponding to their event yields. The background model (teal histogram) is shown with its total post-fit uncertainty (gray band). The final bin includes overflow. Representative spin-0 signal hypotheses (dashed, dotted, and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. The bottom panel shows the difference between the $4b$ distribution and the background model, relative to the background model. No significant excess of data relative to the SM background is observed.

Corrected $m(HH)$ distribution in the resolved $4b$ signal region (dots), after the fit under the background-only hypothesis. The error bars on the $4b$ points represent the Poisson uncertainties corresponding to their event yields. The background model (teal histogram) is shown with its total post-fit uncertainty (gray band). The final bin includes overflow. Representative spin-2 signal hypotheses (dashed, dotted, and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. The bottom panel shows the difference between the $4b$ distribution and the background model, relative to the background model. No significant excess of data relative to the SM background is observed.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The steps up to the $b$-tag categorization are shown for the spin-0.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The steps up to the $b$-tag categorization are shown for the spin-2.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The efficiencies of the three b-tag categories are shown for the spin-0 scenario; this efficiency is obtained after the other selection steps including the SR definition. The signal efficiency in the 4b region has a maximum around 1.5 TeV. Above that value the track-jets start to merge together, and for the highest resonance masses the 2b category becomes the most efficient.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The efficiencies of the three b-tag categories are shown for the spin-2 scenario; this efficiency is obtained after the other selection steps including the SR definition. The signal efficiency in the 4b region has a maximum around 1.5 TeV. Above that value the track-jets start to merge together, and for the highest resonance masses the 2b category becomes the most efficient.

Comparison of the background model (stacked histograms) with data (dots) in the $2b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

Comparison of the background model (stacked histograms) with data (dots) in the $3b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

Comparison of the background model (stacked histograms) with data (dots) in the $4b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

The $m(HH)$ distributions in the boosted $2b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $2b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $3b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $3b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $4b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $4b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

Expected (dashed black lines) and observed (solid black lines) 95% CL upper limits on the cross-section of resonant $HH$ production in the spin-0 signal models. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits (colored bands) are shown. Expected limits using each of the resolved and boosted channels individually (dashed colored lines) are shown. The nominal $H\rightarrow b\bar{b}$ branching ratio is taken as 0.582.

Expected (dashed black lines) and observed (solid black lines) 95% CL upper limits on the cross-section of resonant $HH$ production in the spin-2 signal models. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits (colored bands) are shown. Expected limits using each of the resolved and boosted channels individually (dashed colored lines) are shown. The theoretical prediction for the bulk RS model with $k/\bar{M}_{\text{Pl}} = 1$ (solid red line) is shown; the decrease below 350 GeV is due to a sharp reduction in the $G^{*}_{\text{KK}} \rightarrow HH$ branching ratio. The nominal $H\rightarrow b\bar{b}$ branching ratio is taken as 0.582.

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the High-pT jet SR for observed data events

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the High-pT jet SR for expected signal events in the strong gluino pair pair production model with m(gluino)=1.8 TeV, m(chi0)=0.2 TeV, tau(chi0)=0.1 ns

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the Trackless jet SR for observed data events

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the Trackless jet SR for expected signal events in the electroweak pair production model

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Exclusion limits at 95% CL on the production cross section in the electroweak pair production model.

Exclusion limits at 95% CL on the production cross section in the strong gluino pair production models and m(gluino)=2.4 TeV

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Expected exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Observed exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Expected exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Observed exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Exclusion limits at 95% CL on the production cross section in the strong gluino pair production models and m($ ilde{\chi}^0_1$)=1.25 TeV

Acceptance cutflow for the High-pT SR for representative points in the strong gluino pair production model. See additional resources for more information.

Acceptance cutflow for the Trackless SR for representative points in the electroweak pair production model. See additional resources for more information.

Acceptance cutflow for the Trackless SR for representative points in the electroweak pair production model with heavy-flavor quarks final state. See additional resources for more information.

Acceptance cutflow for the High-pT SR for representative points in the electroweak pair production model with heavy-flavor quarks final state. See additional resources for more information.

Reinterpretation Material: Event-level Efficiency for HighPt SR selections, R < 1150 mm

Reinterpretation Material: Event-level Efficiency for HighPt SR selections, R [1150, 3870] mm

Reinterpretation Material: Event-level Efficiency for HighPt SR selections, R > 3870 mm

Reinterpretation Material: Event-level Efficiency for Trackless SR selections, R < 1150 mm

Reinterpretation Material: Event-level Efficiency for Trackless SR selections, R [1150, 3870] mm

Reinterpretation Material: Event-level Efficiency for Trackless SR selections, R > 3870 mm

Reinterpretation Material: Vertex-level Efficiency for R < 22 mm

Reinterpretation Material: Vertex-level Efficiency for R [22, 25] mm

Reinterpretation Material: Vertex-level Efficiency for R [25, 29] mm

Reinterpretation Material: Vertex-level Efficiency for R [29, 38] mm

Reinterpretation Material: Vertex-level Efficiency for R [38, 46] mm

Reinterpretation Material: Vertex-level Efficiency for R [46, 73] mm

Reinterpretation Material: Vertex-level Efficiency for R [73, 84] mm

Reinterpretation Material: Vertex-level Efficiency for R [84, 111] mm

Reinterpretation Material: Vertex-level Efficiency for R [111, 120] mm

Reinterpretation Material: Vertex-level Efficiency for R [120, 145] mm

Reinterpretation Material: Vertex-level Efficiency for R [145, 180] mm

Reinterpretation Material: Vertex-level Efficiency for R [180, 300] mm

Cutflow (acceptance x efficiency) for the High-pT SR for representative points in the strong gluino pair production model. See additional resources for more information.

Cutflow (acceptance x efficiency) for the Trackless SR for representative points in the electroweak pair production model. See additional resources for more information.

Cutflow (acceptance x efficiency) for the Trackless SR for representative points in the electroweak pair production model with heavy-flavor quarks. See additional resources for more information.

Cutflow (acceptance x efficiency) for the High-pT SR for representative points in the electroweak pair production model with heavy-flavor quarks. See additional resources for more information.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Expected and observed 95% CL upper limits on the branching ratio of the Higgs boson decay to $Za$ times the branching ratio $a $->$ \gamma \gamma$ as a function of the $a$ particle mass in the merged ($m_a \le 2$ GeV) and the resolved ($m_a > 2$ GeV) categories. The figure uses $pp$ collision data at $\sqrt{s}=13$ TeV corresponding to 139 fb$^{-1}$.

ATLAS observed 95% CL exclusion contours limits in terms of the ALP mass and its effective coupling to photons, $|C_{\gamma\gamma}|/\Lambda$, for different values of the Higgs coupling to $Za$, $|C_{Zh}|/\Lambda$. The figure uses $pp$ collision data at $\sqrt{s}=13$ TeV corresponding to 139 fb$^{-1}$. The value 100 in the tables is presentational, denoting the upper edge of the plot.

Bayesian upper exclusion limits on the ratio $g_{W'}/g_{W}$ for an SSM-like W$\prime$ boson are shown. The unity coupling ratio (blue dotted curve) corresponds to the SSM common benchmark. The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian lower exclusion limits on the NUGIM G(221) mixing angle $\cot(\theta_{E})$ are shown as a function of the W$\prime$ boson mass. The theoretically excluded region is shaded in grey. The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian upper exclusion limits at 95% CL on the product of the production cross section and branching fraction of a QBH in an associated $ au$ lepton and neutrino final state. Minimum threshold masses $m_{th}$ of up to 6.6 TeV are excluded at 95% CL. The observed limit (solid line) is obtained from the intersection with the LO QBH cross section (blue dotted curve). The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian upper limits at 95% CL on the cross section of the process $pp\rightarrow\tau\nu$ mediated via LQ exchange in the t-channel. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The predicted LQ cross section at LO in the three coupling benchmark scenarios is depicted in different colors for $g_{U}=1$. The uncertainty bandy correspond to the sum in quadrature of PDF and scale variations. The first benchmark scenario considers only couplings to left-handed SM fermions (i.e. $\beta_{\text{R}}^{ij}=0$) and is referred to as "best fit LH". The second benchmark, referred to as "best fit LH+RH", considers $|\beta_{\text{R}}^{\text{b}\tau}|=1$ and all other $\beta_{\text{R}}^{ij}=0$. In the third "democratic" benchmark, equal couplings only to LH fermions are assumed, i.e. $\beta_{\text{L}}^{ij}=1$ and $\beta_{\text{R}}^{ij}=0$ for all $i$ and $j$.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the LH scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the LH+RH scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the democratic scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Bayesian upper exclusion limits at 95% CL on each of the Wilson coefficients described by the EFT model based on 2016-2018 data. The three different coupling types represent a left-handed vector coupling ($\epsilon^{cb}_{L}$), tensor-like coupling ($\epsilon^{cb}_{T}$), and scalar-tensor-like coupling ($\epsilon^{cb}_{S_{L}}$). The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Summary of exclusion limits (expected and observed) calculated at 95% CL for full Run-2 CMS data.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "year_binnumber", following the binning from Figure 4.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "year_binnumber", following the binning used in Figure 4.

Predicted signal yields for the 2017 data-taking period, corresponding to 41.3 fb$^{-1}$, after the application of each search requirement (cumulative) for various signal hypothesis. The requirements listed are presented as total efficiencies w.r.t. the previous selection step.

Expected exclusion limits at 95% CL from Fig 7(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Expected exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Expected exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRGGWZ-H.

N-1 distribution for $E_{\mathrm{T}}^{\mathrm{miss}}$of observed data and expected background in SRGGSlep-M.

N-1 distribution for $\sum{p_{\mathrm{T}}^\mathrm{jet}}$of observed data and expected background in SRUDD-ge2b.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRLQD.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRSSWZ-H.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRSSSlep-H(loose).

Signal acceptance for SRGGWZ-H signal region from Fig 10(c) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-H signal region from Fig 15(c) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGWZ-M signal region from Fig 10(b) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-M signal region from Fig 15(b) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGWZ-L signal region from Fig 10(a) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-L signal region from Fig 15(a) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-L signal region from Fig 12(a) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-L signal region from Fig 17(a) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-M signal region from Fig 12(b) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-M signal region from Fig 17(b) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-H signal region from Fig 12(c) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-H signal region from Fig 17(c) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRUDD-1b signal region from Fig 14(b) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-1b signal region from Fig 19(b) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-2b signal region from Fig 14(c) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-2b signal region from Fig 19(c) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-ge2b signal region from Fig 14(d) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-ge2b signal region from Fig 19(d) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-ge3b signal region from Fig 14(e) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-ge3b signal region from Fig 19(e) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRLQD signal region from Fig 14(a) in a SUSY scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Signal efficiency for SRLQD signal region from Fig 19(a) in a SUSY scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Signal acceptance for SRSSWZ-L signal region from Fig 11(a) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-L signal region from Fig 16(a) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-ML signal region from Fig 11(b) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-ML signal region from Fig 16(b) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-MH signal region from Fig 11(c) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-MH signal region from Fig 16(c) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-H signal region from Fig 11(d) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-H signal region from Fig 16(d) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-H signal region from Fig 13(d) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-H signal region from Fig 18(d) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-MH signal region from Fig 13(c) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-MH signal region from Fig 18(c) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-L signal region from Fig 13(a) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-L signal region from Fig 18(a) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-ML signal region from Fig 13(b) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-ML signal region from Fig 18(b) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-H(loose) signal region from Fig 13(e) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-H(loose) signal region from Fig 18(e) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-H in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-M in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-L in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGSlep-L in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGSlep-M in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGSlep-H in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-1b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-2b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-ge2b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-ge3b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRLQD in a susy scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2200 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1870 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-L in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-ML in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-MH in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-H in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-H in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-MH in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-L in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-ML in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-H(loose) in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Cross-section upper limits at 95% CL from Fig1(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Cross-section upper limits at 95% CL from Fig1(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Cross-section upper limits at 95% CL from Fig1(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

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.

SM background observation/prediction in the WH signal region

SM background prediction vs. observation in the W signal region

SM background observation/prediction in the W signal region

SM background prediction vs. observation in the H signal region

SM background observation/prediction in the H signal region

Observed exclusion limits assuming the approximate-NLO+NLL cross sections

Expected exclusion limits assuming the approximate-NLO+NLL cross sections

The 95% CL observed upper limits on the production cross sections for $\widetilde{\chi}^\pm_1$ $\widetilde{\chi}^\mp_1$ assuming that each $\widetilde{\chi}^\pm_1$ decays to a W boson and $\widetilde{\chi}^0_1$

Observed exclusion limits assuming the approximate-NLO+NLL cross sections

Expected exclusion limits assuming the approximate-NLO+NLL cross sections

The 95% CL observed upper limits on the production cross sections for $\widetilde{\chi}^\pm_1$ $\widetilde{\chi}^0_2$ assuming that each $\widetilde{\chi}^\pm_1$ decays to a W boson and $\widetilde{\chi}^0_1$ and the $\widetilde{\chi}^0_2$ decays to a Z boson and $\widetilde{\chi}^0_1$

Observed exclusion limits assuming the approximate-NLO+NLL cross sections

Expected exclusion limits assuming the approximate-NLO+NLL cross sections

The 95% CL observed upper limits on the production cross sections for $\widetilde{\chi}^\pm_1$ $\widetilde{\chi}^0_2$ assuming that each $\widetilde{\chi}^\pm_1$ decays to a W boson and $\widetilde{\chi}^0_1$ and the $\widetilde{\chi}^0_2$ decays to a H boson and $\widetilde{\chi}^0_1$

Observed exclusion limits assuming the approximate-NLO+NLL cross sections

Expected exclusion limits assuming the approximate-NLO+NLL cross sections

Observed exclusion limits assuming the approximate-NLO+NLL cross sections

Expected exclusion limits assuming the approximate-NLO+NLL cross sections

Observed exclusion limits assuming the approximate-NLO+NLL cross sections

Expected exclusion limits assuming the approximate-NLO+NLL cross sections

The 95% CL observed upper limits on the production cross sections for mass-degenerate higgsino-like$\widetilde{\chi}^\pm_1$ $\widetilde{\chi}^\mp_1$, $\widetilde{\chi}^\pm_1$ $\widetilde{\chi}^0_2$, $\widetilde{\chi}^\pm_1$ $\widetilde{\chi}^0_3$ and $\widetilde{\chi}^0_2$ $\widetilde{\chi}^0_3$ as functions of the NLSP and LSP masses.

Efficiency of bb-tagger for H(bb), Z(bb) and Z(cc) decays.

Efficiency of W- and V-tagger for W(qq) and Z(qq) decays.

Acceptance times efficiency values with statistical uncertainties for TChiWW in the b-Veto region.

Acceptance times efficiency values with statistical uncertainties for TChiWZ in the b-Veto region.

Acceptance times efficiency values with statistical uncertainties for TChiWH in the b-Veto region.

Acceptance times efficiency values with statistical uncertainties for TChiHZ in the b-Veto region.

Acceptance times efficiency values with statistical uncertainties for TChiWW in the WHSR region.

Acceptance times efficiency values with statistical uncertainties for TChiWW in the WSR region.

Acceptance times efficiency values with statistical uncertainties for TChiWW in the HSR region.

Acceptance times efficiency values with statistical uncertainties for TChiWZ in the WHSR region.

Acceptance times efficiency values with statistical uncertainties for TChiWZ in the WSR region.

Acceptance times efficiency values with statistical uncertainties for TChiWZ in the HSR region.

Acceptance times efficiency values with statistical uncertainties for TChiWH in the WHSR region.

Acceptance times efficiency values with statistical uncertainties for TChiWH in the WSR region.

Acceptance times efficiency values with statistical uncertainties for TChiWH in the HSR region.

Acceptance times efficiency values with uncertainties for TChiHZ in the WHSR region.

Acceptance times efficiency values with uncertainties for TChiHZ in the WSR region.

Acceptance times efficiency values with statistical uncertainties for TChiHZ in the HSR region.

Covariance matrix for the signal regions, derived from a fit to the control regions only under the background-only hypothesis.

Correlation matrix for the signal regions, derived from a fit to the control regions only under the background-only hypothesis.

Post-fit distribution of the $M(ll)$ variable for the ultrahigh-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '2l soft' signal region of the '2/3l soft' analysis.

Post-fit distribution of the $M(ll)$ variable for the low-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '3l soft' signal region of the the '2/3l soft' analysis.

Post-fit distribution of the $M(ll)$ variable for the medium-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '3l soft' signal region of the the '2/3l soft' analysis.

Post-fit distribution of the $m_{\mathrm{T2}}(ll)$ variable for low-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '2l soft' signal region of the '2/3l soft' analysis.

Post-fit distribution of the $m_{\mathrm{T2}}(ll)$ variable for medium-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '2l soft' signal region of the '2/3l soft' analysis.

Post-fit distribution of the $m_{\mathrm{T2}}(ll)$ variable for high-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '2l soft' signal region of the '2/3l soft' analysis.

Post-fit distribution of the $m_{\mathrm{T2}}(ll)$ variable for ultrahigh-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '2l soft' signal region of the '2/3l soft' analysis.

2SS $\ell/{\geq}\,3\ell$ search: observed and expected yields across the SRs in category A, events with three light leptons of which at least two form an OSSF pair, after the requirement that the leading-lepton $p_{\mathrm{T}}$ be greater than 30 GeV is applied.

2SS $\ell/{\geq}\,3\ell$ search: observed and expected yields across the SRs of the '${\geq}\ 3\ell$' search in category B, events with three light leptons and no OSSF pair, after the requirement that the leading-lepton $p_{\mathrm{T}}$ be greater than 30 GeV is applied.

Wino-bino model: cross section limits in the model parameter space, for wino-like chargino-neutralino production in the WZ topology for the full parameter space.

Wino-bino model: cross section limits in the model parameter space, for wino-like chargino-neutralino production in the WZ topology for the compressed space.

Wino-bino model: cross section limits in the model parameter space, for wino-like chargino-neutralino production in the WH topology for the full parameter space.

Wino-bino model: cross section limits in the model parameter space, for wino-like chargino-neutralino production with mixed topology with equal branching fraction to WZ and WH.

Wino-bino model: exclusion contours from the individual and combined analyses targeting WZ topology for the full parameter space. For visualization of the exclusion contours, linear interpolation is employed to account for the limited granularity of the available signal samples.

Wino-bino model: exclusion contours from the individual and combined analyses targeting the corresponding compressed region. For visualization of the exclusion contours, linear interpolation is employed to account for the limited granularity of the available signal samples.

Wino-bino model: exclusion contours from the individual and combined analyses targeting the WH topology for the full parameter space. For visualization of the exclusion contours, linear interpolation is employed to account for the limited granularity of the available signal samples.

Wino-bino model: exclusion contours from the individual and combined analyses targeting combined contours for these two topologies. For visualization of the exclusion contours, linear interpolation is employed to account for the limited granularity of the available signal samples.

GMSB model: expected and observed cross section limits for the neutralino-neutralino production for the ZZ topology.

GMSB model: expected and observed cross section limits for the neutralino-neutralino production for the HH topology.

GMSB model: expected and observed cross section limits for the neutralino-neutralino production for the mixed topology with equal branching fraction to H and Z.

GMSB model: cross section limits for neutralino-neutralino production as a function of the NSLP mass and the branching fraction to the H boson for the combination of the searches.

GMSB model: exclusion limit for neutralino-neutralino production as a function of the NSLP mass and the branching fraction to the H boson for the combination of the searches along with the input searches. For visualization of the exclusion contours, linear interpolation is employed to account for the limited granularity of the available signal samples.

Cross section upper limit(s) in the mass plane of NLSP and LSP masses for the higgsino-bino model.

Mass plane cross section upper limit for direct slepton pair production, with observed and expected exclusion limits in the full mass plane from the combination.

Mass plane cross section upper limit for direct slepton pair production, with observed and expected exclusion limits in the compressed region from '2/3l' soft search.

Distribution of an input to the BDT classifier in the $3\ell$ category: The linear discrimminant $p_{T}^{miss,LD}$. The discriminant is defined by the relation $p_{T}^{miss,LD} = 0.6 p_{T}^{miss} + 0.4 H_{T}^{miss}$, where $H_{T}^{miss}$ corresponds to the magnitude of the vector $p_{T}$ sum of all electrons, muons, taus and small radius jets passing the selection critertia, and $p_{T}^{miss}$ represents the missing transverse momentum vector computed as the negative vector $p_{T}$ sum of all the particles reconstructed by the PF algorithm in the event.

Distributions of the transverse mass $m_{T}$ of the lepton that is not originating from the Z boson decay and the missing transverse momentum in the $3\ell$/WZ control region. The distributions expected for the WZ and ZZ as well as for other background processes are shown for the values of nuisance parameters obtained from the ML fit used in the signal extraction.

Distributions of the four lepton invariant mass $m_{4\ell}$ in the $4\ell$/ZZ control region. The distributions expected for the ZZ as well as for other background processes are shown for the values of nuisance parameters obtained from the ML fit used in the signal extraction.

Distributions of the transverse mass $m_{T}$ of the leading lepton and the missing transverse momentum in the $2\ell$(ss) control region. The distributions expected for the misidentified $\ell$ background as well as for other background processes is shown for the values of nuisance parameters obtained from the background-only ML fit, in which the HH signal is constrained to be zero.

Distributions of the reconstructed HH mass $m_{HH}$ in the $2\ell+2 au_{h}$ control region. The distributions expected for the misidentified $\ell/ au_{h}$ background as well as for other background processes is shown for the values of nuisance parameters obtained from the background-only ML fit, in which the HH signal is constrained to be zero.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $2\ell$(ss) category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $3\ell$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $4\ell$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $3\ell+1 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $2\ell+2 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $1\ell+3 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $4 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Observed and expected $95\%$ CL upper limits on the SM HH production cross section, obtained for both individual search categories and from a simultaneous fit of all seven categories combined.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the Higgs boson self-coupling strength modifier $\kappa_\lambda$ for the combination of all seven categories. All Higgs boson couplings other than $\lambda$ are assumed to have the values predicted in the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the Higgs boson self-coupling strength modifier $\kappa_\lambda$ for the seven different categories as well as their combination. All Higgs boson couplings other than $\lambda$ are assumed to have the values predicted in the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section for the twelve benchmark scenarios from doi:10.1007/JHEP04(2016)126, the additional benchmark scenario 8a from doi:10.1007/JHEP09(2018)057, the seven benchmark scenarios from doi:10.1007/JHEP03(2020)091, and for the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section for the twelve benchmark scenarios from doi:10.1007/JHEP04(2016)126, the additional benchmark scenario 8a from doi:10.1007/JHEP09(2018)057, the seven benchmark scenarios from doi:10.1007/JHEP03(2020)091, and for the SM. The limits are shown for all seven search categories as well as for their combination.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the effective coupling $c_{2}$ for the combination of all seven categories. All Higgs boson couplings other than $c_{2}$ are assumed to have the values predicted in the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the effective coupling $c_{2}$ and the top Yukawa coupling modifier $\kappa_{t}$ for the combination of all seven categories. All Higgs boson couplings other than $c_{2}$ and $\kappa_{t}$ are assumed to have the values predicted in the SM. The position of the SM in the ($c_{2}-\kappa_{t}$) plane, as well as the best fit value of $(c_{2},\kappa_{t})=(1.05, 1.74)$ together with contours for the theory cross section are shown as well.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the Higgs boson self-coupling modifier $\kappa_{\lambda}$ and the top Yukawa coupling modifier $\kappa_{t}$ for the combination of all seven categories. All Higgs boson couplings other than $\kappa_{\lambda}$ and $\kappa_{t}$ are assumed to have the values predicted in the SM. The position of the SM in the ($\kappa_{\lambda}-\kappa_{t}$) plane, as well as the best fit value of $(\kappa_{\lambda},\kappa_{t})=(-3.54, -1.70)$ together with contours for the theory cross section are shown as well.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the effective coupling $c_{2}$ and the Higgs boson self-coupling modifier $\kappa_{\lambda}$ for the combination of all seven categories. All Higgs boson couplings other than $c_{2}$ and $\kappa_{\lambda}$ are assumed to have the values predicted in the SM. The position of the SM in the ($c_{2}-\kappa_{\lambda}$) plane, as well as the best fit value of $(c_{2},\kappa_{\lambda})=(-0.66, 5.33)$ together with contours for the theory cross section are shown as well.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $0$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $0$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs. The Limit is shown for the seven different search categories as well as their combination.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $2$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $2$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs. The Limit is shown for the seven different search categories as well as their combination.

Distribution in the output of the BDT trained for resonances of spin 2 and mass $750\,$GeV production for the $2\ell$(ss) category. The resonant HH signal is shown for a cross section amounting to $1\,$pb. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for resonances of spin 2 and mass $750\,$GeV production for the $3\ell$ category. The resonant HH signal is shown for a cross section amounting to $1\,$pb. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Predicted and observed yields as a function of $m_{4\ell}$ in the VBS-Suppressed region. Overflow events are included in the last bin of the distribution.

Relative systematic uncertainties in the unfolded differential cross section for 4ljj production in the VBS-Enhanced region as a function of m_{jj}

Relative systematic uncertainties in the unfolded differential cross section for 4ljj production in the VBS-Enhanced region as a function of m_{4\ell}

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $m_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $m_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $m_{4\ell}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $m_{4\ell}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $ p_{T,4l}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $ p_{T,4l}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $\Delta\phi_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $\Delta\phi_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $|\Delta y_{jj}|$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $|\Delta y_{jj}|$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $cos(\theta^{*}_{12})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $cos(\theta^{*}_{12})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $cos(\theta^{*}_{34})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $cos(\theta^{*}_{34})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $p_{T,jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $p_{T,jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $p_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $p_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $S_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $S_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Statistical correlation between different bins of the unfolded cross-section for all the measured observables in the VBS-Enhanced region.

Statistical correlation between different bins of the unfolded cross-section for all the measured observables in the VBS-Suppressed region.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 5-jet bHbH channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 5-jet bHbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 5-jet bZbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 5-jet bHtW channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 5-jet bZtW channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 6-jet bHbH channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 6-jet bHbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 6-jet bZbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 6-jet bHtW channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 6-jet bZtW channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the semileptonic 3-jet bHbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the semileptonic 3-jet bZbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the semileptonic 4-jet bHbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the semileptonic 4-jet bZbZ channel.

The limit at 95% CL on the cross section for VLQ pair production for the branching fraction hypothesis 0% $\mathcal{B}(B \to bH)$, 100% $\mathcal{B}(B \to bH)$, and 0% $\mathcal{B}(B \to bH)$

The limit at 95% CL on the cross section for VLQ pair production for the branching fraction hypothesis 25% $\mathcal{B}(B \to bH)$, 25% $\mathcal{B}(B \to bH)$, and 50% $\mathcal{B}(B \to bH)$

The limit at 95% CL on the cross section for VLQ pair production for the branching fraction hypothesis 50% $\mathcal{B}(B \to bH)$, 50% $\mathcal{B}(B \to bH)$, and 0% $\mathcal{B}(B \to bH)$

The limit at 95% CL on the cross section for VLQ pair production for the branching fraction hypothesis 100% $\mathcal{B}(B \to bH)$, 0% $\mathcal{B}(B \to bH)$, and 0% $\mathcal{B}(B \to bH)$

Median expected exclusion limits on the VLQ mass at 95% CL as a function of the branching fractions $\mathcal{B}(B \to bH)$ and $\mathcal{B}(B \to tW)$, with $\mathcal{B}(B \to tW) = 1 - \mathcal{B}(B \to bH) - \mathcal{B}(B \to bZ)$. The grey area corresponds to the region where the exclusion limit is less than 1000 GeV.

Median observed exclusion limits on the VLQ mass at 95% CL as a function of the branching fractions $\mathcal{B}(B \to bH)$ and $\mathcal{B}(B \to tW)$, with $\mathcal{B}(B \to tW) = 1 - \mathcal{B}(B \to bH) - \mathcal{B}(B \to bZ)$. The grey area corresponds to the region where the exclusion limit is less than 1000 GeV.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2016 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 2.3%.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2017 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 3.7%.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2018 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 1.8%.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the high-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the high-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Results of the background-only fit in the low-mass channel discovery region SR_LM_150. Both pre-fit and post-fit values are shown.

Results of the background-only fit in the low-mass channel discovery region SR_LM_300. Both pre-fit and post-fit values are shown.

The experimental efficiency of the low-mass channel for the exclusion and discovery signal regions as a function of higgsino mass. The experimental efficiency is defined as the number of events passing the detector-level event selections divided by the number of events passing the event selections for a perfect detector. The denominator is obtained by implementing particle-level event selections that emulate the detector-level selections. This treats the lack of availability of $b$-jet triggers as an inefficiency.

The particle-level acceptance for the low-mass exclusion and discovery signal regions, shown as a function of higgsino mass. The acceptance is defined as the fraction of signal events passing the particle-level event selection that emulates the detector-level selection. The acceptance calculation considers only those signal events where both higgsinos decay to Higgs bosons.

The experimental efficiency of the high-mass channel discovery regions as a function of higgsino mass. For each higgsino mass, the efficiency is shown for the SR-1 region corresponding to the mass. For masses above 1100 GeV, SR-1-1100 is used. The experimental efficiency is defined as the number of events passing the detector-level event selections divided by the number of events passing the event selections for a perfect detector. The denominator is obtained by implementing particle-level event selections that emulate the detector-level selections. The efficiency calculation considers only those signal events where both higgsinos decay to Higgs bosons.

The particle-level acceptance for the high-mass signal regions, shown as a function of higgsino mass. For each higgsino mass, the acceptance is shown for the SR-1 region corresponding to the mass. For masses above 1100 GeV, SR-1-1100 is used. The acceptance is defined as the fraction of signal events passing the particle-level event selection that emulates the detector-level selection. The acceptance calculation considers only those signal events where both higgsinos decay to Higgs bosons.

Cutflow for the low-mass channel for a representative 130 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 150 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 250 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 600 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 700 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 800 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 900 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 1000 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 1100 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 250 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 600 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 700 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 800 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 900 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1000 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1100 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Correlation matrix for the differential branching fraction extraction between different $q^2$ bins in the simultaneous fit.

The product of acceptance and efficiency ($A\epsilon$) of the $\mathcal{B}(\mathrm{B}^+\to\mathrm{K}^+\mu^+\mu^-)$ channel, as a function of the muon pair $q^2$, as measured in simulated signal events, after all the corrections applied. In the figure, regions corresponding to resonances are displayed with red markers.

Same as Table 5, but for the bins containing the J$\psi$ and $\psi$(2S) resonances.

Integrated branching fraction $\mathcal{B}(\mathrm{B}^\pm \to \mathrm{K}^\pm\mu^+\mu^-)$ in the low-$q^2$ region.

Ratio of branching fractions $\mathcal{B}(\text{B}^\pm \to \text{K}^\pm\mu^+\mu^-)$ to $\mathcal{B}(\text{B}^\pm \to \text{K}^\pm\text{e}^+\text{e}^-)$, determined from the measured double ratio $R$(K) of these decays to the respective branching fractions of the decay chains $\text{B}^\pm\to\text{J/}\psi\text{K}^\pm$ with $\text{J/}\psi\to\mu^+\mu^-$ and $\text{e}^+\text{e}^-$.

Negative log likelihood function from the fit profiled as a function of the inverse ratio $R(\mathrm{K})^{-1}$.

Parametrization of the $A_{X}$ factor, defined as the fraction of diphoton resonances satisfying the fiducial acceptance at the particle level, as function of $m_{X}$ obtained from the narrow width simulated signal samples produced in gluon fusion.

The correction factor, $C_{X}$, defined as the ratio between the number of reconstructed signal events passing the analysis cuts and the number of signal events at the particle level generated within the fiducial volume, and acceptance correction factor, $A_{X}$, defined as the fraction of diphoton resonances satisfying the fiducial acceptance at the particle level. Both are computed for NWA spin-0 models as a function of $m_{X}$.

Effect of event selections on a scalar MC signal sample generated for $m_{X}$ = 15 GeV and on the data. For the MC sample, the efficiencies are shown after applying event weights and a truth level filter that requires two photons with $p^{\gamma\gamma}_{T}>40$ GeV; for the data, the absolute yields are shown. The initial yields for data include a trigger preselection that is the OR of a list of single photon and diphoton triggers. The "2 $loose$ photons" step includes the kinematic acceptance cuts.

Parameterization of the Double Sided Crystal Ball function parameters describing the scalar mass resolution model as a function of $m_{X}$ [GeV].

Background-only fit to the observed numbers of events (black points) as a function of the diphoton invariant mass $m_{\gamma\gamma}$.

Expected and observed 95% CL upper limits on the signal fiducial cross section $\sigma_{\textrm{fid}}$, assuming 100% branching ratio for ALP decay into two photons, as a function of the hypothetical ALP mass $m_{\textrm{X}}$. The 1$\sigma$ and 2$\sigma$ confidence intervals are shown by the coloured bands.

Expected and observed 95% CL upper limits on the the ALP coupling constant, assuming 100% branching ratio for ALP decay into two photons, as a function of the hypothetical ALP mass $m_{\textrm{X}}$. The 1$\sigma$ and 2$\sigma$ confidence intervals are shown by the coloured bands. Contours of the ALP natural width $\Gamma$ are illustrated by the smooth blue solid lines.

Fiducial region acceptance of signal events for coupling constant f$^{-1}=0.05$ TeV$^{-1}$. EL, SD, and DD stand for the exclusive, single-dissociative, and double-dissociative processes. The acceptances are parameterised as $p_0 + p^{i}_{1} e^{p_{2}^{i} m_{\textrm{X}}} + p_{3}^{i} e^{p_{4}^{i} m_{\textrm{X}}} $ , where $i$ specifies the process type, exclusive, single-dissociative, or double-dissociative.

Correction factors $A_{\textrm{X}}/A_{0}$ and $C_{\textrm{X}}$ for the exclusive signal events with ALP coupling constant f$^{-1}$=0.05 TeV$^{-1}$, where acceptance $A_{\textrm{X}}$ and efficiency $C_{\textrm{X}}$ are defined in arXiv:2102.13405. $A_{\textrm{X}}$/$A_{0}$ and $C_{\textrm{X}}$ are parameterised as $p_0 + p^{i}_{1} e^{p_{2}^{i} m_{\textrm{X}}} + p_{3}^{i} e^{p_{4}^{i} m_{\textrm{X}}}$, where $i$ specifies the process type.

Simulated ALP signal cutflow for $m_{\rm{X}}$=300 GeV and $m_{\rm{X}}$=1200 GeV. The yields are normalized to a luminosity of 14.6 fb$^{-1}$. Selections are applied for each side after the selection on acoplanarity. The union of the sets of selected events is taken finally.

Post-fit distributon of BDT discriminants for $\rho_{tc}=1.0$ with $m_A$ = 350 GeV interfered with H.($m_A - m_H$ = 50 GeV)

Observed and expected 95% CL upper limit on the signal strength as a function of $m_A$ for the g2HDM signal model using different coupling assumptions: $\rho_{tu}$ = 0.1, 0.4, 1.0.

Observed and expected 95% CL upper limit on the signal strength as a function of $m_A$ for the g2HDM signal model using different coupling assumptions: $\rho_{tc}$ = 0.1, 0.4, 1.0.

Observed and expected 95% CL upper limit on the signal strength as a function of $m_A$ for the g2HDM signal model using different coupling assumptions: $\rho_{tu}$ = 0.1, 0.4, 1.0 with A–H interference assuming $m_A$ − $m_H$ = 50 GeV.

Observed and expected 95% CL upper limit on the signal strength as a function of $m_A$ for the g2HDM signal model using different coupling assumptions: $\rho_{tc}$ = 0.1, 0.4, 1.0 with A–H interference assuming $m_A$ − $m_H$ = 50 GeV.

Observed 95% CL upper limit on the signal strength as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model.

Expected exclusion as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model

Observed exclusion as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model

Observed 95% CL upper limit on the signal strength as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model.

Expected exclusion as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model

Observed exclusion as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model

Observed 95% CL upper limit on the signal strength as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV.

Expected exclusion as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV

Observed exclusion as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV

Observed 95% CL upper limit on the signal strength as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV.

Expected exclusion as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV

Observed exclusion as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=300 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The m(tt) distribution in the Lhi3_Zin region of the lltt channel. <br><br><a href="?table=overview">return to overview</a>

The m(bb) distribution in the 2 b-tag 0L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The m(bb) distribution in the 3+ b-tag 0L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(gg->A)*B(A->ZH)*B(H->tt) in the (mA,mH) plane. The limits are shown for tanbeta=1. The tanbeta value is relevant only for the choice of the A boson width. <br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->tt)$ in the (mA,mH) plane. The limits are shown for tanbeta=10. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(gg->A)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=1. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=10. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(gg->A)*B(A->ZH)*B(H->tt) in he (mA,mH) plane. The limits are shown for tanbeta=0.5. The tanbeta value is relevant only for the choice of the A boson width. <br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(gg->A)*B(A->ZH)*B(H->tt) in the (mA,mH) plane. The limits are shown for tanbeta=5. The tanbeta value is relevant only for the choice of the A boson width. <br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->tt) in the (mA,mH) plane. The limits are shown for tanbeta=1. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->tt)$ in the (mA,mH) plane. The limits are shown for tanbeta=5. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(gg->A)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=0.5. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(gg->A)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=5. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=1. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=5. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=20. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=450 GeV hypothesis with the bbA signal shown for comparison. <br><br><a href="?table=overview">return to overview</a>

The m(tt) distribution in the L3hi_Zin region of the lltt channel for the mH=450 GeV hypothesis with the bbA signal shown for comparison. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=350 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=400 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=500 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=550 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=600 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=700 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=800 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=130 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=150 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=200 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=250 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=350 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=400 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=450 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=500 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=600 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=700 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=800 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH i+ the 3p b-tag signal region of the vvbb channel for the mH=130 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=150 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=200 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=250 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=350 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=400 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=450 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=500 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=600 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=700 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=800 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag 2L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag em region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag 2L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag em region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

Local significance for lltt ggF signals for tanbeta=1. <br><br><a href="?table=overview">return to overview</a>

Local significance for lltt bbA signals for tanbeta=20. <br><br><a href="?table=overview">return to overview</a>

Local significance for vvbb ggF signals for tanbeta=1. <br><br><a href="?table=overview">return to overview</a>

Local significance for vvbb bbA signals for tanbeta=20. <br><br><a href="?table=overview">return to overview</a>

Signal acceptance times efficiency for lltt ggF signals. <br><br><a href="?table=overview">return to overview</a>

Signal acceptance times efficiency for lltt bbA signals. <br><br><a href="?table=overview">return to overview</a>

Signal acceptance times efficiency for vvbb ggF signals. <br><br><a href="?table=overview">return to overview</a>

Signal acceptance times efficiency for vvbb bbA signals. <br><br><a href="?table=overview">return to overview</a>

The pT distribution of the third lepton in the L3hi_Zin region of the lltt channel. <br><br><a href="?table=overview">return to overview</a>

The eta distribution of the H boson candidate in the rest frame of the A boson candidate in the L3hi_Zin region of the lltt channel. <br><br><a href="?table=overview">return to overview</a>

The ETmiss distribution in the 2 b-tag 0L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The m_top_near distribution in the 2 b-tag 0L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The ETmiss distribution in the 3+ b-tag 0L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The m_top_near distribution in the 3+ b-tag 0L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

Post-fit distributions of $M_{\ell v jj}$ in the $R2B_A$ subregion for electrons. All process yields and nuisance parameters are set to the values obtained from the background plus signal fit. The signal considered for the fit corresponds to the purely right-handed production of a W' with $m_{W'}$ of 3.6 TeV and a relative width of 1$\%$ of the $m_{W'}$, and is represented by the solid red line. The lower panels show the data minus the expected number of events, normalized to the statistical uncertainty of the data. The orange band represents the systematic uncertainties, also normalized to the statistical uncertainty of the data.

Post-fit distributions of $M_{\ell v jj}$ in the $RW'_A$ subregion for muons. All process yields and nuisance parameters are set to the values obtained from the background plus signal fit. The signal considered for the fit corresponds to the purely right-handed production of a W' with $m_{W'}$ of 3.6 TeV and a relative width of 1$\%$ of the $m_{W'}$, and is represented by the solid red line. The lower panels show the data minus the expected number of events, normalized to the statistical uncertainty of the data. The orange band represents the systematic uncertainties, also normalized to the statistical uncertainty of the data.

Post-fit distributions of $M_{\ell v jj}$ in the $RW'_A$ subregion for electrons. All process yields and nuisance parameters are set to the values obtained from the background plus signal fit. The signal considered for the fit corresponds to the purely right-handed production of a W' with $m_{W'}$ of 3.6 TeV and a relative width of 1$\%$ of the $m_{W'}$, and is represented by the solid red line. The lower panels show the data minus the expected number of events, normalized to the statistical uncertainty of the data. The orange band represents the systematic uncertainties, also normalized to the statistical uncertainty of the data.

Post-fit distributions of $M_{\ell v jj}$ in the $RT_A$ subregion for muons. All process yields and nuisance parameters are set to the values obtained from the background plus signal fit. The signal considered for the fit corresponds to the purely right-handed production of a W' with $m_{W'}$ of 3.6 TeV and a relative width of 1$\%$ of the $m_{W'}$, and is represented by the solid red line. The lower panels show the data minus the expected number of events, normalized to the statistical uncertainty of the data. The orange band represents the systematic uncertainties, also normalized to the statistical uncertainty of the data.

Post-fit distributions of $M_{\ell v jj}$ in the $RT_A$ subregion for electrons. All process yields and nuisance parameters are set to the values obtained from the background plus signal fit. The signal considered for the fit corresponds to the purely right-handed production of a W' with $m_{W'}$ of 3.6 TeV and a relative width of 1$\%$ of the $m_{W'}$, and is represented by the solid red line. The lower panels show the data minus the expected number of events, normalized to the statistical uncertainty of the data. The orange band represents the systematic uncertainties, also normalized to the statistical uncertainty of the data.

Observed and expected 95% CL upper limits on the product of the production cross section for a left-handed production of a tb quark pair in the $s$-channel, mediated by a W or W' bosons and including interference terms, given as functions of $m_{W'}$ for a relative width of 1%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a left-handed production of a tb quark pair in the $s$-channel, mediated by a W or W' bosons and including interference terms, given as functions of $m_{W'}$ for a relative width of 10%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a left-handed production of a tb quark pair in the $s$-channel, mediated by a W or W' bosons and including interference terms, given as functions of $m_{W'}$ for a relative width of 20%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a left-handed production of a tb quark pair in the $s$-channel, mediated by a W or W' bosons and including interference terms, given as functions of $m_{W'}$ for a relative width of 30%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of $m_{W'}$ for a relative width of 1%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of $m_{W'}$ for a relative width of 10%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of $m_{W'}$ for a relative width of 20%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of $m_{W'}$ for a relative width of 30%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed 95% CL upper limit on the production cross section for a left-handed W' boson in the tb final state, as functions of $m_{W'}$ and relative width $\Gamma/m_{W'}$. Numbers in red, written diagonally, represent values of the excluded cross sections that are lower than the theoretical ones for the analyzed model.

Observed 95% CL upper limit on the production cross section for a right-handed W' boson in the tb final state, as functions of $m_{W'}$ and relative width $\Gamma/m_{W'}$. Numbers in red, written diagonally, represent values of the excluded cross sections that are lower than the theoretical ones for the analyzed model.

Observed 95% CL upper limit on the production cross section for a generalized left-right coupling of the W' boson to a t and a b quark for a mass of the W' boson of 2 TeV.

Observed 95% CL upper limit on the production cross section for a generalized left-right coupling of the W' boson to a t and a b quark for a mass of the W' boson of 2.8 TeV.

Observed 95% CL upper limit on the production cross section for a generalized left-right coupling of the W' boson to a t and a b quark for a mass of the W' boson of 3.6 TeV.

Observed 95% CL upper limit on the production cross section for a generalized left-right coupling of the W' boson to a t and a b quark for a mass of the W' boson of 4.4 TeV.

Observed 95% CL upper limit on the production cross section for a generalized left-right coupling of the W' boson to a t and a b quark for a mass of the W' boson of 5.2 TeV.

Observed 95% CL upper limit on the production cross section for a generalized left-right coupling of the W' boson to a t and a b quark for a mass of the W' boson of 6 TeV.

Expected 95% CL lower limit on $m_{W'}$ for a generalized left-right coupling of the W' boson to a t and a b quark.

Observed 95% CL lower limit on $m_{W'}$ for a generalized left-right coupling of the W' boson to a t and a b quark.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of the $m_{W'}$ for a relative width of 10%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid curves show the theoretical expectation at LO in the case that the couplings, and thus the partial widths, are varied together with the total width. In this interpretation the branching fraction of the W' --> tb decays is the same for each value of the width of the W' boson.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of the $m_{W'}$ for a relative width of 20%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid curves show the theoretical expectation at LO in the case that the couplings, and thus the partial widths, are varied together with the total width. In this interpretation the branching fraction of the W' --> tb decays is the same for each value of the width of the W' boson.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of the $m_{W'}$ for a relative width of 30%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid curves show the theoretical expectation at LO in the case that the couplings, and thus the partial widths, are varied together with the total width. In this interpretation the branching fraction of the W' --> tb decays is the same for each value of the width of the W' boson.

Measured min-$d_{xy}$ distribution in the 2-match category, after requiring the $d_{xy}$ muon to pass the isolation requirement $I_{\mathrm{rel}}^{\mathrm{PF}} <0.25$ (i. e., the B and D bins of the ABCD plane). Overlaid with a red histogram is the background predicted from the region of the ABCD plane failing the same requirement (the A and C bins), as well as three signal benchmark hypotheses (as defined in the legends), assuming $\alpha_D = \alpha_{\mathrm{EM}}$ (the fine-structure constant). The red hatched bands correspond to the background prediction uncertainty. The last bin includes the overflow.

Two-dimensional exclusion surfaces for $\Delta = 0.1 \, m_1$ as a function of the DM mass $m_1$ and the signal strength $y$, with $m_{A'} = 3 \, m_1$. Filled histograms denote observed limits on $\sigma(\mathrm{pp} \to A' \to \chi_2 \chi_1) \, \mathcal{B}(\chi_2 \to \chi_1 \mu^+ \mu^-)$. Solid (dashed) curves denote the observed (expected) exclusion limits at 95% CL, with 68% CL uncertainty bands around the expectation. Regions above the curves are excluded, depending on the $\alpha_D$ hypothesis: $\alpha_{\mathrm{D}} = \alpha_{\mathrm{EM}}$ (dark blue) or 0.1 (light magenta).

Two-dimensional exclusion surfaces for $\Delta = 0.4 \, m_1$ as a function of the DM mass $m_1$ and the signal strength $y$, with $m_{A'} = 3 \, m_1$. Filled histograms denote observed limits on $\sigma(\mathrm{pp} \to A' \to \chi_2 \chi_1) \, \mathcal{B}(\chi_2 \to \chi_1 \mu^+ \mu^-)$. Solid (dashed) curves denote the observed (expected) exclusion limits at 95% CL, with 68% CL uncertainty bands around the expectation. Regions above the curves are excluded, depending on the $\alpha_D$ hypothesis: $\alpha_{\mathrm{D}} = \alpha_{\mathrm{EM}}$ (dark blue) or 0.1 (light magenta).

Signal region detector-level distribution for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$.

Signal region detector-level distribution for the observable $m_{e\mu}$.

Signal region detector-level distribution for the observable $p_{\mathrm{T}}^{e\mu}$.

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

Observed excluded cross-section at 95% CL for the WZ-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from Fig 8(aux).

Expected exclusion limits at 95% CL for the WZ-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Observed exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from from Fig 13(a) and from Fig 7 and Fig 10(aux).

Observed excluded cross-section at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from Fig 7(aux) and Fig 10(aux).

positive one $\sigma$ observed exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from from Fig 13(a) and from Fig 7 and Fig 10(aux).

negative one $\sigma$ observed exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from from Fig 13(a) and from Fig 7 and Fig 10(aux).

Expected exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Expected exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Expected exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{bRPV}_{2l-SS}$. in a susy scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV, tan$\beta$=5. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{bRPV}_{3l}$. in a susy scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV, tan$\beta$=5. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{WZ}_{high-m_{T2}}$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 150 GeV, $m(\tilde{\chi}^{0} _{1})$ = 50 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{WZ}_{low-m_{T2}}$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 150 GeV, $m(\tilde{\chi}^{0} _{1})$ = 50 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the low mass $SR^{RPV}_{2l1b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the medium mass $SR^{RPV}_{2l1b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the low mass $SR^{RPV}_{2l2b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the medium mass $SR^{RPV}_{2l2b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the high mass $SR^{RPV}_{2l2b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the low mass $SR^{RPV}_{2l3b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the medium mass $SR^{RPV}_{2l3b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the high mass $SR^{RPV}_{2l3b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the $SR^{Wh}_{low-m_{T2} }$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Higgs bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 300 GeV, $m(\tilde{\chi}^{0} _{1})$ = 100 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the $SR^{Wh}_{high-m_{T2} }$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Higgs bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 300 GeV, $m(\tilde{\chi}^{0} _{1})$ = 100 GeV. Only statistical uncertainties are shown.

Signal Hepdataeptance for $SR^{bRPV}_{2l-SS}$ signal region from Fig 13(a)(aux) in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal Hepdataeptance for $SR^{bRPV}_{3l}$ signal region from Fig 13(b)(aux) in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal acceptance for $SR^{WZ}_{high-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal acceptance for $SR^{WZ}_{low-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal acceptance for $SR^{RPV}_{2l1b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l1b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l2b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l2b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l2b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l3b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l3b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l3b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{bRPV}_{2l-SS}$ signal region in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal efficiency for $SR^{bRPV}_{3l}$ signal region in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal efficiency for $SR^{WZ}_{high-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal efficiency for $SR^{WZ}_{low-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal efficiency for $SR^{RPV}_{2l1b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l1b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l2b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l2b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l2b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l3b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l3b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l3b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{Wh}_{high-m_{T2} }$ signal region from Fig 11(a)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Signal acceptance for $SR^{Wh}_{low-m_{T2} }$ signal region from Fig 11(b)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Signal efficiency for $SR^{Wh}_{high-m_{T2} }$ signal region from Fig 15(a)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Signal efficiency for $SR^{Wh}_{low-m_{T2} }$ signal region from Fig 15(b)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Observed 95% X-section upper limits as a function of higgsino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{1}/\tilde{\chi}^{0}_{2}$ mass in the bilinear RPV model from Fig 14.

Observed 95% X-section upper limits as a function of higgsino $\tilde{\chi}^{0}_{1}/\tilde{\chi}^{0}_{2}$ mass in the UDD RPV model from Fig 18.

Observed 95% X-section upper limits as a function of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ mass in the WZ-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from Fig 9(aux).

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{WZ}_{high-m_{T2}}$ from publication's Figure 11(a) . The last bin is inclusive.

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{WZ}_{low-m_{T2}}$ from publication's Figure 11(b) . The last bin is inclusive.

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{bRPV}_{2l-SS}$ from publication's Figure 11(c) . The last bin is inclusive.

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{bRPV}_{3l}$ from publication's Figure 11(d) . The last bin is inclusive.

N-1 distributions for $\sum p^{b-jet}_{T}/\sum p^{jet}_{T}$ of observed data and expected background towards $SR^{RPV}_{2l1b}-L$ from publication's Figure 16(a) . The last bin is inclusive.

N-1 distributions for $\sum p^{b-jet}_{T}/\sum p^{jet}_{T}$ of observed data and expected background towards $SR^{RPV}_{2l2b}-M$ from publication's Figure 16(b) . The last bin is inclusive.

N-1 distributions for $\sum p^{b-jet}_{T}/\sum p^{jet}_{T}$ of observed data and expected background towards $SR^{RPV}_{2l3b}-H$ from publication's Figure 16(c) . The last bin is inclusive.

N-1 distribution for $E_{T}^{miss}$ in $SR^{Wh}_{high-m_{T2} }$ in ee channel

N-1 distribution for $E_{T}^{miss}$ in $SR^{Wh}_{high-m_{T2} }$ in e$\mu$ channel

N-1 distribution for $E_{T}^{miss}$ in $SR^{Wh}_{high-m_{T2} }$ in $\mu\mu$ channel

N-1 distribution for $\mathcal{S}(E_{T}^{miss})$ in $SR^{Wh}_{low-m_{T2} }$ in ee channel

N-1 distribution for $\mathcal{S}(E_{T}^{miss})$ in $SR^{Wh}_{low-m_{T2} }$ in e$\mu$ channel

N-1 distribution for $\mathcal{S}(E_{T}^{miss})$ in $SR^{Wh}_{low-m_{T2} }$ in $\mu\mu$ channel