Search for Heavy Neutral Leptons in Events with Three Charged Leptons in Proton-Proton Collisions at  <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mrow><mml:msqrt><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msqrt><mml:mo>=</mml:mo><mml:mn>13</mml:mn><mml:mtext> </mml:mtext><mml:mtext> </mml:mtext><mml:mrow><mml:mi>TeV</mml:mi></mml:mrow></mml:mrow></mml:math>

A search for a heavy neutral lepton N of Majorana nature decaying into a W boson and a charged lepton is performed using the CMS detector at the LHC. The targeted signature consists of three prompt charged leptons in any flavor combination of electrons and muons. The data were collected in proton-proton collisions at a center-of-mass energy of 13 TeV, with an integrated luminosity of 35.9 fb−1. The search is performed in the N mass range between 1 GeV and 1.2 TeV. The data are found to be consistent with the expected standard model background. Upper limits are set on the values of |VeN|2 and |VμN|2, where VℓN is the matrix element describing the mixing of N with the standard model neutrino of flavor ℓ. These are the first direct limits for N masses above 500 GeV and the first limits obtained at a hadron collider for N masses below 40 GeV.

section_content: The standard model (SM) of particle physics has been successful in describing many phenomena, however several observations remain unexplained, including the nature of dark matter (DM), the baryon asymmetry of the universe, and the smallness of neutrino masses. The latter can be naturally accommodated by the so-called "seesaw" mechanism [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11], in which a new heavy neutral Majorana lepton (or heavy neutrino) N is postulated. 

A model that incorporates the seesaw mechanism, while also providing a DM candidate and giving a possible explanation for the baryon asymmetry, is known as the neutrino minimal standard model (νMSM). In this model, three right-handed heavy sterile neutrinos are added to the SM [12, 13]. The lightest neutrino, N 1, can explain the DM in the universe, while the heavier neutrinos, N 2 and N 3, can be responsible for the baryon asymmetry through leptogenesis [14] or neutrino oscillations [15]. 

In this Letter, a search for a heavy neutrino decaying into a charged lepton and a W boson (either an on-shell W or off-shell W * ) is performed using the CMS detector at the CERN LHC [16]. The heavy neutrinos are produced through mixing with the SM neutrinos, governed by the parameters V eN, V µN, and V τN, where V N is the matrix element describing the mixing of N with the SM neutrino of flavor. The production cross section, decay width, and lifetime of N depend on |V N | 2 and its mass m N. In the νMSM, N 1 is expected to be too light and long-lived to produce an unambiguous signal in the CMS detector, but N 2 and N 3 could decay to W, Zν, or Hν, and are therefore potentially detectable. 

Earlier searches for heavy Majorana neutrinos at the LHC have been undertaken by the AT-LAS and CMS Collaborations at √ s = 7 and 8 TeV [17] [18] [19] [20] [21] [22] [23] [24], employing a signature of same-sign dileptons and jets, exploring the mass range 40 < m N < 500 GeV. Other experiments have searched in the mass region m N < 40 GeV [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] and precision electroweak measurements provide limits on the mixing parameters independent of m N [39] [40] [41] [42] [43]. A recent review of constraints can be found in Ref. [44]. This analysis targets N production in leptonic W ( * ) boson decays, W ( * ) → N ( = e, µ), with subsequent prompt decays of N to W ( * ), where the vector boson decays to ν [45] [46] [47] [48] [49] [50] [51]. The event signature consists of three charged leptons in any combination of electrons and muons, excluding those events containing three leptons of the same charge. Because of the presence of a SM ν escaping detection, a mass peak of N cannot be reconstructed. Therefore the search exploits kinematic properties of the three leptons to discriminate between the signal and SM backgrounds. These backgrounds consist of events containing leptons from hadron decays, leptons from conversions, and SM sources of multiple leptons such as diboson production or top quark (pair) production in association with a boson. Exploiting the trilepton topology allows the mass range 1 GeV < m N < 1. 2 TeV to be explored using pp collision data collected by the CMS experiment at √ s = 13 TeV, corresponding to an integrated luminosity of 35. 9 fb -1. 

The central feature of the CMS apparatus [52] is a superconducting solenoid of 6 m diameter, providing a magnetic field of 3. 8 T. A silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter, each composed of a barrel and two endcap sections, reside within the solenoid. Forward calorimeters extend the pseudorapidity (η) coverage. Muons are detected in gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid. Events of interest are recorded with several triggers [53], requiring the presence of one, two, or three light leptons (e or µ), leading to very high efficiency, nearing 100% in most kinematic regions of the search. 

Samples of simulated events are used to estimate the background from some of the SM processes and to determine the heavy neutrino signal acceptance. The SM background samples are produced using the Monte Carlo (MC) MADGRAPH5 aMC@NLO 2. 2. 2 or 2. 3. 3 generator [54] at leading order (LO) or next-to-leading order (NLO) in perturbative quantum chromodynamics, with the exception of gg → ZZ, which is simulated at LO with MCFM 7. 0 [55], and all other diboson production processes, which are generated at NLO with the POWHEG v2 [56, 57] generator. 

The NNPDF3. 0 [58] LO (NLO) parton distribution functions (PDFs) are used for the simulated samples generated at LO (NLO). Parton showering and hadronization are described using the PYTHIA 8. 212 generator [59] with the CUETP8M1 underlying event tune [60, 61]. Double counting of the partons generated with MADGRAPH5 aMC@NLO and PYTHIA is removed using the MLM [62] and FXFX [63] matching schemes in the LO and NLO samples, respectively. 

Signal samples are generated with MADGRAPH5 aMC@NLO 2. 4. 2 at NLO precision, following the implementation of Ref. [64]. They include processes leading to N production via the charged-current Drell-Yan (DY) process, gluon fusion, and Wγ fusion. The latter production mechanism is important for m N > 600 GeV. The first two production mechanisms employ the NNPDF3. 0 NLO PDF set [58], while the last uses the LUXqed plus PDF4LHC15 nnlo 100 PDF set [65]. Parton showering and hadronization are simulated with PYTHIA. Only the final states with three leptons (electrons or muons) and a neutrino are generated. 

The effects of additional pp interactions in the same or adjacent pp bunch crossings (pileup) are accounted for in the simulations. The MC generated events include the full simulation of the CMS detector based on GEANT4 [66] and are reconstructed using the same CMS software as used for data. 

Information from all subdetectors is combined offline by the CMS particle-flow algorithm [67] used to reconstruct and identify individual particles and to provide a global description of the event. The particles are classified into charged hadrons, neutral hadrons, photons, electrons, and muons. 

Jets are reconstructed using the anti-k T clustering algorithm [68] with a distance parameter of 0. 4, as implemented in the FASTJET package [69, 70]. Jet energies are corrected for residual nonuniformity and nonlinearity of the detector response using simulated and collision data event samples [71] [72] [73]. 

To identify jets originating from b quarks, the CSVv2 algorithm [74, 75] is used. Jets with p T > 25 GeV and |η| < 2. 4 are considered b quark jets ("b jets") if they satisfy the loose working point requirements [75] of this algorithm. Events with one or more identified b jets are vetoed in the analysis to reduce the tt background. 

The missing transverse momentum p miss T is defined as the magnitude of the negative vector sum p miss T of the transverse momenta of all reconstructed particles in the event, and is adjusted to account for jet energy corrections applied to the event [76]. 

The primary pp interaction vertex (PV) is taken to be the reconstructed vertex with the largest value of summed physics-object p 2 T. The physics objects are the jets, clustered using the jet finding algorithm [68, 69] with the tracks assigned to the vertex as inputs, and the associated missing transverse momentum, taken as the negative vector sum of the p T of those jets. 

The analysis depends crucially on the clean and efficient identification of electrons and muons. Electrons are reconstructed by combining the information from the ECAL and the tracker [77]. Electrons are required to have |η| < 2. 5, ensuring that they are within the tracking system volume, and a minimum p T of 10 GeV. Electron identification is performed using a multivariate discriminant that is sensitive to the shower shape and track quality. To reject electrons origi-nating from photon conversions in detector material, electrons must have measurements in all innermost layers of the tracking system and must not be matched to any secondary vertex [77]. 

Muons are reconstructed by combining the information from the tracker and muon spectrometer in a global fit [78]. The geometrical matching between measurements made separately in the two systems is used in the further selection of muons. Only muons within the muon system acceptance of |η| < 2. 4 and with p T > 5 GeV are considered. Electrons within a cone of ∆R < 0. 05 of a muon are discarded as those likely coming from radiation. 

To ensure that electron and muon candidates are consistent with originating from the PV, the transverse (longitudinal) impact parameter of the leptons with respect to this vertex must not exceed 0. 5 (1. 0) mm, and the displacement divided by its uncertainty must not exceed 4. 

Leptons originating from decays of heavy particles, such as electroweak bosons or N, are referred to as "prompt", while leptons produced in hadron decays are called "nonprompt". For convenience, we also include in the nonprompt-lepton classification, particles that result from misidentified hadrons or jets. A powerful discriminator between these two types of leptons is the isolation variable I rel. It is defined as the pileup-corrected scalar p T sum of particles within a cone of ∆R < 0. 3 around the lepton candidate's direction at the vertex, divided by the lepton candidate p T. The summation comprises the reconstructed charged hadrons originating from the PV, neutral hadrons, and photons. 

Electrons and muons that pass all the aforementioned requirements and satisfy I rel < 0. 6 are referred to as "loose leptons". Leptons that additionally satisfy I rel < 0. 1 and, in the case of electrons, pass a more stringent requirement on the multivariate discriminant, chosen to maximize the signal over background ratio, are referred to as "tight leptons". Events containing exactly three loose leptons, not all having the same charge, are retained in the analysis. 

To distinguish between SM background and N production, the three leptons are required to pass the tight selection, and the following variables are used: the flavor, charge, and p T of the leptons in the event; the invariant mass of the trilepton system M 3 ; the minimum invariant mass of all opposite-sign lepton pairs in the event M min 2 OS ; and the transverse mass

, where p T is the transverse momentum of the lepton that is not used in the M min 2 OS calculation, and ∆φ is the azimuthal angle between p T and p miss T. 

To address the kinematically distinct cases of N masses below and above that of the W boson, two search regions are defined, referred to as the low-and high-mass regions. 

In the low-mass region (targeting m N < m W ), N is produced in the decay of an on-shell W boson, leading to the decay W → ν via an intermediate N. To reflect the targeted m N range in this region, and to suppress the background from Z → + -production with an accompanying high-p T lepton from an asymmetric photon conversion, the requirement M 3 < 80 GeV is imposed. The background from Wγ * events is reduced by rejecting events that contain an opposite-sign same-flavor (OSSF) lepton pair. The effectiveness of this requirement relies on the fact that N is a Majorana particle and can decay to a lepton of equal or opposite charge to that of its parent W boson. 

Events in the low-mass region are required to have p miss T < 75 GeV to suppress tt background. The highest p T (leading) lepton must satisfy p T > 15 GeV, while the next-to-highest p T (subleading) lepton must have p T > 10 GeV. The third (trailing) lepton must have p T > 10 (5) GeV if it is an electron (muon). In addition the following conditions are imposed to avoid trigger threshold effects: in eµµ events, if a trailing electron has 10 < p T < 15 GeV, the leading muon is required to have p T > 23 GeV, and if a trailing muon has 5 < p T < 8 GeV, the leading and subleading electrons must satisfy p T > 25 and 15 GeV; in eeµ events, if a trailing muon has p T > 8 GeV, either the leading electron must have p T > 23 GeV, or the subleading electron must have p T > 15 GeV. These requirements lead to a signal selection efficiency of 5-7% for a trilepton final state depending on the N mass. 

The events are subdivided into two ranges of leading lepton p T : 15 < p leading T < 30 GeV and 30 < p leading T < 55 GeV. The lower range has higher signal efficiencies for m N close to m W, leading to three leptons with similar p T spectra. The higher range targets very low m N down to 1 GeV, with one energetic and two soft leptons in the event. 

Finally, the variable M min 2 OS, which is correlated with m N, is used to further subdivide the events into four bins (<10, 10-20, 20-30, and >30 GeV) giving a total of eight search regions, as shown in Fig. 1. 

In the high-mass region (targeting m N > m W ), N is produced in the decay of a high-mass off-shell W boson, leading to three relatively energetic leptons, and more sizable p miss T. In this region, the three selected leptons must satisfy p T > 55, 15, 10 GeV. With these requirements, the background from Wγ * production is negligible, and events containing an OSSF lepton pair are therefore retained, but with the invariant mass of any OSSF lepton pair required to satisfy M > 5 GeV. The backgrounds from WZ and Zγ ( * ) production are respectively suppressed by vetoing events having an OSSF lepton pair with |M (or M 3 ) -M Z | < 15 GeV. Signal selection efficiency for trilepton final state reaches up to 50%. 

In order to improve the discrimination of signal from background, the high-mass region is divided into two event categories: events containing an OSSF lepton pair and events with no such pair. Both categories are divided into bins of M 3 and M min 2 OS, each further subdivided according to M T. This results in a total of 25 search regions, as shown in Fig. 1. 

To extract values of |V eN | 2 and |V µN | 2 separately, the results for both the low-and high-mass regions are obtained for events with ≥2 electrons or ≥2 muons, respectively. 

We now consider the most important sources of background and their associated systematic uncertainties. The tt and DY processes, with an additional nonprompt lepton, constitute the main background for events in the low-and high-mass regions with no OSSF lepton pair. It is estimated by using the tight-to-loose ratio method described in Ref. [79]. The probability for a loose nonprompt lepton to pass the tight selection criteria is measured in a multijet sample in data enriched in nonprompt leptons. This probability is applied to events that pass the full signal selection, but contain at least one lepton that fails the tight selection, while satisfying the loose selection requirements. The method is validated using simulation and data in control regions enriched in tt or DY+jets events. Agreement between the predicted and observed yields in the various control regions is found to be within the overall systematic uncertainty of 30% assigned to this background estimate. 

The background from WZ and Wγ * production is dominant in the high-mass region containing an OSSF lepton pair. It is obtained from simulation, with the simulated yield normalized to data in a control region formed by selecting three tight leptons with p T > 25, 15, 10 GeV, and requiring an OSSF lepton pair with invariant mass M consistent with a Z boson: |M -M Z | < 15 GeV. In addition, events are required to have p miss T > 50 GeV. The ratio of the predicted to observed WZ yield in this control region is found to be 1. 08 ± 0. 09. This ratio is used to normalize the MC generated event samples, and its associated uncertainty is propagated to the result. 

Production of ZZ events with both Z bosons decaying leptonically, and one lepton not identi-fied, results in a trilepton signature. This contribution is estimated from simulation, and the simulated yield is normalized in a control region containing four leptons that form two OSSF lepton pairs with invariant masses consistent with a Z boson. The ratio of data to simulation in the control region is found to be 1. 03 ± 0. 10. An additional uncertainty of 25% is assigned to the prediction of events with M T > 75 GeV, based on a comparison of the observed and predicted event yields in the control region. 

External and internal photon conversions (Xγ ( * ) ) contribute to the trilepton final state when a photon is produced with a Z boson, and this photon undergoes an asymmetric conversion in which one of the leptons has very low p T and fails the lepton selection criteria. This contribution is obtained from simulation and verified in a data control region enriched in conversions from the Z+jets process, with Z → γ ( * ) and γ ( * ) →, where one of the leptons is outside the detector acceptance. The control region is defined by |M -M Z | > 15 GeV and |M 3 -M Z | < 15 GeV. The ratio of data to expected background in the control region is 0. 95 ± 0. 08, and is used to normalize the simulation. Kinematic properties of the events in data are used to set a systematic uncertainty in the photon conversion background of 15%. 

Other SM processes that can yield three or more prompt leptons include triboson production (W, Z, H, or a prompt γ) and single-boson production in association with a single top quark or a tt pair (tt/t+X). Such processes generally have very small production rates and in some cases are further suppressed by the b jet veto. They are estimated from simulation with an uncertainty of 50%, which includes uncertainties due to experimental effects, event simulation, and theoretical calculations of the cross sections. 

The background from the mismeasurement of charge arises from events with an e + e -pair in which the charge of one of the electrons is misreconstructed. It is predicted using simulation, which is validated to agree within 10% of an estimate obtained from data [80]. 

Systematic uncertainties affecting any process whose yield is estimated from simulation are considered, such as those from trigger efficiency, lepton selection efficiencies, jet energy scale, b jet veto efficiency, pileup modeling, and those related to fixed-order calculations in event simulation. The effect of each uncertainty on the event yields is computed and accounted for. 

The uncertainty in the trigger efficiency is obtained by measuring the efficiencies of all trigger components using the tag-and-probe technique [77, 78], and is estimated to be 2% for events with leading lepton p T > 30 GeV and 5% otherwise. Lepton identification efficiencies are also computed using the tag-and-probe technique with an uncertainty of 2% per lepton. 

The impact of the jet energy scale uncertainty is determined by shifting the jet energy correction factors up and down by their estimated uncertainty for each jet, and recalculating all kinematic quantities obtained from jets. This results in an uncertainty in event yields of up to 3%, depending on the search region. Correlation effects due to the migration of events from one search region to another are also taken into account. Similarly, the b jet veto efficiency is corrected for differences between data and simulation, leading to an uncertainty in event yields of 1-5%. The uncertainty in yields due to modeling of pileup is computed by modifying the total inelastic scattering cross section by 5% [81], and is measured to be 1-5%, depending on the search region. The uncertainty in the integrated luminosity is 2. 5% [82]. 

Further uncertainties in background yields estimated from simulation arise from the unknown higher-order effects in the theoretical calculations of cross sections, and from uncertainties in the knowledge of the proton PDFs. Uncertainties in the renormalization and factorization scales affect the signal cross section and acceptance. These are evaluated by independently varying the aforementioned scales up and down by a factor of two relative to their nominal values. The uncertainties from the PDFs are estimated using the replica weights of the nominal PDF sets [83]. 

The limited statistical precision of the available MC samples leads to an additional uncertainty of 1-30%, depending on the process and search region. The expected and observed yields together with the relative contributions of the different background sources in each search region, are shown in Fig. 1. Tabulated results and enlarged ver-sions of Fig. 1, with potential signals superimposed, are provided in Appendix A. We see no evidence for a significant excess in data beyond the expected SM background. We compute 95% confidence level (CL) upper limits on |V eN | 2 and |V µN | 2 separately, while assuming other matrix elements to be 0, using the CL s criterion [84, 85] under the asymptotic approximation for the test statistic [86, 87]. A simultaneous fit of all search regions is performed and all systematic uncertainties are treated as log-normal nuisance parameters in the fit. 

The interpretation of the results is presented in Fig. 2. The N lifetime is inversely proportional to m 5 N |V N | 2 [45, 50]. At low masses this becomes significant, resulting in displaced decays and lower efficiency than if the decays were prompt, illustrated by comparison of the black dotted line in Fig. 2 (prompt assumption) with the final result. This is accounted for by calculating the efficiency vs. N lifetime, and propagating this to the limits on mixing parameter vs. mass. The dashed black curve is the expected upper limit, with one and two standard-deviation bands shown in dark green and light yellow, respectively. The solid black curve is the observed upper limit, while the dotted black curve is the observed limit in the approximation of prompt N decays. Also shown are the best upper limits at 95% CL from other collider searches in L3 [36], DELPHI [33], ATLAS [23], and CMS [22]. 

In summary, a search has been performed for a heavy neutral lepton N of Majorana nature produced in the decays of a W boson, with subsequent prompt decays of N to W, where the vector boson decays to ν. The event signature consists of three charged leptons in any combination of electrons and muons. No statistically significant excess of events over the expected standard model background is observed. 

Upper limits at 95% confidence level are set on the mixing parameters |V eN | 2 and |V µN | 2, ranging between 1. 2 × 10 -5 and 1. 8 for N masses in the range 1 GeV < m N < 1. 2 TeV. These results surpass those obtained in previous searches carried out by the ATLAS [23] and CMS [22, 24] Collaborations, and are the first direct limits for m N > 500 GeV. This search also provides the first probes for low masses (m N < 40 GeV) at the LHC, improving on the limits set previously by the L3 [29] and DELPHI [33] Collaborations. For N masses below 3 GeV, the most stringent limits to date are obtained from the beam-dump experiments: CHARM [26, 31], BEBC [25], FMMF [32], and NuTeV [34]. 

Grid for delivering so effectively the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agencies: BMWFW and FWF (Austria); FNRS and FWO (Belgium); CNPq, CAPES, FAPERJ, and FAPESP (Brazil); MES (Bulgaria search region
section_name: 1

section_name: Acknowledgments

section_content: We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centers and personnel of the Worldwide LHC Computing
section_name: Acknowledgments