Search for Majorana neutrinos in same-sign W W scattering events from pp collisions at √ s = 13 TeV

A search for Majorana neutrinos in same-sign W W scattering events is presented. The analysis uses √ s = 13 TeV proton–proton collision data with an integrated lumi-nosityof140fb − 1 recordedduring2015–2018bytheATLAS detector at the Large Hadron Collider. The analysis targets ﬁnal states including exactly two same-sign muons and at least two hadronic jets well separated in rapidity. The modelling of the main backgrounds, from Standard Model same-sign W W scattering and W Z production, is constrained with data in dedicated signal-depleted control regions. The distribution of the transverse momentum of the second-hardest muon is used to search for signals originating from a heavy Majorana neutrino with a mass between 50 GeV and 20 TeV. No signiﬁcant excess is observed over the background expectation. The results are interpreted in a benchmark scenario of the Phenomenological Type-I Seesaw model. In addition, the sensitivity to the Weinberg operator is investigated. Upper limits at the 95% conﬁdence level are placed on the squared muon-neutrino–heavy-neutrino mass-mixing matrix element | V μ N | 2 as a function of the heavy Majorana neutrino’s mass m N , and on the effective μμ Majorana neutrino mass | m μμ | .


Introduction
The Higgs field's Yukawa couplings explain in a satisfactory way, and in agreement with data, the generation of the charged-fermion masses.However, the precise origin of tiny but non-zero neutrino masses, brought to light by neutrino oscillation experiments [1,2], is not explained by the Standard Model (SM) of particle physics and remains a mystery today.Among the most compelling solutions to this mystery is the hypothesised existence of Majorana masses for neutrinos, which can be explored by high-and low-energy experiments [3][4][5][6].
Majorana neutrino masses are present in beyond-the-SM (BSM) theories that include new, heavy leptons or extended scalar sectors, collectively known as Seesaw models, or in theories with extended gauge sectors (Left-Right Symmetric Model, Grand Unified Theories) [7].Even if these new particles have masses that are much higher than the energies reached at the Large Hadron Collider (LHC), and are therefore kinematically inaccessible, the Majorana nature of active neutrinos can still be probed indirectly using effective field theories, with the lowest-dimensional operator being the so-called dimension  = 5 Weinberg operator [8].
At energies below 100 MeV, the Majorana nature of neutrinos can be probed by searching for neutrinoless double beta decays of certain isotopes, which are characterised by the appearance of two same-sign electrons () with an absence of neutrinos () in the final state [9].Although the production of same-sign leptons (ℓ) involving muons () or -leptons is kinematically forbidden in such isotope decays, the same process can be studied at high energies when mediated at dimension  = 5, 7, or higher [10][11][12], as shown in Figure 1(a).This is the subject of the present work, which searches for signs of heavy neutral leptons in the same-sign dimuon final state with the ATLAS detector.This search also connects to existing LHC searches for heavy Majorana neutrinos in resonant production [13][14][15][16], with the goal of extending the reach for heavy neutrino masses to beyond the TeV scale [10], where the resonant production channels become kinematically inaccessible.A similar search was reported recently by the CMS Collaboration [17].Searches for heavy Majorana neutrinos lighter than the  boson were conducted at LEP [18].
Figure 1: Diagrammatic representation of (a) same-sign  +  + production in  +  + scattering mediated by a Majorana neutrino  in proton-proton collisions and (b) electroweak same-sign  +  + scattering, which is the main background in this search.The corresponding diagrams with negative-charged particles are also relevant to this work but not included here for simplicity.
This search is conducted in the context of two benchmark models.The first is the Phenomenological Type-I Seesaw model [3,19], where heavy Majorana neutrinos couple to SM particles through mass-mixing with light neutrinos, which is parameterised with a complex-valued matrix element.For simplicity, this work investigates only the phenomenology of the lightest heavy-neutrino mass eigenstate  by assuming that all other heavy eigenstates are decoupled.The benchmark scenario is considered without  and  flavour mixing, so the only non-zero entry in the mass-mixing matrix is the element |   | 2 .In this case the signal cross section (   →  ±  ± + ) can be written as a function of the cross section  0 (   →  ±  ± + ) for the flavour-aligned scenario (in which |   | = 1): The two free parameters of the model are |   | 2 and the mass   of the heavy Majorana neutrino.
The second benchmark model is the  = 5 Weinberg operator [8], which extends the SM Lagrangian by: where Λ is the scale at which the particles mediating the BSM signal become relevant degrees of freedom;  ℓℓ ′ 5 is a complex, flavour-dependent Wilson matrix coefficient;   ℓ = ( ℓ , ℓ) is the left-handed lepton doublet; and Φ is the SM Higgs doublet.Considering the Higgs vacuum expectation value  =

√
2⟨Φ⟩ ≈ 246 GeV, it is possible to define the quantity which for ℓℓ ′ =  reduces to the effective Majorana mass usually quoted by neutrinoless double beta decay experiments, which is a measure of the scale at which lepton-number symmetry is broken [9,11].For simplicity, it is assumed here that   5 is the only non-zero Wilson coefficient of this matrix.The signal cross section (   →  ±  ± + ) can then be written as a function of the dimensionless effective cross section  5 (   →  ±  ± + ) for Λ/|  5 | = 1 TeV as follows: The analysis targets final states with two same-sign muons, and also two particle jets with a large angular separation to capture the signature of vector-boson scattering.The main backgrounds are SM  ±  ±   production, shown in Figure 1(b), and   production, where the bosons decay leptonically.They are estimated with simulated events and their modelling is constrained in dedicated signal-depleted control regions.A data-driven technique assisted by simulation is used to model backgrounds which do not originate from prompt muons.The presence of a signal is tested in a region designed to be signal-enriched using the benchmark models discussed above.

The ATLAS detector
The ATLAS detector [20] at the LHC is a multipurpose particle detector with a forward-backward symmetric cylindrical geometry with a nearly 4 coverage in solid angle. 1 It consists of an inner tracking detector (ID) surrounded by a superconducting solenoid, electromagnetic (EM) and hadron calorimeters, and a muon spectrometer (MS).
The ID covers the pseudorapidity range || < 2.5.The high-granularity silicon pixel detector covers the vertex region and typically provides four measurements per track, the first hit normally being in the insertable B-layer installed before Run 2 [21,22].It is surrounded by a silicon microstrip detector and a straw tube transition-radiation tracking detector.The field integral of the toroids ranges between 2.0 and 6.0 T m across most of the detector.The MS includes a system of precision chambers for tracking and fast detectors for triggering.The ATLAS trigger and data acquisition system [23] contains a two-level trigger system in order to select events of interest.The first-level trigger is implemented in hardware and uses a subset of the detector information to accept events at a maximum rate of nearly 100 kHz.The second-level is a software-based trigger that reduces the accepted event rate to 1 kHz, on average, depending on the data-taking conditions [23].An extensive software suite [24] is used in data simulation, in the reconstruction and analysis of real and simulated data, in detector operations, and in the trigger and data acquisition systems of the experiment.

Object reconstruction and event selection
This analysis uses a set of   collision data events collected by the ATLAS detector between 2015 and 2018 at a centre-of-mass energy of √  = 13 TeV.Only events for which all detector subsystems were operational are considered.The data set corresponds to an integrated luminosity of 140 fb −1 [25,26].This analysis is based on events where the detector readout is triggered by the presence of at least one muon, referred to as 'single-muon triggers'.Single-muon triggers require the presence of a single muon with a minimum transverse momentum ( T ) in the range 20-26 GeV depending on the data-taking period [27].
Events are required to have at least one collision vertex reconstructed from at least two ID tracks with  T > 500 MeV.For events with several collision vertices, the one with the largest sum of the squared transverse momenta of the associated tracks is taken as the primary vertex [28].
Three types of muon candidates are used to select events for the signal region, to estimate the background contribution from fake muons, and to veto additional muons in the event selection.These muon candidates have different selection criteria and are referred to as 'signal muons', 'spurious muons', and 'baseline muons', respectively.Baseline muon candidates are reconstructed by combining tracks in the ID with tracks in the MS [29] and are required to have || < 2.5 and  T > 3 GeV.Further selections on the longitudinal and transverse impact parameters are imposed.The transverse impact parameter  0 , measured from the beamline and divided by its estimated uncertainty, must satisfy | 0 |/( 0 ) < 15.The longitudinal impact parameter  0 , measured from the primary vertex, is required to satisfy | 0 sin()| < 1.5 mm.Baseline muon candidates must satisfy the 'Loose' cut-based identification working point defined in Ref. [30], based on the numbers of hits in the different ID and MS subsystems, and on the significance of the charge-to-momentum ratio /.To select signal muons, these selection criteria are tightened to  T > 10 GeV, || < 2.5, | 0 |/( 0 ) < 3 and | 0 sin()| < 0.5 mm, and used together with the 'Medium' cut-based identification working point.In addition, an isolation is required where the scalar sum of the  T of tracks inside a variable-size cone around the muon (excluding its own track) plus 40% of the transverse energy of neutral particle-flow objects in a cone of size Δ = 0.2 around the muon, must be less than 4.5% of the muon  T [30].To select spurious muons, the above isolation requirement is loosened to be less than 16% of the muon  T and | 0 |/( 0 ) < 10.In addition, they must fail the tighter signal-muon isolation or transverse impact parameter requirements.
Electron candidates are reconstructed from an isolated EM calorimeter energy deposit which is matched to a track in the ID [31] and they are required to have  T > 4.5 GeV.The pseudorapidity of the calorimeter energy cluster,  cluster , must satisfy | cluster | < 2.47.Electron candidates must pass a 'Loose' likelihoodbased identification working point [31] that employs calorimeter and tracking information to discriminate between electrons and jets, and their tracks must have impact parameters satisfying | 0 |/( 0 ) < 5 and | 0 sin()| < 0.5 mm.In this analysis, electron candidates are used only to veto events.Jets are reconstructed by using the anti-  algorithm [32,33] with a radius parameter of  = 0.4 after applying a particle-flow algorithm [34] to tracking information and topological clusters of calorimeter-cell energies.Jets are required to have  T > 25 GeV and || < 4.5, and are calibrated as described in Ref. [35].The jet-vertex tagger (JVT) discriminant [36] is used to identify jets originating from the hard-scatter interaction through the use of tracking and vertexing information.Jets with  T < 60 GeV and || < 2.5 are required to satisfy the requirement JVT > 0.50, corresponding to a selection efficiency of about 96% for hard-scatter jets [37].If either of the two leading jets has || > 2.5 and  T < 60 GeV, it must satisfy a threshold for the forward jet-vertex tagger (fJVT) discriminant [38,39] that has an efficiency of about 93% for hard-scatter jets.
In order to suppress contributions from background processes that involve top quarks or leptonic -hadron decays, the DL1r classification algorithm based on recurrent neural networks [40,41] is used to identify jets containing -hadrons (referred to as '-jets').This algorithm identifies -jets amongst a background of light-flavour-and charm-quark-initiated jets using information about the impact parameters of tracks associated with the jet and the topological properties of displaced vertices reconstructed within the jet.This analysis uses a -jet efficiency working point of 85% as measured for jets with  T > 20 GeV and || < 2.5 in simulated  t events [42].The tagging algorithm at that working point gives an expected rejection factor (defined as the inverse of the efficiency) of about 40 for light-flavour jets, and about 2.9 for jets originating from charm quarks [41,43].
To resolve the potential ambiguities due to a single detector response being assigned to two objects by the reconstruction algorithm, a sequential overlap-removal procedure is applied.First, any electron found to share a track with another electron with higher  T is removed.Second, electrons sharing their track with a muon candidate are removed.Third, labelling the rapidity of an object by ,2 any jet within 2 of an electron is removed, unless it is a -jet and the electron  T is below 100 GeV, in which case the electron is removed.Fourth, any electrons within Δ  = 0.4 of a remaining jet are removed.Fifth, any jet that is within Δ  = 0.2 of a muon and has less than three associated tracks is removed, unless it is a -jet, in which case the muon is removed.Finally, remaining muons are removed if their track is within Δ  = 0.4 of a remaining jet.Resolving the overlaps between / and jets in two steps based on Δ  values of 0.2 and 0.4 allows the reconstruction efficiency to be maximised for -jets including semileptonic -and -hadron decays.
The missing transverse momentum, with magnitude  miss T , is defined as the negative vectorial sum of the transverse momenta of all selected and calibrated electrons, muons and jets, and an additional term to account for soft energy in the event that is not associated with any of the selected objects.This soft term is calculated from unassociated ID tracks matched to the primary vertex.While it suffers from an inherent underestimation, which is covered by dedicated uncertainties, it is more resilient to contamination from multiple interactions in the same or neighbouring bunch crossings (pile-up) compared to the calorimeter-based soft term [44,45].The object-based  miss T significance (S) is computed [46] by including the expected resolutions for all the objects used in the  miss T calculation and is used to reject background events with genuine  miss T originating from undetected particles such as neutrinos.
At preselection level, the events are required to have at least two same-sign signal muons with  T > 27 GeV and one of them has to match the muon that fired the trigger.At least two jets are required and one of the jets has to have a transverse momentum larger than 30 GeV to reduce the contamination from pile-up jets.Events with -jets are vetoed.In order to preferentially capture the signature of  ±  ± scattering in   collisions, the invariant mass of the two jets with the highest  T , labelled as    , has to be larger than 300 GeV and their rapidity separation must satisfy |Δ   | > 4.These preselection cuts are applied to every event which is used further in the analysis.

Event simulation
Monte Carlo (MC) samples are used to model the expected signal and background yields, while backgrounds from misidentified objects are estimated using data-driven techniques described in Section 5.All MC samples are normalised to the best theory predictions available.The NNPDF3.0nlo [47] PDF set was used in all matrix element calculations if not stated otherwise.The generated event samples were processed through the full ATLAS detector simulation [48] based on Geant4 [49].Additional simulated   collisions generated using Pythia 8.186 [50] with the A3 set of tuned parameters (tune) [51] and the MSTW2008lo parton distribution function (PDF) set [52] were overlaid to model the effects of pile-up.The distribution of the number of additional   interactions in the MC samples is reweighted to match the one observed in data.All simulated samples are processed through the same reconstruction algorithms and analysis chain as the data.
The Phenomenological Type-I Seesaw model signal samples were simulated following the prescription given in Ref. [10].The samples were simulated using the MadGraph5_aMC@NLO 2.7.2 [53] generator at next-to-leading order (NLO) in QCD after importing the default variant of the HeavyN [54,55] FeynRules UFO libraries [56,57].All HeavyN events were interfaced to Pythia 8.243 [58].Twenty mass points were generated, spanning the range   = 50 GeV to 20 TeV.Since the mixing parameter    only impacts the total cross section, as  ∼ |   | 4 , these mass points were generated with a mixing parameter of unity and are rescaled to other values by using Eq. ( 1).
The signal sample for the Weinberg operator was simulated by using the prescription given in Ref. [11].For the hard parton-level scattering processes, MadGraph5_aMC@NLO 2.9.3 [53] was used at NLO in QCD after importing the default variant of the SMWeinberg [11] FeynRules UFO libraries [56,57].The NNPDF3.1luxQED [47] PDF was used.All SMWeinberg events were interfaced to Pythia 8.245 [58].The sample was generated at an energy scale Λ = 200 TeV and with the Wilson coefficient value  For all samples of simulated signal events, the parton shower, hadronisation, and underlying-event modelling by Pythia used the A14 tune [59] and the NNPDF2.3lo[47] PDF set, and the EvtGen 1.7.0 program [60] was used to model the decays of bottom and charm hadrons.
Electroweak (EW)  ±  ± production and decay in association with two jets was simulated with diagrams including exactly six orders of the EW coupling [61].The simulation of the strong production processes, and the interference between EW and strong contributions, includes diagrams with exactly four and five EW vertices [62], respectively.All SM  ±  ±   contributions were simulated with Mad-Graph5_aMC@NLO 2.6.7 [53], which was interfaced to the Pythia 8.244 [58] parton shower model.The matching of the matrix element simulation to parton showers used the MC@NLO method [63] based on FKS subtraction [64].These samples used the A14 set of tunable parameters [59].More details of the generator settings used to obtain a better description of deep-inelastic scattering data [65] are given in Ref. [61].
Samples of diboson final states (,  = /) were simulated with the Sherpa 2.2.1 or 2.2.2 [66] generator depending on the process.Fully leptonic final states and semileptonic final states, where one boson decays leptonically and the other hadronically, were generated using matrix elements with NLO accuracy in QCD for up to one additional parton and with LO accuracy for up to three additional parton emissions.Samples for the loop-induced processes  →  were generated using LO matrix elements for up to one additional parton emission.The production of triboson () events was simulated with the Sherpa 2.2.2 generator, accurate to NLO in QCD for the inclusive process and to LO for up to two additional parton emissions, using factorised gauge-boson decays.The matrix element calculations for individual  and  processes were matched and merged with the Sherpa parton shower based on Catani-Seymour dipole factorisation [67,68] using the MEPS@NLO prescription [69][70][71][72].The virtual QCD corrections for matrix elements with NLO accuracy were provided by the OpenLoops 1 library [73][74][75].
Additional samples are used to simulate minor backgrounds with prompt muons.The production of  t,  t,  and   events was modelled using the MadGraph5_aMC@NLO 2.3.3 generator at NLO, interfaced to Pythia [58].The Sherpa 2.2.11 generator was used to simulate   final states.For studies of background which does not originate from prompt muons, the production of  t and single-top-quark events was simulated using the Powheg Box v2 [76][77][78][79] generator at NLO in QCD, interfaced to Pythia 8.230 [58], and the production of +jets events was simulated with the Sherpa 2.2.11 generator.

Analysis strategy
The signal considered in this search is characterised by two same-sign muons with large transverse momentum and two jets with a large rapidity separation, which is already exploited in the preselection requirements.The signal region (SR) was optimised for the Phenomenological Type-I Seesaw model with Majorana neutrino masses above 1 TeV, since this is the kinematic region where this search is expected to T ).The distribution has a maximum above 200 GeV for   = 500 GeV, with a long tail towards higher momenta [10], while all SM backgrounds have a steeply falling   2 T distribution.The main background contributions arise from SM  ±  ±   and   production.The SR is defined by requiring the  miss T to have a small significance value S < 4.5, and events with a third baseline muon or electron are vetoed.This suppresses background from SM  ±  ±   events, which have two neutrinos in the final state (cf. Figure 1(b)), and   events, which have a third lepton.In contrast, these two requirements are inverted to define signal-depleted control regions (CR) where the data are used to constrain the normalisation of these simulated background processes.The event selection for the SR and the two CRs is summarised in Table 1.
A same-sign  ±  ±   control region (ssWW-CR) is selected with S > 5.8.In order to suppress possible contamination from signal events, T is restricted to be below 120 GeV.This event selection also ensures that background from misidentified jets is adequately suppressed.The ssWW-CR has a SM  ±  ±   purity of 70% inclusively (64% for 80 <   2 T < 120 GeV), with the main contamination coming from   events, whose modelling is, however, constrained in another dedicated control region.Events which fulfil the ssWW-CR requirement but with a  miss T significance 4.5 < S < 5.8 or with   2 T > 120 GeV could have a non-negligible signal contribution.These two regions are not included in the profile likelihood fit, but are used to validate the modelling of SM  ±  ±   and   events.The final background predictions and the collision data agree well in both validation regions.
The second-largest background in the SR is   production where both bosons decay leptonically and the -boson-decay lepton opposite in sign to the  boson falls outside the detector acceptance (mainly due to the muon spectrometer's rapidity coverage).To select a control region enriched in   production (WZ-CR), the SR is modified such that events are required to have three signal muons, with one opposite-sign (OS) muon pair to be compatible with a  boson decay.Events with an additional baseline muon are vetoed to reduce the contribution from   production.Furthermore, a large three-lepton invariant mass ( ℓℓℓ > 100 GeV) is required in order to reduce the fake-lepton contribution and to make the WZ-CR orthogonal to the regions where fake leptons from photons are studied, as explained later in this section.The WZ-CR has a   purity of 88% inclusively and 83% for   2 T > 120 GeV.Non-prompt muons, defined as originating from bottom-or charm-hadron, kaon or pion decays, constitute a non-negligible source of background in the analysis.Most of these muons arise from +jets events and  t events where one  boson (originating from a top-quark decay in the  t case) decays leptonically.This background contribution is evaluated using the so-called fake-factor method [80], which is based on collision data.Fake factors are derived as a function of muon  T in a dedicated region enriched in non-prompt muons which selects events that include a jet recoiling almost back-to-back against a fake muon candidate.In particular, the presence of +jets events in this sample leads to a non-negligible contamination from prompt muons.This contamination is reduced by requiring the transverse mass formed by the non-prompt muon candidate's  T and the missing transverse momentum to be below 50 GeV.The residual prompt-muon contribution, estimated from MC predictions, is subtracted from the data, and the fake factors are calculated as the ratio of non-prompt 'signal' to 'spurious' muon events.The background from non-prompt muons is then estimated by multiplying the fake factors by the data yields observed in the so-called 'spurious-SR', which has the exact same requirements as the SR except that the object definition of one of the muons is changed from 'signal' to 'spurious'.Non-prompt contributions originating from either the leading or subleading muon candidates are considered.The contamination from prompt muons in the spurious-SR, originating mostly from same-sign  and   events, is subtracted using MC predictions.The final data-driven estimate gives a low non-prompt-muon event yield in the different analysis regions, and hence this contribution to the   2 T distribution is smoothed by fitting it with an exponential function.Alternative functions are used to estimate the systematic uncertainty due to the functional form choice.
Background from  t,  t,  , ,   and   events is estimated from MC simulation and denoted as 'Other' backgrounds in the remainder of the paper, in which   and  constitute the dominant contributions.
Photons originating from final-state radiation splitting into muons ( * → ) were studied in a region selecting leptonic  boson decays produced in association with an extra muon (which is the non-prompt candidate from final-state radiation) with  ℓℓℓ < 100 GeV to suppress the events with initial-state radiation.Their contribution to the SR was found to be negligible.Muons with the wrong reconstructed charge were studied with MC simulation and were also found to be negligible.The probability of misreconstructing two inelastic collisions, each contributing one of a pair of same-sign  bosons, as a single event was evaluated.A Gaussian function with a standard deviation as measured in Ref. [81] was used to model the longitudinal distribution of inelastic collisions and to compute the conditional probability.Less than 0.1 events from this source are expected to pass the analysis selection.Finally, single   interactions where more than one set of partons from each proton interact (double parton scattering) to produce two  bosons with the same sign were considered.Although produced at a rate several times larger than that of coincidental  events, less than one such event is expected to pass the analysis selction, so it is also safe to ignore this process.

Results
This search is conducted using the binned   2 T distributions.A combined profile likelihood fit is used in the SR and all CRs to obtain the final background model, as shown in Figure 2, and make further statistical inferences regarding the presence of a signal.The likelihood function depends on the free-floating signal normalisation, the free-floating SM  ±  ±   and   normalisations, and a set of nuisance parameters encoding the effect of systematic uncertainties.The normalisation corrections obtained by the fit (defined as the ratio between post-fit and pre-fit predictions) are 1.17 ± 0.23 for SM  ±  ±   and 0.92 ± 0.16 for   backgrounds.The nuisance parameters, implemented as Gaussian constraints, allow variations of the expectations for signal and background according to the corresponding systematic uncertainties, whilst penalising the likelihood function for any deviation from the nominal expectations according to the specified constraints.The fit also reduces the impact of systematic uncertainties on the search sensitivity by exploiting the background-dominated CRs.The fit and the subsequent statistical analysis are performed for each signal hypothesis.
Uncertainties affecting the measurement which originate from statistical sources are considered together with systematic uncertainties related to the detector calibration and physics modelling.Uncertainty sources related to muon reconstruction and energy calibration, jet energy calibration, -jet identification, rejection of pile-up jets, and missing transverse momentum reconstruction were studied.The uncertainty in the non-prompt muon background is evaluated by propagating both the systematic and statistical uncertainties in the fake factors and in the spurious-SR non-prompt muon yields.The systematic uncertainty is dominated by the modelling of the contamination from prompt muons in the region used to extract the fake factors.The statistical uncertainty, which dominates the uncertainty in the non-prompt background, is driven by the limited size of the spurious-SR data set.Systematic uncertainties in the physics modelling of SM same-sign  and   events are estimated by varying the factorisation ( f ) and renormalisation ( r ) scales, as well as the strong coupling constant  s and the PDF choice.Factorisation and renormalisation scale variations are used to account for the uncertainty due to truncated higher-order corrections.The effects of a seven-point pairwise variation of  f and  r , where each scale is either halved or doubled (the off-diagonal variations by [0.5, 2] and [2, 0.5] are not included), are computed and then combined into a single uncertainty by forming an envelope of the effects of these seven variations.The total PDF uncertainty is obtained from the standard deviation of NNPDF3.0nlovariations and the envelope of alternative PDFs (CT14 and MMHT2014), while the effect of the  s uncertainty is obtained by varying the nominal  s = 0.118 within its uncertainty of ±0.001.Alternative SM same-sign  ( ) samples generated with Sherpa 2.2.12 and MadGraph5_aMC@NLO interfaced to Herwig 7.2 [82] (MadGraph5_aMC@NLO 2.2.2 interfaced to Pythia 8.186) are used to evaluate the residual generator dependencies.A similar scheme is applied to derive uncertainties from the  f / r scale variations and from the PDF/ s choice for the HeavyN and SMWeinberg signal samples, while other smaller backgrounds are assigned a conservative ±50% uncertainty correlated across analysis bins.The largest systematic uncertainty in the total background estimate is the uncertainty in the non-prompt-muon background estimation (8%).None of the considered systematic uncertainties were pulled appreciably by the fit or found to have an impact on the analysis sensitivity, which is driven by the statistical uncertainty of the data in the tail of the   2 T distribution.The expected limits change at the per-mille level if systematic uncertainties are not considered.
In the absence of any significant excess above the background expectation, upper limits on |   | 2 and |  5 |/Λ are derived with the CL s method [83,84].The asymptotic approximation [85] is used to derive the upper limits, which agree very well with those estimated via background-only pseudo-experiments.The observed and expected 95% confidence level (CL) upper limits on the heavy Majorana neutrino mixing element |   | 2 are shown in Figure 3 as a function of   .The form of the exclusion curve follows what is expected from the bare cross section  0 as defined in Eq. ( 1), which rises moderately with increasing   , reaching a maximum around 700 GeV for 13 TeV   collisions before decreasing again due to a kinematic suppression in the  ±  ± matrix element [10].The strongest exclusion is achieved at   = 500 GeV, with an observed (expected) upper limit of 0.10 (0.08) on |   | 2 .For a heavy neutrino mass of 20 TeV the exclusion curve reaches the mixing element's upper bound of unity.Overall, the expected limits constitute a sensitivity improvement of 15%-20% compared to the results in Ref. [17].Over the entire mass range the observed limit is around one standard deviation weaker than the expected limit, which is due to the slight excess of data in the last two bins of the

Conclusion
A new search for heavy Majorana neutrinos is presented, targeting the same-sign  production mode.Since in this process the heavy neutrinos are produced in the -channel, this analysis extends the reach to higher neutrino masses and complements previous searches in resonant production channels by adding valuable sensitivity for neutrino masses of a few hundred GeV.The search is performed using 140 fb −1 of   collision data at √  = 13 TeV collected by the ATLAS detector at the LHC.To preferentially capture the signature of same-sign  scattering in   collisions, events are required to have two jets with a large rapidity separation and two muons with the same electric charge.Since signal events typically have a significantly higher muon transverse momentum than those arising from SM backgrounds, the shape of the transverse momentum distribution of the subleading muon is used to discriminate between the signal and background.No significant excess was observed above the expected backgrounds, which were modelled in dedicated control regions.Upper limits on the muon-neutrino-heavy-neutrino mass-mixing matrix element |   | 2 at 95% confidence level are computed for signals with Majorana neutrino masses   between 50 GeV and 20 TeV based on a Phenomenological Type-I Seesaw model.The strongest exclusion was achieved at   = 500 GeV, with an observed (expected) limit of 0.10 (0.08) on |   | 2 .Constraints on |   | 2 are extended from previous results at   values of the order of 1 TeV in the resonant production channels to   values of up to 20 TeV in this analysis.In addition, limits are set on the  = 5 Weinberg operator, which translate to an observed (expected) upper limit of 16.7 (13.1)GeV on the effective  Majorana neutrino mass, which cannot be probed in nuclear decays.

𝜇𝜇 5 = 1 .
Both Λ and   5 only impact the total cross section, as  ∼ (  5 /Λ) 2 , so a global rescaling of a single signal sample according to Eq. (3) is used to vary their values.

Figure 2 : 2 T
Figure 2: Comparison of the   2T distributions in data and the post-fit prediction under the background-only hypothesis in the different analysis regions: (a) SR, (b) ssWW-CR, and (c) WZ-CR.The shaded area surrounding the background expectation represents the total uncertainty from the post-fit computation.The ratio of the observed yields to the prediction of the background model post-fit is shown by the points in the bottom plot and the change in the prediction relative to the pre-fit yields is shown by the blue line.The 'Other' contribution consists mainly of   and  events.Two benchmark signal samples are overlaid on top of the total background prediction in (a), and were obtained by using their theoretical cross sections.The last bin in (a) and (c) includes the overflow bin.

𝜇 2 T 5 |
distribution.The observed (expected) 95% CL lower limit on the  = 5 Weinberg operator Λ/|  is 3.635 (4.608) TeV, with a probability of  0 = 9.3% that the observed data are compatible with the background-only hypothesis.Using Eq. (2), this translates into an upper limit of 16.7 (13.1)GeV on the effective  Majorana neutrino mass.

Figure 3 :
Figure 3: Observed and expected 95% CL upper limits on the heavy Majorana neutrino mixing element |   | 2 as a function of   in the Phenomenological Type-I Seesaw model.The one and two standard deviation bands of the expected limit are indicated in green and yellow, respectively.

Table 1 :
Summary of the main kinematic requirements for the different selection regions.The first four requirements belong to the preselection stage and are described in Section 3.