CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-EXO-24-016 ; CERN-EP-2026-132
Search for long-lived particles decaying into muons in proton-proton collisions at $ \sqrt{s} = $ 13.6 TeV using data scouting
Submitted to the Journal of High Energy Physics
Abstract: A search for long-lived particles decaying into muons is performed using proton-proton collisions at $ \sqrt{s}= $ 13.6 TeV, collected by the CMS experiment at the LHC in 2022 and 2023, corresponding to an integrated luminosity of 62.4 fb$ ^{-1} $. The data were collected using dedicated dimuon triggers with low transverse momentum thresholds, recorded with a high-rate data scouting trigger stream. This data stream retains a reduced amount of information at the high-level trigger, to explore otherwise inaccessible phase space at low multimuon invariant mass and nonzero displacement from the primary interaction vertex. No significant excess of events above the standard model prediction is found. Upper limits on branching fractions at 95% confidence level are set for a wide range of mass and lifetime hypotheses in several beyond the standard model frameworks, where the Higgs boson decays into long-lived dark photons or into dark partons that produce showers containing long-lived particles, or where a long-lived scalar resonance is produced from the decay of a b hadron. The resulting constraints improve and extend existing ones in large regions of the parameter space.
Figures & Tables Summary References CMS Publications
Figures

png pdf
Figure 1:
Diagrams illustrating an SM-like Higgs boson (H) decay to four fermions ($ f $) via two intermediate dark photons, $ \mathrm{Z}_{\mathrm{D}} $ [6]: through the hypercharge portal (left) and the Higgs portal (right) via a dark Higgs boson ($ \mathrm{H}_{\mathrm{D}} $). The long-lived particles are highlighted in green.

png pdf
Figure 1-a:
Diagrams illustrating an SM-like Higgs boson (H) decay to four fermions ($ f $) via two intermediate dark photons, $ \mathrm{Z}_{\mathrm{D}} $ [6]: through the hypercharge portal (left) and the Higgs portal (right) via a dark Higgs boson ($ \mathrm{H}_{\mathrm{D}} $). The long-lived particles are highlighted in green.

png pdf
Figure 1-b:
Diagrams illustrating an SM-like Higgs boson (H) decay to four fermions ($ f $) via two intermediate dark photons, $ \mathrm{Z}_{\mathrm{D}} $ [6]: through the hypercharge portal (left) and the Higgs portal (right) via a dark Higgs boson ($ \mathrm{H}_{\mathrm{D}} $). The long-lived particles are highlighted in green.

png pdf
Figure 2:
Diagrams illustrating two dark-shower scenarios [8], where a dark meson $ \tilde{\eta} $ is allowed to decay into dark pions ($ \tilde{\pi} _3 $), which in turn decay into pairs of dark photons ($ A^{\prime} $). Scenario A (left) assumes prompt $ \tilde{\pi} _3 $ decays into a pair of long-lived dark photons, $ A^{\prime} $, each decaying into fermions. Scenario B1 (right) assumes long-lived $ \tilde{\pi} _3 $ decays into a pair of promptly decaying dark photons ($ A^{\prime} $). The long-lived particles are highlighted in green.

png pdf
Figure 2-a:
Diagrams illustrating two dark-shower scenarios [8], where a dark meson $ \tilde{\eta} $ is allowed to decay into dark pions ($ \tilde{\pi} _3 $), which in turn decay into pairs of dark photons ($ A^{\prime} $). Scenario A (left) assumes prompt $ \tilde{\pi} _3 $ decays into a pair of long-lived dark photons, $ A^{\prime} $, each decaying into fermions. Scenario B1 (right) assumes long-lived $ \tilde{\pi} _3 $ decays into a pair of promptly decaying dark photons ($ A^{\prime} $). The long-lived particles are highlighted in green.

png pdf
Figure 2-b:
Diagrams illustrating two dark-shower scenarios [8], where a dark meson $ \tilde{\eta} $ is allowed to decay into dark pions ($ \tilde{\pi} _3 $), which in turn decay into pairs of dark photons ($ A^{\prime} $). Scenario A (left) assumes prompt $ \tilde{\pi} _3 $ decays into a pair of long-lived dark photons, $ A^{\prime} $, each decaying into fermions. Scenario B1 (right) assumes long-lived $ \tilde{\pi} _3 $ decays into a pair of promptly decaying dark photons ($ A^{\prime} $). The long-lived particles are highlighted in green.

png pdf
Figure 3:
Diagram illustrating the production of a scalar resonance $ \phi $ in a b hadron decay, through the mixing with an SM-like Higgs boson (H). The long-lived particle $ \phi $ is highlighted in green.

png pdf
Figure 4:
Overall signal efficiency as measured for the HAHM signal model with $ m_{\mathrm{Z}_{\mathrm{D}}} = 2, 5, $ 12, and 20 GeV, as a function of the mean proper decay length of the dark photon, $ c\tau_{0}^{\mathrm{Z}_{\mathrm{D}}} $. The efficiency is defined as the fraction of events with at least one dark photon decaying into a muon pair, reconstructed as a dimuon candidate within the $ l_\mathrm{xy} $ bins of $ [0, 0.2, 1, 2.4, 3.1, 7, 11, 16, 70] \text{cm} $. The average efficiency over the 2022 and 2023 data-taking periods is shown.

png pdf
Figure 4-a:
Overall signal efficiency as measured for the HAHM signal model with $ m_{\mathrm{Z}_{\mathrm{D}}} = 2, 5, $ 12, and 20 GeV, as a function of the mean proper decay length of the dark photon, $ c\tau_{0}^{\mathrm{Z}_{\mathrm{D}}} $. The efficiency is defined as the fraction of events with at least one dark photon decaying into a muon pair, reconstructed as a dimuon candidate within the $ l_\mathrm{xy} $ bins of $ [0, 0.2, 1, 2.4, 3.1, 7, 11, 16, 70] \text{cm} $. The average efficiency over the 2022 and 2023 data-taking periods is shown.

png pdf
Figure 4-b:
Overall signal efficiency as measured for the HAHM signal model with $ m_{\mathrm{Z}_{\mathrm{D}}} = 2, 5, $ 12, and 20 GeV, as a function of the mean proper decay length of the dark photon, $ c\tau_{0}^{\mathrm{Z}_{\mathrm{D}}} $. The efficiency is defined as the fraction of events with at least one dark photon decaying into a muon pair, reconstructed as a dimuon candidate within the $ l_\mathrm{xy} $ bins of $ [0, 0.2, 1, 2.4, 3.1, 7, 11, 16, 70] \text{cm} $. The average efficiency over the 2022 and 2023 data-taking periods is shown.

png pdf
Figure 4-c:
Overall signal efficiency as measured for the HAHM signal model with $ m_{\mathrm{Z}_{\mathrm{D}}} = 2, 5, $ 12, and 20 GeV, as a function of the mean proper decay length of the dark photon, $ c\tau_{0}^{\mathrm{Z}_{\mathrm{D}}} $. The efficiency is defined as the fraction of events with at least one dark photon decaying into a muon pair, reconstructed as a dimuon candidate within the $ l_\mathrm{xy} $ bins of $ [0, 0.2, 1, 2.4, 3.1, 7, 11, 16, 70] \text{cm} $. The average efficiency over the 2022 and 2023 data-taking periods is shown.

png pdf
Figure 4-d:
Overall signal efficiency as measured for the HAHM signal model with $ m_{\mathrm{Z}_{\mathrm{D}}} = 2, 5, $ 12, and 20 GeV, as a function of the mean proper decay length of the dark photon, $ c\tau_{0}^{\mathrm{Z}_{\mathrm{D}}} $. The efficiency is defined as the fraction of events with at least one dark photon decaying into a muon pair, reconstructed as a dimuon candidate within the $ l_\mathrm{xy} $ bins of $ [0, 0.2, 1, 2.4, 3.1, 7, 11, 16, 70] \text{cm} $. The average efficiency over the 2022 and 2023 data-taking periods is shown.

png pdf
Figure 5:
Four-muon invariant mass distributions of selected events in the multivertex (left) and overlapping vertex (right) categories, for data (black points with uncertainties) and representative signal models (colored histograms).

png pdf
Figure 5-a:
Four-muon invariant mass distributions of selected events in the multivertex (left) and overlapping vertex (right) categories, for data (black points with uncertainties) and representative signal models (colored histograms).

png pdf
Figure 5-b:
Four-muon invariant mass distributions of selected events in the multivertex (left) and overlapping vertex (right) categories, for data (black points with uncertainties) and representative signal models (colored histograms).

png pdf
Figure 6:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} > $ 25 GeV. Distributions are shown separately for each of the six $ l_\mathrm{xy} < 11 \text{cm} $ bins. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 6-a:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} > $ 25 GeV. Distributions are shown separately for each of the six $ l_\mathrm{xy} < 11 \text{cm} $ bins. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 6-b:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} > $ 25 GeV. Distributions are shown separately for each of the six $ l_\mathrm{xy} < 11 \text{cm} $ bins. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 6-c:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} > $ 25 GeV. Distributions are shown separately for each of the six $ l_\mathrm{xy} < 11 \text{cm} $ bins. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 6-d:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} > $ 25 GeV. Distributions are shown separately for each of the six $ l_\mathrm{xy} < 11 \text{cm} $ bins. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 6-e:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} > $ 25 GeV. Distributions are shown separately for each of the six $ l_\mathrm{xy} < 11 \text{cm} $ bins. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 6-f:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} > $ 25 GeV. Distributions are shown separately for each of the six $ l_\mathrm{xy} < 11 \text{cm} $ bins. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 7:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} > $ 25 GeV. Distributions are shown separately for each of the two $ l_\mathrm{xy} > 11 \text{cm} $ bins (i.e.,, for the displacement range inaccessible to the previous results presented in Ref. [18]). The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 7-a:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} > $ 25 GeV. Distributions are shown separately for each of the two $ l_\mathrm{xy} > 11 \text{cm} $ bins (i.e.,, for the displacement range inaccessible to the previous results presented in Ref. [18]). The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 7-b:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} > $ 25 GeV. Distributions are shown separately for each of the two $ l_\mathrm{xy} > 11 \text{cm} $ bins (i.e.,, for the displacement range inaccessible to the previous results presented in Ref. [18]). The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 8:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} < $ 25 GeV. Distributions are shown separately for each of the six $ l_\mathrm{xy} < 11 \text{cm} $ bins. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 8-a:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} < $ 25 GeV. Distributions are shown separately for each of the six $ l_\mathrm{xy} < 11 \text{cm} $ bins. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 8-b:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} < $ 25 GeV. Distributions are shown separately for each of the six $ l_\mathrm{xy} < 11 \text{cm} $ bins. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 8-c:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} < $ 25 GeV. Distributions are shown separately for each of the six $ l_\mathrm{xy} < 11 \text{cm} $ bins. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 8-d:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} < $ 25 GeV. Distributions are shown separately for each of the six $ l_\mathrm{xy} < 11 \text{cm} $ bins. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 8-e:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} < $ 25 GeV. Distributions are shown separately for each of the six $ l_\mathrm{xy} < 11 \text{cm} $ bins. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 8-f:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} < $ 25 GeV. Distributions are shown separately for each of the six $ l_\mathrm{xy} < 11 \text{cm} $ bins. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 9:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} < $ 25 GeV. Distributions are shown separately for each of the two $ l_\mathrm{xy} > 11 \text{cm} $ bins (i.e.,, for the displacement range inaccessible to the previous results presented in Ref. [18]). The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 9-a:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} < $ 25 GeV. Distributions are shown separately for each of the two $ l_\mathrm{xy} > 11 \text{cm} $ bins (i.e.,, for the displacement range inaccessible to the previous results presented in Ref. [18]). The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 9-b:
Dimuon invariant mass distributions of selected events with a single, isolated muon pair with $ p_{\mathrm{T}}^{\mu\mu} < $ 25 GeV. Distributions are shown separately for each of the two $ l_\mathrm{xy} > 11 \text{cm} $ bins (i.e.,, for the displacement range inaccessible to the previous results presented in Ref. [18]). The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. Distributions are shown as obtained in data (black points with uncertainties) and for representative signal models (colored histograms).

png pdf
Figure 10:
The dimuon invariant mass distribution (left) is shown for data collected in 2022 (black markers) in a mass window centered around 5 GeV, in one of the dimuon search bins (0.2 $ < l_\mathrm{xy} < 1 \text{cm} $, $ p_{\mathrm{T}}^{\mu\mu} > $ 25 GeV, with two isolated muons). The results of background+signal fits to the data of three selected functional forms are also shown (different colors, partially overlapping lines), highlighting the Bernstein polynomial as the best fit (solid, purple line) for this mass window. The corresponding pre-fit function (solid blue line) for a representative $ \mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}} $ signal, assuming a cross section of $ \sigma(\mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}}) = $ 100 fb, is shown for visualization purposes. The fit (solid blue line for the total fit function, and dashed red and orange lines for its individual components) to the dimuon invariant mass distribution (right) expected in the same dimuon search bin is shown for a representative $ \mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}} $ signal (black markers) with $ m_{\mathrm{Z}_{\mathrm{D}}}= $ 5 GeV and $ c\tau_{0}^{\mathrm{Z}_{\mathrm{D}}}=1 \text{cm} $.

png pdf
Figure 10-a:
The dimuon invariant mass distribution (left) is shown for data collected in 2022 (black markers) in a mass window centered around 5 GeV, in one of the dimuon search bins (0.2 $ < l_\mathrm{xy} < 1 \text{cm} $, $ p_{\mathrm{T}}^{\mu\mu} > $ 25 GeV, with two isolated muons). The results of background+signal fits to the data of three selected functional forms are also shown (different colors, partially overlapping lines), highlighting the Bernstein polynomial as the best fit (solid, purple line) for this mass window. The corresponding pre-fit function (solid blue line) for a representative $ \mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}} $ signal, assuming a cross section of $ \sigma(\mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}}) = $ 100 fb, is shown for visualization purposes. The fit (solid blue line for the total fit function, and dashed red and orange lines for its individual components) to the dimuon invariant mass distribution (right) expected in the same dimuon search bin is shown for a representative $ \mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}} $ signal (black markers) with $ m_{\mathrm{Z}_{\mathrm{D}}}= $ 5 GeV and $ c\tau_{0}^{\mathrm{Z}_{\mathrm{D}}}=1 \text{cm} $.

png pdf
Figure 10-b:
The dimuon invariant mass distribution (left) is shown for data collected in 2022 (black markers) in a mass window centered around 5 GeV, in one of the dimuon search bins (0.2 $ < l_\mathrm{xy} < 1 \text{cm} $, $ p_{\mathrm{T}}^{\mu\mu} > $ 25 GeV, with two isolated muons). The results of background+signal fits to the data of three selected functional forms are also shown (different colors, partially overlapping lines), highlighting the Bernstein polynomial as the best fit (solid, purple line) for this mass window. The corresponding pre-fit function (solid blue line) for a representative $ \mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}} $ signal, assuming a cross section of $ \sigma(\mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}}) = $ 100 fb, is shown for visualization purposes. The fit (solid blue line for the total fit function, and dashed red and orange lines for its individual components) to the dimuon invariant mass distribution (right) expected in the same dimuon search bin is shown for a representative $ \mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}} $ signal (black markers) with $ m_{\mathrm{Z}_{\mathrm{D}}}= $ 5 GeV and $ c\tau_{0}^{\mathrm{Z}_{\mathrm{D}}}=1 \text{cm} $.

png pdf
Figure 11:
Upper limits at 95% CL on the branching fraction $ \mathcal{B}(\mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}}) $ as functions of the dark-photon mass, $ m_{\mathrm{Z}_{\mathrm{D}}} $, for $ c\tau_{0}^{\mathrm{Z}_{\mathrm{D}}}= $ 0.1 (upper left), 1 (upper right), 10 (lower left), and 100$ \text{cm} $ (lower right). The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all event categories. The dark blue (magenta) solid line represents the observed upper limits previously set in Ref. [18] ( [20]). The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. The blue (magenta) line can overlap with the gray bands, owing to the different mass ranges covered by the results from Ref. [18] ( [20]).

png pdf
Figure 11-a:
Upper limits at 95% CL on the branching fraction $ \mathcal{B}(\mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}}) $ as functions of the dark-photon mass, $ m_{\mathrm{Z}_{\mathrm{D}}} $, for $ c\tau_{0}^{\mathrm{Z}_{\mathrm{D}}}= $ 0.1 (upper left), 1 (upper right), 10 (lower left), and 100$ \text{cm} $ (lower right). The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all event categories. The dark blue (magenta) solid line represents the observed upper limits previously set in Ref. [18] ( [20]). The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. The blue (magenta) line can overlap with the gray bands, owing to the different mass ranges covered by the results from Ref. [18] ( [20]).

png pdf
Figure 11-b:
Upper limits at 95% CL on the branching fraction $ \mathcal{B}(\mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}}) $ as functions of the dark-photon mass, $ m_{\mathrm{Z}_{\mathrm{D}}} $, for $ c\tau_{0}^{\mathrm{Z}_{\mathrm{D}}}= $ 0.1 (upper left), 1 (upper right), 10 (lower left), and 100$ \text{cm} $ (lower right). The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all event categories. The dark blue (magenta) solid line represents the observed upper limits previously set in Ref. [18] ( [20]). The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. The blue (magenta) line can overlap with the gray bands, owing to the different mass ranges covered by the results from Ref. [18] ( [20]).

png pdf
Figure 11-c:
Upper limits at 95% CL on the branching fraction $ \mathcal{B}(\mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}}) $ as functions of the dark-photon mass, $ m_{\mathrm{Z}_{\mathrm{D}}} $, for $ c\tau_{0}^{\mathrm{Z}_{\mathrm{D}}}= $ 0.1 (upper left), 1 (upper right), 10 (lower left), and 100$ \text{cm} $ (lower right). The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all event categories. The dark blue (magenta) solid line represents the observed upper limits previously set in Ref. [18] ( [20]). The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. The blue (magenta) line can overlap with the gray bands, owing to the different mass ranges covered by the results from Ref. [18] ( [20]).

png pdf
Figure 11-d:
Upper limits at 95% CL on the branching fraction $ \mathcal{B}(\mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}}) $ as functions of the dark-photon mass, $ m_{\mathrm{Z}_{\mathrm{D}}} $, for $ c\tau_{0}^{\mathrm{Z}_{\mathrm{D}}}= $ 0.1 (upper left), 1 (upper right), 10 (lower left), and 100$ \text{cm} $ (lower right). The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all event categories. The dark blue (magenta) solid line represents the observed upper limits previously set in Ref. [18] ( [20]). The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. The blue (magenta) line can overlap with the gray bands, owing to the different mass ranges covered by the results from Ref. [18] ( [20]).

png pdf
Figure 12:
Upper limits at 95% CL on the branching fraction $ \mathcal{B}(\mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}}) $ as functions of $ c\tau_{0}^{\mathrm{Z}_{\mathrm{D}}} $, for representative mass hypotheses, $ m_{\mathrm{Z}_{\mathrm{D}}} = $ 2 (upper left), 5 (upper right), 12 (lower left), and 20 (lower right) GeV. The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all event categories. The dark blue (magenta) solid line represents the observed upper limits previously set in Ref. [18] ( [20]).

png pdf
Figure 12-a:
Upper limits at 95% CL on the branching fraction $ \mathcal{B}(\mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}}) $ as functions of $ c\tau_{0}^{\mathrm{Z}_{\mathrm{D}}} $, for representative mass hypotheses, $ m_{\mathrm{Z}_{\mathrm{D}}} = $ 2 (upper left), 5 (upper right), 12 (lower left), and 20 (lower right) GeV. The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all event categories. The dark blue (magenta) solid line represents the observed upper limits previously set in Ref. [18] ( [20]).

png pdf
Figure 12-b:
Upper limits at 95% CL on the branching fraction $ \mathcal{B}(\mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}}) $ as functions of $ c\tau_{0}^{\mathrm{Z}_{\mathrm{D}}} $, for representative mass hypotheses, $ m_{\mathrm{Z}_{\mathrm{D}}} = $ 2 (upper left), 5 (upper right), 12 (lower left), and 20 (lower right) GeV. The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all event categories. The dark blue (magenta) solid line represents the observed upper limits previously set in Ref. [18] ( [20]).

png pdf
Figure 12-c:
Upper limits at 95% CL on the branching fraction $ \mathcal{B}(\mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}}) $ as functions of $ c\tau_{0}^{\mathrm{Z}_{\mathrm{D}}} $, for representative mass hypotheses, $ m_{\mathrm{Z}_{\mathrm{D}}} = $ 2 (upper left), 5 (upper right), 12 (lower left), and 20 (lower right) GeV. The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all event categories. The dark blue (magenta) solid line represents the observed upper limits previously set in Ref. [18] ( [20]).

png pdf
Figure 12-d:
Upper limits at 95% CL on the branching fraction $ \mathcal{B}(\mathrm{H}\to\mathrm{Z}_{\mathrm{D}}\mathrm{Z}_{\mathrm{D}}) $ as functions of $ c\tau_{0}^{\mathrm{Z}_{\mathrm{D}}} $, for representative mass hypotheses, $ m_{\mathrm{Z}_{\mathrm{D}}} = $ 2 (upper left), 5 (upper right), 12 (lower left), and 20 (lower right) GeV. The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all event categories. The dark blue (magenta) solid line represents the observed upper limits previously set in Ref. [18] ( [20]).

png pdf
Figure 13:
Upper limits at 95% CL on the branching fraction $ \mathcal{B}(\mathrm{H}\to {\Psi} \bar{\Psi} ) $, for representative mass hypotheses of the Scenario A dark-shower model. Limits are shown as functions of $ c\tau_{0}^{{A}{\prime} } $, for $ m_{\tilde{\pi} } = $ 2 GeV and $ m_{{A}{\prime} } = $ 0.67 GeV (upper left), $ m_{\tilde{\pi} } = $ 5 GeV and $ m_{{A}{\prime} } = $ 1.67 GeV (upper right), and $ m_{\tilde{\pi} } = $ 7.5 GeV and $ m_{{A}{\prime} } = $ 2.5 GeV (lower). For this model, it is assumed $ \tilde{\Lambda} = m_{\tilde{\eta} } = 4m_{\tilde{\pi} } $, $ \sin\theta = $ 0.1 and $ \mathcal{B}(\tilde{\pi} _3 \to {A}{\prime} {A}{\prime} )= $ 1. The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all event categories. The dark blue solid line represents the observed upper limits previously set in Ref. [21], rescaled to the same Higgs boson production cross section of 59.8 pb [38,39].

png pdf
Figure 13-a:
Upper limits at 95% CL on the branching fraction $ \mathcal{B}(\mathrm{H}\to {\Psi} \bar{\Psi} ) $, for representative mass hypotheses of the Scenario A dark-shower model. Limits are shown as functions of $ c\tau_{0}^{{A}{\prime} } $, for $ m_{\tilde{\pi} } = $ 2 GeV and $ m_{{A}{\prime} } = $ 0.67 GeV (upper left), $ m_{\tilde{\pi} } = $ 5 GeV and $ m_{{A}{\prime} } = $ 1.67 GeV (upper right), and $ m_{\tilde{\pi} } = $ 7.5 GeV and $ m_{{A}{\prime} } = $ 2.5 GeV (lower). For this model, it is assumed $ \tilde{\Lambda} = m_{\tilde{\eta} } = 4m_{\tilde{\pi} } $, $ \sin\theta = $ 0.1 and $ \mathcal{B}(\tilde{\pi} _3 \to {A}{\prime} {A}{\prime} )= $ 1. The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all event categories. The dark blue solid line represents the observed upper limits previously set in Ref. [21], rescaled to the same Higgs boson production cross section of 59.8 pb [38,39].

png pdf
Figure 13-b:
Upper limits at 95% CL on the branching fraction $ \mathcal{B}(\mathrm{H}\to {\Psi} \bar{\Psi} ) $, for representative mass hypotheses of the Scenario A dark-shower model. Limits are shown as functions of $ c\tau_{0}^{{A}{\prime} } $, for $ m_{\tilde{\pi} } = $ 2 GeV and $ m_{{A}{\prime} } = $ 0.67 GeV (upper left), $ m_{\tilde{\pi} } = $ 5 GeV and $ m_{{A}{\prime} } = $ 1.67 GeV (upper right), and $ m_{\tilde{\pi} } = $ 7.5 GeV and $ m_{{A}{\prime} } = $ 2.5 GeV (lower). For this model, it is assumed $ \tilde{\Lambda} = m_{\tilde{\eta} } = 4m_{\tilde{\pi} } $, $ \sin\theta = $ 0.1 and $ \mathcal{B}(\tilde{\pi} _3 \to {A}{\prime} {A}{\prime} )= $ 1. The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all event categories. The dark blue solid line represents the observed upper limits previously set in Ref. [21], rescaled to the same Higgs boson production cross section of 59.8 pb [38,39].

png pdf
Figure 13-c:
Upper limits at 95% CL on the branching fraction $ \mathcal{B}(\mathrm{H}\to {\Psi} \bar{\Psi} ) $, for representative mass hypotheses of the Scenario A dark-shower model. Limits are shown as functions of $ c\tau_{0}^{{A}{\prime} } $, for $ m_{\tilde{\pi} } = $ 2 GeV and $ m_{{A}{\prime} } = $ 0.67 GeV (upper left), $ m_{\tilde{\pi} } = $ 5 GeV and $ m_{{A}{\prime} } = $ 1.67 GeV (upper right), and $ m_{\tilde{\pi} } = $ 7.5 GeV and $ m_{{A}{\prime} } = $ 2.5 GeV (lower). For this model, it is assumed $ \tilde{\Lambda} = m_{\tilde{\eta} } = 4m_{\tilde{\pi} } $, $ \sin\theta = $ 0.1 and $ \mathcal{B}(\tilde{\pi} _3 \to {A}{\prime} {A}{\prime} )= $ 1. The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all event categories. The dark blue solid line represents the observed upper limits previously set in Ref. [21], rescaled to the same Higgs boson production cross section of 59.8 pb [38,39].

png pdf
Figure 14:
Upper limits at 95% CL on the branching fraction $ \mathcal{B}(\mathrm{H}\to {\Psi} \bar{\Psi} ) $, for representative mass hypotheses of the Scenario B1 dark-shower model. Limits are shown as functions of $ c\tau_{0}^{{A}{\prime} } $, for $ m_{\tilde{\pi} } = $ 2 GeV and $ m_{{A}{\prime} } = $ 0.67 GeV (upper left), $ m_{\tilde{\pi} } = $ 5 GeV and $ m_{{A}{\prime} } = $ 1.67 GeV (upper right), and $ m_{\tilde{\pi} } = $ 7.5 GeV and $ m_{{A}{\prime} } = $ 2.5 GeV (lower). For this model, it is assumed $ \tilde{\Lambda} = m_{\tilde{\eta} } = 4m_{\tilde{\pi} } $, $ \sin\theta = $ 0.1 and $ \mathcal{B}(\tilde{\pi} _3 \to {A}{\prime} {A}{\prime} )= $ 1. The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all event categories. The dark blue solid line represents the observed upper limits previously set in Ref. [21], rescaled to the same Higgs boson production cross section of 59.8 pb [38,39].

png pdf
Figure 14-a:
Upper limits at 95% CL on the branching fraction $ \mathcal{B}(\mathrm{H}\to {\Psi} \bar{\Psi} ) $, for representative mass hypotheses of the Scenario B1 dark-shower model. Limits are shown as functions of $ c\tau_{0}^{{A}{\prime} } $, for $ m_{\tilde{\pi} } = $ 2 GeV and $ m_{{A}{\prime} } = $ 0.67 GeV (upper left), $ m_{\tilde{\pi} } = $ 5 GeV and $ m_{{A}{\prime} } = $ 1.67 GeV (upper right), and $ m_{\tilde{\pi} } = $ 7.5 GeV and $ m_{{A}{\prime} } = $ 2.5 GeV (lower). For this model, it is assumed $ \tilde{\Lambda} = m_{\tilde{\eta} } = 4m_{\tilde{\pi} } $, $ \sin\theta = $ 0.1 and $ \mathcal{B}(\tilde{\pi} _3 \to {A}{\prime} {A}{\prime} )= $ 1. The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all event categories. The dark blue solid line represents the observed upper limits previously set in Ref. [21], rescaled to the same Higgs boson production cross section of 59.8 pb [38,39].

png pdf
Figure 14-b:
Upper limits at 95% CL on the branching fraction $ \mathcal{B}(\mathrm{H}\to {\Psi} \bar{\Psi} ) $, for representative mass hypotheses of the Scenario B1 dark-shower model. Limits are shown as functions of $ c\tau_{0}^{{A}{\prime} } $, for $ m_{\tilde{\pi} } = $ 2 GeV and $ m_{{A}{\prime} } = $ 0.67 GeV (upper left), $ m_{\tilde{\pi} } = $ 5 GeV and $ m_{{A}{\prime} } = $ 1.67 GeV (upper right), and $ m_{\tilde{\pi} } = $ 7.5 GeV and $ m_{{A}{\prime} } = $ 2.5 GeV (lower). For this model, it is assumed $ \tilde{\Lambda} = m_{\tilde{\eta} } = 4m_{\tilde{\pi} } $, $ \sin\theta = $ 0.1 and $ \mathcal{B}(\tilde{\pi} _3 \to {A}{\prime} {A}{\prime} )= $ 1. The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all event categories. The dark blue solid line represents the observed upper limits previously set in Ref. [21], rescaled to the same Higgs boson production cross section of 59.8 pb [38,39].

png pdf
Figure 14-c:
Upper limits at 95% CL on the branching fraction $ \mathcal{B}(\mathrm{H}\to {\Psi} \bar{\Psi} ) $, for representative mass hypotheses of the Scenario B1 dark-shower model. Limits are shown as functions of $ c\tau_{0}^{{A}{\prime} } $, for $ m_{\tilde{\pi} } = $ 2 GeV and $ m_{{A}{\prime} } = $ 0.67 GeV (upper left), $ m_{\tilde{\pi} } = $ 5 GeV and $ m_{{A}{\prime} } = $ 1.67 GeV (upper right), and $ m_{\tilde{\pi} } = $ 7.5 GeV and $ m_{{A}{\prime} } = $ 2.5 GeV (lower). For this model, it is assumed $ \tilde{\Lambda} = m_{\tilde{\eta} } = 4m_{\tilde{\pi} } $, $ \sin\theta = $ 0.1 and $ \mathcal{B}(\tilde{\pi} _3 \to {A}{\prime} {A}{\prime} )= $ 1. The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all event categories. The dark blue solid line represents the observed upper limits previously set in Ref. [21], rescaled to the same Higgs boson production cross section of 59.8 pb [38,39].

png pdf
Figure 15:
Upper limits at 95% CL on the branching fraction product $ \mathcal{B}(\mathrm{h}_\mathrm{b} \to \phi X)\mathcal{B}(\phi\to\mu\mu) $ as functions of the long-lived scalar resonance mass, $ m_{\phi} $, for $ c\tau_{0}^{\phi}= $ 0.1 (upper left), 1 (upper right), and 10$ \text{cm} $ (lower). The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all dimuon event categories. The dark blue solid line represents the exclusion limits previously set in Ref. [18]. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. The blue line can overlap with the gray bands, owing to the different mass ranges covered by the results from Ref. [18].

png pdf
Figure 15-a:
Upper limits at 95% CL on the branching fraction product $ \mathcal{B}(\mathrm{h}_\mathrm{b} \to \phi X)\mathcal{B}(\phi\to\mu\mu) $ as functions of the long-lived scalar resonance mass, $ m_{\phi} $, for $ c\tau_{0}^{\phi}= $ 0.1 (upper left), 1 (upper right), and 10$ \text{cm} $ (lower). The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all dimuon event categories. The dark blue solid line represents the exclusion limits previously set in Ref. [18]. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. The blue line can overlap with the gray bands, owing to the different mass ranges covered by the results from Ref. [18].

png pdf
Figure 15-b:
Upper limits at 95% CL on the branching fraction product $ \mathcal{B}(\mathrm{h}_\mathrm{b} \to \phi X)\mathcal{B}(\phi\to\mu\mu) $ as functions of the long-lived scalar resonance mass, $ m_{\phi} $, for $ c\tau_{0}^{\phi}= $ 0.1 (upper left), 1 (upper right), and 10$ \text{cm} $ (lower). The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all dimuon event categories. The dark blue solid line represents the exclusion limits previously set in Ref. [18]. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. The blue line can overlap with the gray bands, owing to the different mass ranges covered by the results from Ref. [18].

png pdf
Figure 15-c:
Upper limits at 95% CL on the branching fraction product $ \mathcal{B}(\mathrm{h}_\mathrm{b} \to \phi X)\mathcal{B}(\phi\to\mu\mu) $ as functions of the long-lived scalar resonance mass, $ m_{\phi} $, for $ c\tau_{0}^{\phi}= $ 0.1 (upper left), 1 (upper right), and 10$ \text{cm} $ (lower). The solid black (dashed black) line represents the observed (median expected) exclusion. The inner blue (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The limits are obtained using the combination of all dimuon event categories. The dark blue solid line represents the exclusion limits previously set in Ref. [18]. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. The blue line can overlap with the gray bands, owing to the different mass ranges covered by the results from Ref. [18].
Tables

png pdf
Table 1:
List of known SM resonances and corresponding masked mass windows, equal to $ \pm5\sigma_{\text{mass}} $ around the mean mass, where mean and resolution ($ \sigma_{\text{mass}} $) are determined from a fit to data.
Summary
A search for displaced multimuon resonances has been performed using proton-proton collisions at a center-of-mass energy of 13.6 TeV, collected by the CMS experiment at the LHC in 2022 and 2023, corresponding to an integrated luminosity of 62.4 fb$^{-1}$. The data were collected using a dedicated dimuon trigger stream with low transverse momentum thresholds, recorded at high rate by retaining a reduced amount of information at the high-level trigger, to explore otherwise inaccessible phase space at low multimuon invariant mass and nonzero displacement from the primary interaction vertex. No significant excess above the standard model prediction is found. The data are used to set constraints on a wide range of mass and lifetime hypotheses for various scenarios of physics beyond the standard model. Three sets of models are considered: the hidden Abelian Higgs model, where the Higgs boson decays into a pair of long-lived dark photons; dark-shower models, where dark showers originate from the Higgs boson decays into dark quarks; and models where a long-lived scalar resonance is produced from the decay of a b hadron. In models with a b hadron, the improved trigger and analysis strategy compensate for the smaller data set, leading to constraints comparable with the previous results. For other models, the constraints obtained are instead the most stringent to date in a large fraction of the explored parameter space. In the hidden Abelian Higgs model, the most stringent limits to date are obtained for dark-photon proper decay lengths greater than 1$ \text{cm} $ and dark-photon masses less than 5 GeV. In the dark-shower models, the most stringent limits to date are set for proper decay lengths larger than 10$ \text{cm} $.
References
1 G. Bertone and J. Silk Particle dark matter: Observations, models and searches Cambridge Univ. Press, Cambridge,, ISBN 978-1-107-65392-4, 2010
link
2 J. L. Feng Dark matter candidates from particle physics and methods of detection Ann. Rev. Astron. Astrophys. 48 (2010) 495 1003.0904
3 T. A. Porter, R. P. Johnson, and P. W. Graham Dark matter searches with astroparticle data Ann. Rev. Astron. Astrophys. 49 (2011) 155 1104.2836
4 Planck Collaboration Planck 2015 results. XIII. Cosmological parameters Astron. Astrophys. 594 (2016) A13 1502.01589
5 R. Essig et al. Working group report: New light weakly coupled particles in: Snowmass on the Mississippi, 2013
Proc. Community Summer Study 201 (2013) 3
1311.0029
6 D. Curtin, R. Essig, S. Gori, and J. Shelton Illuminating dark photons with high-energy colliders JHEP 02 (2015) 157 1412.0018
7 CMS Collaboration Dark sector searches with the CMS experiment Phys. Rept. 1115 (2025) 448 CMS-EXO-23-005
2405.13778
8 S. Born, R. Karur, S. Knapen, and J. Shelton Scouting for dark showers at CMS and LHCb PRD 108 (2023) 035034 2303.04167
9 A. Konaka et al. Search for neutral particles in electron beam dump experiment PRL 57 (1986) 659
10 APEX Collaboration Search for a new gauge boson in electron-nucleus fixed-target scattering by the APEX experiment PRL 107 (2011) 191804 1108.2750
11 BaBar Collaboration Search for dimuon decays of a light scalar boson in radiative transitions $ \Upsilon\rightarrow\gamma A^{0} $ PRL 103 (2009) 081803 0905.4539
12 SINDRUM I Collaboration Search for weakly interacting neutral bosons produced in $ {\pi^{-}\mathrm{p}} $ interactions at rest and decaying into $ {\mathrm{e}^+\mathrm{e}^-} $ pairs PRL 68 (1992) 3845
13 LHCb Collaboration Proposed inclusive dark photon search at LHCb PRL 116 (2016) 251803 1603.08926
14 LHCb Collaboration Search for dark photons produced in 13 TeV pp collisions PRL 120 (2018) 061801 1710.02867
15 LHCb Collaboration Search for $ \mathrm{A}^{\prime}\rightarrow\mu^{+}\mu^{-} $ decays PRL 124 (2020) 041801 1910.06926
16 ATLAS Collaboration Search for light long-lived neutral particles that decay to collimated pairs of leptons or light hadrons in pp collisions at $ \sqrt{s} = $ 13 TeV with the ATLAS detector JHEP 06 (2023) 153 2206.12181
17 CMS Collaboration Search for a narrow resonance lighter than 200 GeV decaying to a pair of muons in proton-proton collisions at $ \sqrt{s}= $ 13 TeV PRL 124 (2020) 131802 CMS-EXO-19-018
1912.04776
18 CMS Collaboration Search for long-lived particles decaying into muon pairs in proton-proton collisions at $ \sqrt{s}= $ 13 TeV collected with a dedicated high-rate data stream JHEP 04 (2022) 062 CMS-EXO-20-014
2112.13769
19 CMS Collaboration Search for direct production of GeV-scale resonances decaying to a pair of muons in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JHEP 12 (2023) 070 CMS-EXO-21-005
2309.16003
20 CMS Collaboration Search for long-lived particles decaying to final states with a pair of muons in proton-proton collisions at $ \sqrt{s} = $ 13.6 TeV JHEP 05 (2024) 047 CMS-EXO-23-014
2402.14491
21 CMS Collaboration Search for low-mass hidden-valley dark showers with non-prompt muon pairs in proton-proton collisions at $ \sqrt{s}= $ 13 TeV JHEP 03 (2026) 189 CMS-EXO-24-008
2511.11888
22 F. Bezrukov and D. Gorbunov Light inflaton after LHC8 and WMAP9 results JHEP 07 (2013) 140 1303.4395
23 J. A. Evans, A. Gandrakota, S. Knapen, and H. Routray Searching for exotic $ {\mathrm{B}} $ meson decays with the CMS L1 track trigger PRD 103 (2021) 015026 2008.06918
24 CHARM Collaboration Search for axion like particle production in 400 GeV proton-copper interactions PLB 157 (1985) 458
25 LHCb Collaboration Search for hidden-sector bosons in $ {\mathrm{B}^0}\to\mathrm{K}^{\ast0}\mu^{+}\mu^{-} $ decays PRL 115 (2015) 161802 1508.04094
26 LHCb Collaboration Search for long-lived scalar particles in $ {\mathrm{B}^{+}}\to\mathrm{K^+}\chi\left(\mu^{+}\mu^{-}\right) $ decays PRD 95 (2017) 071101 1612.07818
27 CMS Collaboration Enriching the physics program of the CMS experiment via data scouting and data parking Phys. Rept. 1115 678, 2025 CMS-EXO-23-007
2403.16134
28 CMS Collaboration HEPData record for this analysis link
29 CMS Collaboration Particle-flow reconstruction and global event description with the CMS detector JINST 12 (2017) P10003 CMS-PRF-14-001
1706.04965
30 Tracker Group of the CMS Collaboration The CMS phase-1 pixel detector upgrade JINST 16 (2021) P02027 2012.14304
31 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004
32 CMS Collaboration Development of the CMS detector for the CERN LHC Run 3 JINST 19 (2024) P05064 CMS-PRF-21-001
2309.05466
33 CMS Collaboration Performance of the CMS Level-1 trigger in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JINST 15 (2020) P10017 CMS-TRG-17-001
2006.10165
34 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
35 CMS Collaboration Performance of the CMS high-level trigger during LHC Run 2 JINST 19 (2024) P11021 CMS-TRG-19-001
2410.17038
36 CMS Collaboration Strategy and performance of the CMS long-lived particle trigger program in proton-proton collisions at $ \sqrt{s} = $ 13.6 TeV Submitted to Phys. Rept, 2026 CMS-EXO-23-016
2601.17544
37 J. Alwall et al. The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations JHEP 07 (2014) 079 1405.0301
38 M. Cepeda et al. Report from working group 2: Higgs physics at the HL-LHC and HE-LHC CERN Yellow Rep. Monogr. 7 (2019) 221 1902.00134
39 A. Karlberg et al. Ad interim recommendations for the Higgs boson production cross sections at $ \sqrt{s} = $ 13.6 TeV 2402.09955
40 C. Bierlich et al. A comprehensive guide to the physics and usage of PYTHIA 8.3 SciPost Phys. Codeb. 2022 (2022) 8 2203.11601
41 CMS Collaboration Extraction and validation of a new set of CMS PYTHIA 8 tunes from underlying-event measurements EPJC 80 (2020) 4 CMS-GEN-17-001
1903.12179
42 NNPDF Collaboration Parton distributions from high-precision collider data EPJC 77 (2017) 663 1706.00428
43 \GEANTfour Collaboration GEANT 4---a simulation toolkit NIM A 506 (2003) 250
44 M. J. Oreglia A study of the reactions $ \psi^\prime \rightarrow \gamma\gamma \psi $ PhD thesis, Stanford University, SLAC Report SLAC-R-236, see Appendix D, 1980
link
45 J. E. Gaiser Charmonium spectroscopy from radiative decays of the $ \mathrm{J}/\psi $ and $ \psi^\prime $ PhD thesis, Stanford University, SLAC Report SLAC-R-255, 1982
link
46 R. A. Fisher On the interpretation of $ \chi^{2} $ from contingency tables, and the calculation of P J. R. Stat. Soc 85 (1922) 87
47 E. Gross and O. Vitells Trial factors for the look elsewhere effect in high energy physics EPJC 70 (2010) 525 1005.1891
48 P. D. Dauncey, M. Kenzie, N. Wardle, and G. J. Davies Handling uncertainties in background shapes: the discrete profiling method JINST 10 (2015) P04015 1408.6865
49 T. Junk Confidence level computation for combining searches with small statistics NIM A 434 (1999) 435 hep-ex/9902006
50 A. L. Read Presentation of search results: The $ \text{CL}_\text{s} $ technique JPG 28 (2002) 2693
51 ATLAS and CMS Collaborations Procedure for the LHC Higgs boson search combination in summer 2011 ATL-PHYS-PUB-2011-011, CMS NOTE-2011/005, 2011
link
52 CMS Collaboration The CMS statistical analysis and combination tool: Combine Comput. Softw. Big Sci. 8 (2024) 19 CMS-CAT-23-001
2404.06614
53 W. Verkerke and D. Kirkby The RooFit toolkit for data modeling in the International Conference on Computing in High Energy and Nuclear Physics (CHEP ): La Jolla CA, United States, March 24--28,, 2003
Proc. 1 (2003) 3
physics/0306116
54 L. Moneta et al. The RooStats project in the International Workshop on Advanced Computing and Analysis Techniques in Physics Research (ACAT ): Jaipur, India, February 22--27,, 2010
Proc. 1 (2010) 3
1009.1003
55 CMS Collaboration Precision luminosity measurement in proton-proton collisions at $ \sqrt{s}= $ 13 TeV in 2015 and 2016 at CMS EPJC 81 (2021) 800 CMS-LUM-17-003
2104.01927
56 CMS Collaboration Luminosity measurement in proton-proton collisions at 13.6 TeV in 2022 at CMS CMS Physics Analysis Summary,, 2024
CMS-PAS-LUM-22-001
CMS-PAS-LUM-22-001
57 CMS Collaboration Measurement of the offline integrated luminosity for the CMS proton-proton collision dataset recorded in 2023 CMS Detector Performance Note, CMS-DP-2024-068, 2024
CDS
Compact Muon Solenoid
LHC, CERN