Inclusive Displaced Vertex Searches for Heavy Neutral Leptons at the LHC

The inclusion of heavy neutral leptons to the Standard Model particle content could provide solutions to many open questions in particle physics and cosmology. The modification of the charged and neutral currents from active-sterile mixing of neutral leptons can provide novel signatures in Standard Model processes. We revisit the displaced vertex signature that could occur in collisions at the LHC via the decay of heavy neutral leptons with masses of a few GeV emphasizing the implications of flavor, kinematics, inclusive production and number of these extra neutral fermions. We study in particular the implication on the parameter space sensitivity when all mixings to active flavors are taken into account. We also discuss alternative cases where the new particles are produced in a boosted regime.


Introduction
A natural and simple extension of the Standard Model (SM) accommodating theoretical and observational challenges calls upon the existence of neutral leptons which are sterile from the SM gauge point of view (see Ref. [1] for a comprehensive review). Their unique interaction with the SM particles is via their mixing with the active neutrinos (Yukawa interaction). There is a priori no limit on the mass scale associated to these extra fermions and when their masses are below the electroweak scale, they are named 'heavy neutral leptons', HNL.
In this work we focus on collider experimental tests of their existence, in particular in the mass range up to a few tens of GeV, without invoking a link of the HNL to a specific neutrino mass generation mechanism. Many experimental searches have focused on the high mass regime (including masses above the W boson mass) where the HNL are produced directly or in some prompt decay channels [2][3][4][5][6][7][8], with numerous dedicated analysis, for instance in Refs. [9][10][11][12][13][14][15][16][17][18][19][20][21][22][23][24]. Having relatively light heavy neutral fermions that do not decouple since they could have sufficiently large mixings with active SM neutrinos may lead to important consequences, a major one being the modification of the charged and neutral currents with a leptonic mixing matrix (encoding the PMNS mixing matrix [25,26] and the active-HNL mixings). Moreover, and if sufficiently light, the HNL can be produced as final states. Both of these two features impact several observables, leading at the same time to abundant constraints on parameter space, i.e. mixing angles and masses of the HNL (see Refs. [12,27] and references therein).
An exciting possibility occurs in a region of parameter space for which the HNL are long-lived particles that can decay with a sizable displacement in the LHC detectors, and such displaced vertices (DV) would be a distinctive signature of their existence [28]. Recently, several studies have pushed for this kind of dedicated LHC searches for HNL with associated charged leptons [29][30][31][32][33][34], from Higgs decays [35][36][37], for LHCb [38] or at future detectors proposed for searching for long-lived particles [39][40][41]. There is also a potential to search for such signatures at DUNE [42], future lepton colliders [43,44] and SHiP [45,46], the latter expected to highly improve sensitivity to HNL below the charm quark mass, with the HNL being abundantly produced by meson decays [47].
At the LHC, the cross sections of the single production of W ± , Z and H bosons have been measured [48][49][50][51] and given the expected integrated luminosity, it would be possible to have a large number of events in which these bosons decay into a HNL. In addition, the LHC has successfully measured a large number of SM processes such as diboson, tt, W/Z+ jets, etc. As mentioned above, the modification of the charged and neutral currents due to the active-sterile mixings opens up additional decay channels for these SM processes occurring at the LHC, in which a W ± , Z and the H boson will decay.
In this work we address the DV signature in two types of scenarios: i) large number of events like the Drell-Yan processes, ii) other SM processes with a smaller cross section mostly relevant for the high luminosity phase of the LHC. In both cases we look at the DV signature possibility focusing on the production mechanisms of the HNL when they are long-lived and decay in the detector.
In the first case, we consider the inclusive production from Drell-Yan, W ± , Z and H decay processes for which we do not apply selection cuts on the primary vertex producing the HNL, but instead focus on the signature of the DV from the HNL decay. This allows a larger number of total events. The experiments at the LHC will be sensitive for a given region of the parameter space (mixings and masses of the HNL) that we identify in this work.
We also exemplify in this work the flavor dependence for a particular production and decay channel given that the decay width is proportional to the sum of the square of the mixings, implying that all these channels are complementary and necessary to probe the whole parameter space (masses of the HNL and their mixings to the SM active neutrinos); for this we consider different cases of the active-sterile mixing pattern. This is in contrast with what has been done for instance in Refs. [29][30][31][32][33], where only one flavor was considered.
In the second case, we consider other SM processes in which the gauge bosons are produced in combination with other particles. In this case, the gauge bosons decay to boosted HNLs and the kinematic of these decays might be significantly different from that of the Drell-Yan processes. In particular, for those events in which the HNL is strongly boosted, distinctive DV signatures can take place. We illustrate in this work the potential of this DV signature for a specific process.
Throughout the study, we do not consider a specific neutrino mass generation mechanism but rather consider a bottom-up approach extending the SM with N sterile fermions. Most of the analysis is conducted for N = 1 and extended to N = 2 to show (new) distinctive features of DV signatures compared to the case with N = 1.
This work is organized as follows: in Section 2 we present the framework we adopt that is the simple extension of the SM via sterile fermions, 3 + N and provide details of the numerical simulations we perform. Section 3 is devoted to the determination of the HNL decay width discussing the different HNL decay channels and highlighting the role of flavor. We present and discuss the results of DV arising from inclusive HNL production at the LHC in Section 4, considering both flavored and flavor blind cases for the production of HNL. We also extend the study to the 3 + 2 model (N = 2). Finally, we discuss in Section 5 the results we obtain in a specific SM process in which the HNL is strongly boosted (W/Z+ jets). Our results are summarized in Section 6. The parametrization we use for the 3 + N model (N = 1, 2) is given in Appendix A.

Theoretical Framework
In this work we are interested in studying collider phenomenology of HNLs without assuming any specific underlying model or mechanism of light neutrino mass generation. In order to do this, we follow a bottom-up approach where the SM is extended by ad-hoc masses for the 3 active neutrinos, as required by oscillation phenomena, and by N additional sterile fermions (the HNL). We refer to this kind of scenarios as 3 + N models, and they are useful to understand the general phenomenology of a broad class of models where the SM is enlarged only by sterile fermions, such as the type-I seesaw model and its variants.
In most of our forthcoming analysis, we will focus on the simplest case of having only one HNL (N = 1), although we will discuss the generalization to the case with more HNLs in Section 4.4. We will also assume that neutrinos are Majorana fermions. The 3 + 1 model consists of four neutrino masses and a unitary 4 × 4 leptonic mixing matrix, where the 3 × 3 matrixŨ PMNS is the usual PMNS matrix up to non-unitarity corrections due to the presence of light-sterile V N mixings, with = e, µ, τ . These mixings, together with the HNL mass m N , define the interaction strength of the HNL via charged currents, as well as the neutral currents to both Z and H bosons. Therefore, they will be the relevant parameters for our study. 1 In the case of an extension of the SM with N HNL, the relevant terms in the Lagrangian (in the Feynman-'t Hooft gauge) are given below: where g is the SU (2) L gauge coupling, U αi ≡ U 3+N ν αi are the lepton mixing matrix components, m i are the mass eigenvalues of the neutrinos and m α are the charged lepton masses. The indices α and i run as α = e, µ, τ and i = 1. . . 3 + N .
The existence of a HNL has been tested in different observables, depending on its mass, and at present there are strong bounds on its mixings (see Ref. [27] for a recent update). In the range of masses we will be interested on, i.e. the few GeV regime, the strongest upper bounds are those from DV searches by DELPHI [2], which constrains the sum of the squared mixings to be below 2 × 10 −5 . We will consider these bounds in our numerical analysis and discuss how the LHC could improve them searching also for DV signatures.
For our forthcoming numerical simulations, we have implemented the 3+1 and 3+2 models in FeynRules [52,53] to generate the UFO model file [54], taking also into account the effective Higgs coupling to gluons for Higgs production. Then, we use MadGraph5 aMC@NLO [55] to generate the HNL production events from proton-proton collisions and Pythia 8.2 [56] for its subsequent decay and treatment of the DV. Finally, Les Houches Event files are obtained using MadAnalysis5 [57]. The HNL lifetimes are included using the time of flight option in MadGraph5.

Decays of the Heavy Neutral Lepton
When exploring signatures from DV, the total width of the decaying particle is a crucial parameter. Generally speaking, its decay length must be of the same order of the size of the detector, or more specifically of the tracker system, which in the case of detectors such as ATLAS or CMS cover transverse displacements between roughly 1 mm and 1 m from the interaction point. In this section we explore the decays of the HNL and discuss the parameter space that could be explored by DV searches at the LHC. In the following, we will consider the Majorana HNL case, whose total width is twice the one in the Dirac case, since it can decay to both CP conjugated final states. The decays of the HNL can be divided in different regimes, depending on its mass [12,27,58]. When the mass is below the GeV scale, it mainly decays via off-shell W or Z bosons, leading to three body leptonic or two body semileptonic final states, i.e. N → ν, N → ννν or N → M , M being a light meson. Above the GeV but still below the W mass, the semileptonic decays are better estimated by the three body decay to quarks, N → qq , which also accounts for the possible hadron multiproduction. Finally, if the HNL is very massive, m N > m W , it tends to decay very fast via on-shell W , Z and H bosons, and, thus, one needs to search for prompt decay signatures.
We show in Fig. 1 the total neutrino width Γ N (left) and the decay length τ N c (right) as a function of the HNL mass m N and mixing |V N | in a "democratic" scenario of |V eN | = |V µN | = |V τ N |. The vertical gray line is only to illustrate the transition between the semileptonic decay to mesons, N → M , and to quarks, N → qq . We see how the increase with m N and |V N | in the left panel is translated into the diagonal contour lines in the right one. The black lines in the latter show contours for τ N c = 10 3 , 1, 10 −3 and 10 −6 m, while shaded areas are experimentally excluded [27]. Interestingly, the region with decay lengths relevant for DV searches at the LHC, say τ N c ∈ [1 mm, 1 m], lies in the few GeV region, where present experimental constraints on the mixing angles are weaker. In a general case of three different mixing angles, we have checked that the general behavior of the total width can be approximated as Distributions of the DV coming from a HNL decay in the (z DV , DV ) plane. The production channel is pp → e ± N , with p e T > 25GeV and |η e | < 2.5. We fix m N = 15 GeV, |V eN | 2 = 10 −7 and V µN = 0 in all the cases, while |V τ N | 2 = 0 (red), 10 −6 (yellow) and 10 −5 (purple). As |V τ N | 2 increases, the HNL becomes more prompt and, therefore, insensitive to the DV searches. The white region is the potential DV area.
which works very well within this area of interest for the DV at LHC, especially above the tau mass threshold.
At this point, a remark about the role of flavor is in order. In many collider searches for HNL, one often considers simplified hypothesis where the HNL mixes only to one flavor at a time. This simplification is well justified in prompt decayed HNL searches with dilepton or trilepton final states, since the number of events with charged leptons of a flavor depends only on the corresponding mixing V N . Nevertheless, this does not apply to the DV searches. In this case, the total decay width plays a crucial role defining where the HNL will decay and, since it depends on all the mixings V N , = e, µ, τ , one cannot conclude independently of each of the mixings in DV searches.
In order to illustrate this effect, we display in Fig. 2 the distributions for the HNL decay position in the (z DV , DV ) plane, where z DV is the displacement along the beam axis and DV in the transverse plane. We have generated pp → e ± N events with fixed values of m N = 15 GeV, |V eN | 2 = 10 −7 and V µN = 0, and take the |V τ N | 2 mixing equal to zero (red), 10 −6 (yellow) and 10 −5 (purple). Since the mass and electron mixing are fixed, the same number of HNL are produced in all cases. Nevertheless, the sensitivity of the DV searches to the HNL depends on the amount of events within the DV fiducial volume, DV ∈ [1 mm, 1 m] and z DV < 300 mm. In this example, this kind of searches are very sensitive to the single mixing case, as most of the red points are within this area. However, increasing values of V τ N enhances the total width without affecting the production cross section, shifting the distribution towards lower displacements and thus reducing the efficiency of the DV searches. Alternatively, some of the cases escaping the detector in the single mixing scenario could lead to DV signatures when all mixings are taken into account. Consequently, this kind of searches cannot explore independently the mixing of each flavor and they should take into account the flavor combination entering in the total width in Eq. (3.1). We will discuss the impact of this effect on the final sensitivity estimates in Section 4.        Before concluding this section, we focus on the branching ratios for the different HNL decay channels, which are important to understand the kind of signature we can expect in the DV. As we already said, in the few GeV mass region that we are interested in, the main decay channels are N → ν, N → ννν and N → qq . Since all of them are three body decays, they will have the same dependence with m N , implying that the branching ratios will be almost mass independent. The most relevant dependence will be that from the mixings |V N |, or more precisely from the ratio of mixings |V eN | 2 :|V µN | 2 :|V τ N | 2 , as the total width normalizes their global size. In this situation, one can display the different decay channels in triangular plots, like those given 2 in Fig. 3.
The largest branching ratio is for the semileptonic N → qq decay, with a total of approximately 60% divided between the three flavors of the lepton. The concrete value of each channel depends on |V eN | 2 :|V µN | 2 :|V τ N | 2 , and therefore the ratio between different final state flavors would provide information on the mixing ratios. Unfortunately, these semileptonic channels are more difficult to explore in a hadronic collider. On the other hand, purely leptonic channels have a clearer signature, however their branching ratios are smaller, especially those with leptons of the same flavor. Interestingly, channels such as N → e ± µ ∓ ν have a larger ratio and could be promising to explore, since one could take advantage of having two leptons of a different flavor from a DV. We will discuss further the potential of each decay channel is Section 4.

Displaced Vertices from Inclusive HNL Production at the LHC
In this work we consider the production of a HNL in the regime in which there could be a detectable DV, i.e. m N in the GeV ballpark. We restrict our analysis to the case in which the displaced decay occurs in the inner tracker of ATLAS or CMS, although it is possible that the HNL decay occurs in the calorimeters.
The dominant HNL production channels involve mesons for low HNL masses, below the b-quark mass [58], however in this case i) the decay products of such a light HNL are typically very soft and difficult to study at the LHC, and ii) the HNL tends to be long-lived and to decay outside the detector. These production channels are important for beam-dumped experiments such as SHiP [45,46] or for proposed future detectors dedicated to long-lived particles searches [39][40][41], among others.
For the mass region we are interested in, the relevant production mechanism is via W ± , Z and H bosons that decay into a HNL. Some other (subdominant) channels including additional particles are also interesting to explore, since their different kinematics gives access to alternative regions of the parameter space via DV signatures, as will be discussed in Section 5. Consequently, we focus on the production from W/Z/H boson decays, which are indeed the main channels for the mass region we are interested in. Some important features should be noted: • The decays from W ± and Z bosons are largely independent of m N , for m N m W and m Z .
• The decays from W ± bosons are flavor dependent, proportional to a single flavor mixing.
• The decays from Z and H bosons are flavor independent, proportional to the sum of the square of all mixings, which is the same combination that enters in the HNL total decay width, Eq. (3.1).
• The Higgs boson decay has an additional suppression proportional to the square of the ratio of m N /m W .
The branching ratios can then be expressed as: where for Z, H boson decays we have summed over the three light neutrinos, as one cannot distinguish them at the LHC. Though the contribution from Higgs bosons is subdominant for small values of the HNL mass, this production channel is relevant to infer information on the neutrino mass generation mechanism. We show in Fig. 4 (left) the explicit values of these branching ratios as a function of HNL mass for an example of the mixing pattern: |V eN | 2 = |V µN | 2 = 10 −6 and |V τ N | 2 = 4 × 10 −6 . Full analytical expressions for these decays can be found e.g. in Refs. [35,59], and are therefore not reported here. As a large number of W ± , Z and H bosons are produced at the LHC, we can set parameter space sensitivity limits for this inclusive DV analysis. We show in Fig. 4 (right) the number of HNL produced in pp collisions at √ s = 13 TeV, with an integrated luminosity L = 300 fb −1 . We have generated the pp → N and pp → νN processes and we have checked that they are indeed dominated by the production and decay of on-shell W ± and Z bosons, respectively. We emphasize however that the relative importance of each W ± channel depends on the relative size of each mixing, unlike the case of the Z boson decays. In Fig. 4 (right) we show as well the events corresponding to the (also) flavor blind production channels pp → νN j and pp → νN γ. These channels could be useful as the initial state radiated jet  Ev(pp → N X) or photon could be used as trigger for the interaction point in the Z channel, or to access different kinematic regimes of the HNL. In order to estimate the impact of the chosen HNL flavor hypothesis on DV searches, we define three illustrative benchmark scenarios for the forthcoming analysis: Based on our results, we classify the different production channels into flavor dependent and independent categories, showing how the sensitivity predictions of DV searches vary in the former case but not in the latter. This motivates us to focus on the inclusive production analysis for the DV with both leptonic and semileptonic HNL decays.

Flavor Dependent Production in pp → N
The details of the leptonic flavor structure are prominent only in this decay channel dominated by W ± decays, thus in this case any experimental analysis will be flavor dependent. This channel is interesting since the charged lepton, if detected, could be used as a trigger for the primary vertex. This production mechanism has been explored in the literature [29,30] with the conclusion that the LHC could probe mixings up to |V N | 2 ∼ 10 −7 , with = e, µ, after collecting L = 300 fb −1 of data, implying that the LHC could explore this mass range beyond present constraints. In all our figures experimental bounds correspond to the shaded green areas of the plots. Nevertheless, these analyses assumed that the HNL mixes only to one flavor, which is a condition that is typically not realized in most of the BSM models, and the deviations from this simplified hypothesis may change the conclusions.
The impact of the mixing to different flavors is exemplified in the left panel of Fig. 5, where we consider DV events from the HNL production in association with an electron or a positron for the three benchmark scenarios described in the previous section. We have imposed cuts on the electron and positron (p e T > 25 GeV and |η e | < 2.5), and on the DV (1 mm < DV < 1 m and z DV < 300 mm). The solid (dot-dashed) lines are contour lines for 10 2 (10 4 ) events. From this figure we can see that the number of events, and therefore the sensitivity, is different in each scenario, the most optimistic numbers being those of the simplified scenario with only one non-vanishing mixing. The differences come from the fact that the total width, defining the area where DV may occur, depends on the sum of all mixings, see Eq. (3.1), while the production rate of pp → N is only sensitive to |V N |. This is the same effect we already discussed in Fig. 2. Therefore, when considering mixings also to other flavors, the relative importance of |V eN | decreases and so does the sensitivity via the pp → e ± N channel. Moreover, notice that the number of events in Fig. 5 corresponds to the total number of produced HNL that would decay, to any channel, at a DV. Therefore the full process including the HNL decay to a given channel will have an even stronger dependence on the flavor hypothesis.

Flavor Independent Production in pp → ν N
We have performed the same exercise for the pp → ν N production channel, where ν stands for the sum over the three light neutrinos. The three considered scenarios lead to the same number of DV events, with small differences at low values of m N m τ , and are thus shown as a single color in the right panel of Fig. 5. As we explained before, both the total width and the production rate depend on the same combination |V eN | 2 + |V µN | 2 + |V τ N | 2 . Therefore, this flavor blind production mechanism is closely related to DV searches, as it was already pointed out when exploring HNL via H decays [35].
The idea of the flavor independent production can be extended to an inclusive search, where we focus only in detecting the DV. In this case, the production cross section and the decay width have the same flavor dependence, and thus one can set bounds without making further assumptions. Of course, a flavor dependence will enter if the HNL decays to charged leptons, nevertheless it can be explored by searching for different final state channels, as discussed in the next section.

Inclusive HNL Production and Displaced Vertices
In order to focus on an inclusive HNL production, we consider all possible contributions pp → XN (referring to pp → ν N and pp → N with = e, µ, τ ) followed by the decay of HNL into different channels, as discussed in Section 3. This inclusive production is flavor blind and it depends on the sum of squared mixings, similar to the HNL total width. Therefore, we show the sensitivity plots in terms of the combination |V eN | 2 + |V µN | 2 + |V τ N | 2 , the proper variable to be explored in DV analysis. We do not place any cuts on the primary vertex decay products. However, in order to identify the DV and select the events, we have imposed the following cuts: • 1 mm < DV < 1 m and z DV < 300 mm, • |η | < 2.5 and |η j | < 2.5 for tracks from the DV, • p e T > 10 GeV, p µ T > 5 GeV and p j T > 10 GeV for tracks from the DV, • m , m jj < m N when there are two charged leptons or one charged lepton and two quarks from the DV, • ∆R , ∆R jj < 1 when there are two charged leptons/quarks from the DV. Fig. 6 shows the number of events for the inclusive HNL production and different decay channels: N → e − e + ν, µ − µ + ν, µ ± e ∓ , e ± qq and µ ± qq . Different flavor hypothesis are explored in different colors, |V eN | 2 : |V µN | 2 : |V τ N | 2 =1:0:0 (black), 1:1:0 (blue) and 1:1:1 (green). The solid (dot-dashed) lines are contours for 5 (100) events with a DV at the LHC 13 TeV with L = 300 fb −1 . The shaded green areas are the experimental bounds. 3 Let us first note that contours in Fig. 6 follow the ones of Fig. 5 but with a smaller number of events, due to the suppression of the branching fractions and to the selection cuts. For the channel pp → X N → X e − e + ν (upper left panel of Fig. 6), the maximum sensitivity is reached in the case where the HNL mixes only to electrons (black lines). However, if the HNL is allowed to mix with the muons (blue lines) or the muons and the taus (green lines), other decay channels open up, reducing the branching ratio of the HNL into e − e + ν and therefore decreasing the sensitivity to that particular channel. A similar behavior is featured in the pp → X µ ± e ∓ ν channel (central panel), although the flavor dependence is milder. In this case there are more events due to the larger branching ratios to this channel, and the experimental signature is more interesting because of the different flavor of the final charged leptons. On the contrary, for the pp → X µ − µ + ν channel (upper right panel), the sensitivity is minimal (but not zero) for the case when V µN = 0. In fact, even if the HNL cannot decay via the W ± bosons, it can always via the exchange of a Z boson. The sensitivity to this channel grows with the mixing |V µN | and, in this sense, it is complementary to the purely electronic channel.
The case for the semileptonic decays (lower panels) resembles the pure leptonic one. However, some comments are in order. First, the pp → X µ ± qq channel (lower right panel) vanishes in the scenario where V µN = 0, because the corresponding decay can only occur by  the mediation of a W ± . Second, the semileptonic branching ratios are typically larger, and so is the sensitivity. Finally, the possible backgrounds are also more important in the fully leptonic channels. Nevertheless, one would like to finally explore all the possible final states, as they provide complementary information about the flavor structure of the HNL.

Two Heavy Neutral Leptons
We close this section by discussing the possibility of having more than one HNL in our region of interest. In general, the 3+1 model already captures most of the collider implications of the existence of a HNL. Nevertheless, in seesaw-like models at least two heavy neutrinos are needed in order to accommodate neutrino oscillation data and, consequently, some remarks should be made in this latter case with two HNL, which we effectively describe by a 3+2 model. Under the lack of a positive experimental signal, one might consider the possibility of reinterpreting the bounds on the 3+1 model in terms of the 3+2 parameters. The production cross section in the 3+2 would be the sum of the 2 HNL contributions and, if they do not have a large separation in mass, one can set bounds on the sum |V N 1 | 2 + |V N 2 | 2 . However, the decay length of each HNL depends on each mixing separately, and therefore the interpretations in DV searches are more involved.
On the other hand, if a positive signal is found, one would be interested in disentangling the N HNL hypothesis from the single one. We show in Fig. 7 two differential distributions that could help to differentiate the N = 2 from the N = 1 case. This particular example is for the process pp → e + N i with N i → e − µ + ν and shows the invariant mass m DV (left panel) and ∆R DV (right panel) of the two charged leptons originating at the DV. Dashed lines represent the contribution of each HNL to the total red distribution that would be observed. It is particularly interesting in the case of m DV , since the individual distributions show edges at m N i , which translates to a kink-like distortion in the total spectrum. The distribution of ∆R DV is as expected in this case, with a small deviation from the single HNL hypothesis at large values for ∆R DV .

Displaced Vertices in Other Standard Model Processes
The success of the LHC and detectors like ATLAS, CMS and LHCb have allowed new and more precise measurements of SM processes. The embedding of a HNL into the SM allows for all typical SM processes that produce Z, W ± and H bosons to have additional decay channels. In the previous sections we discussed the possibilities when a single Z, W ± or H boson is produced. We now explore additional processes such as: • Dibosons, • t, tt with one the W ± bosons decaying via the HNL, considering the case in which a single gauge boson is decaying via a HNL. As before, depending on the relevant variables m N and the mixing angle the HNL can be either prompt or long-lived in these SM channels. Nevertheless, the different kinematics of these processes may modify the parameter space area where displaced vertices happen.
In particular, for these processes the gauge boson can have a large transverse momenta p T and for the mass ranges we are considering the corresponding HNL will be boosted. The smoking gun signature would require a single DV as in the analysis presented in the previous section. This implies that the whole effective displacement region shown in Figs. 5 and 6 would be shifted to the right. To clarify, this will allow a specific point in the (m N , |V N | 2 ) plane for which the N decay is prompt when it is produced via a Drell-Yan process to be long-lived when produced via one of the above mentioned SM processes. Therefore the precise transition between the prompt and non-prompt decay of the HNL is slightly different in these processes with respect to the W ± and Z decays we explored in the previous section.
We illustrate this effect for the case of the process pp → e + N j, where we can induce a boosted topology for the HNL by asking for a high-p T jet. If we select events that have a large transverse momenta of the HNL, then the decay length becomes much larger and significant. Fig. 8 (left) shows the shift in the transverse displacement from the boosted HNL. We see that for a non-negligible shift can be obtained, especially when a very hard jet is present. Moreover, the decay products of the HNL become more collimated in this case, see Fig. 8 (right), which can help to reduce backgrounds when an appropriate cut is implemented. Nevertheless, these strong cuts severely reduce the number of events and therefore this boosted HNL scenario could be probed only at the high luminosity LHC.
A more dedicated analysis is needed to study the potential of these SM processes, exploring the different channels mentioned above and considering more energetic hadron and leptonic colliders. Although experimentally challenging, these new processes would aid in exploring areas of parameter space that could be difficult to access via other channels.

Conclusions
We have revisited collision phenomenology when heavy neutral leptons (HNL) are added to the particle content of the SM. We do not consider a specific mechanism or model in which these additional neutral fermions are embedded. We have especially focused on the region of parameter space of mass and mixing of the heavy neutral leptons for which the distinctive signature of a displaced vertex (DV) in the inner tracker of a detector is produced when this neutral fermion decays at the LHC.
We have further analyzed the implications on the decay width and the corresponding decay length when different flavor configurations are considered in the production channel via W ± and Z gauge bosons. While Z and H bosons decays are flavor blind, the W ± boson decays are not and the decay length is modified when the HNL couples to more than one flavor with similar (or larger) strength.
Our next step has been to consider an inclusive and flavor blind production in order to determine the precise values of the mass and mixing angles of the HNL in which a sufficient number of events with measurable displaced vertices can occur. We find that a dedicated experimental analysis could improve the bounds by at least one order of magnitude in the sum of the squared mixings for values of m N 10−15 GeV, and for an integrated luminosity of 300 fb −1 .
We also show how the distribution for the invariant mass of the leptons (or quarks) produced at the DV can help distinguishing the number of additional heavy neutral leptons that are important in the leptonic mixing matrix.
For the mass ranges of the HNL we are considering for all SM processes in which a W ± , Z or H boson is produced, their decays are modified by the additional channels that contain the HNL and its subsequent decays. However, the kinematics associated to the HNL varies significantly according to the specific SM process in which it is produced. We have shown the impact that a boosted HNL has on the decay length and state that the analysis of the DV signatures is relevant for these SM processes in the high luminosity phase of the LHC.
In summary, we have highlighted in this work the importance of multiple production channels in order to improve the constraints on the region of parameter space for which DV signatures occur when the HNL decays. We clarified the role of flavor and kinematics showing how different SM processes and channels are complementary to establish the existence of a HNL with a mass on the order of a few tens of GeV at the LHC.
Grant FPA2017-84543-P. ML thanks LPT-Orsay and CEA-Saclay for hospitality during the completion of this work. In addition to the software packages cited above, this research made use of IPython [60], Matplotlib [61], SciPy [62] and python-ternary [63].

A Parametrization of the 3 + N Minimal SM Extensions
The extension of the SM with one HNL state reflects into three new mixing angles (θ 14  where R ij is the rotation matrix between i and j. For instance, the rotation matrix R 45 is explicitly given by and likewise for the other matrices R ij (in terms of θ ij and δ ij ).