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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Correlation matrix associated to the statistical covariance matrix of the data and MC for the primary Lund jet plane density for AK8 jets in a one-dimensional representation with bin indices. The mapping between the bin indices and the physical binning can be imported from the XML file attached to this HepData record using the TUnfoldBinningXML class of ROOT (qualitatively, it corresponds to slicing the Lund plane horizontally from low kT to high kT).

Two-dimensional exclusion limits for Dirac HNL pure electron coupling scenario

Two-dimensional exclusion limits for Majorana HNL pure muon coupling scenario

Two-dimensional exclusion limits for Dirac HNL pure muon coupling scenario

Two-dimensional exclusion limits for Majorana HNL pure tau lepton coupling scenario

Two-dimensional exclusion limits for Dirac HNL pure tau lepton coupling scenario

Two-dimensional exclusion limits for Majorana HNL mixed electron-muon coupling scenario

DTwo-dimensional exclusion limits for Dirac HNL mixed electron-muon coupling scenario

Two-dimensional exclusion limits for Majorana HNL mixed electron-tau lepton coupling scenario

Two-dimensional exclusion limits for Dirac HNL mixed electron-tau lepton coupling scenario

Two-dimensional exclusion limits for Majorana HNL mixed muon-tau lepton coupling scenario

Two-dimensional exclusion limits for Dirac HNL mixed muon-tau lepton coupling scenario

Two-dimensional exclusion limits for Majorana HNL democratic lepton coupling scenario

Two-dimensional exclusion limits for Dirac HNL democratic lepton coupling scenario

Lower limits on Majorana HNL mass for fixed proper decay length of $c\tau=0.1\,\text{mm}$

Lower limits on Dirac HNL mass for fixed proper decay length of $c\tau=0.1\,\text{mm}$

Lower limits on Majorana HNL mass for fixed proper decay length of $c\tau=1\,\text{mm}$

Lower limits on Dirac HNL mass for fixed proper decay length of $c\tau=1\,\text{mm}$

Lower limits on Majorana HNL proper decay length for fixed mass of 4.5 GeV

Lower limits on Dirac HNL proper decay length for fixed mass of 4.5 GeV

Lower limits on Majorana HNL proper decay length for fixed mass of 8 GeV

Lower limits on Dirac HNL proper decay length for fixed mass of 8 GeV

The ratios between the hydrodynamic probes and their Taylor expansions. The ratios are plotted as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The ratios between the hydrodynamic probes and their Taylor expansions. The ratios are plotted as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The magnitudes of the measured skewness $\gamma_1^\textrm{exp}$ and "corrected" skewness $\gamma_{1,\textrm{corr}}^\textrm{exp}$ as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The magnitudes of the measured kurtosis $\gamma_2^\textrm{exp}$ and "corrected" kurtosis $\gamma_{2,\textrm{corr}}^\textrm{exp}$ as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The magnitudes of the measured superskewness $\gamma_3^\textrm{exp}$ and "corrected" superskewness $\gamma_{3,\textrm{corr}}^\textrm{exp}$ as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

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.

Transverse dilepton-$p_{T}^{miss}$ mass pre-fit distribution for selected events in SR2 of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton mass pre-fit distribution for selected events in SR3 of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton $p_T$ pre-fit distribution for selected events in SR3 of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dijet $p_T$ pre-fit distribution for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Lepton-dijet $p_T$ pre-fit distribution for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Maximum $p_T$ of the visible particles divided by the missing $p_T$, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between jets, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and missing $p_T$, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and the dijet system, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta R$ distance between the leptons, for selected events in the non-prompt validation region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the di-lepton system and the missing $p_T$, for selected events in the non-prompt validation region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and the di-jet system, for selected events in the non-prompt validation region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton-dijet system and the missing $p_T$, for selected events in the non-prompt validation region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Leading lepton $p_T$ pre-fit distribution for selected events in the top control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Trailing lepton $p_T$ pre-fit distribution for selected events in the top control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton mass pre-fit distribution for selected events in the Drell-Yan control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Missing $p_T$ pre-fit distribution for selected events in the Drell-Yan control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Missing $p_T$ pre-fit distribution for selected events in the $WW$ control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton $p_T$ pre-fit distribution for selected events in the $WW$ control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Minimum $p_T$ of the visible particles divided by the missing $p_T$, for selected events in the top control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between jets, for selected events in the top control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and missing $p_T$, for selected events in the $W$+jets control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and the dijet system, for selected events in the $W$+jets control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Unrolled $m_{ll}-m_{T}^{l min, p_{T}^{miss}}$ post-fit distributions in the di-leptonic channel in SR1 for the full data set. The histogram bins are spaced uniformly. Each group of five bins (from left to right) corresponds to $m_{T}^{l min, p_{T}^{miss}}$ distribution in a $m_{ll}$ region, placed in ascending order. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$.

Unrolled $m_{ll}-m_{T}^{l min, p_{T}^{miss}}$ post-fit distributions in the di-leptonic channel in SR2 for the full data set. The histogram bins are spaced uniformly. Each group of five bins (from left to right) corresponds to $m_{T}^{l min, p_{T}^{miss}}$ distribution in a $m_{ll}$ region, placed in ascending order. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$.

Unrolled $m_{ll}-m_{T}^{l min, p_{T}^{miss}}$ post-fit distributions in the di-leptonic channel in SR3 for the full data set. The histogram bins are spaced uniformly. Each group of five bins (from left to right) corresponds to $m_{T}^{l min, p_{T}^{miss}}$ distribution in a $m_{ll}$ region, placed in ascending order. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$.

Post-fit BDT distributions in the semi-leptonic channel for the full data set in the top control region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Post-fit BDT distributions in the semi-leptonic channel for the full data set in the $W$+jets control region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Post-fit BDT distributions in the semi-leptonic channel in the 2016 signal region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Post-fit BDT distributions in the semi-leptonic channel in the 2017+2018 signal region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Observed upper limits at 90% CL on the square of the kinetic mixing coefficient $\epsilon$ in the minimal model of a dark photon

Observed upper limits at 90% CL on the mixing angle $\theta_H$ for the type IV 2HDM+S scenario with fixed $tan\beta=0.5$

The combined efficiency of the dimuon scouting trigger and the MVA muon selection, averaged between 2017 and 2018, weighted by the integrated luminosity of each year. The trigger and selection efficiencies are also included in the table.

Theory cross-section times branching fraction to muons times acceptance for the dark photon and 2HDM+S models.

Compatibility of normalized differential cross section measurements with modeling predictions. The compatibility is quantified with z scores for each of the theoretical predictions, given the unfolded normalized differential cross sections and their covariances. A lower value indicates a better agreement between prediction and measurement. A value of z = 2 indicates a p-value of 5%. In the calculation of the z score only the measurement uncertainties and the statistical uncertainties of the modeling predictions are taken into account

Normalized differential cross section of $|\eta(\mathrm{b}^{\mathrm{add.}}_{1})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}^{\mathrm{add.}}_{1})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{1})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{1})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}^{\mathrm{add.}}_{2})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}^{\mathrm{add.}}_{2})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{2})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{2})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{add.}})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{add.}})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\Delta\mathrm{R}_{\mathrm{b}\mathrm{b}}^{\mathrm{avg}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\Delta\mathrm{R}_{\mathrm{b}\mathrm{b}}^{\mathrm{avg}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}_{3})|$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}_{3})|$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}_{3})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}_{3})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}_{3})$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}_{3})$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}_{3})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}_{3})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}_{4})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}_{4})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}_{4})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}_{4})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{b}}_{\mathrm{T}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{b}}_{\mathrm{T}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{b}}_{\mathrm{T}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{b}}_{\mathrm{T}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $|\eta(\mathrm{b}^{\mathrm{extra}}_{1})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}^{\mathrm{extra}}_{1})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{1})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{1})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}^{\mathrm{extra}}_{2})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}^{\mathrm{extra}}_{2})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{2})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{2})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{extra}})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{extra}})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $H^{\mathrm{j}}_{\mathrm{T}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{j}}_{\mathrm{T}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{j}}_{\mathrm{T}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{j}}_{\mathrm{T}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\mathrm{m}_{\mathrm{b}\mathrm{b}}^{\mathrm{max}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\mathrm{m}_{\mathrm{b}\mathrm{b}}^{\mathrm{max}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{light}}_{\mathrm{T}}$ in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $H^{\mathrm{light}}_{\mathrm{T}}$ in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $H^{\mathrm{light}}_{\mathrm{T}}$ in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $H^{\mathrm{light}}_{\mathrm{T}}$ in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $N_{\mathrm{jets}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $N_{\mathrm{jets}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $N_{\mathrm{jets}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $N_{\mathrm{jets}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $N_{\mathrm{b}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $N_{\mathrm{b}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{add.}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{add.}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\Delta\mathrm{R}_{\mathrm{b}\mathrm{b}}^{\mathrm{avg}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\Delta\mathrm{R}_{\mathrm{b}\mathrm{b}}^{\mathrm{avg}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}_{4})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}_{4})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}_{4})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}_{4})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{extra}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{extra}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\mathrm{m}_{\mathrm{b}\mathrm{b}}^{\mathrm{max}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\mathrm{m}_{\mathrm{b}\mathrm{b}}^{\mathrm{max}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $N_{\mathrm{jets}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $N_{\mathrm{jets}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $N_{\mathrm{jets}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $N_{\mathrm{jets}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $N_{\mathrm{b}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $N_{\mathrm{b}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{add.}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{add.}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\Delta\mathrm{R}_{\mathrm{b}\mathrm{b}}^{\mathrm{avg}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\Delta\mathrm{R}_{\mathrm{b}\mathrm{b}}^{\mathrm{avg}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}_{4})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}_{4})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}_{4})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}_{4})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{extra}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{extra}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\mathrm{m}_{\mathrm{b}\mathrm{b}}^{\mathrm{max}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\mathrm{m}_{\mathrm{b}\mathrm{b}}^{\mathrm{max}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $N_{\mathrm{jets}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $N_{\mathrm{jets}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $N_{\mathrm{jets}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $N_{\mathrm{jets}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $N_{\mathrm{b}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $N_{\mathrm{b}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Fiducial cross sections from the measurements of all observables, compared to predictions from different ttbb simulation approaches. For each of the normalized differential measurements the fiducial cross section in the respective phase space is also determined. In the paper only one representative observable is quoted for each fiducial phase space, while here the measured cross section with the uncertainties from the fit to the respective observable is summarized.