Figure 3 top left. Expected and observed 95% CL upper limits on the product of the production cross section and branching fraction as a function of the RS graviton mass $m_{G}$ for the full Run 2 data set are shown. Expected $1\sigma$ and $2\sigma$ limit bands are shown in green and yellow, respectively

Figure 3 top right. Expected and observed 95% CL upper limits on the product of the production cross section and branching fraction as a function of the heavy higgs mass $m_{S}$ for the full Run 2 data set are shown. Expected $1\sigma$ and $2\sigma$ limit bands are shown in green and yellow, respectively

Figure 3 mid left. Expected and observed 95% CL upper limits on the product of the production cross section and branching fraction as a function of the RS graviton mass $m_{G}$ for the full Run 2 data set are shown. Expected $1\sigma$ and $2\sigma$ limit bands are shown in green and yellow, respectively

Figure 3 mid right. Expected and observed 95% CL upper limits on the product of the production cross section and branching fraction as a function of the heavy higgs mass $m_{S}$ for the full Run 2 data set are shown. Expected $1\sigma$ and $2\sigma$ limit bands are shown in green and yellow, respectively

Figure 3 bottom left. Expected and observed 95% CL upper limits on the product of the production cross section and branching fraction as a function of the RS graviton mass $m_{G}$ for the full Run 2 data set are shown. Expected $1\sigma$ and $2\sigma$ limit bands are shown in green and yellow, respectively

Figure 3 bottom right. Expected and observed 95% CL upper limits on the product of the production cross section and branching fraction as a function of the heavy higgs mass $m_{S}$ for the full Run 2 data set are shown. Expected $1\sigma$ and $2\sigma$ limit bands are shown in green and yellow, respectively

Figure 4 top left. Expected 95% CL upper limits on the product of the cross section and branching fraction as a function of the $m_{G}$ versus the resonance width for the full Run 2 data set.

Figure 4 top right. Observed 95% CL upper limits on the product of the cross section and branching fraction as a function of the $m_{G}$ versus the resonance width for the full Run 2 data set.

Figure 4 bottom left. Expected 95% CL upper limits on the product of the cross section and branching fraction as a function of the $m_{S}$ versus the resonance width for the full Run 2 data set.

Figure 4 bottom right. Observed 95% CL upper limits on the product of the cross section and branching fraction as a function of the $m_{S}$ versus the resonance width for the full Run 2 data set.

Figure 6. The $m_{\gamma\gamma}$ spectra and background prediction after nuisance paramter marginalization (post-fit) due to SM diphoton production ($\gamma\gamma$) and fake photon production (j$\gamma$, jj) for the EBEB, combining the 2016, 2017, and 2018 data sets. The prediction with a large extra dimensions signal (GRW convention with $M_{S}$ = 9 TeV) is also shown. The pull distributions, normalized to the statistical uncertainty only, are shown underneath the spectra. The last bin contains the overflow of events beyond $m_{\gamma\gamma}$ > 3.5 TeV.

Figure 6. The $m_{\gamma\gamma}$ spectra and background prediction after nuisance paramter marginalization (post-fit) due to SM diphoton production ($\gamma\gamma$) and fake photon production (j$\gamma$, jj) for the EBEE, combining the 2016, 2017, and 2018 data sets. The prediction with a large extra dimensions signal (GRW convention with $M_{S}$ = 9 TeV) is also shown. The pull distributions, normalized to the statistical uncertainty only, are shown underneath the spectra. The last bin contains the overflow of events beyond $m_{\gamma\gamma}$ > 3.5 TeV.

Figure 7. The exclusion limit for the clockwork framework over the k--M5 parameter space. The shaded region denotes where the theory becomes nonperturbative. The region below and to the left of the solid line constitutes the excluded region. Expected 1std and 2std limit bands are shown in green and yellow, respectively.

Electroweak corrections for the $R_{\Delta\phi}(p_\mathrm{T})$ distribution. The multiplicative combination is considered as the main electroweak correction factor and the additive combination as an uncertainty estimate.

Example description

Example description

Example description

Observed and expected 95% CL upper limits on the nonresonant HH production via vector boson fusion cross section obtained for both single-lepton and dilepton channels, and from their combination

Observed and expected 95% CL upper limits on the nonresonant HH production cross section for two different benchmark scenarios “JHEP04(2016)01” and “JHEP03(2020)91”

Observed and expected 95% CL upper limits on the nonresonant HH production cross section as a function of the Higgs boson self-coupling strength modifier κλ All Higgs boson couplings other than λ are assumed to have the values predicted by the SM

Observed and expected 95% CL upper limits on the nonresonant HH production via VBF cross section as a function of the effective coupling κ2V. The ggF contribution in this case is set to the SM expectation. All other Higgs boson couplings are assumed to have the values predicted by the SM.

Observed and expected 95% CL upper limits on the nonresonant HH production cross section as a function of the effective coupling c2. All other Higgs boson couplings are assumed to have the values predicted in the SM

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.33, 0.33, 0.33) and in the Majorana scenario.

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.0, 1.0, 0.0) and in the Dirac-like scenario.

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.0, 0.5, 0.5) and in the Dirac-like scenario.

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.5, 0.5, 0.0) and in the Dirac-like scenario.

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.33, 0.33, 0.33) and in the Dirac-like scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 1.0 GeV in the Majorana scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 1.5 GeV in the Majorana scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 2.0 GeV in the Majorana scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 1.0 GeV in the Dirac-like scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 1.5 GeV in the Dirac-like scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 2.0 GeV in the Dirac-like scenario.

Observed profiled likelihood on $f_{\Lambda1}$ using Approach 1. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{a3}$ using Approach 1. The signal strength modifiers are treated as free parameters.

Observed profiled likelihood on $f_{a3}$ using Approach 1. The signal strength modifiers are treated as free parameters.

Expected profiled likelihood on $f_{\Lambda1}^{Z\gamma}$ using Approach 1. The signal strength modifiers are treated as free parameters.

Observed profiled likelihood on $f_{\Lambda1}^{Z\gamma}$ using Approach 1. The signal strength modifiers are treated as free parameters.

Expected profiled likelihood on $f_{a2}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Observed profiled likelihood on $f_{a2}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{\Lambda1}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Observed profiled likelihood on $f_{\Lambda1}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{a3}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters.

Observed profiled likelihood on $f_{a3}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters.

Expected profiled likelihood on $f_{a2}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Observed profiled likelihood on $f_{a2}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{\Lambda1}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Observed profiled likelihood on $f_{\Lambda1}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{a3}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters.

Observed profiled likelihood on $f_{a3}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters.

Expected profiled likelihood on the $\delta c_{Z}$ coupling of the SMEFT Higgs basis.

Observed profiled likelihood on the $\delta c_{Z}$ coupling of the SMEFT Higgs basis.

Expected profiled likelihood on the $c_{Z\Box}$ coupling of the SMEFT Higgs basis.

Observed profiled likelihood on the $c_{Z\Box}$ coupling of the SMEFT Higgs basis.

Expected profiled likelihood on the $c_{ZZ}$ coupling of the SMEFT Higgs basis.

Observed profiled likelihood on the $c_{ZZ}$ coupling of the SMEFT Higgs basis.

Expected profiled likelihood on the $\tilde{c}_{ZZ}$ coupling of the SMEFT Higgs basis.

Observed profiled likelihood on the $\tilde{c}_{ZZ}$ coupling of the SMEFT Higgs basis.

Expected profiled likelihood on $f_{a3}^{ggH}$. The signal strength modifiers and the CP-odd HVV anomalous coupling cross section fraction are treated as free parameters.

Observed profiled likelihood on $f_{a3}^{ggH}$. The signal strength modifiers and the CP-odd HVV anomalous coupling cross section fraction are treated as free parameters.

The 95% CL limits on $|V_{\mu N}|^2$ as a function of the HNL mass for a Majorana HNL.

The 95% CL limits on $|V_{\tau N}|^2$ as a function of the HNL mass for a Majorana HNL.

The 95% CL limits on $|V_{\tau N}|^2$ as a function of the HNL mass for a Majorana HNL.

The 95% CL limits on $|V_{eN}|^2$ as a function of the HNL mass for a Dirac HNL.

The 95% CL limits on $|V_{eN}|^2$ as a function of the HNL mass for a Dirac HNL.

The 95% CL limits on $|V_{\mu N}|^2$ as a function of the HNL mass for a Dirac HNL.

The 95% CL limits on $|V_{\mu N}|^2$ as a function of the HNL mass for a Dirac HNL.

The 95% CL limits on $|V_{\tau N}|^2$ as a function of the HNL mass for a Dirac HNL.

The 95% CL limits on $|V_{\tau N}|^2$ as a function of the HNL mass for a Dirac HNL.

Comparison of the number of predicted background events and the observed events in search regions La. The predicted background yields are shown with the values of the normalizations and nuisance parameters obtained in the background-only fit applied.

Comparison of the number of predicted background events and the observed events in search regions La. The predicted background yields are shown with the values of the normalizations and nuisance parameters obtained in the background-only fit applied.

Signal yields predicted in search regions La as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and an electron neutrino.

Signal yields predicted in search regions La as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and an electron neutrino.

Signal yields predicted in search regions La as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a muon neutrino.

Signal yields predicted in search regions La as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a muon neutrino.

Signal yields predicted in search regions La as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a tau neutrino.

Signal yields predicted in search regions La as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a tau neutrino.

Comparison of the number of predicted background events and the observed events in search regions Lb. The predicted background yields are shown with the values of the normalizations and nuisance parameters obtained in the background-only fit applied.

Comparison of the number of predicted background events and the observed events in search regions Lb. The predicted background yields are shown with the values of the normalizations and nuisance parameters obtained in the background-only fit applied.

Signal yields predicted in search regions Lb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and an electron neutrino.

Signal yields predicted in search regions Lb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and an electron neutrino.

Signal yields predicted in search regions Lb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a muon neutrino.

Signal yields predicted in search regions Lb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a muon neutrino.

Signal yields predicted in search regions Lb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a tau neutrino.

Signal yields predicted in search regions Lb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a tau neutrino.

Comparison of the number of predicted background events and the observed events in search regions Ha. The predicted background yields are shown with the values of the normalizations and nuisance parameters obtained in the background-only fit applied.

Comparison of the number of predicted background events and the observed events in search regions Ha. The predicted background yields are shown with the values of the normalizations and nuisance parameters obtained in the background-only fit applied.

Signal yields predicted in search regions Ha as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and an electron neutrino.

Signal yields predicted in search regions Ha as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and an electron neutrino.

Signal yields predicted in search regions Ha as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a muon neutrino.

Signal yields predicted in search regions Ha as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a muon neutrino.

Signal yields predicted in search regions Ha as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a tau neutrino.

Signal yields predicted in search regions Ha as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a tau neutrino.

Comparison of the number of predicted background events and the observed events in search regions Hb. The predicted background yields are shown with the values of the normalizations and nuisance parameters obtained in the background-only fit applied. Search regions Hb2 and Hb3 are merged for final states including a hadronically decayed tau lepton. This is represented by the value of -1 in the Hb2 entry.

Comparison of the number of predicted background events and the observed events in search regions Hb. The predicted background yields are shown with the values of the normalizations and nuisance parameters obtained in the background-only fit applied. Search regions Hb2 and Hb3 are merged for final states including a hadronically decayed tau lepton. This is represented by the value of -1 in the Hb2 entry.

Signal yields predicted in search regions Hb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and an electron neutrino. Search regions Hb2 and Hb3 are merged for final states including a hadronically decayed tau lepton. This is represented by the value of -1 in the Hb2 entry.

Signal yields predicted in search regions Hb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and an electron neutrino. Search regions Hb2 and Hb3 are merged for final states including a hadronically decayed tau lepton. This is represented by the value of -1 in the Hb2 entry.

Signal yields predicted in search regions Hb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a muon neutrino. Search regions Hb2 and Hb3 are merged for final states including a hadronically decayed tau lepton. This is represented by the value of -1 in the Hb2 entry.

Signal yields predicted in search regions Hb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a muon neutrino. Search regions Hb2 and Hb3 are merged for final states including a hadronically decayed tau lepton. This is represented by the value of -1 in the Hb2 entry.

Signal yields predicted in search regions Hb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a tau neutrino. Search regions Hb2 and Hb3 are merged for final states including a hadronically decayed tau lepton. This is represented by the value of -1 in the Hb2 entry.

Signal yields predicted in search regions Hb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a tau neutrino. Search regions Hb2 and Hb3 are merged for final states including a hadronically decayed tau lepton. This is represented by the value of -1 in the Hb2 entry.

Covariance matrix for the postfit background prediction in signal region bins La.

Covariance matrix for the postfit background prediction in signal region bins La.

Covariance matrix for the postfit background prediction in signal region bins Lb.

Covariance matrix for the postfit background prediction in signal region bins Lb.

Covariance matrix for the postfit background prediction in signal region bins Ha.

Covariance matrix for the postfit background prediction in signal region bins Ha.

Covariance matrix for the postfit background prediction in signal region bins Hb.

Covariance matrix for the postfit background prediction in signal region bins Hb.

Simulated signal yields passing each stage of the low mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and an electron neutrino.

Simulated signal yields passing each stage of the low mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and an electron neutrino.

Simulated signal yields passing each stage of the low mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and a muon neutrino.

Simulated signal yields passing each stage of the low mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and a muon neutrino.

Simulated signal yields passing each stage of the low mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and a tau neutrino.

Simulated signal yields passing each stage of the low mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and a tau neutrino.

Simulated signal yields passing each stage of the high mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and an electron neutrino.

Simulated signal yields passing each stage of the high mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and an electron neutrino.

Simulated signal yields passing each stage of the high mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and a muon neutrino.

Simulated signal yields passing each stage of the high mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and a muon neutrino.

Simulated signal yields passing each stage of the high mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and a tau neutrino.

Simulated signal yields passing each stage of the high mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and a tau neutrino.

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.