CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-PAS-EXO-20-014
Search for long-lived particles decaying into two muons in proton-proton collisions at $\sqrt{s}= $ 13 TeV using data collected with high rate triggers
Abstract: A search for displaced dimuon resonances is performed using proton-proton collisions at a center-of-mass energy of 13 TeV, collected by the CMS experiment at the LHC in 2017--2018, corresponding to an integrated luminosity of 101 fb$^{-1}$. The data sets used in this search were collected using a dedicated dimuon trigger stream with low transverse momentum thresholds, recorded at high rate by retaining a reduced amount of trigger-level information, in order to explore otherwise inaccessible phase space at low dimuon mass and non-zero displacement from the interaction point. We find no significant excess, and use the data to set stringent constraints on a wide range of mass and lifetime hypotheses for models of physics beyond the standard model where a Higgs boson decays to a pair of long-lived dark photons, or where a long-lived scalar resonance arises from the decay of a B hadron.
Figures Summary Additional Figures & Tables References CMS Publications
Figures

png pdf
Figure 1:
Diagrams illustrating a SM-like Higgs boson (h) decay to four leptons via one or two intermediate $\mathrm{Z} _{\mathrm {D}}$ [6]: (left) $\mathrm{h} \rightarrow \mathrm{Z} \mathrm{Z} _{\mathrm {D}}\rightarrow 4\ell $, through the hypercharge portal; (right) $\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}}\rightarrow 4\ell $, through the Higgs portal.

png pdf
Figure 1-a:
Diagrams illustrating a SM-like Higgs boson (h) decay to four leptons via one or two intermediate $\mathrm{Z} _{\mathrm {D}}$ [6]: (left) $\mathrm{h} \rightarrow \mathrm{Z} \mathrm{Z} _{\mathrm {D}}\rightarrow 4\ell $, through the hypercharge portal; (right) $\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}}\rightarrow 4\ell $, through the Higgs portal.

png pdf
Figure 1-b:
Diagrams illustrating a SM-like Higgs boson (h) decay to four leptons via one or two intermediate $\mathrm{Z} _{\mathrm {D}}$ [6]: (left) $\mathrm{h} \rightarrow \mathrm{Z} \mathrm{Z} _{\mathrm {D}}\rightarrow 4\ell $, through the hypercharge portal; (right) $\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}}\rightarrow 4\ell $, through the Higgs portal.

png pdf
Figure 2:
Diagram illustrating the production of a scalar resonance $\phi $ via a B hadron decay.

png pdf
Figure 3:
The dimuon invariant mass distribution is shown in bins of ${l_\mathrm {xy}}$ as obtained from selected dimuon events in data where both muons are isolated with $ {{p_{\mathrm {T}}} ^{\mu \mu}} < $ 25 GeV: (upper left) 0.0 $ < {l_\mathrm {xy}} < $ 0.2 cm; (upper right) 0.2 $ < {l_\mathrm {xy}} < $ 1.0 cm; (middle left) 1.0 $ < {l_\mathrm {xy}} < $ 2.4 cm; (middle right) 2.4 $ < {l_\mathrm {xy}} < $ 3.1 cm; (lower left) 3.1 $ < {l_\mathrm {xy}} < $ 7.0 cm; (lower right) 7.0 $ < {l_\mathrm {xy}} < $ 1.0 cm. The distribution expected for representative signal models is overlaid. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Figure 4:
The dimuon invariant mass distribution is shown in bins of ${l_\mathrm {xy}}$ as obtained from selected dimuon events in data where both muons are isolated with $ {{p_{\mathrm {T}}} ^{\mu \mu}} \geq $ 25 GeV: (upper left) 0.0 $ < {l_\mathrm {xy}} < $ 0.2 cm; (upper right) 0.2 $ < {l_\mathrm {xy}} < $ 1.0 cm; (middle left) 1.0 $ < {l_\mathrm {xy}} < $ 2.4 cm; (middle right) 2.4 $ < {l_\mathrm {xy}} < $ 3.1 cm; (lower left) 3.1 $ < {l_\mathrm {xy}} < $ 7.0 cm; (lower right) 7.0 $ < {l_\mathrm {xy}} < $ 1.0 cm. The distribution expected for representative signal models is overlaid. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Figure 5:
The dimuon invariant mass distribution is shown in a mass window around 5 GeV, in one of the dimuon search bins (0.2 $ < {l_\mathrm {xy}} < $ 1.0 cm, $ {{p_{\mathrm {T}}} ^{\mu \mu}} < $ 25 GeV, with two isolated muons). The result of the background-only fit to data is also shown together with the dimuon invariant mass distribution expected for a representative signal model (green), with $m_{\mathrm {LLP}} = $ 5 GeV and $c\tau _{0}^{\mathrm {LLP}}=$ 1 mm.

png pdf
Figure 6:
Exclusion limits at 95% CL on the branching fraction $\mathcal {B}(\mathrm{B} \rightarrow \phi X) \times \mathcal {B}(\phi \rightarrow \mu \mu)$ as a function of the signal mass ($m_{\phi}$) hypothesis for (upper left) $c\tau _{0}^{\phi}=$ 1 mm, (upper right) $c\tau _{0}^{\phi}=$ 10 mm and (lower) $c\tau _{0}^{\phi}=$ 100 mm, for the $\mathrm{B} \rightarrow \phi X$ signal model. The solid (dashed) black line represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Figure 6-a:
Exclusion limits at 95% CL on the branching fraction $\mathcal {B}(\mathrm{B} \rightarrow \phi X) \times \mathcal {B}(\phi \rightarrow \mu \mu)$ as a function of the signal mass ($m_{\phi}$) hypothesis for (upper left) $c\tau _{0}^{\phi}=$ 1 mm, (upper right) $c\tau _{0}^{\phi}=$ 10 mm and (lower) $c\tau _{0}^{\phi}=$ 100 mm, for the $\mathrm{B} \rightarrow \phi X$ signal model. The solid (dashed) black line represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Figure 6-b:
Exclusion limits at 95% CL on the branching fraction $\mathcal {B}(\mathrm{B} \rightarrow \phi X) \times \mathcal {B}(\phi \rightarrow \mu \mu)$ as a function of the signal mass ($m_{\phi}$) hypothesis for (upper left) $c\tau _{0}^{\phi}=$ 1 mm, (upper right) $c\tau _{0}^{\phi}=$ 10 mm and (lower) $c\tau _{0}^{\phi}=$ 100 mm, for the $\mathrm{B} \rightarrow \phi X$ signal model. The solid (dashed) black line represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Figure 6-c:
Exclusion limits at 95% CL on the branching fraction $\mathcal {B}(\mathrm{B} \rightarrow \phi X) \times \mathcal {B}(\phi \rightarrow \mu \mu)$ as a function of the signal mass ($m_{\phi}$) hypothesis for (upper left) $c\tau _{0}^{\phi}=$ 1 mm, (upper right) $c\tau _{0}^{\phi}=$ 10 mm and (lower) $c\tau _{0}^{\phi}=$ 100 mm, for the $\mathrm{B} \rightarrow \phi X$ signal model. The solid (dashed) black line represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Figure 7:
Exclusion limits at 95% CL on the branching fraction $\mathcal {B}(\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}}) \times \mathcal {B}(\mathrm{Z} _{\mathrm {D}}\rightarrow \mu \mu)$ as a function of the signal mass ($m_{\mathrm{Z} _{\mathrm {D}}}$) hypothesis for (upper left) $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 1 mm, (upper right) $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 10 mm and $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 100 mm, for the $\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}}$ signal model. The solid (dashed) black line represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Figure 7-a:
Exclusion limits at 95% CL on the branching fraction $\mathcal {B}(\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}}) \times \mathcal {B}(\mathrm{Z} _{\mathrm {D}}\rightarrow \mu \mu)$ as a function of the signal mass ($m_{\mathrm{Z} _{\mathrm {D}}}$) hypothesis for (upper left) $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 1 mm, (upper right) $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 10 mm and $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 100 mm, for the $\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}}$ signal model. The solid (dashed) black line represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Figure 7-b:
Exclusion limits at 95% CL on the branching fraction $\mathcal {B}(\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}}) \times \mathcal {B}(\mathrm{Z} _{\mathrm {D}}\rightarrow \mu \mu)$ as a function of the signal mass ($m_{\mathrm{Z} _{\mathrm {D}}}$) hypothesis for (upper left) $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 1 mm, (upper right) $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 10 mm and $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 100 mm, for the $\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}}$ signal model. The solid (dashed) black line represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Figure 7-c:
Exclusion limits at 95% CL on the branching fraction $\mathcal {B}(\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}}) \times \mathcal {B}(\mathrm{Z} _{\mathrm {D}}\rightarrow \mu \mu)$ as a function of the signal mass ($m_{\mathrm{Z} _{\mathrm {D}}}$) hypothesis for (upper left) $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 1 mm, (upper right) $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 10 mm and $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 100 mm, for the $\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}}$ signal model. The solid (dashed) black line represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Figure 8:
Exclusion limits at 95% CL on the branching fraction $\mathcal {B}(\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}})$ as a function of the signal mass ($m_{\mathrm{Z} _{\mathrm {D}}}$) hypothesis for (upper left) $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 1 mm, (upper right) $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 10 mm and $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 100 mm, for the $\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}}$ signal model. The solid (dashed) black line represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Figure 8-a:
Exclusion limits at 95% CL on the branching fraction $\mathcal {B}(\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}})$ as a function of the signal mass ($m_{\mathrm{Z} _{\mathrm {D}}}$) hypothesis for (upper left) $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 1 mm, (upper right) $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 10 mm and $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 100 mm, for the $\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}}$ signal model. The solid (dashed) black line represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Figure 8-b:
Exclusion limits at 95% CL on the branching fraction $\mathcal {B}(\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}})$ as a function of the signal mass ($m_{\mathrm{Z} _{\mathrm {D}}}$) hypothesis for (upper left) $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 1 mm, (upper right) $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 10 mm and $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 100 mm, for the $\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}}$ signal model. The solid (dashed) black line represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Figure 8-c:
Exclusion limits at 95% CL on the branching fraction $\mathcal {B}(\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}})$ as a function of the signal mass ($m_{\mathrm{Z} _{\mathrm {D}}}$) hypothesis for (upper left) $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 1 mm, (upper right) $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 10 mm and $c\tau _{0}^{\mathrm{Z} _{\mathrm {D}}}=$ 100 mm, for the $\mathrm{h} \rightarrow \mathrm{Z} _{\mathrm {D}}\mathrm{Z} _{\mathrm {D}}$ signal model. The solid (dashed) black line represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing 68 (95)% of the distribution of limits expected under the background-only hypothesis. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.
Summary
We presented a search for displaced dimuon resonances using proton-proton collision at a center-of-mass energy of 13 TeV, collected by the CMS experiment at the LHC in 2017--2018, corresponding to an integrated luminosity of 101 fb$^{-1}$. The data sets used in this search were collected using a dedicated dimuon trigger stream with low transverse momentum thresholds, recorded at high rate by retaining a reduced amount of trigger-level information, in order to explore otherwise inaccessible phase space at low dimuon mass and non-zero displacement from the interaction point. We found no significant excess, and used the data to set constraints on a wide range of mass and lifetime hypotheses for models of physics beyond the standard model where a Higgs boson decays to a pair of long-lived dark photons, or where a long-lived scalar resonance arises from the decay of a B hadron.
Additional Figures

png pdf
Additional Figure 1:
Distribution of ${{p_{\mathrm {T}}} ^{{\mu} {\mu}}}$ for data events and illustrative benchmark signal models, after applying the full event selection.

png pdf
Additional Figure 2:
Distribution of $\Delta (\mu\mu, \vec{\mathrm{DV}})$ for data events and illustrative benchmark signal models, after applying the trigger selections and muon and DV quality criteria.

png pdf
Additional Figure 3:
Distribution of ${\log_{10}(|\Delta \eta _{\mu \mu}|/|\Delta \phi _{\mu \mu}|)}$ in data events and illustrative benchmark signal models. The cut value is indicated by the red vertical dashed line. All selections except those involving muon isolation, ${d_\mathrm {xy}}$, and number of excess hits in the pixel tracker for each muon are applied.

png pdf
Additional Figure 4:
Distribution of ${{| {d_\mathrm {xy}} |}/({l_\mathrm {xy}} {m_{{\mu} {\mu}}} / {{p_{\mathrm {T}}} ^{{\mu} {\mu}}})}$ in data events and illustrative benchmark signal models. The cut value is indicated by the red vertical dashed line. All selections except those involving ${d_\mathrm {xy}}$ are applied.

png pdf
Additional Figure 5:
Distribution of the number of excess pixel hits (number of observed hits minus expected hits) in data and illustrative benchmark signal models for the trailing muon. The cut value is indicated by the red vertical dashed line. All selections except those involving muon isolation, ${d_\mathrm {xy}}$, and ${\log_{10}(|\Delta \eta _{\mu \mu}|/|\Delta \phi _{\mu \mu}|)}$ are applied.

png pdf
Additional Figure 6:
Distribution of DV position in the $x-y$ plane in data events passing trigger selections before a pixel tracker material veto is applied.

png pdf
Additional Figure 7:
Distribution of DV position in the $x-y$ plane in data events passing trigger selections after a pixel tracker material veto is applied.

png pdf
Additional Figure 8:
Distribution of ${l_\mathrm {xy}}$ for data events passing trigger selections before (black) and after (red) a pixel tracker material veto is applied.

png pdf
Additional Figure 9:
The dimuon invariant mass distribution is shown in bins of ${l_\mathrm {xy}}$ as obtained from all selected dimuon events.

png pdf
Additional Figure 10:
The dimuon invariant mass distribution is shown in bins of ${l_\mathrm {xy}}$ as obtained from selected dimuon events in data where none of the two muons are isolated with $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} < $ 25 GeV: (upper left) 0.0 $ < {l_\mathrm {xy}} < 0.2 cm $; (upper right) 0.2 $ < {l_\mathrm {xy}} < 1.0 cm $; (middle left) 1.0 $ < {l_\mathrm {xy}} < 2.4 cm $; (middle right) 2.4 $ < {l_\mathrm {xy}} < 3.1 cm $; (lower left) 3.1 $ < {l_\mathrm {xy}} < 7.0 cm $; (lower right) 7.0 $ < {l_\mathrm {xy}} < $ 11.0 cm. The distribution expected for representative signal models is overlaid. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 11:
The dimuon invariant mass distribution is shown in bins of ${l_\mathrm {xy}}$ as obtained from selected dimuon events in data where none of the two muons are isolated with $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 25 GeV: (upper left) 0.0 $ < {l_\mathrm {xy}} < 0.2 cm $; (upper right) 0.2 $ < {l_\mathrm {xy}} < 1.0 cm $; (middle left) 1.0 $ < {l_\mathrm {xy}} < 2.4 cm $; (middle right) 2.4 $ < {l_\mathrm {xy}} < 3.1 cm $; (lower left) 3.1 $ < {l_\mathrm {xy}} < 7.0 cm $; (lower right) 7.0 $ < {l_\mathrm {xy}} < $ 11.0 cm. The distribution expected for representative signal models is overlaid. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 12:
The dimuon invariant mass distribution is shown in bins of ${l_\mathrm {xy}}$ as obtained from selected dimuon events in data where only one muon is isolated with $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} < $ 25 GeV: (upper left) 0.0 $ < {l_\mathrm {xy}} < 0.2 cm $; (upper right) 0.2 $ < {l_\mathrm {xy}} < 1.0 cm $; (middle left) 1.0 $ < {l_\mathrm {xy}} < 2.4 cm $; (middle right) 2.4 $ < {l_\mathrm {xy}} < 3.1 cm $; (lower left) 3.1 $ < {l_\mathrm {xy}} < 7.0 cm $; (lower right) 7.0 $ < {l_\mathrm {xy}} < $ 11.0 cm. The distribution expected for representative signal models is overlaid. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 13:
The dimuon invariant mass distribution is shown in bins of ${l_\mathrm {xy}}$ as obtained from selected dimuon events in data where only one muon is isolated with $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 25 GeV: (upper left) 0.0 $ < {l_\mathrm {xy}} < 0.2 cm $; (upper right) 0.2 $ < {l_\mathrm {xy}} < 1.0 cm $; (middle left) 1.0 $ < {l_\mathrm {xy}} < 2.4 cm $; (middle right) 2.4 $ < {l_\mathrm {xy}} < 3.1 cm $; (lower left) 3.1 $ < {l_\mathrm {xy}} < 7.0 cm $; (lower right) 7.0 $ < {l_\mathrm {xy}} < $ 11.0 cm. The distribution expected for representative signal models is overlaid. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 14:
Observed exclusion limits at 95% CL on the branching fraction $\mathcal {B}({\mathrm {B}} \rightarrow \phi X) \times \mathcal {B}(\phi \rightarrow \mu \mu)$ as a function of $c\tau _{0}^{\phi}$ for various benchmark $m_{\phi}$ hypotheses, for the $ {\mathrm {B}} \rightarrow \phi X$ signal model. The upper limits are obtained using the combination of all dimuon event categories.

png pdf
Additional Figure 15:
Observed exclusion limits at 95% CL on the branching fraction $\mathcal {B}({\mathrm {h}} \rightarrow {\mathrm {Z}} _{\mathrm {D}} {\mathrm {Z}} _{\mathrm {D}})$ as a function of $c\tau _{0}^{{\mathrm {Z}} _{\mathrm {D}}}$ for various benchmark $m_{{\mathrm {Z}} _{\mathrm {D}}}$ hypotheses, for the $ {\mathrm {h}} \rightarrow {\mathrm {Z}} _{\mathrm {D}} {\mathrm {Z}} _{\mathrm {D}}$ signal model. The upper limits are obtained using the combination of all dimuon event categories and the four-muon event category.

png pdf
Additional Figure 16:
Observed exclusion limits at 95% CL on the branching fraction $\mathcal {B}({\mathrm {B}} \rightarrow \phi X) \times \mathcal {B}(\phi \rightarrow \mu \mu)$ as a function of $m_{\phi}$ for $c\tau _{0}^{\phi}=1 mm $, for the $ {\mathrm {B}} \rightarrow \phi X$ signal model. The black line represents the observed exclusion when the inclusive $ {\mathrm {B}} $ hadron production process is considered. The blue (magenta) line represents the observed exclusion after scaling the signal acceptance by the fraction of $ {\mathrm {B}} ^{\pm}$ ($ {\mathrm {B}} ^{0}$) hadrons at generator level. The exclusion limits are compared to recent LHCb results on the corresponding exclusive production and decay channels, represented by the triangular markers. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. The upper limits are obtained using the combination of all dimuon event categories.

png pdf
Additional Figure 17:
Observed exclusion limits at 95% CL on the branching fraction $\mathcal {B}({\mathrm {B}} \rightarrow \phi X) \times \mathcal {B}(\phi \rightarrow \mu \mu)$ as a function of $m_{\phi}$ for $c\tau _{0}^{\phi}=10 mm $, for the $ {\mathrm {B}} \rightarrow \phi X$ signal model. The black line represents the observed exclusion when the inclusive $ {\mathrm {B}} $ hadron production process is considered. The blue (magenta) line represents the observed exclusion after scaling the signal acceptance by the fraction of $ {\mathrm {B}} ^{\pm}$ ($ {\mathrm {B}} ^{0}$) hadrons at generator level. The exclusion limits are compared to recent LHCb results on the corresponding exclusive production and decay channels, represented by the triangular markers. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. The upper limits are obtained using the combination of all dimuon event categories.

png pdf
Additional Figure 18:
Observed exclusion limits at 95% CL on the branching fraction $\mathcal {B}({\mathrm {B}} \rightarrow \phi X) \times \mathcal {B}(\phi \rightarrow \mu \mu)$ as a function of $m_{\phi}$ for $c\tau _{0}^{\phi}=100 mm $, for the $ {\mathrm {B}} \rightarrow \phi X$ signal model. The black line represents the observed exclusion when the inclusive $ {\mathrm {B}} $ hadron production process is considered. The blue (magenta) line represents the observed exclusion after scaling the signal acceptance by the fraction of $ {\mathrm {B}} ^{\pm}$ ($ {\mathrm {B}} ^{0}$) hadrons at generator level. The exclusion limits are compared to recent LHCb results on the corresponding exclusive production and decay channels, represented by the triangular markers. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. The upper limits are obtained using the combination of all dimuon event categories.

png pdf
Additional Figure 19:
Observed exclusion limits at 95% CL on the branching fraction $\mathcal {B}({\mathrm {B}} \rightarrow \phi X) \times \mathcal {B}(\phi \rightarrow \mu \mu)$ as a function of $c\tau _{0}^{\phi}$ and $m_{\phi}$, for the $ {\mathrm {B}} \rightarrow \phi X$ signal model. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. The upper limits are obtained using the combination of all dimuon event categories.

png pdf
Additional Figure 20:
Observed exclusion limits at 95% CL on the branching fraction $\mathcal {B}({\mathrm {h}} \rightarrow {\mathrm {Z}} _{\mathrm {D}} {\mathrm {Z}} _{\mathrm {D}})$ as a function of $c\tau _{0}^{{\mathrm {Z}} _{\mathrm {D}}}$ and $m_{{\mathrm {Z}} _{\mathrm {D}}}$, for the $ {\mathrm {h}} \rightarrow {\mathrm {Z}} _{\mathrm {D}} {\mathrm {Z}} _{\mathrm {D}}$ signal model. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. The upper limits are obtained using the combination of all dimuon event categories and the four-muon event category.

png pdf
Additional Figure 21:
Observed exclusion limits at 95% CL on the branching fraction $\mathcal {B}({\mathrm {h}} \rightarrow {\mathrm {Z}} _{\mathrm {D}} {\mathrm {Z}} _{\mathrm {D}})$ as a function of the kinetic mixing parameter $\epsilon $ and $m_{{\mathrm {Z}} _{\mathrm {D}}}$, for the $ {\mathrm {h}} \rightarrow {\mathrm {Z}} _{\mathrm {D}} {\mathrm {Z}} _{\mathrm {D}}$ signal model. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search. The upper limits are obtained using the combination of all dimuon event categories and the four-muon event category.

png pdf
Additional Figure 22:
Trigger efficiency as a function of ${l_\mathrm {xy}}$ and trailing muon ${p_{\mathrm {T}}}$.

png pdf
Additional Figure 23:
Selection efficiency after trigger requirements as a function of generated ${l_\mathrm {xy}}$ and generated ${{p_{\mathrm {T}}} ^{{\mu} {\mu}}}$ for muon pairs with dimuon mass between 0 and 3 GeV.

png pdf
Additional Figure 24:
Selection efficiency after trigger requirements as a function of generated ${l_\mathrm {xy}}$ and generated ${{p_{\mathrm {T}}} ^{{\mu} {\mu}}}$ for muon pairs with dimuon mass between 3 and 10 GeV.

png pdf
Additional Figure 25:
Selection efficiency after trigger requirements as a function of generated ${l_\mathrm {xy}}$ and generated ${{p_{\mathrm {T}}} ^{{\mu} {\mu}}}$ for muon pairs with dimuon mass above 10 GeV.

png pdf
Additional Figure 26:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 0.2 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 0 GeV and no requirement on muon isolation. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 27:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 0.2 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 0 GeV and two isolated muons. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 28:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 0.2 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 25 GeV and no requirement on muon isolation. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 29:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 0.2 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 25 GeV and two isolated muons. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 30:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 1.0 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 0 GeV and no requirement on muon isolation. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 31:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 1.0 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 0 GeV and two isolated muons. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 32:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 1.0 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 25 GeV and no requirement on muon isolation. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 33:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 1.0 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 25 GeV and two isolated muons. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 34:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 2.4 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 0 GeV and no requirement on muon isolation. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 35:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 2.4 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 0 GeV and two isolated muons. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 36:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 2.4 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 25 GeV and no requirement on muon isolation. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 37:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 2.4 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 25 GeV and two isolated muons. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 38:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 3.1 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 0 GeV and no requirement on muon isolation. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 39:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 3.1 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 0 GeV and two isolated muons. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 40:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 3.1 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 25 GeV and no requirement on muon isolation. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 41:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 3.1 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 25 GeV and two isolated muons. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 42:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 7.0 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 0 GeV and no requirement on muon isolation. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 43:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 7.0 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 0 GeV and two isolated muons. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 44:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 7.0 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 25 GeV and no requirement on muon isolation. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.

png pdf
Additional Figure 45:
Observed upper limit at 95% CL on number of signal events in aggregate signal region with 7.0 $ < {l_\mathrm {xy}} < $ 11.0 cm, $ {{p_{\mathrm {T}}} ^{{\mu} {\mu}}} \geq $ 25 GeV and two isolated muons. The vertical gray bands indicate mass ranges containing known SM resonances, which are masked for the purpose of this search.
Additional Tables

png pdf
Additional Table 1:
List of known resonances and corresponding masked mass windows, equal to $ \pm $5$ \sigma $ around the mean mass, where mean and resolution ($\sigma $) are determined from a fit to data.

png pdf
Additional Table 2:
List of aggregate signal regions for reinterpretation of the results. The aggregate signal regions are not mutually exclusive.
References
1 G. Bertone and J. Silk Particle Dark Matter: Observations, Models and Searches Cambridge Univ. Press, Cambridge, 2010 , ISBN 978-1-107-65392-4
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 Community Summer Study 2013: Snowmass on the Mississippi 10 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 A. Konaka et al. Search for Neutral Particles in Electron Beam Dump Experiment PRL 57 (1986) 659
8 APEX Collaboration Search for a New Gauge Boson in Electron-Nucleus Fixed-Target Scattering by the APEX Experiment PRL 107 (2011) 191804 1108.2750
9 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
10 SINDRUM I Collaboration Search for weakly interacting neutral bosons produced in $ \pi^{-}p $ interactions at rest and decaying into $ e^{+}e^{-} $ pairs. PRL 68 (1992) 3845
11 LHCb Collaboration Proposed Inclusive Dark Photon Search at LHCb PRL 116 (2016) 251803 1603.08926
12 LHCb Collaboration Search for Dark Photons Produced in 13 $ TeV pp $ Collisions PRL 120 (2018) 061801 1710.02867
13 LHCb Collaboration Search for $ A'\to\mu^+\mu^- $ Decays PRL 124 (2020) 041801 1910.06926
14 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
15 J. A. Evans, A. Gandrakota, S. Knapen, and H. Routray Searching for exotic B meson decays with the CMS L1 track trigger PRD 103 (2021) 015026 2008.06918
16 CHARM Collaboration Search for Axion Like Particle Production in 400-$ GeV\ $ Proton - Copper Interactions PLB 157 (1985) 458
17 LHCb Collaboration Search for hidden-sector bosons in $ B^0 \!\to K^{*0}\mu^+\mu^- $ decays PRL 115 (2015) 161802 1508.04094
18 LHCb Collaboration Search for long-lived scalar particles in $ B^+ \to K^+ \chi (\mu^+\mu^-) $ decays PRD 95 (2017) 071101 1612.07818
19 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004 CMS-00-001
20 CMS Tracker Group Collaboration The CMS Phase-1 Pixel Detector Upgrade JINST 16 (2021) P02027 2012.14304
21 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
22 CMS Collaboration Search for narrow resonances in dijet final states at $ \sqrt(s)= $ 8 TeV with the novel CMS technique of data scouting PRL 117 (2016) 031802 CMS-EXO-14-005
1604.08907
23 M. Cacciari, G. P. Salam, and G. Soyez The anti-$ {k_{\mathrm{T}}} $ jet clustering algorithm JHEP 04 (2008) 063 0802.1189
24 M. Cacciari, G. P. Salam, and G. Soyez FastJet user manual EPJC 72 (2012) 1896 1111.6097
25 P. Nason A new method for combining NLO QCD with shower Monte Carlo algorithms JHEP 11 (2004) 040 hep-ph/0409146
26 S. Frixione, P. Nason, and C. Oleari Matching NLO QCD computations with parton shower simulations: the POWHEG method JHEP 11 (2007) 070 0709.2092
27 S. Alioli, P. Nason, C. Oleari, and E. Re NLO single-top production matched with shower in POWHEG: $ s $- and $ t $-channel contributions JHEP 09 (2009) 111 0907.4076
28 Y. Gao et al. Spin Determination of Single-Produced Resonances at Hadron Colliders PRD 81 (2010) 075022 1001.3396
29 S. Bolognesi et al. On the spin and parity of a single-produced resonance at the LHC PRD 86 (2012) 095031 1208.4018
30 I. Anderson et al. Constraining Anomalous HVV Interactions at Proton and Lepton Colliders PRD 89 (2014) 035007 1309.4819
31 A. V. Gritsan, R. Rontsch, M. Schulze, and M. Xiao Constraining anomalous Higgs boson couplings to the heavy flavor fermions using matrix element techniques PRD 94 (2016) 055023 1606.03107
32 A. V. Gritsan et al. New features in the JHU generator framework: constraining Higgs boson properties from on-shell and off-shell production PRD 102 (2020), no. 5, 056022 2002.09888
33 T. Sjostrand et al. An Introduction to PYTHIA 8.2 CPC 191 (2015) 159 1410.3012
34 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
35 M. Cacciari, M. Greco, and P. Nason The P(T) spectrum in heavy flavor hadroproduction JHEP 05 (1998) 007 hep-ph/9803400
36 M. Cacciari, S. Frixione, and P. Nason The p(T) spectrum in heavy flavor photoproduction JHEP 03 (2001) 006 hep-ph/0102134
37 M. Cacciari et al. Theoretical predictions for charm and bottom production at the LHC JHEP 10 (2012) 137 1205.6344
38 M. Cacciari, M. L. Mangano, and P. Nason Gluon PDF constraints from the ratio of forward heavy-quark production at the LHC at $ \sqrt{S}= $ 7 and 13 TeV EPJC 75 (2015) 610 1507.06197
39 CMS Collaboration Extraction and validation of a new set of CMS PYTHIA8 tunes from underlying-event measurements EPJC 80 (2020) 4 CMS-GEN-17-001
1903.12179
40 NNPDF Collaboration Parton distributions from high-precision collider data EPJC 77 (2017) 663 1706.00428
41 GEANT4 Collaboration GEANT4--a simulation toolkit NIMA 506 (2003) 250
42 J. E. Gaiser Charmonium Spectroscopy From Radiative Decays of the $ J/\psi $ and $ \psi^\prime $ PhD thesis
43 R. A. Fisher On the Interpretation of $ \chi^{2} $ from Contingency Tables, and the Calculation of P J. R. Stat. Soc 85 (1922) 87
44 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
45 T. Junk Confidence level computation for combining searches with small statistics NIMA 434 (1999) 435 hep-ex/9902006
46 A. L. Read Presentation of search results: The $ CL_{\mathrm{s}} $ technique JPG 28 (2002) 2693
47 G. Cowan, K. Cranmer, E. Gross, and O. Vitells Asymptotic formulae for likelihood-based tests of new physics EPJC 71 (2011) 1554 1007.1727
48 ATLAS and CMS Collaborations Procedure for the LHC Higgs boson search combination in summer 2011 ATL-PHYS-PUB-2011-011, CMS NOTE-2011/005
49 CMS Collaboration CMS luminosity measurement for the 2017 data-taking period at $ \sqrt{s} = $ 13 TeV
50 CMS Collaboration CMS luminosity measurement for the 2018 data-taking period at $ \sqrt{s} = $ 13 TeV
Compact Muon Solenoid
LHC, CERN