Selection efficiencies for the barrel region at object level. Definitions used in this table: Track-Match efficiency (barrel / endcap) is the fraction of reconstructed superclusters matched to at least two tracks (within a ΔR < 0.15 cone) with $p_T$>25GeV, normalized to the number of generated electron pairs. Merged-electron-pair (MEP) reconstruction efficiency (barrel / endcap) is the fraction of events in which an MEP is successfully reconstructed, using the same normalization. Iso-cuts efficiency (barrel / endcap) is the fraction of reconstructed MEPs that pass the isolation requirement. BDT efficiency (barrel / endcap) is the fraction of MEPs that pass the BDT selection among those passing the isolation requirement. Merged-electron-pair ID efficiency (barrel / endcap) is the fraction of MEPs that satisfy all signal-selection criteria, again normalized to the number of generated electron pairs.

Selection efficiencies for the endcap region at object level. Definitions as above.

Trigger-level and overall event-level efficiency.Definitions used in this table: Trigger efficiency is defined as the number of events passing both the trigger and the HAA4e offline selection, divided by the number passing the offline selection alone. Number of normalized events passing all selections gives the luminosity-scaled event yield after applying all cuts. Total event efficiency is the fraction of generated events in which a four-electron invariant mass is successfully reconstructed within the 80-250 GeV range, with event weights applied.

Selection efficiencies for the barrel region at object level. Definitions used in this table: Track-Match efficiency (barrel / endcap) is the fraction of reconstructed superclusters matched to at least two tracks (within a ΔR < 0.15 cone) with $p_T$>25GeV, normalized to the number of generated electron pairs. Merged-electron-pair (MEP) reconstruction efficiency (barrel / endcap) is the fraction of events in which an MEP is successfully reconstructed, using the same normalization. Iso-cuts efficiency (barrel / endcap) is the fraction of reconstructed MEPs that pass the isolation requirement. BDT efficiency (barrel / endcap) is the fraction of MEPs that pass the BDT selection among those passing the isolation requirement. Merged-electron-pair ID efficiency (barrel / endcap) is the fraction of MEPs that satisfy all signal-selection criteria, again normalized to the number of generated electron pairs.

Selection efficiencies for the endcap region at object level. Definitions as above.

Trigger-level and overall event-level efficiency.Definitions used in this table: Trigger efficiency is defined as the number of events passing both the trigger and the HAA4e offline selection, divided by the number passing the offline selection alone. Number of normalized events passing all selections gives the luminosity-scaled event yield after applying all cuts. Total event efficiency is the fraction of generated events in which a four-electron invariant mass is successfully reconstructed within the 80-250 GeV range, with event weights applied.

Selection efficiencies for the barrel region at object level. Definitions used in this table: Track-Match efficiency (barrel / endcap) is the fraction of reconstructed superclusters matched to at least two tracks (within a ΔR < 0.15 cone) with $p_T$>25GeV, normalized to the number of generated electron pairs. Merged-electron-pair (MEP) reconstruction efficiency (barrel / endcap) is the fraction of events in which an MEP is successfully reconstructed, using the same normalization. Iso-cuts efficiency (barrel / endcap) is the fraction of reconstructed MEPs that pass the isolation requirement. BDT efficiency (barrel / endcap) is the fraction of MEPs that pass the BDT selection among those passing the isolation requirement. Merged-electron-pair ID efficiency (barrel / endcap) is the fraction of MEPs that satisfy all signal-selection criteria, again normalized to the number of generated electron pairs.

Selection efficiencies for the endcap region at object level. Definitions as above.

Trigger-level and overall event-level efficiency.Definitions used in this table: Trigger efficiency is defined as the number of events passing both the trigger and the HAA4e offline selection, divided by the number passing the offline selection alone. Number of normalized events passing all selections gives the luminosity-scaled event yield after applying all cuts. Total event efficiency is the fraction of generated events in which a four-electron invariant mass is successfully reconstructed within the 80-250 GeV range, with event weights applied.

Selection efficiencies for the barrel region at object level. Definitions used in this table: Track-Match efficiency (barrel / endcap) is the fraction of reconstructed superclusters matched to at least two tracks (within a ΔR < 0.15 cone) with $p_T$>25GeV, normalized to the number of generated electron pairs. Merged-electron-pair (MEP) reconstruction efficiency (barrel / endcap) is the fraction of events in which an MEP is successfully reconstructed, using the same normalization. Iso-cuts efficiency (barrel / endcap) is the fraction of reconstructed MEPs that pass the isolation requirement. BDT efficiency (barrel / endcap) is the fraction of MEPs that pass the BDT selection among those passing the isolation requirement. Merged-electron-pair ID efficiency (barrel / endcap) is the fraction of MEPs that satisfy all signal-selection criteria, again normalized to the number of generated electron pairs.

Selection efficiencies for the endcap region at object level. Definitions as above.

Trigger-level and overall event-level efficiency.Definitions used in this table: Trigger efficiency is defined as the number of events passing both the trigger and the HAA4e offline selection, divided by the number passing the offline selection alone. Number of normalized events passing all selections gives the luminosity-scaled event yield after applying all cuts. Total event efficiency is the fraction of generated events in which a four-electron invariant mass is successfully reconstructed within the 80-250 GeV range, with event weights applied.

$M_{\tau}$ distribution in signal region, (OS$\mu$, Belle II)

$M_{\tau}$ distribution in signal region, (SS$e$, Belle)

$M_{\tau}$ distribution in signal region, (SS$e$, Belle II)

$M_{\tau}$ distribution in signal region, (SS$\mu$, Belle)

$M_{\tau}$ distribution in signal region, (SS$\mu$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$ and $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (OS$e$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$ and $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (OS$e$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$ and $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (OS$\mu$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$ and $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (OS$\mu$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$ and $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (SS$e$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$ and $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (SS$e$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$ and $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (SS$\mu$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$ and $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (SS$\mu$, Belle II)

Efficiency as function the helicity angle of the kaon $\cos\theta_{K}$, and the helicity angle of the lepton $\cos\theta_{\ell}$ in bin e the four-momentum-transfer $q^{2}$, (OS$e$, Belle)

Efficiency as function the helicity angle of the kaon $\cos\theta_{K}$, and the helicity angle of the lepton $\cos\theta_{\ell}$ in bin e the four-momentum-transfer $q^{2}$, (OS$e$, Belle II)

Efficiency as function the helicity angle of the kaon $\cos\theta_{K}$, and the helicity angle of the lepton $\cos\theta_{\ell}$ in bin e the four-momentum-transfer $q^{2}$, (OS$\mu$, Belle)

Efficiency as function the helicity angle of the kaon $\cos\theta_{K}$, and the helicity angle of the lepton $\cos\theta_{\ell}$ in bin e the four-momentum-transfer $q^{2}$, (OS$\mu$, Belle II)

Efficiency as function the helicity angle of the kaon $\cos\theta_{K}$, and the helicity angle of the lepton $\cos\theta_{\ell}$ in bin e the four-momentum-transfer $q^{2}$, (SS$e$, Belle)

Efficiency as function the helicity angle of the kaon $\cos\theta_{K}$, and the helicity angle of the lepton $\cos\theta_{\ell}$ in bin e the four-momentum-transfer $q^{2}$, (SS$e$, Belle II)

Efficiency as function the helicity angle of the kaon $\cos\theta_{K}$, and the helicity angle of the lepton $\cos\theta_{\ell}$ in bin e the four-momentum-transfer $q^{2}$, (SS$\mu$, Belle)

Efficiency as function the helicity angle of the kaon $\cos\theta_{K}$, and the helicity angle of the lepton $\cos\theta_{\ell}$ in bin e the four-momentum-transfer $q^{2}$, (SS$\mu$, Belle II)

Efficiency as function of the helicity angle of the kaon $\cos\theta_{K}$, (OS$e$, Belle)

Efficiency as function of the helicity angle of the kaon $\cos\theta_{K}$, (OS$e$, Belle II)

Efficiency as function of the helicity angle of the kaon $\cos\theta_{K}$, (OS$\mu$, Belle)

Efficiency as function of the helicity angle of the kaon $\cos\theta_{K}$, (OS$\mu$, Belle II)

Efficiency as function of the helicity angle of the kaon $\cos\theta_{K}$, (SS$e$, Belle)

Efficiency as function of the helicity angle of the kaon $\cos\theta_{K}$, (SS$e$, Belle II)

Efficiency as function of the helicity angle of the kaon $\cos\theta_{K}$, (SS$\mu$, Belle)

Efficiency as function of the helicity angle of the kaon $\cos\theta_{K}$, (SS$\mu$, Belle II)

Efficiency as function of the helicity angle of the lepton $\cos\theta_{\ell}$, (OS$e$, Belle)

Efficiency as function of the helicity angle of the lepton $\cos\theta_{\ell}$, (OS$e$, Belle II)

Efficiency as function of the helicity angle of the lepton $\cos\theta_{\ell}$, (OS$\mu$, Belle)

Efficiency as function of the helicity angle of the lepton $\cos\theta_{\ell}$, (OS$\mu$, Belle II)

Efficiency as function of the helicity angle of the lepton $\cos\theta_{\ell}$, (SS$e$, Belle)

Efficiency as function of the helicity angle of the lepton $\cos\theta_{\ell}$, (SS$e$, Belle II)

Efficiency as function of the helicity angle of the lepton $\cos\theta_{\ell}$, (SS$\mu$, Belle)

Efficiency as function of the helicity angle of the lepton $\cos\theta_{\ell}$, (SS$\mu$, Belle II)

Efficiency as function of the acoplanarity angle $\phi$, (OS$e$, Belle)

Efficiency as function of the acoplanarity angle $\phi$, (OS$e$, Belle II)

Efficiency as function of the acoplanarity angle $\phi$, (OS$\mu$, Belle)

Efficiency as function of the acoplanarity angle $\phi$, (OS$\mu$, Belle II)

Efficiency as function of the acoplanarity angle $\phi$, (SS$e$, Belle)

Efficiency as function of the acoplanarity angle $\phi$, (SS$e$, Belle II)

Efficiency as function of the acoplanarity angle $\phi$, (SS$\mu$, Belle)

Efficiency as function of the acoplanarity angle $\phi$, (SS$\mu$, Belle II)

Efficiency as function of the four-momentum-transfer $q^{2}$, (OS$e$, Belle)

Efficiency as function of the four-momentum-transfer $q^{2}$, (OS$e$, Belle II)

Efficiency as function of the four-momentum-transfer $q^{2}$, (OS$\mu$, Belle)

Efficiency as function of the four-momentum-transfer $q^{2}$, (OS$\mu$, Belle II)

Efficiency as function of the four-momentum-transfer $q^{2}$, (SS$e$, Belle)

Efficiency as function of the four-momentum-transfer $q^{2}$, (SS$e$, Belle II)

Efficiency as function of the four-momentum-transfer $q^{2}$, (SS$\mu$, Belle)

Efficiency as function of the four-momentum-transfer $q^{2}$, (SS$\mu$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$, (OS$e$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$, (OS$e$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$, (OS$\mu$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$, (OS$\mu$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$, (SS$e$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$, (SS$e$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$, (SS$\mu$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$, (SS$\mu$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (OS$e$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (OS$e$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (OS$\mu$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (OS$\mu$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (SS$e$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (SS$e$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (SS$\mu$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (SS$\mu$, Belle II)

Distributions of the number of hits in the cluster (Nhits) for the DT category in the signal region (SR). The last histogram bin contains all overflow events.

Distributions of the number of hits in the cluster (Nhits) for the DT category in the out-of-time (OOT) region. The last histogram bin contains all overflow events.

Distributions of the number of hits in the cluster (Nhits) for the CSC category in the signal region (SR). The last histogram bin contains all overflow events.

Distributions of the number of hits in the cluster (Nhits) for the CSC category in the out-of-time (OOT) region. The last histogram bin contains all overflow events.

Distributions of the number of hits in the cluster (Nhits) for the DT category in the out-of-time (OOT) region. The last histogram bin contains all overflow events.

The 95% CL observed and expected upper limits on the VLL production cross section as a function of the VLL mass for the pseudoscalar mean proper decay length c$\tau$ = 0.025 m. The pseudoscalar mass is 2 GeV.

Distributions of the number of hits in the cluster (Nhits) for the CSC category in the out-of-time (OOT) region. The last histogram bin contains all overflow events.

The 95% CL observed and expected upper limits on the VLL production cross section as a function of the pseudoscalar mean proper decay length c$\tau$ for VLL Mass of 700 GeV. The pseudoscalar mass is 2 GeV.

The 95% CL observed and expected upper limits on the VLL production cross section as a function of the VLL mass for the pseudoscalar mean proper decay length c$\tau$ = 0.025 m. The pseudoscalar mass is 2 GeV.

The 95% CL observed upper limits on the VLL production cross section as a function of the VLL mass and the pseudoscalar mean proper decay length c$\tau$. The pseudoscalar mass is 2 GeV.

The 95% CL observed and expected upper limits on the VLL production cross section as a function of the pseudoscalar mean proper decay length c$\tau$ for VLL Mass of 700 GeV. The pseudoscalar mass is 2 GeV.

The 95% CL observed exclusion contour on the VLL production cross section as a function of the VLL mass and the pseudoscalar mean proper decay length c$\tau$. The pseudoscalar mass is 2 GeV.

The 95% CL observed upper limits on the VLL production cross section as a function of the VLL mass and the pseudoscalar mean proper decay length c$\tau$. The pseudoscalar mass is 2 GeV.

Selection efficiencies in the CSC category signal region (SR) for different signal hypothesis. The pseudoscalar mass is 2 GeV.

The 95% CL observed exclusion contour on the VLL production cross section as a function of the VLL mass and the pseudoscalar mean proper decay length c$\tau$. The pseudoscalar mass is 2 GeV.

Selection efficiencies in the DT category signal region (SR) for different signal hypothesis. The pseudoscalar mass is 2 GeV.

Selection efficiencies in the CSC category signal region (SR) for different signal hypothesis. The pseudoscalar mass is 2 GeV.

Selection efficiencies in the DT category signal region (SR) for different signal hypothesis. The pseudoscalar mass is 2 GeV.

The cluster reconstruction efficiency, including both DT and CSC clusters, as a function of the simulated r and |z| decay positions of the pseudoscalar into photons in events with MET > 200 GeV, for a VLL mass of 300 GeV and a pseudoscalar mass of 2 GeV, and a range of ctau values uniformly distributed between 0.01 and 0.1 m.

The cluster reconstruction efficiency, including both DT and CSC clusters, as a function of the simulated r and |z| decay positions of the pseudoscalar into photons in events with MET > 200 GeV, for a VLL mass of 500 GeV and a pseudoscalar mass of 2 GeV, and a range of ctau values uniformly distributed between 0.01 and 0.1 m.

Upsilon 1S+2S+3S integrated cross section.

Upsilon 1S+2S+3S invariant cross section.

Upsilon 2S and 3S invariant cross section.

Upsilon 1S+2S+3S invariant cross section vs. rapidity.

Upsilon 1S invariant cross section.

Upsilon 2S invariant cross section.

Upsilon 3S invariant cross section.

Upsilon 1S invariant cross section vs. rapidity.

Upsilon 2S invariant cross section vs. rapidity.

Upsilon 3S invariant cross section vs. rapidity.

Upsilon x_{T} dependence

Upsilon x_{T} dependence.

Upsilon cross section ratios vs. collision energy.

Upsilon cross section ratios vs. charged particle multiplicity.

Upsilon cross section ratios vs. charged particle multiplicity.

Upsilon cross section ratios vs. charged particle multiplicity.

Upsilon cross section ratios vs. charged particle multiplicity.

Upsilon cross section ratios vs. charged particle multiplicity.

Self-normalized Upsilon ratio vs. charged particle multiplicity

Self-normalized Upsilon ratio vs. charged particle multiplicity

Measured single differential cross-sections in bins of $y^{Z}$. The first uncertainty is statistical, the second systematic, and the third is from the uncertainty on the integrated luminosity.

Measured single differential cross-sections in bins of $p_{T}^{Z}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured single differential cross-sections in bins of $\phi_{\eta}^{*}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Statistical correlation matrix for the one-dimensional $y^{Z}$ measurement.

Statistical correlation matrix for the one-dimensional $p_{T}^{Z}$ measurement.

Statistical correlation matrix for the one-dimensional $\phi_{\eta}^{*}$ measurement.

Correlation matrix for the efficiency uncertainty of the one-dimensional $y^{Z}$ measurement.

Correlation matrix for the efficiency uncertainty of the one-dimensional $p_{T}^{Z}$ measurement.

Correlation matrix for the efficiency uncertainty of the one-dimensional $\phi_{\eta}^{*}$ measurement.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Total efficiencies for the selection of signal events after categorization and application of the tighter photon $E_{\mathrm{T}}^{\gamma}$ selection used to optimize the signal significance spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$. In addition to the total efficiency, contributions to the signal selection from each of the separate event categories are shown. The efficiencies calculated from MC samples with $W/Z$ hadronic decays are shown as the points on each curve. The line presents interpolated results.

Total efficiencies for the selection of signal events after categorization and application of the tighter photon $E_{\mathrm{T}}^{\gamma}$ selection used to optimize the signal significance spin-2 $q\bar{q}\to X^0 \to Z\gamma$. In addition to the total efficiency, contributions to the signal selection from each of the separate event categories are shown. The efficiencies calculated from MC samples with $W/Z$ hadronic decays are shown as the points on each curve. The line presents interpolated results.

Total efficiencies for the selection of signal events after categorization and application of the tighter photon $E_{\mathrm{T}}^{\gamma}$ selection used to optimize the signal significance spin-2 $gg\to X^0 \to Z\gamma$. In addition to the total efficiency, contributions to the signal selection from each of the separate event categories are shown. The efficiencies calculated from MC samples with $W/Z$ hadronic decays are shown as the points on each curve. The line presents interpolated results.

The $m_{J\gamma}$ distributions of data events selected for the spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$ search in the WMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-0 $gg \to X^{0} \to Z\gamma$ search in the BTAG category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-0 $gg \to X^{0} \to Z\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-0 $gg \to X^{0} \to Z\gamma$ search in the ZMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $gg \to X^{0} \to Z\gamma$ search in the BTAG category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $gg \to X^{0} \to Z\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $gg \to X^{0} \to Z\gamma$ search in the ZMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $q\bar{q} \to X^{0} \to Z\gamma$ search in the BTAG category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $q\bar{q} \to X^{0} \to Z\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $q\bar{q} \to X^{0} \to Z\gamma$ search in the ZMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-0 $gg \to X^{0} \to Z\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-2 $gg \to X^{0} \to Z\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-2 $q\bar{q} \to X^{0} \to Z\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The jet mass distribution of large-$R$ jets originating from the hadronic decay of $W$ and $Z$ bosons produced from the decay of BSM bosons with mass $m_X = 2000$ GeV. The decays simulated are for the production models $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$ with a spin-1 resonance $X^{\pm}$ and $gg\to X^0 \to Z\gamma$ with a spin-0 resonance $X^{0}$.

The photon $E_{T}$ selection criteria as a function of resonance mass for D2 categories. Selections for four different signal channels are shown together.

The photon $E_{T}$ selection criteria as a function of resonance mass for ZMASS or WMASS categories. Selections for four different signal channels are shown together.

Level-1 muon trigger efficiency in cosmic-ray muon data (blue) and signal simulation (red) as a function of $d_0$, for the Level-1 trigger $p_T$ threshold used in the 2018 analysis triggers. The denominator in the efficiency calculation is the number of STA muons with $|\eta| < 1.2$ and $p_T > 28$ GeV.

Fractions of signal events with zero (green), one (blue), and two (red) STA muons matched to TMS muons by the STA-to-TMS muon association procedure, as a function of true $L_{xy}$, in all simulated $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal samples combined. The fractions are computed relative to the number of signal events passing the trigger and containing two STA muons with more than 12 muon detector hits and $p_T > 10$ GeV matched to generated muons from $X \rightarrow \mu \mu$ decays.

Fractions of signal events with zero (green), one (blue), and two (red) STA muons matched to TMS muons by the STA-to-TMS muon association procedure, as a function of true $L_{xy}$, in all simulated $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal samples combined. The fractions are computed relative to the number of signal events passing the trigger and containing two STA muons with more than 12 muon detector hits and $p_T > 10$ GeV matched to generated muons from $X \rightarrow \mu \mu$ decays.

Comparison of the number of events observed in 2016 data in the STA-STA dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legends also include the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{\mu \mu}$ intervals combined.

Comparison of the number of events observed in 2016 data in the STA-STA dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legends also include the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{\mu \mu}$ intervals combined.

Comparison of the number of events observed in 2018 data in the STA-STA dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legends also include the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{\mu \mu}$ intervals combined.

Comparison of the number of events observed in 2018 data in the STA-STA dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legends also include the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{\mu \mu}$ intervals combined.

Comparison of the number of events observed in 2016 data in the STA-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 30 and 60 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legends also include the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{\mu \mu}$ intervals combined.

Comparison of the number of events observed in 2016 data in the STA-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 30 and 60 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legends also include the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{\mu \mu}$ intervals combined.

Comparison of the number of events observed in 2018 data in the STA-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 30 and 60 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legends also include the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{\mu \mu}$ intervals combined.

Comparison of the number of events observed in 2018 data in the STA-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 30 and 60 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legends also include the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{\mu \mu}$ intervals combined.

Comparison of the number of events observed in 2016 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $6 < min(d_0 / \sigma_{d_0}) \leq 10$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2016 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $6 < min(d_0 / \sigma_{d_0}) \leq 10$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2016 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $10 < min(d_0 / \sigma_{d_0}) \leq 20$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2016 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $10 < min(d_0 / \sigma_{d_0}) \leq 20$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2016 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $min(d_0 / \sigma_{d_0}) > 20$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2016 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $min(d_0 / \sigma_{d_0}) > 20$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2018 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $6 < min(d_0 / \sigma_{d_0}) \leq 10$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2018 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $6 < min(d_0 / \sigma_{d_0}) \leq 10$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2018 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $10 < min(d_0 / \sigma_{d_0}) \leq 20$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2018 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $10 < min(d_0 / \sigma_{d_0}) \leq 20$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2018 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $min(d_0 / \sigma_{d_0}) > 20$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2018 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $min(d_0 / \sigma_{d_0}) > 20$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2016 data in the TMS-TMS dimuon category with the expected number of background events, as a function of the smaller of the two $d_0 / \sigma_{d_0}$ values in the TMS-TMS dimuon. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility.

Comparison of the number of events observed in 2016 data in the TMS-TMS dimuon category with the expected number of background events, as a function of the smaller of the two $d_0 / \sigma_{d_0}$ values in the TMS-TMS dimuon. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility.

Comparison of the number of events observed in 2018 data in the TMS-TMS dimuon category with the expected number of background events, as a function of the smaller of the two $d_0 / \sigma_{d_0}$ values in the TMS-TMS dimuon. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility.

Comparison of the number of events observed in 2018 data in the TMS-TMS dimuon category with the expected number of background events, as a function of the smaller of the two $d_0 / \sigma_{d_0}$ values in the TMS-TMS dimuon. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 125\ GeV$ and $m(X) = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 125\ GeV$ and $m(X) = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 125\ GeV$ and $m(X) = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 125\ GeV$ and $m(X) = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 200\ GeV$ and $m(X) = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 200\ GeV$ and $m(X) = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 200\ GeV$ and $m(X) = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 200\ GeV$ and $m(X) = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 400\ GeV$ and $m(X) = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 400\ GeV$ and $m(X) = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 400\ GeV$ and $m(X) = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 400\ GeV$ and $m(X) = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 400\ GeV$ and $m(X) = 150\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 400\ GeV$ and $m(X) = 150\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 1000\ GeV$ and $m(X) = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 1000\ GeV$ and $m(X) = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 1000\ GeV$ and $m(X) = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 1000\ GeV$ and $m(X) = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 1000\ GeV$ and $m(X) = 150\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 1000\ GeV$ and $m(X) = 150\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 1000\ GeV$ and $m(X) = 350\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 1000\ GeV$ and $m(X) = 350\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 10\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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 horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 10\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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 horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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 horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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 horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 30\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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 horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 30\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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 horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 40\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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 horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 40\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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 horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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 horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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 horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 60\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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 horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 60\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; 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 horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

Observed 95% CL exclusion contours in the HAHM model, in the ($m(Z_D)$, $c\tau(Z_D)$) plane. The contours correspond to several representative values of $B(H \rightarrow Z_DZ_D$) ranging from 0.005 to 1%.

Observed 95% CL exclusion contours in the HAHM model, in the ($m(Z_D)$, $c\tau(Z_D)$) plane. The contours correspond to several representative values of $B(H \rightarrow Z_DZ_D$) ranging from 0.005 to 1%.

Observed 95% CL exclusion contours in the HAHM model, in the ($m(Z_D)$, $\epsilon$) plane. The contours correspond to several representative values of $B(H \rightarrow Z_DZ_D$) ranging from 0.005 to 1%.

Observed 95% CL exclusion contours in the HAHM model, in the ($m(Z_D)$, $\epsilon$) plane. The contours correspond to several representative values of $B(H \rightarrow Z_DZ_D$) ranging from 0.005 to 1%.

Background estimation and observed number of events in the STA-STA dimuon category in 2016 and 2018 data. For each probed LLP mass, the chosen mass interval is shown. The mass interval is followed by the estimated and observed counts for the given year. The quoted uncertainties are statistical only.

Background estimations and observed numbers of events in the STA-STA dimuon category in 2016 and 2018 data. For each probed LLP mass, the chosen mass interval is shown, followed by the predicted background yield $N^\text{est}_\text{bkg}$ and the observed number of events $N^\text{obs}$ for the given year. The quoted uncertainties are statistical only.

Background estimation and observed number of events in the TMS-TMS dimuon category in 2016 data. The mass interval is followed by the estimated and observed counts within each $min(d_0 / \sigma_{d_0})$ bin in this mass interval. The quoted uncertainties are statistical only.

Background estimations and observed numbers of events in the TMS-TMS dimuon category in 2016 data. For each mass interval, the table shows the predicted background yield $N^\text{est}_\text{bkg}$ and the observed number of events $N^\text{obs}$ in each of the three $\text{min}(d_0 / \sigma_{d_0})$ bins. The quoted uncertainties are statistical only

Background estimation and observed number of events in the TMS-TMS dimuon category in 2018 data. The mass interval is followed by the estimated and observed counts within each $min(d_0 / \sigma_{d_0})$ bin in this mass interval. The quoted uncertainties are statistical only.

Background estimations and observed numbers of events in the TMS-TMS dimuon category in 2016 data. For each mass interval, the table shows the predicted background yield $N^\text{est}_\text{bkg}$ and the observed number of events $N^\text{obs}$ in each of the three $\text{min}(d_0 / \sigma_{d_0})$ bins. The quoted uncertainties are statistical only

Correspondence between the mass intervals in the TMS-TMS category and the parameters of the simulated signal samples.

Correspondence between the probed LLP masses and the chosen mass intervals in the TMS-TMS category.

Background estimation and observed number of events in the STA-TMS dimuon category in 2016 and 2018 data. For each probed LLP mass, the chosen mass interval is shown. The mass interval is followed by the estimated and observed counts for the given year. The quoted uncertainties are statistical only.

Background estimations and observed numbers of events in the STA-TMS dimuon category in 2016 and 2018 data. For each probed LLP mass, the chosen mass interval is shown, followed by the predicted background yield $N^\text{est}_\text{bkg}$ and the observed number of events $N^\text{obs}$ for the given year. The quoted uncertainties are statistical only.

Number of events passing consecutive sets of selection criteria for 2018 collision data and the signal process $\Phi(125) \rightarrow XX(20\ GeV, c\tau = 13\ cm) \rightarrow \mu\mu$. Each row introduces a new criterion that is applied in addition to the selection of the previous row. In addition to the total number of events, N(events), the event yields of the individual dimuon vertex categories, STA-STA, TMS-TMS, and STA-TMS, are shown in separate columns for each data set. In these columns, events containing selected dimuons of different categories are independently counted for each category.

Number of events passing consecutive sets of selection criteria, in 2018 data and in a sample of simulated $\Phi \rightarrow XX \rightarrow \mu\mu$ signal events with $m(H) = 125\ GeV$, $m(X) = 20\ GeV$, and $c\tau = 13\ cm$. Each row introduces a new criterion that is applied in addition to the selection of the previous row. In addition to the total number of events $N(\text{total})$, the event yields in the individual dimuon categories, STA-STA, TMS-TMS, and STA-TMS, are shown in separate columns for each data set. In these columns, events containing selected dimuons of different categories are counted independently for each category.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 125\ GeV$ and $m(X) = 20\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 125\ GeV$ and $m(X) = 20\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 125\ GeV$ and $m(X) = 50\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 125\ GeV$ and $m(X) = 50\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 200\ GeV$ and $m(X) = 20\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 200\ GeV$ and $m(X) = 20\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 200\ GeV$ and $m(X) = 50\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 200\ GeV$ and $m(X) = 50\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 20\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 20\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 50\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 50\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 150\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 150\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 1000\ GeV$ and $m(X) = 20\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 1\ TeV$ and $m(X) = 20\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 1000\ GeV$ and $m(X) = 50\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 1\ TeV$ and $m(X) = 50\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 1000\ GeV$ and $m(X) = 150\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 1\ TeV$ and $m(X) = 150\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 1000\ GeV$ and $m(X) = 350\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 1\ TeV$ and $m(X) = 350\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 125\ GeV$ and $m(X) = 20\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 125\ GeV$ and $m(X) = 20\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 125\ GeV$ and $m(X) = 50\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 125\ GeV$ and $m(X) = 50\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 200\ GeV$ and $m(X) = 20\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 200\ GeV$ and $m(X) = 20\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 200\ GeV$ and $m(X) = 50\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 200\ GeV$ and $m(X) = 50\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 20\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 20\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 50\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 50\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 150\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 150\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 1000\ GeV$ and $m(X) = 20\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 1\ TeV$ and $m(X) = 20\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 1000\ GeV$ and $m(X) = 50\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 1\ TeV$ and $m(X) = 50\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 1000\ GeV$ and $m(X) = 150\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 1\ TeV$ and $m(X) = 150\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 1000\ GeV$ and $m(X) = 350\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 1\ TeV$ and $m(X) = 350\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 10\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 10\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 20\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 20\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 30\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 30\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 40\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 40\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 50\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 50\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 60\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 60\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well the combined efficiency (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. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Signal efficiencies 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\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2016 samples. 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 applied in the STA-STA dimuon category 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 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\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2018 samples. 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 applied in the STA-STA dimuon category 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 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\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2016 samples. 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 applied in the STA-TMS dimuon category 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 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\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2018 samples. 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 applied in the STA-TMS dimuon category 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 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\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2016 samples. 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 applied in the TMS-TMS dimuon category 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 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\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2018 samples. 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 applied in the TMS-TMS dimuon category 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 as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $20\ cm < L_{xy}^\mathrm{true} < 70\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2016 samples. 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 applied in the STA-STA dimuon category 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 as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $20\ cm < L_{xy}^\mathrm{true} < 70\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2018 samples. 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 applied in the STA-STA dimuon category 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 as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $20\ cm < L_{xy}^\mathrm{true} < 70\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2016 samples. 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 applied in the STA-TMS dimuon category 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 as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $20\ cm < L_{xy}^\mathrm{true} < 70\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2018 samples. 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 applied in the STA-TMS dimuon category 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 as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $20\ cm < L_{xy}^\mathrm{true} < 70\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2016 samples. 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 applied in the TMS-TMS dimuon category 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 as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $20\ cm < L_{xy}^\mathrm{true} < 70\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2018 samples. 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 applied in the TMS-TMS dimuon category 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 as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $70\ cm < L_{xy}^\mathrm{true} < 320\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2016 samples. 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 applied in the STA-STA dimuon category 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. Efficiencies for dimuons with $70\ cm < L_{xy}^\mathrm{true} < 320\ cm$ in the STA-TMS and TMS-TMS dimuon categories are equal to zero.

Signal efficiencies as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $70\ cm < L_{xy}^\mathrm{true} < 320\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2018 samples. 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 applied in the STA-STA dimuon category 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. Efficiencies for dimuons with $70\ cm < L_{xy}^\mathrm{true} < 320\ cm$ in the STA-TMS and TMS-TMS dimuon categories are equal to zero.

Selection acceptance times efficiency as a function of WPRIME mass at truth level for left- and right-handed WPRIME MC. The total efficiency curves correspond to the sum of the efficiencies of the one b-tag and two b-tag categories.

Cutflow (efficiency with respect to total number of events in %) for several WPRIME masses in the left-handed and in the right-handed model for hadronic top-quark decays.

Overview of the signal acceptance times efficiency (A x Eff) for the hadronic top-quark decay for both categories, the predicted cross section times branching ratio to TOP BOTTOM, as well as the observed 95% CL limit on the cross section times branching ratio for several WPRIME masses in the left-handed and in the right-handed model. The expected signal event yields in the two categories for 20.3 fb-1 of 8 TeV data are also shown.

The observed 95% CL limits on the minimal GMSB model parameters Lambda and tan(beta) over the tan(beta) region 0.8 to 6.8.

The observed 95% CL limits on the minimal GMSB model parameters Lambda and tan(beta) over the tan(beta) region 6.8 to 12.8.

The observed 95% CL limits on the minimal GMSB model parameters Lambda and tan(beta) over the tan(beta) region 12.8 to 18.8.

The observed 95% CL limits on the minimal GMSB model parameters Lambda and tan(beta) over the tan(beta) region 18.8 to 24.8.

The observed 95% CL limits on the minimal GMSB model parameters Lambda and tan(beta) over the tan(beta) region 24.8 to 30.8.

The observed 95% CL limits on the minimal GMSB model parameters Lambda and tan(beta) over the tan(beta) region 30.8 to 36.8.

The observed 95% CL limits on the minimal GMSB model parameters Lambda and tan(beta) over the tan(beta) region 36.8 to 42.8.

The observed 95% CL limits on the minimal GMSB model parameters Lambda and tan(beta) over the tan(beta) region 42.8 to 48.8.

The observed 95% CL limits on the minimal GMSB model parameters Lambda and tan(beta) over the tan(beta) region 48.8 to 51.2.

The observed and expected 95% CL limits on the minimal GMSB model parameters Lambda and tan(beta).

The acceptance, efficiency and their product in the Lambda-Tan(Beta)-plane for M_mess=250 TeV, N_5=3 and C_grav=1.

The relative systematic uncertainty from the selection, theoretical and statistical uncertainities and their combined values.

The observed 95% CL upper limits on the visible and production cross sections in the Lambda-Tan(Beta) plane for stau and slepton NLSP production in the minimal GMSB model.