Time-Delayed Electrons from Higgs Decays to Right-Handed Neutrinos

We consider the physics reach at the high luminosity LHC (HL-LHC) using the timing capability of a minimum ionizing particle (MIP) timing layer and the electromagnetic calorimeter in the CMS experiment for a simplified right-handed neutrino effective theory. In the simplified model, the lightest right-handed neutrinos are produced through higgs decay with mean decay lengths of order $\mathcal{O}(0.1~{\rm m}-100~{\rm m}$) and then subsequently decay into final states that include electrons. We consider the effect of several different kinematic cuts for a timing based analysis of this model. We demonstrate that if timing information can be successfully incorporated into the event trigger, it is possible to probe the higgs branching ratio to two right-handed neutrinos to the order $\mathcal{O}(10^{-5}) $ for a range of different right-handed neutrino lifetimes.


Introduction and overview
The origin of neutrino mass remains an un-determined aspect of the Standard Model. It is well-known that a type I seesaw mechanism [1,2] can efficiently generate the hierarchically small neutrino masses (m ν ≤ 0.1 eV [3]) in a theory that is renormalizable up to the electroweak scale and beyond (see [4] for a review). If the the right-handed neutrinos have an electroweak scale majorana mass, they can be kinematically acessible at colliders [5][6][7][8][9][10]. It is natural in such models for the right-handed neutrinos to have mean decay lengths of order O(0.1 m − 100 m) due to small coupling to standard model particles. Since the right-handed neutrinos in such models typically decay to final states with leptons, recent investigations of such scenarios have focussed on the displaced-vertex/track nature of the final state leptons [11][12][13][14][15][16][17][18][19][20][21] and have proven to be an efficient probe of this class of models.
In general, if long-lived particles produce a displaced-vertex/track type signature at colliders, they can also be associated with calorimeter energy deposition that is "timedelayed" with respect to the energy deposition of massless states produced at the primary collision vertex. Such a time-based signature is not new but a classic signature of lowscale gauge mediation (GMSB) [22]. The proposal to include a minimum ionizing particle (MIP) detector on the inner barrel and endcap regions of the electromagnetic calorimeter (Ecal) as well as planned improvements to the timing resolution of the Ecal itself at the CMS experiment promises to increase time resolution for charged particles to 30 ps at the HL-LHC [23]. The utility of this MIP/Ecal timing capability has been recently analyzed JHEP07(2019)089 in [24] to investigate the reach of the timing capability to GMSB and Higgs decays to long-lived particles where peak sensitivity is achieved at cτ = 1 m. Since a timing based analysis need not depend crucially on tracking, the time-delayed signature can arise from a signal event so long as it occurs in the main inner detector volume. Contrast this with the displaced-vertex/track searches which typically reach peak sensitivity at cτ = 1−10 cm for final state electrons. The distinction is due to the deterioration of tracking capability as electron final states produced from long-lived particles that decay farther than a distance of 50 cm transverse to the beamline [25]. Additionally, the volume of the region of efficient tracking is roughly only 10% of the volume enclosed by the inner detector. This opens up the possibility that a timing based analysis may provide a complimentary probe of models having displaced-vertex/tracking signatures, in particular targeting models with longer lifetimes where a majority of the delayed decays occur outside of the region where tracking is most efficient.
The purpose of this paper is to analyze the physics reach of a timing based analysis of right-handed neutrino models where the Higgs boson can decay to "long-lived" righthanded neutrinos which subsequently decay to final states with electrons. We shall work within the context of the type I seesaw effective theory described in [11,12] using the timing analysis techniques of [24]. This model is a well-motivated member of a much larger class of models recently reviewed in [26], in which the Higgs decays to long-lived particles. In order to fully utilize MIP/Ecal timing capability we will not rely on track reconstruction of the electrons in the signal. It will be critical for our purposes that the MIP/Ecal timing capability be incorporated into the event trigger at the CMS experiment and we shall assume for the remainder of the paper that this will be realized at the HL-LHC. Since the performance of such a timing trigger is not finalized, we will need to make reasonable assumptions regarding the efficiency of such a trigger, and we shall clearly state them below. The paper is organized as follows. In section 2 we describe the electroweak scale seesaw effective theory and then in section 3 focus on a particular region of parameter space where only one right-handed neutrino is phenomenologically relevant, parameterizing a simplified model in terms of three parameters. In section 4, we describe the general kinematic features of new Higgs decay channels that are present in this model. Based on these features, in section 5 we outline search criteria, describe our event simulations, estimate backgrounds, and present results in terms of exclusions on Higgs branching ratios to right-handed neutrinos as a function of lifetime for different masses. We conclude in section 6.

Electroweak scale type I seesaw effective theory
The field content of the model is the standard model plus three majorana right-handed neutrinos: N I , where I = 1, 2, 3 labels the three right-handed neutrino flavors. In addition to the particle content, the model also has the following renormalizable interactions:

JHEP07(2019)089
In this notationH = iσ 2 H where H transforms under SU(2) L × U(1) Y as a (2, +1), gauge indices are not shown and i = 1, 2, 3 (I, J = 1, 2, 3) are light (heavy) neutrino flavor indices and N I R,L = P R,L N I . It is well known that at energies below M IJ , these interactions generate a dimension-five Weinberg operator and yield Majorana neutrino masses, m ν , given by where v = 174 GeV is the vacuum expectation value of the Higgs field and M N (y) is the scale of the majorana mass matrix M IJ (y iJ ). If M N is of order the electroweak-scale, neutrino masses of order O(0.1 eV) will be generated for yukawa couplings y ∼ 5 × 10 −7 . Due to the mass mixing induced in the full neutrino sector of this theory, the neutral and charged current interactions include W and Z interactions with the right handed-neutrinos.
Here V andṼ are 3 × 3 matrices that depend on the M IJ , (y) iJ and rotation matrices that diagonalize the neutrino and lepton sector. The yukawa interaction of (2.1) gives rise to Higgs mediated decays of the right-handed neutrinos: We briefly review the derivation of these terms and the precise definitions of V andṼ in appendix A. Parametrically, the couplings V andṼ are of order O( yv M ). We shall focus on a right-handed neutrino mass, 20 GeV < M N 1 < 60 GeV, in which case the dimensionless couplings in V andṼ are O(10 −6 ). The phenomenological effect of such a small coupling is to give the right-handed neutrino a relatively large lifetime of order O(0.1 − 100 ns). We graph the parametric relation between the scale of the left handed neutrino mass, the right handed neutrino mass and the mean decay length of the right-handed neutrino in figure 1.
Interestingly, there is considerable overlap of the parameter space for the right-handed neutrino mean decay length and the size of the CMS inner detector, where the shortest distance from the interaction point (IP) to the inner radius of the electromagnetic calorimeter (Ecal) is 1.29 meters. Note also that significant portions of parameter space for this model have cτ > 1m. This opens up the opportunity to probe the nature of the right-handed neutrinos using precision timing. As the neutrino model presented above presently stands, all couplings are too small to give any appreciable cross-section for production of right-handed neutrinos. Of course in order to probe the nature of the right-handed neutrinos, there must be a mechanism for them to be produced. We now turn to motivate and parameterize such a mechanism.
The model so far includes renormalizable interactions. However, if we view the above theory as an effective theory, we should include higher dimensional operators suppressed by a mass scale Λ where this effective theory breaks down. If Λ is not too much greater than the electroweak scale, then these higher dimension operators can affect phenomenology. Presumably this would be an adequate description of some neutrino models that extend the most basic type I seesaw with other particles. However, matching existing models to this effective theory is beyond the scope of this paper. Ignoring flavor indices, there are only two linearly independent gauge invariant dimension-five operators involving the right-handed neutrino, N I R [11] where B µν is the U(1) Y field strength. In this paper we shall restrict our attention to a regime of parameter space where the effects of O IJ B are negligible. Specifically, we shall work with only one phenomenologically accessible lightest right-handed majorana neutrino where N c 1R σ µν N 1R = 0 due to the majorana nature of the right-handed neutrino. Turning our attention to O IJ h , one main phenomenological consequence results from including it in our description of the effective theory; when the Higgs field gets a vev, O IJ h can lead to a significant production mode of right-handed neutrinos through Higgs boson decay if the mass of the right-handed neutrino mass eigenstate is less than m h /2. Approximately, If we take Λ = 10 TeV then BR(h → N 1 N 1 ) < 0.035 for c < 0.68. The Higgs vev also leads to a contribution to the right-handed neutrino mass, however that contribution will always be sub-dominant compared to the already present M -term in (2.1) for the parameters studied in this paper.
To summarize, the completely general right-handed neutrino effective theory lagrangian considered in this paper is give by:

JHEP07(2019)089
where L SM is the standard model lagrangian which includes L W/Z and L h when all fields are rotated to their mass eigenbasis and L N from 2.1. The parameter space of this model is quite large, depending on y iJ , M IJ and c IJ . In the next section we will restrict the parameter space of this model with some basic simplifying assumptions in order to more easily facilitate a collider study of the model.

A simplified model
We shall now restrict our attention to a particular region of parameters of the right-handed neutrino effective theory that is of phenomenological interest. The region of interest is chosen as follows: 1. We assume the mass of the lightest right-handed neutrino (Ñ 1 ), MÑ 1 < 62.5 GeV and the other right-handed neutrinos have masses above 62.5 GeV and are thus inaccessible or otherwise produced with such low rates that they are not phenomenologically relevant.
2. We assume that the lightest right-handed neutrino (Ñ 1 ) couples dominantly to e ± or ν e . We will impose this by making V 11 =Ṽ 11 = O(10 −6 ) and all other V iJ = V iJ = 0. Our choice for the value of this coupling comes from the requirement that the neutrino mass is given by the relation: 3. Even with the above constraints, there remains enough freedom in the choice of y iJ and M IJ that the neutrino mass mixings can be reproduced without spoiling the above choices for couplings.
The above restrictions applied to 2.7 result in the following interaction lagrangian between the lightest kinematically accessible right-handed neutrino and standard model leptons, gauge bosons, and the Higgs: where we have now lowered the flavor index on the right-handed neutrino and removed the "L,R" labels as well and the "∼" for neatness. The lagrangian (3.1) depends effectively on only three parameters (V 11 , M N 1 , cv Λ ). However, if the see-saw mechanism is responsible for neutrino masses, then there is a constraint between two of the parameters in terms of the light neutrino masses (m ν ): GeV and m ν ∼ 0.1 eV. Restricting our attention to masses in the range: 20 GeV < M N 1 < 62.5 GeV, this constraint combined with the above interactions will lead to three-body decays of the right-handed neutrinos through off-shell W/Z/h to a number of final states with decay lengths (cτ ) that range between: 0.1 m < cτ < 100 m (see figure 1). Note that due to range of m ν that may be physically realized in nature, a number of different values of V 11 (and thus cτ ) are possible given a choice of mass M N 1 . The third parameter, cv Λ , Figure 2. An example of a Higgs decay process through off-shell W-bosons into final states with electrons.
controls the branching ratio of the Higgs boson to a pair of right-handed neutrinos. For cv Λ in the range 10 −4 < cv Λ < 10 −2 , the branching ratio, BR(h → N 1 N 1 ), ranges between 2.5 × 10 −6 < BR(h → N 1 N 1 ) < 2.5 × 10 −2 . Since this acts as an additional decay mode for the Higgs, the indirect bound on h → invisible [27] constrains the BR(h → N 1 N 1 ) < 0.27 in this model. We shall exchange the parameters that appear in the lagrangian for the corresponding phenomenological parameters: where BR = BR(h → N 1 N 1 ) and describe this simplified model in terms of the phenomenological parameters for the remainder of this work. The only appreciable production mode open to the right-handed neutrinos is through the intermediate production and decay of the Higgs boson. After being produced in pairs by some fraction of Higgs decays, the lightest right-handed neutrinos have a variety of decay modes open to them. For the assumptions listed above, the decays will always be threebody decays through off-shell Z, W or h. A sample decay is shown in figure 2. Roughly 75% of the time there will be at least one e ± and jets or missing energy associated with the final state. See [12,28,29] for a detailed description of final states. interaction point. The efficiency for track reconstruction falls quickly as the true location of the origin of the electron track becomes greater than a distance of R = 50 cm transverse to the beamline [25,32]. Thus a significant fraction of the decay events will produce either poorly re-constructed electrons or events that resemble photons in that the electrons deposit energy in the Ecal but are not associated with tracks after being passed through the multilevel track reconstruction algorthim. These issues are carefully considered in searches for models that produce leptons with tracks having large transverse impact parameters [32,33]. Closely related to the displaced nature of the decays, and most important for our purposes, is the "delayed" time at which electrons from right-handed neutrino decays deposit energy in the Ecal. For instance, after being produced, the right-handed neutrinos will travel some time t N before decaying into electrons and jets. The electrons then travel an additional distance, L, from their production vertex location to an Ecal cell located at some (η, φ), in a time, t e = L/c. Thus the arrival time is given by: t arrival = t N + t e . Had the electrons originated promptly from the primary vertex, they would have a prompt arrival time of t prompt = 1.29m 2c sin(tan −1 (e −η )) for an Ecal cell located at (η, φ) within the central barrel. With these parameters, a time delay: can be associated with each Ecal energy deposit in an event. A distribution of ∆t is shown in figure 4 for signal events of a 50 GeV right-handed neutrino with a 1 meter mean decay length. For heavier mass right-handed neutrinos, which tend to be produced with lower velocities, the time delay is largely due to the longer lifetime of the particle. For lighter mass right-handed neutrinos, which are produced with higher boosts, the time delay is largely due to the longer "kinked" path taken by the right-handed neutrino and subsequently produced electron. The right-handed neutrino decays we are discussing here are very similar JHEP07(2019)089 to the classic signature for low-scale gauge mediated supersymmetry breaking models [22], although the soft kinematics are certainly unique to the delayed neutrino scenario.

Search proposal using MIP timing at CMS
We now propose a strategy to search for events predicted by the simplified model described in section 3: Higgs decays to two "long-lived" right-handed neutrinos which in turn decay to electrons, jets and neutrinos through interactions with charged and neutral currents. Previous work on this model has focused on the displaced nature of the electron and muon tracks [16]. Here we focus on the "delayed" nature of the events and follow very closely the analysis described in [24]. We first define our search criteria, describe our event simulation then we estimate backgrounds and give out results.

Search criteria
We now consider ways that the HL-LHC can capitalize on the unique kinematics of delayed right-handed neutrino decays and be sensitive to these models using precision timing information. A key component of our analysis will be the proposed MIP timing layer [23] and Ecal timing resolution for which the final design, calibration and installation is anticipated at the HL-LHC. The Ecal timing resolution is projected to measure the time of energy deposition at the inner radius of the Ecal with a 30 ps time resolution when p T ≥ 50 GeV. The MIP timing layer will have similar timing resolution for charged particles with much smallet p T . The timing layer's main purpose is to veto events that are out of time with the hard interaction as a means to reduce pile-up. However, the spread of pile-up from one bunch crossing leads to a 190 ps time spread for all of the energy deposition in the Ecal for that crossing. With 25 ns between bunch crossings, this leaves roughly 12.5 ns of time to remain sensitive to Ecal energy deposition before the detector must be prepared for the next bunch crossing. This opens up the possibility for the MIP/Ecal timing information JHEP07(2019)089 to be incorporated into the trigger such that if a event shows significant out of time energy deposition, greater than 1 ns for instance, then such an event can be triggered on rather than vetoed. Furthermore, we shall include objects identified as a photon or as an e ± in our analysis since we do not want counting of signal events to rely solely on objects that are associated with well reconstructed tracks. In the following, we shall assume a timing trigger will be implemented and let the current level-1 trigger menu [31] act as a guide for how we should model the timing trigger performance and efficiency in combination with the p T of the object. We will model the timing trigger as 50% efficient for events with out of time Ecal energy deposit greater than 1 ns (∆t > 1 ns and p T roughly equal to the level-1 trigger threshold). This will be our standard timing trigger benchmark. However, since the timing trigger's performance has not been finalized, we shall vary the trigger threshold around this benchmark value to demonstrate the search's sensitivity to different trigger thresholds. In order to register a reference time for each event, we require a jet produced from initial state radiation with p T > 30 GeV. Based on the location in η of this jet, its time of arrival, and the location of the primary vertex an expected time of arrival for any promptly produced massless particle, t prompt , can be assigned to any Ecal energy cell given it's value of η. Given the above considerations, we define our benchmark search criteria as: Event Selection Criteria (e ± or γ) 1. A leading e ± or γ with: p T > 20 GeV, |η| < 2.5 and ∆t > 1 ns 2. At least one ISR jet with p T > 30 GeV These criteria also meet or exceed our assumptions regarding the trigger threshold as discussed above.

Event simulation
The effective right-handed neutrino model is implemented using FeynRules [34] based on a modification of the model file built in [35,36] and used for [37]. Events are generated using MadGraph [38] and Pythia 8.2 [39]. Detector simulation is handled with Delphes [40] using the CMS detector card. Events are produced in Higgs production through gg fusion from pp collisions at 13 TeV center of mass energy with a scaled production cross-section of 44 pb.
We will define signal acceptance as the percentage of signal events that pass the cuts listed in the previous section after being run through the standard Delphes detector simulation using the default CMS card. Due to η and p T -dependent efficiencies on electrons and photons, detector efficiencies are also incorporated into our signal acceptance values. The signal acceptance for the search criteria is plotted as the darker of the lines in the left plot of figure 5. This left figure also shows the acceptance as the p T trigger threshold is varied and the right shows the same variations for a ∆t threshold of 2 ns. As expected, the acceptance decreases significantly as the p T threshold is increased. Increasing the P T threshold from 20 GeV to 30 GeV (40 GeV) decreases the acceptance by a factor of three (ten). Including tracking information to identify the electrons results in an order of magnitude less acceptance for the same kinematic cuts. The reason for this is the requirement JHEP07(2019)089  that the right-handed neutrino decays within R < 50 cm restricts the decays to fall within only 10% of the volume of the inner detector significantly reducing sensitivity to model parameter space with long lifetimes. Furthermore, large time delays of ∆t > 1 ns are less likely when the right-handed neutrinos decay so close to the interaction point. In these cases, we expect that an order of magnitude fewer signal events would pass cuts that utilize a timing based analysis in combination with a displaced-vertex/track type of analysis since they are somewhat mutually exclusive in this way an do not consider a tracking based analysis in this work. Figure 6 shows the signal acceptance as a function of right-handed neutrino mass for fixed mean decay length (cτ = 1 m) where acceptance is larger for heavier masses. Lighter right-handed neutrino masses are produced with larger boosts from Higgs decay compared to heavier right-handed neutrinos. This results in less time delay due to the geometric similarity of this decay process to that of promptly produced particles.
Generating events at leading order in Madgraph and then scaling up the cross-section with using the approriate k-factor, as we have done above, overpopulates (under-populates) the Higgs' low p T (higher but less than 100 GeV p T ) distribution. Explicitly including the gg → h + j process in an NLO simulation would have two main results; it would harden the p T distributions above and reduce the "kinked" path contribution to time delay due to the additional boost of provided to the Higgs. We expect that these effects would somewhat compensate each other in the acceptances above resulting in a change of no more than 10% in our results [41].

Backgrounds
Ordinary Standard Model processes in pp collisions are generically prompt compared to nano-second time scales and do not produce events with significant values for the time delay (∆t). However, jet and photon production cross-sections are so large for p T ∼ 20 GeV that very unlikely timing mis-measurements can be realized. A main source of background is time of arrival mis-measurement due to timing resolution and pile-up due to other pp collisions that occur during the same bunch crossing. We shall model these backgrounds with gaussian distributions closely following [24] and review and apply their estimates with slight modifications below. To estimate background arising from jets we utilize the low-p T data of [42] which gives a jet cross-section of σ j = 1.6× 10 8 pb for p T > 30 GeV. The crosssection for inclusive photon production are σ γ = 47 × 10 3 pb [43]. Other important values for estimating backgrounds rates are the fraction of jets that are track-less, f j = 10 −3 , the rate for a jet to fake a photon, f γ = 10 −4 [44], and a full integrated luminosity of L = 3 × 10 6 pb −1 . Production of a photon and jet can happen through actual photon production in association with a jet or through dijet production where one of the jets is mis-identified as a photon; we estimate the numbers of such events as N γ = σ γ L and N j = σ j Lf γ , respectively. Our estimate is N bg = σ γ L + σ j Lf γ = 2 × 10 11 . In fact this is conservative due to the fact that photon cross-section used in inclusive and would be less if a p T cut was put on an associated jet. If we model a δ = 30 ps time resolution with a gaussian smear ( dP dτ = 1 √ 2πδ e −τ 2 /2δ 2 ), the number of background events estimated to be mis-measured as having a time delay of greater than 1 ns is negligible since 1 ns is more than 30 standard deviations outside this distribution. Anticipating less resolution in timing for lower p T objects since the p T of signal electrons are only 20 GeV, we note that even a 60 ps resolution in timing would result in negligible background. A more relevant background however arises from pile-up. A signal-like event can arise when in addition to the collision JHEP07(2019)089 that occurs to trigger the event, there is another collision in the same bunch crossing that does not produce tracks and remains unidentified. In this case, the trackless jets associated with this second collision can be interpreted as out of time energy deposition in the ECAL. An estimate of the number of these events can be given by: N The factor of f γ is the rate for a jet to fake a photon and the other factor of f j is the rate for the recoiling jet to be trackless. This results in N (P U ) j = 10 7 . Applying a gaussian smear of 190 ps resolution to model the timing distribution, where the value of 190 ps is determined by the longitudinal spread of the bunch crossing, the background estimate for ∆t > 1 ns (∆t > 2 ns) is 0.7 (zero) from pileup.
Additional backgrounds that can lead to out of time signatures are from beam halo, satellite bunch collisions and cosmic rays. Techniques for dealing with these backgrounds have been recently developed in the analysis of [45] where the focus was on searching for BSM models with a non-prompt jet signature. Assuming that these techniques can can be applied to the delayed-electron search, we shall roughly estimate the background from these three effects based on the results of that work. Utilizing vetos based on correlations between the energy measured by cathode strip chambers (CSCs) and the ECAL as well as their respective azimuthal location, the background from beam halo effects was found to be 0.02 for an integrated luminosity of 137 f b −1 . Naively scaling this to the 3 ab −1 luminosity of the HL-LHC, we can anticipate a background from beam halo of about half an event. Satellite bunches with 10 −5 protons compared to the main bunch can be expected to undergo hard collisions at a rate of 10 −10 compared to the rate of the main bunch crossing. Thus we would expect roughly O(10) background events after scaling down the pre-timing cut background estimate of 10 11 above. However, with a total forward crosssection of σ inc = 80 mb, 100 pp collisions per bunch crossing and 10 8 pb cross-section for jets with p T > 30 GeV, only about 10% of main bunch crossings will have a jet with p T > 30 GeV. In order for satellite bunch crossings to be background they must occur immediately after a main bunch produces a jet with p T > 30 GeV. This reduces our background estimate due to satellite bunches to O(1) event. Finally, since the HL-LHC will achieve its 3 ab −1 of integrated luminosity from an increased luminosity upgrade rather than increased running time, we estimate the background from cosmic ray muons to be the same as estimated in [45] and be of order one event. Overall, we estimate that the background to our signal events from beam halo, satellite bunch collisions and cosmic rays would be a few events. Ultimately, a proper characterization of background due beam halo, satellite bunches and cosmics at the HL-LHC would require data driven approach (e.g. the use of validation regions) as was accomplished in [45]. Finally, because the particular signal described in this paper does not produce muons and the electrons are not correlated in the azimuthial directions, the vetos required to reduce background have negligible effect in signal acceptance. Due to the preliminary and rough nature of our estimates for these three background sources in particular, we choose not to include them directly in our exclusion limit plots below, but we will comment on them in the text.

Results
Using the signal acceptances discussed above, we can characterize the experimental reach of the HL-LHC in terms of an exclusion bound on the Higgs branching ratio to two right-  handed neutrinos for a given right-handed neutrino mass and lifetime. The 95% C.L. exclusion is defined for models that predict at least four (three) events at the 13 TeV HL-LHC with the full integrated luminosity of 3 ab −1 for a timing cut of ∆t > 1 ns (∆t > 2 ns); this corresponds roughly to a 95% C.L. exclusion where the ∆t > 1 ns (∆t > 2 ns) has a background of 0.7 (zero). In what follows we will consider the case with a timing trigger for events with p T > 20 GeV and ∆t > 1 ns and assume a flat trigger efficiency of trig = 0.5. In figure 7 we plot the value of the Higgs Branching Ratio (BR hN 1 N 1 ) that produces an average of four events that pass the kinematic cuts: p T > 20 GeV , |η| < 2.5, ∆t > 1 ns in the top plot and p T > 20 GeV , |η| < 2.5, ∆t > 2 ns in the bottom plot at the HL-LHC at 13 TeV with and integrated luminosity of 3 ab −1 . Since the physical neutrino mass parametrically depends on the Higgs-lepton-neutrino yukawa coupling as m ν ∼ (yv) 2 /M N , then as the right-handed neutrino lifetime is decreased for a given right-handed neutrino mass, the physical neutrino mass increases. At some lifetime for each mass, this parametric JHEP07(2019)089 dependence dictates that m ν ∼ 0.1 eV. This is indicated in each plot at the transition from solid to dashed lines. While this bound is somewhat model dependent, physical neutrino masses much larger than 0.1 eV are disfavored by standard cosmology [46]. As expected, timing sensitivity peaks near cτ = 1 m. For either the ∆t > 1 ns or ∆t > 2 ns cuts, the sensitivity is best for a Higgs branching ratio of BR hN 1 N 1 = 2 × 10 −6 . The slightly lower background of the ∆t > 2 ns cut compensates the lower acceptance, resulting in similar limits. This demonstrates that the signal is generally robust against moderate increases of the timing cut. Looking back at the discussion around figure 5, one can expect that the sensitivity will decrease by a factor of three (ten) as the p T cut on the photons is increased to 30 GeV (40 GeV). Finally we mention that for cτ in the range: 0.1 m < cτ < 10 m, exclusions were derived directly from simulation. Outside of this range however, we have extrapolated the bounds to reduce computational time. We use a simple model of exponentially small likelihood of the number (N) of long decay lengths for cτ < 0.1 m: dN dt ∼ e −t/(γτ ) and simple geometric dependence for the number (N) of short decay lengths: N ∼ (cτ ) −1 when cτ > 10 m. The boost factor for lighter right-handed neutrinos leads to a more slowly decreasing exclusion for smaller mean decay lengths. If we include an additional background of a total of three events from beam halo, satellite bunch collisions and cosmic rays, the 95% C.L. exclusion limit on branching ratio would reduce the sensitivity by a factor of two compared to what appears in figure 7.
A variation of the effective right-handed neutrino model discussed here was previously considered in [16] where a displaced track analysis reached peak sensitivity at cτ ∼ 0.1 m and became weaker for longer lifetimes. That analysis showed a gap in sensitivity in the 1 m < cτ < 10 m range before the MATHUSLA detector [47] sensitivity could potentially take over for larger mean decay lengths. This gap remains after scaling up the analysis of [16] to the anticipated integrated luminosity at HL-LHC. Based on the sensitivities shown in figure 7, we have shown that a timing-based analysis is a promising avenue to fill that gap in sensitivity to better understand and constrain the right-handed effective neutrino model discussed in this work.

Conclusion
We have shown that there exists an opportunity to probe a simplified right-handed neutrino effective theory by utilizing MIP/Ecal timing at the CMS experiment at HL-LHC. Generally, the final states of Higgs decay can include soft electrons in the final state whose time of arrival at the MIP layer and Ecal is delayed with respect to the particles produced at the primary vertex. If MIP/Ecal timing can be incorporated in such a way as to trigger on events with greater than a nano-second time delay and with p T near the existing level-1 trigger threshold, then the CMS experiment can reach sensitivities of Higgs branching ratios of order O(10 −5 ) for a range of physically motivated right-handed neutrino masses and lifetimes for the simplified model described above. For mean decay lengths of cτ > 10 m, other experiments (e.g. the MATHUSLA detector [47]) are likely to provide better sensitivity to this model. For decay lengths of cτ < 1 m we expect displaced vertex/tracking would provide better sensitivity. However for the intermediate range 1 m < cτ < 10 m,

JHEP07(2019)089
we have shown that MIP/Ecal timing information can lead to considerable sensitivity to the type I seesaw effective theory. Depending on the ability to incorporate timing into the event trigger, timing may provide the best limit for for this range of mean decay lengths.

Acknowledgments
I would like to thank Samuel Alipour-Fard, Matthew Citron, Nathaniel Craig, Seth Koren and Dave Sutherland for useful conversations regarding this work and for comments on the manuscript. I am especially indebted to Physics Department at UCSB for their hospitality in hosting me while on sabbatical from Western State Colorado University.
A Right-handed neutrino notation In this appendix we briefly describe the charged and neutral current interactions arising in the type I seesaw model using the notation of [37]. The particle content is that of the Standard Model plus three right-handed neutrinos. The fermions in the lepton sector are: In this notation the charged and neutral current interactions are: In addition to the other standard model interactions, the model has the following additional interactions at the renormalizable level: whereν andÑ are mass eignestates and L is given by: Writing the charged and neutral current interactions in terms of neutrino mass eigenstates then introduces the following interactions: where O L is the matrix used to diagonalize the charged lepton sector: O † L m l O R = diag(m e , m µ , m τ ). These interactions lead to decay modes of heavy right-handed neutrinos governed by the couplings (V † O L ) Ji . In the limit where the yukawa matrix (m d ) iJ M IJ the approximate form of V is given by V iJ = (m d ) iI (M −1 ) IJ . In order to streamline our discussion in the text, we shall define V ≡ (m d )(M −1 ) † O L andṼ ≡ U † (m d )(M −1 ). In this case the interactions read: Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.