Seeking for sterile neutrinos with displaced leptons at the LHC

# Seeking for sterile neutrinos with displaced leptons at the LHC

Jia Liu Enrico Fermi Institute, University of Chicago, Chicago, IL 60637, USA    Zhen Liu Maryland Center for Fundamental Physics, Department of Physics, University of Maryland, College Park, MD 20742, USA    Lian-Tao Wang Enrico Fermi Institute, University of Chicago, Chicago, IL 60637, USA Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA    Xiao-Ping Wang High Energy Physics Division, Argonne National Laboratory, Argonne, IL 60439, USA
July 2, 2019
###### Abstract

We study the signal of long-lived sterile neutrino at the LHC produced through the decay of the boson. It decays into charged lepton and jets. The characteristic signature is a hard prompt lepton and a lepton from the displaced decay of the sterile neutrino, which leads to a bundle of displaced tracks with large transverse impact parameter. Different from other studies, we neither reconstruct the displaced vertex nor place requirement on its invariant mass to maintain sensitivity for low sterile neutrino masses. Instead, we focus on the displaced track from the lepton. A difficulty for low mass sterile neutrino study is that the displaced lepton is usually non-isolated. Therefore, leptons from heavy flavor quark is the major source of background. We closely follow a search for displaced electron plus muon search at CMS and study their control regions, which is related to our signal regions, in great detail to develop a robust estimation of the background for our signals. After further optimization on the signal limiting the number of jets, low and large lepton displacement to suppress SM background, we reach an exclusion sensitivity of about () for the mixing angle square at 10 (2) GeV sterile neutrino mass respectively. The strategy we propose can cover the light sterile masses complimentary to beam dump and forward detector experiments.

preprint: EFI-19-6

## I Introduction

The origin of the neutrino masses is a puzzle of the Standard Model Fukuda:1998mi (); Fukuda:1998ah (); Ahmad:2002jz (); Eguchi:2002dm (). The parameters in the neutrino sector, such as the mass differences and mixing angles, have been measured with higher and higher precision Tanabashi:2018oca (). At the same time, cosmological observations set an upper bound on the sum of neutrino masses to be smaller than eV and the effective extra relativistic degrees of freedom to be from Planck 2018 data Aghanim:2018eyx (). Many mechanisms have been introduced to incorporate the presence of the tiny neutrino masses. Among them, the seesaw type of solutions to the small neutrino mass is probably the most plausible. A heavier, often sterile, neutrino to realize the seesaw mechanism of a certain type would be a smoking gun signal for this class of models.

Many active experimental programs have been developed to search for sterile neutrinos, through oscillation via light (eV-keV) ones at current or future short and long-baseline neutrino facilities Kopp:2013vaa (); Giunti:2015wnd (); Capozzi:2016vac (); Gariazzo:2017fdh (); Dentler:2017tkw (); Dentler:2018sju (), neutrinoless double-beta decay (MeV) for intermediate ones Benes:2005hn (); Atre:2009rg (); Blennow:2010th (); Rodejohann:2011mu (); Barea:2015zfa () and at LHC through same-sign dileptons for heavy ones (100 GeV or above) Han:2006ip (); delAguila:2007qnc (); Atre:2009rg (). Furthermore, in the mass range of GeV to 100 GeV, the sterile neutrino would be metastable or long-lived at the detector scale, and can be probed at beam-dump types of experiments Badier:1985wg (); Bergsma:1985is (); CooperSarkar:1985nh (); Bernardi:1987ek (); Baranov:1992vq (); Vilain:1994vg (); Vaitaitis:1999wq (); Astier:2001ck (); Orloff:2002de (); Adams:2013qkq (); Anelli:2015pba (); Drewes:2018gkc () (see also a review Deppisch:2015qwa ()), Very recently, the searches for sterile neutrino at the LHC started to develop actively as part of the long-lived particles searches covering the GeV scale Ahmad:2002jz (); Graesser:2007pc (); Graesser:2007yj (); Helo:2013esa (); Maiezza:2015lza (); Batell:2016zod (); Antusch:2016vyf (); Nemevsek:2016enw (); Antusch:2017hhu (); Caputo:2017pit (); Cottin:2018kmq (); Helo:2018qej (); Kling:2018wct (); Curtin:2018mvb (); Dib:2019ztn (); Drewes:2019fou (); Bondarenko:2019tss ().

In this work, we focus on sterile neutrino production from boson decay, which leads to a signature of a prompt lepton and a displaced but non-isolated lepton. A central challenge for new phenomenological studies on long-lived particles at the LHC is how to estimate the corresponding background. We overcome, at least partially, this challenge by extracting background behavior information from two very similar control regions measured and validated by an experimental search targeting different signatures at CMS CMS-PAS-EXO-16-022 (). We show the LHC sensitivity would be improved significantly in the regime of sterile neutrino mass around 1-20 GeV with a mixing angle squared between to .

The paper is organized as follows. In Sec.II, we provide a brief review of the seesaw models. In Sec.III, we discuss the properties of the sterile neutrino relevant for this study. In Sec. IV, we present the proposed analysis with the corresponding background estimation and the resulting for model parameter coverage. Finally, we conclude and discuss future directions in Sec. V.

## Ii Sterile neutrino models

In this section, we briefly review a few classes of seesaw models to motivate the parameter space we focus on. We begin with the original seesaw models Minkowski:1977sc (); GellMann:1980vs (); Mohapatra:1979ia (); Yanagida:1980xy (); Schechter:1980gr () with the interaction Lagrangian of the new sterile neutrino sector

 ΔLν=−λν¯L~HN−mN2¯NcN+h.c., (1)

where . The mass matrix in the flavor basis is

 Mν=(0mDmDmN), (2)

where with the vacuum expectation value of the Higgs field GeV. The mass of light and heavy neutrino are

 mν≡m1≃m2DmN,m2≃mN+m2DmN≃mN (3)

respectively, in the heavy Majorana mass limit. The mixing angle is , which yields a relation between the mixing angle, light (active) neutrino mass and heavy neutrino mass:

 sin2θ≃mνmN. (4)

As we shall see later, the mixing angle controls the phenomenology of the searches for the sterile neutrino . Beyond the original (and the most basic) seesaw model, a number of extensions and variations have been proposed, with new parameters and deviations beyond the simple relation shown in Eq. 4. We mention two such examples, namely, inverse seesaw model Mohapatra:1986aw (); Mohapatra:1986bd () or linear seesaw model Wyler:1982dd (); Akhmedov:1995ip (); Akhmedov:1995vm (). These variations of seesaw models can be parameterized similarly, and they are more accessible (and testable) with a larger mixing angle while capable of reproducing the light neutrino mass. The phenomenology is controlled by the same set of phenomenologically relevant parameters and but with a different relation to Lagrangian parameters. In inverse seesaw, we have

 sin2θ=mνμ(inverse seesaw), (5)

where is a free parameter with the relation . Here becomes a Dirac mass term of new sterile neutrinos and is the Majorana mass of one of the new species. For linear seesaw, the relation is

 sinθ=mνmψ(linear seesaw), (6)

where is an additional Dirac mass term small than , and its presence will violate the lepton number. With the additional parameters introduced, we can attribute neutrino mass to these new parameters. Thereafter, we only focus on and which are both free parameters.

In the later sections, we will constrain the mixing angle as a function of sterile neutrino mass with experimental studies. To simplify this discussion and focus on demonstrating our strategy, we will make the simple assumption that all three generations of the active neutrinos mix with equally. In the context of more general seesaw models, making this assumption does not commit us to a particular neutrino mass hierarchy. The collider signal studied in this paper mainly rely on final states with electrons and muons, while the mixing with neutrino enters the decay width of the sterile neutrino if . If we relax the assumption of universal mixing, the reach can be obtained from our result by rescaling relevant rates.

In the basic seesaw model, the mixing angle is not a free parameter for a particular sterile neutrino mass. It is fixed to be around

 sin2θ≃10−12(mν0.01 eV)(10 GeVmN), (7)

which is very difficult to be probed at the colliders. Still, in the context of more extended models, larger mixing angle is allowed. It is these non-minimal models which would be the main target for LHC searches.

## Iii Sterile neutrino properties

In this section, we discuss the features and properties of sterile neutrino from the extended seesaw models, with mixing angle and mass as free parameters. Due to the mixing between SM neutrino and sterile neutrino , the relevant interactions are

 L=gsinθ√2(Wμ¯ℓLγμN+h.c.)−gcosθsinθ2cosθwZμ(¯νLγμN+¯NγμνL)+gsin2θ2cosθwZμ¯NγμPLN, (8)

where is the weak mixing angle.

In a large region of parameter spaces of interests, the right-handed neutrino has a macroscopic lifetime

 cτ≃12 km×(10−12sin2θ)(10 GeVmN)5, (9)

with the details of partial decay widths following Helo:2010cw (); Helo:2013esa (), see also Bondarenko:2018ptm (). The reference mixing angle squared is chosen for SM neutrino mass to  eV in the basic seesaw model. For an intermediate value of the mixing angle squared , a 10 GeV sterile neutrino will still be long-lived at collider level with a proper lifetime of 12 mm. In Fig. 1, we show the partial widths, branching ratios and lifetime for heavy sterile neutrino decay, with the assumption mixing angles between the three SM flavor neutrino are democratic .

At colliders, the production of the sterile neutrino is the same as SM neutrinos, with additional suppression from mixing angle and kinematics. The production through the on-shell boson is of particular interest here, given that the associated charge leptons can be used to trigger the signal processes. The long-lived sterile neutrino can then be analyzed without trigger penalties. The expected total number of sterile neutrino at the HL-LHC for one generation of democratic sterile neutrino is approximately

 L×σ(pp→W±)Br(W±→ℓ±N)≃1.8×105(sin2θ10−6), (10)

where is the integrated luminosity at the HL-LHC.

It is worth noting that the existence of two other SM venues for heavy sterile neutrino production, namely the boson decay and the Higgs boson decay. These two channels will become competitive if displaced-track triggers becoming available Gershtein:2017tsv (), which can provide a similar amount of sterile neutrino in addition to the boson decay. Furthermore, the Higgs channel has two exciting features, namely Higgs specific trigger options and large branching fractions. There are many sub-leading Higgs production channels which can be triggered on, especially the vector boson fusion channel and weak boson associated production channel. Although the Higgs boson production rate is three to four orders of magnitude smaller than the and boson, the branching fraction of the Higgs boson to can be five orders of magnitude larger, given its small total width. In this study, we do not include these channels and save these interesting new production modes for future studies in association with displaced triggers.

On the top panel of Fig. 2, we show the transverse momentum and pseudo-rapidity distribution for sterile neutrino from production . The process is generated by MadGraph5_aMCNLO Alwall:2014hca (), and the parton shower is performed by Pythia8 Sjostrand:2006za (); Sjostrand:2007gs (). The distribution of is symmetric and dominantly within the range between . Moreover, its transverse momentum peaks around 3040 GeV, which is dictated by the maximum momentum it can obtain from decay in the center of mass frame, . For the events with larger than this value, the initial state radiation provides additional transverse momentum of the boson system. In the bottom panel of Fig. 2, we show the distribution for both the prompt lepton from decay and displaced lepton from decay. The distribution of prompt lepton is similar to that of , with the difference that the maximum momentum in center of mass frame becomes . The displaced lepton transverse momentum, , is dominantly distributed within GeV, because the other two particles in the decay take away approximately two thirds of the available energy.

## Iv New searches strategy for sterile neutrinos

There are many challenges to overcome to reach good sensitivities in the search for sterile neutrinos at the LHC. A major challenge is how to achieve a good signal selection efficiency with effective background suppression. Typical neutrino mass models point to tiny couplings between the sterile neutrino and the SM electroweak gauge bosons, leading to displaced decays of the sterile neutrinos. Therefore, a very effective strategy is necessary to pick out these displaced events at the LHC. Another aspect is triggering. The sterile neutrino can be produced in electroweak processes, which typically give rise to soft objects. At the same time, due to the low signal rates, the triggering needs to be as efficient as possible.

We consider the process of , with a subsequent displaced decay of the sterile neutrino . Since the prompt lepton is hard, one can in principal trigger the event using single lepton trigger. According to CMS Phase-2 upgrade of the level-1 trigger Collaboration:2283192 (), the trigger thresholds (with track trigger) on a single lepton is GeV for isolated (non-isolated) , and GeV for . As shown in the prompt lepton distribution in Fig. 2, the signal events will be reduced substantially by the single lepton trigger. Therefore, we consider double lepton trigger with the leading lepton GeV and the sub-leading lepton with GeV, also benefits from the track trigger Collaboration:2283192 (). In reality, the combination of triggers will be used in the experiment. Therefore, it is expected that the trigger efficiency will be better than using a single trigger category. Therefore, our approach on the triggering in this study is conservative.

For the off-line signal selection, the major advantage is the presence of a bundle of tracks with large (transverse) impact parameter, . We denote the decay position of sterile neutrino as with the production point being the origin, and the momentum of a daughter particle from decay in x-y plane is denoted as . We can define the transverse impact parameter as111For simplicity, we ignore the bending of tracks in this study as the displaced leptons is sufficiently large. Furthermore, the hard prompt lepton enable us to reconstruct the three-dimensional position of the primary vertex, additional requirement on the three-dimensional impact parameter can compensate this simplification and improve our results.

 d0= ⎷x2+y2−(xpx+ypy)2p2x+p2y. (11)

We require the sterile neutrino to decay within the detector region m and m, where is the radial distance to the beam line. These numbers are chosen to guarantee that the lepton from the decay of the sterile neutrino would pass at least four Outer Tracker layers in either central or forward region, based on the Phase-2 upgrade of the CMS tracker Collaboration:2272264 (). Moreover, each layer of Outer Tracker consists of two closely spaced silicon sensors, which are called modules to reject the low tracks. The threshold is about 2 GeV which is easy to satisfy for our signal, based on the distribution of displaced lepton in Fig. 2 and in particular for the events passing our trigger selection. The silicon sensors have a good granularity to provide sufficient spatial resolution and the module has a resolution of . If the electromagnetism calorimetry also has an excellent pointing resolution, using one layer of the Outer Tracker could be good enough to identify the direction of the track and measure the , the detector region can be enlarged to m and m. In this case, the volume of the detector is increased by a factor of 8, leading to a significant improvement of the reach. We leave this possibility to future experimental studies and keep our conservative requirement on the number of track layers to hit.

The signal is selected first by imposing the dilepton trigger with and requiring the long-lived sterile neutrino decays in the detector region m and m. We further require the displaced lepton in the event has a large transverse impact parameter . The survival possibility is the number of such events divided by the total signal events, which we denoted as . We plot its distributions with different sterile neutrino mass and proper lifetime in Fig. 3. For a given proper lifetime, a smaller sterile neutrino mass leads to a lower . This is a result of the kinematics of the resonance decay. Smaller sterile neutrino mass has larger boost factor, reducing the decay probability of sterile neutrino inside the detector region and reducing the displaced lepton opening angles to give a smaller . For m, a significant portion of the sterile neutrino produced decay within the region of the detector we focused on. A notable exception is when the sterile neutrino is light. For example, the is significantly smaller for GeV than other masses since it would often decay outside the region due to the large boost. For m and m, the survival probabilities for different s are roughly proportional to due to the boost factor of . For the same , the probabilities of m and m differ by a factor of 100, exactly proportional to as expected in the long lifetime limit. In summary, the survival probabilities in Fig. 3 are mostly determined by the probability for decay inside the required region, which is proportional to distance parameter . Here, is the expected decay distance in the laboratory frame, where is the boost factor and is the velocity of .

The signal contains one prompt lepton and another displaced lepton. The estimation of the Standard Model (SM) background of a newer search for a long-lived particle is always a challenge at the LHC. Fortunately, many important features of the corresponding background have been effectively explored as the control regions in a search for displaced electron plus muon search at CMS CMS-PAS-EXO-16-022 (). This search was performed with of 13 TeV LHC data, and it looks for a pair of displaced isolated leptons with different flavors with minimum of 38 GeV each, targeting signals from pair produced top squarks decaying leptonically. The control regions (CR) of the background for this search helps to identify the dominant background of our new signals for long-lived sterile neutrinos. Specifically, the CR-III and CR-IV of this analysis requires one prompt lepton with transverse impact parameter smaller than 200 m and a displaced lepton with transverse impact parameter greater than 100 m. The CR-III looks for a heavy flavor jet plus a displaced electron and CR-IV looks for a heavy flavor jet plus a displaced muon. CR-III is triggered by a singlet electron trigger (with electron  GeV), while CR-IV is triggered by a muon plus jet trigger (with a muon  GeV and jet GeV), both with low threshold for the leptons Khachatryan:2016bia ().

The detailed studies of these control regions in Ref. CMS-PAS-EXO-16-022 () confirm the following crucial fact about the displaced background distribution. First of all, the displaced leptons are dominantly from heavy flavor jets (b-jets) by using the “tag and probe” method where the jet recoil against the displaced lepton is tagged as a b-jet. The subleading background is from which is smaller by more than one order of magnitude, agreeing well with our simulation. Furthermore, for the displaced lepton, the (normalized) differential distribution as a function of the transverse impact parameter is shown to be the same for isolated and non-isolated leptons222CR-III and CR-IV are dominated by heavy flavor plus lepton background. When describing the properties, we use them interchangeably..

In Fig. 2 of Ref. CMS-PAS-EXO-16-022 (), the background events are dominated by heavy flavor. Summing up all the bins, the corresponding cross-sections are 16.6 pb for HF (CR-III) and 259 pb for HF (CR-IV) respectively. We denote them as and respectively. The HF background is much larger than HF background because the muon plus jet trigger requires a softer lepton than the single electron trigger. We used these results to validate our own simulation. In particular, we generated events and require or to show up in the event after hadron fragmentation using Pythia8. We apply the CR-III and CR-IV triggers separately and further require at least one b-tagged jet. For HF, the isolation requirement for the electron is and the additional sum within the isolation cone should be less than () of the electron’s in the barrel (endcap) region. For HF, the isolation requirement for the muon is and the additional sum within the isolation cone should be less than of the muon . We implement these requirements by modifying the Delphes3 deFavereau:2013fsa (). We found the corresponding cross-sections to be 21.9 pb and 241.2 pb for HF and HF backgrounds, which are consistent with our extraction from Ref. CMS-PAS-EXO-16-022 (). We denote them as and which are the cross-sections after cuts with isolation requirement on leptons, denoted as “icut” in the superscript. Such consistency provides confidence when we calculate the background to our new search333We thank Bingxuan Liu on the analysis of Ref. CMS-PAS-EXO-16-022 () for clarifying the details of the control region in the analysis and confirming our scaling method..

The displaced lepton in the signal is often non-isolated. Therefore, a major background comes from events with a displaced non-isolated lepton from heavy flavor jets and a prompt lepton. Consequently, the leading backgrounds for our signal could be with decaying leptonically, and with one of the top quarks decay leptonically. The number of the background can be calculated as

 Nbkg=σCMSHF+e+σCMSHF+μσicutb¯b(e)+σicutb¯b(μ)(σncutW+b¯b,W→ℓν×ϵW+b¯bopt+σncutt¯t→b¯b+ℓ+X×ϵt¯topt)×LHL−LHC, (12)

where the “ncut” in the upper script means requiring jets with GeV while having one non-isolated lepton in the final states. In the event generation, we require b quark GeV at parton level to ensure an energetic non-isolated lepton from its hadronic fragmentation. The lepton represent both and . Further optimization selections are applied to both the signal and background, whose selection efficiencies are denoted as in Eq. 12. The cross-sections and are found to be 136 pb and 3.8 pb respectively, after applying the “ncut”.

To further reduce the SM background, especially those from the process, additional cuts on the hadronic activities can help. In Fig. 4, we show the normalized distribution for the number of b-jets with GeV denoted as , the number of jets with GeV denoted as , the number of jets with GeV denoted as and , where is the scalar sum of for all visible objects, including hadronic jets and leptons. We choose the optimization condition , , and GeV to suppress the SM background without a significant reduction of the signal events. In Fig. 4, the above cuts have been applied after its distribution have been shown in sub-figures (a), (b), (c), and (d). We also include the efficiency information for the dilepton trigger requirement GeV and GeV. The effects of various selection cuts and the total efficiencies are given in Table. 1.

After combining all the cuts and optimizations, the backgrounds from and are 0.017 pb and 0.022 pb respectively, which are comparable in size. Multiplying by the integrated luminosity from HL-LHC, the corresponding number of SM backgrounds are about 51000 and 65100 respectively. We assume for the SM background, the distribution for the lepton does not depend strongly on the of the lepton and other optimization cuts on jets. Together with the fact that isolated and non-isolated lepton has the same distribution, we can use the normalized differential distribution from Ref. CMS-PAS-EXO-16-022 () to determine the background.

In Fig. 5, we plot the total number of SM background after applying with all the cuts and optimizations applied 444Since the differential distribution of from CMS study CMS-PAS-EXO-16-022 () stops at 0.5 cm, we linearly extrapolate their data beyond that for the case of a larger cut of 2 cm.. The SM background decreases linearly with cut. After comparing with the accumulative distribution for signal in Fig. 3, we adopt two benchmarks for the displacement, requiring 0.2 cm and 2 cm, respectively. From Fig. 5, the total number of background event is around 2000 for cm and 100 for cm. With these, we derive the C.L. sensitivity for sterile neutrino in red shaded region in Fig. 6, where the left and right panel are for  cm and  cm, respectively. The upper edge of the sensitivity region corresponds to the shorter lifetime. It is driven by the lower cut on the transverse impact parameter and we see clear advantage for smaller cut on , even though the number of SM background is larger. On the other hand, the lower edge of the red shaded region corresponds to the longer lifetime, which is insensitive to the cut.

In Fig. 6, we also show the projected sensitivity from a proposed displaced vertex search at the LHC in dot-dashed curves Cottin:2018nms () and various fix target or satellite detector experiments Adams:2013qkq (); Alekhin:2015byh (); Anelli:2015pba (); Curtin:2018mvb (); Kling:2018wct () in dashed curves. It is clear that our more inclusive search allows us to explore parameter spaces of the sterile neutrino with lighter masses. Furthermore, when compared with other fixed target experiment and satellite detector experiments, our proposal covers new regions of sterile neutrinos whose lifetime is too short to have sufficient flux for these experiments. In the long run, if any discoveries were made, the regions overlapping with different searches and experiments will provide valuable information for the underlying new physics model.

## V Conclusion and outlook

In this paper, we study the signal of long-lived sterile neutrino from -boson decay, , with the subsequent decay of . The characteristic feature is a hard prompt lepton with a displaced lepton with large transverse impact parameter . We neither reconstruct the displaced vertex nor cut on its invariant mass, therefore it can be sensitive for very low sterile neutrino mass. However, there is a crucial subtlety that with such small masses, the displaced lepton is usually non-isolated from the other two jets in the same decay. To estimate the background, we have used the information from a search for displaced electron plus muon search at CMS CMS-PAS-EXO-16-022 () which studied relevant background in its control regions. It shows that for non-isolated lepton, those from heavy flavor quarks are the dominant SM background. Moreover, it demonstrates the important fact that the normalized differential distribution has the same shape for isolated and non-isolated leptons. Therefore, we can use their distribution for the non-isolated lepton from heavy flavor quark background. We recast their control region selection and found a good agreement with their observations. This ensures our background estimations are robust. After proposing the optimization cuts for the signal, we obtain the result for the sensitivity on parameters for the sterile neutrino. Our results can explore the region with light and relative large mixing angle, which is not covered by the displaced vertex searches at LHC, beam dump and far detectors experiments.

There are many new exciting opportunities at the LHC, such as displaced tracker triggers proposed to be implemented at low-level at CMS by Ref. CMS-PAS-FTR-18-018 (). This trigger can be handy as a general trigger for long-lived particles in a broad class of theories Alimena:2019zri (), especially for physics signals hard to be triggered on using traditional trigger. Although we have applied two lepton trigger in this study, adding the displaced tracker trigger into the trigger menu can effectively lower the requirement on the leptons. Moreover, the Phase-2 upgrade of the tracker system of CMS contains the Outer Tracker layer which is consisted of two closely spaced silicon sensors. It effectively doubles the hits of the track; therefore one can lower the requirement on the number of penetrating layers, which can increase the volume of the long-lived particle decay. It can significantly increase the sensitivity but needs further studies.

Note added. While we were finalizing this paper, Ref. Dib:2019ztn (), Drewes:2019fou () and Bondarenko:2019tss () appeared on the arXiv studying similar topics. The main focus of Ref. Dib:2019ztn () is looking for sterile neutrino decay to lepton plus pions, using displaced vertex associated with low muon. And their mass coverage for sterile neutrino is from 5 GeV to 20 GeV. The second paper Drewes:2019fou () studies long-lived sterile neutrino in displaced vertex searches and also using muon chambers to detect the muons from such vertex. It uses invariant mass cut and displacement cut to suppress the SM background. The third paper Bondarenko:2019tss () requires displaced vertex reconstructed by two muons in the CMS muon detector.

Acknowledgements: We thank Felix Kling, Jared Evans and Bingxuan Liu for helpful discussion. JL acknowledges support by an Oehme Fellowship. ZL is supported in part by the NSF under Grant No. PHY1620074 and by the Maryland Center for Fundamental Physics. XPW is supported by the U.S. Department of Energy under Contract No. DE-AC02-06CH11357. LTW is supported by the DOE grant DE-SC0013642. ZL thanks the Aspen Center for Physics, which is supported by National Science Foundation grant PHY-1607611.

## References

You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters