Fatjet signatures of heavy neutrinos and heavy leptons in a left-right model with universal seesaw at the HL-LHC

We perform a collider search for fatjet signals originating from boosted heavy neutral and charged leptons with masses between a few hundred GeV to a TeV. These heavy leptons originate from the decay of heavy gauge bosons with masses above 4 TeV in a left-right symmetric extension of the Standard Model (SM), which considers a universal seesaw mechanism for the generation of all the SM fermion masses. The fatjet signals arise naturally in this model due to the presence of heavy seesaw partners of the SM fermions which decay to SM gauge bosons carrying large boosts. We employ substructure based variables lepton sub-jet fraction ($LSF$) and lepton mass drop ($LMD$) together with kinematic variables of fatjets to look for fatjet signals associated with non-isolated leptons. These variables help in reducing the SM backgrounds while retaining enough statistics for signal events, which leads to a robust discovery potential at the high-luminosity Large Hadron Collider (HL-LHC).


Introduction
Non zero neutrino mass observed in experiments [1][2][3][4][5] is one of the key hints to look for physics beyond the Standard Model (SM).Some other aspects unexplained in the SM are dark matter, hierarchy of masses in the three known fermion generations, gauge hierarchy problem, Baryogenesis, CP-violation, a email:atridey@hri.res.inb email:rafiqulrahaman@hri.res.inc email:skrai@hri.res.inetc. Left-right symmetric extension of the SM (LRSM) [6][7][8][9] is one of the well motivated models of new physics, which offers an explanation for some of the experimental as well as theoretical limitations of the SM mentioned above.Leftright symmetric models resolve the issue of maximal parity violation in the weak sector and naturally accommodate righthanded neutrinos in the framework which leads to the popular seesaw mechanism for neutrino mass generation [10][11][12][13].The seesaw framework gives rise to Majorana masses to the neutrinos and the heavy neutrinos have lepton-number violating interactions.These heavy states can decay via their CP-violating Yukawa interactions to generate a lepton asymmetry.This lepton asymmetry is then partially converted to a baryon asymmetry through the SM sphaleron processes that can explain the matter-antimatter asymmetry [14][15][16][17].The LRSM models also account for CP-violation and resolve the strong CP problem [18][19][20][21][22][23].
The minimal left-right models include scalar triplets and a scalar bi-doublet along with right-handed SM fermions arranged in the SU(2) R as doublets [24].This leads to the SM neutrinos getting their eventual mass via the seesaw mechanism, while the rest of the SM fermions get mass through their Yukawa interactions with the bidoublet scalar in the usual manner.We adopt an LRSM framework that suggests that all SM fermions, including the neutrinos, get their mass from a seesaw mechanism similar to that of the neutrino [25].This can be achieved by modifying the scalar sector of LRSM with four SU(2) doublet scalars, one each for the lepton and quark doublets in the left sector as well as in the right sector [25,26].

arXiv:2207.06857v3 [hep-ph] 8 Feb 2024
A big advantage of generating the fermion masses in this manner is to prevent a highly hierarchical Yukawa structure, like in the SM.The fermion sector, however, also gets modified and heavy singlet charged fermions along with heavy singlet Majorana neutrinos need to be included to achieve a universal seesaw mechanism for the generation of all the fermion masses.
In such models, even the strong CP problem is resolved without an axion if a discrete parity symmetry is imposed [27].
This framework offers rich and interesting phenomenology whose signals can be observed at current and future colliders.An interesting signal for the model will be the observation of fatjet signatures of heavy neutrinos (ν i ) and exotic heavy charged leptons (E i ) with masses around a few hundred GeV to TeV.This signal originates from the decay of the heavy right-handed charged gauge boson (W R ) with mass above a few TeV.We note that fatjet signatures of heavy neutrinos have been studied in the literature for LRSM models with minimal scalar sector that generate seesaw masses for neutrinos [28][29][30] at the Large Hadron Collider (LHC).However, in the universal seesaw model, the presence of heavy charged fermions along with the heavy Majorana neutrinos * lead to an even more interesting picture of final states with multiple fatjets and different lepton charge multiplicities.The signal originates from the production of the W R via pp collision, which then decays to a heavy neutrino and a heavy charged lepton.
The heavy lepton further decays to a heavy neutrino and a pair of jets via an off-shell W R .Each of the heavy neutrinos (dominantly right-handed) further decay to a lepton and a pair of jets through off-shell W R .The heavy neutrinos being much lighter than the W R are produced with a substantial boost, and their decay products (pair of jets including a charged lepton) form boosted fatjets.The emerging signal becomes a topology of three fatjets where two of them include a lepton and one is characterized by two sub-jets.This three-fatjet signal results from the following sub-process in the model, This LRSM framework also offers two-fatjet signals including a lepton, originating from the decay of the heavy right-handed * We shall refer to these heavy Majorana neutrinos as "heavy neutrinos" in the remainder of the text.
neutral gauge boson (Z R ): Interestingly, a possibility of a more exotic signal including four fatjets also exists because of the production of the heavy charged leptons E i produced via A complete account of the possible signals is given in section 4. We note that the presence of leptons in the fatjet can be a crucial identification since the leptons and jets produced from the decay of the boosted heavy neutrinos are mostly nonisolated, following the standard isolation criteria used by the CMS and ATLAS [31,32] Collaborations.We, therefore, use substructure-based variables employed by the experimental collaborations called lepton sub-jet fraction (LSF) and lepton mass drop (LMD) to distinguish our fatjet signal from the non-reducible quantum chromodynamics (QCD) background [33].These variables (LSF, LMD) will be defined later in section 4.
The rest of the paper is organized as follows.We give a brief description of the LRSM model with a universal seesaw mechanism in section 2. In section 3, we discuss the constraints on the model parameters and phenomenologically allowed benchmark points for our analysis.In section 4, we present our collider analysis and results using the variables LSF and LMD.We conclude with our findings in section 5.

The Model
A brief description of the model [25] is given below.The leftright symmetric model is based on the gauge group try to incorporate a lepton-specific scenario.The charge of a particle is defined as The particle spectrum of the model along with their quantum number and Z 2 charges are shown in Table 1.The matter structure consists of three families of SU(2) L and SU(2) R quark and lepton doublets, Table 1 Particle spectrum along with their quantum numbers.
where the numbers in the parentheses denote the quantum numbers under SU(3) C ×SU(2) L ×SU(2) R ×U(1) B−L gauge groups, respectively.For the generation of SM quark and lepton masses through a universal seesaw, the model includes heavy singlet quarks ); and singlet heavy neutrinos The Higgs sector consists of four SU(2) doublets: Under the Z 2 symmetry, the singlet leptons , where v EW = 246 GeV is the electroweak (EW) vacuum expectation value.The hierarchy of the VEVs is arranged as The VEVs of the SU(2) R scalar doublets, v RQ and v Rl , are responsible for generating the mass of W R and Z R gauge bosons, while the VEVs to the SU(2) L scalar doublets,v LQ and v Ll , break the EW symmetry and generate the SM gauge boson (W and Z) masses.
The full Lagrangian of this model can be written as where L kinetic contains the kinetic terms for the gauge boson, scalars and fermions in the model, while the Yukawa La-grangian, L Y and the scalar potential, V (H) are given later in Eq. ( 14) and Eq. ( 21), respectively.
The covariant derivatives appearing in L kinetic can be written as, where V µ is the U(1) B−L gauge boson and g V its gauge coupling, while W L , W R and g L , g R are the gauge bosons and gauge couplings corresponding to the SU(2) L and SU(2) R gauge groups, respectively.The mass of the charged gauge bosons W R and W L are given by The masses of the two neutral massive gauge bosons are given by in the limit v EW << v RQ , v Rl .Here, g Y is the effective SM U(1) Y gauge coupling given as In this model, the Z R is always heavier than W R , which implies that a strong limit on W R mass indirectly puts a stronger bound on Z R mass.We discuss the generation of masses for the fermions through the seesaw mechanism in the the following subsection.

Seesaw and fermion masses
The quarks and leptons obtain their masses through a universal seesaw with the help of heavy singlet fermionic states.The Lagrangian in this model contains the following gauge invariant Yukawa terms: Here, Y iA 's are the Yukawa coupling matrices, and M X 's are the heavy singlet mass terms.The H L/R is the conjugate scalar defined as All the charged fermions, i.e., quarks and charged leptons, acquire their masses by diagonalizing 6 × 6 matrices through a universal seesaw mechanism.Three light eigenstates among the six eigenstates are identified as the three SM states.
The up quarks acquire their masses from the following mass matrix in the (u,U) basis.This matrix is diagonalized by the following bi-unitary transformation: where U uL and U uR are 6 × 6 unitary matrices transforming the left-handed and right-handed fermions from the gauge basis to their mass basis.The down quarks also acquire their masses through the seesaw mechanism, where the corresponding matrices are U dL and U dR .For simplicity, we choose the Y uL/R and Y dR to be diagonal and keep only Y dL to be off-diagonal so as to generate the correct Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix for the SM quarks.The 6 × 6 mixing matrix contains the SM CKM matrix in the top-left (bottom-right) 3 × 3 block if mass eigenvalues are arranged in ascending (descending) order.The mixing between the SM light quarks and the heavy quarks is negligible except for the top quark as it is heavy.
The charged leptons acquire their masses from the following 6 × 6 mass matrix: The Yukawa coupling matrices (Y eL/R ) and the heavy lepton mass matrix (M E ) are chosen diagonal to prevent any charged lepton flavor violating interactions.The hierarchical structure of the block matrices that give the correct masses to the SM charged leptons (e, µ, τ) allow negligible mixing between the SM charged leptons and heavy charged leptons.
The neutrinos obtain their masses from the following 12 × 12 mass matrix: formed in the basis . Similar to the quark sector, all the Yukawa couplings and mass matrices are taken diagonal except Y νL that generates the Pontecorvo-Maka-Nakagawa-Sakata (PMNS) mixing matrix [4].

Scalar masses
The scalars obtain their masses by minimizing the following gauge-invariant potential where The terms µ 12 and µ 34 break the discrete Z 2 symmetry softly and prevent the formation of domain walls which could otherwise destabilize the model [34,35].
The Higgs boson spectrum consists of four CP-even states, two CP-odd states, and two charged Higgs bosons.A major distinctive feature from the more widely studied minimal LRSM model which includes scalar triplets, is the absence of a doubly charged Higgs boson in the particle spectrum of our model.Two charged goldstone bosons are eaten up by the W L and W R gauge bosons to give them mass, while two neutral goldstone states provide mass to the Z and Z R .

Phenomenological Benchmark points
We now highlight some of the constraints on the parameter space in our model before discussing the choice of our benchmark points for the collider analysis.We note that the most relevant constraints arise from the experimental bound on heavy gauge boson masses.We also mention a few theoretical constraints that affect our parameter choices.

Theoretical and phenomenological constraints
The theoretical constraints on some of the couplings in the model come from the requirement of perturbativity and unitarity.The perturbativity condition requires that at least at the electroweak symmetry breaking (EWSB) scale, all quartic couplings satisfy whereas the Yukawa and gauge couplings need to be less than √ 4π [36].Furthermore, tree-level unitarity in the scattering of Higgs bosons and longitudinal components of EW gauge bosons necessitates that the eigenvalues of the scattering matrices must be less than 16π [37,38].
Additional constraints appear as a result of electroweak precision measurements, particularly the oblique parameters [39].
The addition of extra Higgs doublets has little effect on the oblique parameter in general, particularly the T -parameter because the custodial SU(2) stays intact at the tree level.It can be broken at the loop level by the effect of other states.However, the additional states in our model primarily belong to the right-handed sector and are rather heavy.Although the mixing between the left and right sectors is tiny, it cannot cause significant changes in the EW W or Z boson masses.The consistency of our chosen benchmarks has been checked against existing limits [40].

Experimental constraints
The non-observation of any direct signal at the LHC has put for the electron (e ± ) and muon (µ ± ) channels, respectively.
On the other hand, ATLAS [42] di-jet search puts a lower bound of 4 TeV on W R mass with SM-like couplings (i.e., 50% branching in di-jet mode).These limits differ depending on the branchings of W R in the di-jet and the l ± N channels.These dependencies are summarized in Fig. 6 of Ref. [41] and Fig. 4c of Ref. [42], and one can extract the bound on M W R depending on the value of production cross sections times the branching ratios.Both the CMS [41] and In the di-jet mode, one notes that the branching for W R is 66% in our model, and thus the limit on W R mass increases (see right-panel of Fig. 1).We have a common limit of 4.3 TeV on W R mass coming either from l ± N or di-jet search for our parameter choices.
The heavy neutral gauge boson (Z R ) has a rather weak limit from the direct di-jet search [47] excluding m Z R < 2.9 TeV having SM-like gauge couplings with the SM fermions.However, in LRSM the Z R is heavier than W R , and thus a strong limit on W R mass indirectly puts a stronger limit on Z R mass (see Eqs. ( 11) and ( 12)).The 4.3 TeV lower limit on W R mass, as discussed above, puts a lower limit of 5 TeV on Z R mass.
In addition, we have also considered heavy neutrino search results from the LHC which put limits on m N [44, 48, 49], and these are included in our benchmark selection.TeV.As our non-standard neutral and charged scalar states are considered to be very heavy they do not affect our analysis.

Chosen parameter space
Motivated by the signal of fatjet searches arising from heavy vector-like leptons and heavy neutrinos, we choose the heavy charged lepton mass to be ∼ 0.   2 for the three benchmark points.In Table 3, we show fermion masses chosen for our phenomenological study.
All the Yukawa matrices in Eq. ( 14    in Table 3.The Yukawa couplings are chosen as We require Y 33 dR ≃ 0.092(0.078) in BP1 (BP2 and BP3); Y 33 eL ≃ 0.146(0.132) in BP1 and BP2 (BP3) to keep the third generation heavy fermion masses in few TeV range.We set the offdiagonal elements in Y dL matrix by fitting them to get the experimentally measured CKM mixing matrix [60]

Collider Analysis
In this section, we study the collider signatures of our model at the high luminosity run of the LHC with the chosen benchmark parameters discussed above.Even though the exotic charged and neutral leptons are very heavy and have very limited production rates through SM gauge bosons, they can be resonantly produced through the heavy SU(2) R gauge bosons with comparatively large rates.We therefore focus on the production of heavy neutrinos and heavy charged leptons through the s-channel processes mediated by the heavy right-handed gauge bosons.Once produced, these heavy charged leptons (E i 's) and neutrinos (ν i , where i ≥ 4) finally decay to SM leptons and jets with almost 100% branching ratios, as shown in Table 4 for The branching ratios for the other benchmark points, viz.In our work, we therefore look for fatjets containing charged leptons and analyze these fatjets by studying their substructure.We consider signals with different fatjet multiplicities, that include inclusive two, three, and four fatjets.Note that our primary production channel involves particles (E i and ν i ) that carry a lepton number which leads to some of the fatjets having a non-isolated charged lepton (e, µ) within the fatjet radius.We use the jet-substructure techniques [33] to identify the sub-jet associated with these boosted charged leptons when they cannot be isolated at the detector.
Fig. 3 Representative Feynman diagram for three fatjet signal.

Signal and background
Here we describe the fatjet signals we are interested in, followed by the possible SM background in detail.We have generated events for all our signal and background processes using MadGraph5_aMC@NLO v2.7.3 [61] at leading order (LO) in QCD with a dynamic choice of factorization scale given by , where M i is the transverse mass of final state particles.We use nn23lo1 [62] for the parton distribution functions (PDFs).Events are passed to PYTHIA8 [63] for showering and hadronization followed by a fast detector simulation in Delphes v3.4.2 [64], with added pile-up events embedded in Delphes v3.4.2 for the high luminosity LHC (HL-LHC).The final state hadrons with non-isolated high P T leptons are clustered using the Cambridge-Aachen algorithm of FastJet [65,66] with a jet radius R 0 = 0.8 (AK8 fatjet).We groom [67] the fatjets using soft-drop method [68,69] to remove soft and wide angle radiations after some pre-selection cuts.We use the soft radiation fraction parameter z cut = 0.1 and the angular exponent parameter β = 0 for the grooming [41].The groomed AK8 fatjets are required to have a threshold on the soft-drop mass (m SD ) as m SD > 40 GeV.
To identify the sub-jets associated with leptons and leading hadrons we employ the N-subjettiness technique, where we use 'OnePass General ET General KT Axes' with p = 0.6 (for kt and Cambridge-Aachen axis choice of p = 1 and p = 0, respectively) [70,71].To use this method for finding sub-jets out of a fatjet, we use Fastjet Contrib [72].The t t background Table 4 The relevant branching ratios are listed for the decay channel of the heavy gauge bosons W R and Z R , heavy leptons E i , and lightest three heavy neutrinos ν i , (i = 4, 5, 6) in BP1.

Fatjet Signals
We consider three types of signals containing fatjets for our analysis which are described below.Note that all the signal fat- iii.
then decays into e/µ with additional two jets).
Four-fatjet inclusive searches: In this signal, the final state contains at least four high P T fatjets with the leading two fatjets inclusive of only one high P T lepton within their jet radius.
For this case we find that only pp → Z R , Z R → E 4 Ē4 , E 4 → ν 4 j j contributes dominantly where ν 4 as before decays into e/µ + j j.

All individual contributions in the respective final states
are summed together and identified as the signal.

Backgrounds
For the above signals, we consider the following SM background.
QCD multi-jet: QCD multi-jet events will be a major source of background due to its large production cross section.We primarily focus on the 4-jet final state, which can lead to our two, three, or four fatjet signals where some jet could also fake a lepton as a part of the fatjet.We however find that a large P T cut on the fatjet and a cut on the jet invariant mass helps to reduce the QCD background by a significant amount.

Top pair production (t t):
The other dominant background for our signal comes from the SM production of t t where the leptonic decay mode of t t can give us fatjets with two high P T lepton as part of the fatjets.On the other hand, the semi-leptonic and fully hadronic decay mode of t t can also be a dominant background where one or two hadrons can fake as leptons or have leptonic decays in the detector.
W /Z+jets and tW : The weakly produced ZW and tW processes can also be possible sources of fatjet background.However, the production cross section of both these processes are very small compared to the QCD and t t background.In addition, for the tW case the probability of getting at least two heavy and highly boosted fatjets is much smaller than t t due to the lower mass of W .The other SM subprocesses that can give fatjet signatures are W /Z+jets.Although the production cross section for W /Z+jets is comparatively much larger than tW , it has less probability of giving at least two high-energy leptons in two separate fatjets, when compared to the t t leptonic decay case.We therefore, neglect the above backgrounds in our analysis.

Variables
To achieve good signal significance over the aforementioned backgrounds, we identify some kinematic variables which have characteristically different distributions for the signal compared to the SM background.After analyzing various distributions, we find the relevant variables to be the P T and invariant mass of the fatjets, P T of the sub-jets associated with leptons, and two interesting substructure variables LSF and LMD.
The P T -ordered jets and the jet-mass are denoted by P T j i and m j i where i = 1, 2, 3, 4. The invariant mass of any two, three, and four fatjets are also important variables and are called m j 1 ... j n .In addition, we include two jet substructure variables that have been used for fatjet analyses in CMS searches [41], viz. the lepton sub-jet fraction (LSF) and lepton mass drop (LMD) of the sub-jet associated with charged leptons [33].
To compute the variables LSF and LMD, we cluster all final state particles in the event (including leptons) into a fatjet.
For each fatjet we then use jet-substructure (JSS) algorithms to cluster its constituents in three sub-jets using N-subjettiness techniques [72].We then check for sub-jets associated with all high P T leptons.Thereafter we calculate the lepton sub-jet fraction of each lepton (defined as the ratio between the lepton P T and P T of the sub-jet associated with it) which is given by where n = 1, 2. LSF 1 is associated with the leading lepton while LSF 2 is for the sub-leading one.The lepton mass drop parameter is constructed in a similar way and given by where m s j represents the invariant mass of the sub-jet associated with the lepton and m s j−l n represents the invariant mass of the same sub-jet with the n th lepton subtracted out.
We now concentrate on the distributions of the aforemen-     values of P T because of the large mass gap between the decay products and their mother particles.In contrast, the background events, say from the t t production where the jets come from the decay of SM W bosons (including the b-jet), have relatively softer P T beyond the pre-selection cut.Note that for the top quarks to be highly boosted, they would have to be produced at large √ ŝ which would lead to a smaller cross section for the t t production.The QCD multijet background also falls off rapidly for large P T .The middle-panel of Fig. 4 shows the jet mass (m j ) distributions of the two leading fatjets, while the bottom-panel shows the distribution of the invariant mass of the two leading fatjets.For the signal events, one finds that the distributions in both m j 1 and m j 2 show two peaks that appear at the mass values of the heavy neutrinos (especially ν 4 and ν 5 ), whose highly boosted decay products correspond to the reconstructed fatjets.The jet-mass variable therefore helps us reconstruct the intermediate particles for the signal, which in turn also helps us to reject a large part of the SM background, where no such peak is observed at large values of m j 's.The four-fatjet inclusive search is expected to be almost without any SM background and we focus only on the event rates for this final state.

Results
We now use the kinematic variables discussed in the previous section and exploit their characteristic distributions for signal and background to define the event selection cuts we use to analyze the three different fatjet signals.We have used a preselection criterion for signal and background in all our subsequent analyses given by: Pre-selection criteria : For the pre-selection of events, we demand at least two charged leptons (e, µ) with P T > 25 GeV at the event generation level along with two fat jets with radius R = 0.8 and minimum P T > 200 GeV which must lie within the rapidity gap of |η| < 4 for the ungroomed jet.As pointed out earlier, the strong pre-selection cut on the fatjet P T helps suppress the background events, dominantly coming from the QCD multijet processes.As the QCD background has a huge cross section, it is quite challenging to generate uniformly populated events in our simulations.The pre-selection cuts at the generation level allow us to lower the large QCD cross sec-Table 5 Cut summary for two-fatjet inclusive searches.

Two-fatjet inclusive searches:
The two-fatjet signal gives the largest cross section among all the fatjet signals.The signal and background distributions shown in Fig. 4 and Fig. 5 are used to define the event selection criterion to help achieve a good signal significance.We classify the kinematic selections for the signal and background events into two categories, which we call Cut-A and Cut-B.
Cut-A : i.The leading groomed fatjet is required to have P T j 1 > 750 GeV, while the sub-leading groomed fatjet should have ii.The jet mass of both the leading and sub-leading groomed fatjet should have a minimum value of m j i > 125 GeV (i = 1, 2) while their invariant mass must satisfy m j 1 j 2 > 1.5 TeV.
iii.The sub-jets associated with the leading charged lepton must have P T s j l 1 > 400 GeV and sub-jets associated with the sub-leading charged lepton are required to have P T s j l 2 > 250 GeV.iv.The variables LSF 1 and LSF 2 need to be > 0.9, while LMD 1 and LMD 2 has an upper limit satisfying < 0.05 after applying the grooming method.

Cut-B :
i.The leading groomed fatjet is required to have P T j 1 > 750 GeV, while the sub-leading groomed fatjet must have ii.The minimum value for the variables LSF 1 and LSF 2 is > 0.9, while the maximum value for LMD 1 and LMD 2 is < 0.05 after applying grooming method.
A summary of all cuts mentioned here is also described in Tables 5.
Note that the main difference between Cut-A and Cut-B is simply on the way the lepton inclusive fatjets are treated, along with how we utilize the jet mass in the analysis which carries a bias of the mass of the new exotics.We list the number of signal and background events expected with an integrated luminosity of L = 3000 fb −1 for BP1 in Table 6.We also list the cross sections for different signal sub-processes contributing to the two-fatjet inclusive search.The events with Cut-A and Cut-B are shown in Tables 13 and 14 for BP2, and BP3 in Appendix C.
As the signal events are peaked at very high P T (Fig. 4) the strong cut on their P T removes all of the background events coming from the QCD multijet and the hadronic decay of the t t production.In addition, the choice on the jet mass, invariant mass of the di-fatjet system and the additional cuts on the jet substructure variables are used to exploit the presence of the highly boosted charged lepton and are helpful to remove the background events coming from the leptonic decay modes of the t t production (Fig. 5).For all the BPs, we conclude that it is possible to remove the SM background with moderate efficiency using the pre-selection cuts, which is significantly im-Table 6 Signal and background events surviving after applying the pre-selection criteria as well as selection cuts at LHC with center-of-mass energy, √ s = 14 TeV and L = 3000 fb −1 in the inclusive two-fatjet analysis for BP1.proved after using Cut-A and Cut-B.We also note that Cut-A is a much stronger selection criterion compared to Cut-B as it requires an additional condition on the jet mass and invariant mass of the fatjet pairs, which reduces the signal events too.Therefore, Cut-A will yield better significance where the production cross section of the signal is relatively large, while

Data sets
Cut-B will perform better where the signal cross section is small, and we do not want a signification reduction in the event rates for the signal.
We calculate the signal significance [74] using the formula where s and b stand for the total number of signal and background events surviving after cuts and L is the integrated luminosity.The signal significances for all the BPs in the inclusive two-fatjet final state are summarized in Table 7, with integrated luminosities of 3000 fb −1 , 600 fb −1 and 300 fb −1 at the LHC with a center-of-mass energy √ s = 14 TeV.For BP1, we are able to achieve a significance of more than 5σ at L = 600 fb −1 , while for BP2, which represents a point with a heavier W R , we need 3000 fb −1 luminosity to achieve 5σ significance.This is due to the large mass of W R and Z R in BP2 compared to BP1, which reduces the production cross section of the con-tributing sub-processes for BP2, although the cut efficiency for the signal events is better in the case of BP2.In the case of BP3, the signal significance is better compared to BP2 because of the larger production cross section and higher boost for the heavy leptons.For BP3, an integrated luminosity of 600 fb −1 is required for a 5σ discovery.We note that Cut-A and Cut-B provide roughly the same amount of significance for BP1.
While for BP2 and BP3, Cut-A gives better significance as the jet mass and P T of sub-jets play an important role in reducing the background further.Cut-A : i.The leading groomed fatjet is required to have P T j 1 > 500 GeV, while the sub-leading groomed fatjet should have P T j 2 > 300 GeV.ii.The jet mass of both the leading and sub-leading groomed fatjet should have a minimum value of m j i > 150 GeV (i = 1, 2), while their invariant mass must satisfy m j 1 j 2 > 800 GeV.Additionally, the next-to-sub-leading groomed fatjet has jet mass of m j 3 > 100 GeV.
iii.The sub-jets associated with leading charged lepton must have P T s j l 1 > 100 GeV.

Cut-B :
i.The leading groomed fatjet is required to have P T j 1 > 750 GeV, while the groomed sub-leading fatjet should have P T j 2 > 500 GeV.
ii.The variables LSF 1 and LSF 2 have a minimum value of 0.9, while LMD 1 and LMD 2 are bounded from above by < 0.05 after applying grooming techniques.

Cut-C :
i.The leading groomed fatjet is required to have P T j 1 > 750 GeV, while the sub-leading groomed fatjet should have P T j 2 > 500 GeV.
ii.The jet mass of both the leading and sub-leading groomed fatjet should have a minimum value of m j i > 150 GeV (i = 1, 2), while their invariant mass must satisfy m j 1 j 2 > 1000 GeV.Additionally, the next-to-sub-leading groomed fatjet has jet mass of m j 3 > 100 GeV.
iii.The sub-jets associated with leading charged lepton must have P T s j l 1 > 100 GeV.
Note that the role of the cuts exploiting the P T , jet mass, invariant mass of fatjet pairs, and the substructure variables, LSF and LMD remain very similar to the two-fatjet analysis, albeit with differences in the threshold choices for the event selection.The effect of the three sets of cuts on our signal and SM background is tabulated in Table 15, 16, and 17 in Appendix C for BP1, BP2, and BP3, respectively.Here, Cut-C represents the strongest selection criteria for the events as it includes the additional requirements on the jet mass and a stronger cut on the di-fatjet invariant mass.We note that Cut-B which is the only set that includes the cut on the leptonic sub-jet variables, is useful in suppressing the the SM background most and also retains the maximum signal events.This is in agreement with our earlier observation in the two-fatjet analysis where we found that the jet-mass and invariant mass selections also suppress the signal events considerably.The signal significances for all the BPs are shown in Table 9 for L = 3000 fb −1 , L = 600 fb −1 , and 300 fb −1 .We find that Cut-B that includes the LSF and LMD substructure variables, gives the best signal significance as it retains the maximum signal events for all the BPs.Cut-A and Cut-C give similar but lower significance.In the three-fatjet inclusive analysis, we require L = 600 fb −1 for BP1 and L ≊ 3000 fb −1 for BP2 and BP3 to achieve 5σ discovery.
Four-fatjet inclusive searches: An interesting and more exotic signal for the model would be in the form of four fatjets.However, for such a final state, one can expect a very low signal cross section for our benchmark choices as only the p p → Z R → E 4 Ē4 process contributes and is suppressed due to the large Z R mass (already bounded from below due to the W R mass limits in an LRSM framework).A part of the parameter space can be tuned to generate a conducive set of masses to improve some decay branching fractions, but we are still constrained by the luminosity reach of the LHC.The SM background is expected to be quite negligible for such a final state and this signal may prove useful to consider at more high-energy machines such as the future 100 TeV hadron collider [75].As the LRSM gauge boson masses are excluded to be heavier and go beyond the reach of the LHC, the future hadron collider with √ s = 100 TeV will be able to pro- • N l ≥ 2 with P T > 25 GeV • P T j 1 ( j 2 ) > 500(300) GeV • P T j 1 ( j 2 ) > 750(500) GeV • P T j 1 ( j 2 ) > 750(500) GeV • Two fatjets with radius R = 0.8.
• m j i > 150 GeV (i = 1, 2), • LSF 1 > 0.9, LSF 2 > 0.9 • m j i > 150 GeV (i = 1, 2), duce them easily and will give significantly large cross sections [76].The large production rates will help study the fourfatjet signal very easily.For the sake of comparison, we go ahead and analyze our four-fatjet signal for the previously considered benchmark points to show what sensitivity could be expected at the LHC with a center-of-mass energy √ s = 14 TeV.
Placing similar pre-selection criteria as in the previous cases, we choose a single set of cuts to show our result: Cut-A : i.The P T of leading fatjet, P T j 1 > 500.0GeV and sub-leading fatjet, P T j 2 > 300.0 GeV.ii.Jet mass of leading jet, m j 1 > 150.0 GeV and sub-leading jet, m j 2 > 120.0 GeV while invariant mass of leading and subleading jet, m j 1 j 2 > 800.0 GeV.Additionally, the invariant mass of all four fatjets, m j 1 j 2 j 3 j 4 > 1000.0GeV.
iii.The P T of sub jets associated with leading charged lepton satisfies P T s j l 1 > 200.0 GeV.
We summarize our cuts for four-fatjet inclusive searches in Tables 10.
As expected, we find that the four-fatjet analysis is not feasible with the available cross section for the signal and, therefore, will not be accessible at the LHC, even with L = 3000 fb −1 as shown in Table 11.This significance is further diminished if the LSF and LMD variables are invoked in our analysis.

Conclusion
In this work, we perform a collider study to look for fatjets

Cut summary
Pre-selection cut Cut-A • N l ≥ 2 with P T > 25 GeV • P T j 1 > 500 GeV,P T j 2 > 300 GeV • Two fatjets with radius R = 0.8 • m j i > 150 GeV (i = 1, 2), m j 2 > 120 GeV (i = 1, 2) • m j 1 j 2 > 800 GeV, m j 1 j 2 j 3 j 4 > 1000 GeV    Table 13 Signal and background events surviving after applying the pre-selection criteria as well as selection cuts at √ s = 14 TeV and L = 3000 fb −1 for two-fatjet searches in BP2.The individual contributions of SM subprocesses to the background are as shown in Table 6.
stringent constraints on the mass of the heavy right-handed gauge bosons Z R , W R as well as the heavy Majorana neutrinos in LRSM.The different constraints leading to limits on the model are briefly summarized below.Searches for heavy resonances: The search for heavy resonances decaying to di-lepton/di-jet final state at the LHC put strong constraints on their production cross section.The major experimental constraint in our model comes from the W R search in the l ± N channel (N represents the heavy right-handed neutrino) where the final state is either one charged lepton and a fatjet containing a high P T lepton [43] or two same-sign leptons and two jets [44, 45], depending on the mass difference of W R and N. The most recent CMS search [41] puts a lower limit on W R which excludes its masses below 4.7 TeV and 5.0 TeV ATLAS [43] collaborations have used the minimal LRSM (MLRSM) model in their search for the heavy W R boson, where the W R → eN decay branching ratio is about 10% for large values of the W R mass [46].In our case, the branching ratio is relatively smaller (< 5%) due to new modes of decay available for the W R boson into heavy leptons (seesaw partners) and heavy neutrinos, as shown in Table 4.As the decay branching ratio of W R → eN is suppressed, limits on the W R mass are relaxed (see left-panel of Fig. 1) when compared to the CMS [41] 4.7 TeV bound in the electron channel.
Higgs searches: We avoid flavorviolating interactions of charged leptons with neutral scalars by choosing diagonal configurations for the corresponding Yukawa matrices involving the SM and heavy exotic leptons.There are no flavor-changing neutral Higgs (FCNH) interactions in our model at the tree-level because of the diagonal couplings for the neutral scalars.Heavy non-standard neutral scalars, pseudoscalars and the charged Higgses have strong bounds [50, 51] from low energy flavour-changing neutral current (FCNC) effects, especially from the K-meson and B-meson mass mixing.These bounds have been studied for the MLRSM in Ref. [52].As the FCNH effects are absent, these bounds do not apply on our neutral states.The W R mass limits coming from the meson mixings are much weaker than the direct collider bounds.The signal intensity measurements of the 125 GeV Higgs in several final states, including ZZ, WW , bb, ττ, and γγ final states [53, 54], provide additional constraints, which have been considered while choosing our benchmark points.The gauge and Yukawa couplings of the 125 GeV scalar seem to lie extremely close to their SM value, as suggested by the experimental results.We thus limit our study to the alignment limit, i.e., |y V h | ∼ 1, which is the modification factor to the hVV (V = W /Z) couplings due to new physics.The parameter space of interest is further restricted by collider searches for non-standard neutral scalar states and charged scalar states.These searches are carried out at the LHC [55-59] in a variety of SM final states limiting their masses in the range of 1 to 2

72 − 1
TeV while three heavy neutrino masses in the range of ∼ 180 -700 GeV and m W R ≈

Fig. 1
Fig.1Limit for the cross section times the branching ratio of W R as a function of its mass at 95% C.L. are shown from CMS heavy neutrino search[41] (left-panel) and ATLAS dijet search[42] (right-panel).Theoretical expectations for the cross sections are also shown with the branching ratios used in the experimental search (10% in left-panel and 50% in right panel ) and in our analysis (5% in left-panel and 66% in right-panel ).Please refer to the text for more details.

m U = 11 19 − 33 νLFig. 2
Fig.2The CKM satisfied points for the Y dL , and the U PMNS and ∆ m 2 i j satisfied points for the Y νL are shown in the left-panel and right-panel, respectively for BP1.
within 4σ error, with Y 21 dL = Y 32 dL = Y 31 dL = 0.The values of Y 12 dL , Y 13 dL , and Y 23 dL are shown in Fig. 2 in left-panel for different scanning configurations (horizontal axis) in BP1.The values are |Y 12 dL | ≃ 10 −3 , |Y 13 dL | ≃ 3 × 10 −3 , and |Y 23 dL | ≃ 4 × 10 −2 for all scanning configuration.The values of Y 12 dL , Y 13 dL , and Y 23 dL are of the same order of magnitude as above, for all our benchmark points.Similar to the quark sector Yukawa, we keep the Y νR diagonal while Y νL is chosen non-diagonal in order to obtain the neutrino mixing matrix U PMNS .We choose Y νR = Diag (0.2, 0.3, 0.4) and scan the parameter space for Y i j νL which is shown in Fig. 2 in the right-panel with different scanned configurations (horizontal axis) which satisfy 6.82 × 10 −5 eV 2 < ∆ m 2 21 < 8.04 × 10 −5 eV 2 , 2.43 × 10 −3 eV 2 < ∆ m 2 31 < 2.60 × 10 −3 eV 2 and U PMNS mixing elements within 3σ error [4, 5].The relative values of Y i j νL are below 4 × 10 −5 keeping the relative order in-tact among themselves.We have chosen the first point in the scanned configuration in our analysis.

BP2
and BP3 are shown in Table12in the Appendix A. The decay of these heavy fermions (produced from the decay of very heavy SU(2) R gauge bosons) leads to the interesting possibility of final states with boosted objects.What we finally observe are some clusters of highly boosted and collinear particles, which can only be detected as a large radius jet in the detector.The presence of these fatjets play an important role in achieving good signal significance over the estimated SM background.In addition, the fatjet from the decay of heavy charged lepton also contains a non-isolated SM charged lepton along with two sub-jets.The fatjet originating from the heavy neutrino decay contains two sub-jets, as shown in Eq. 1.
jets have P T > 200 GeV before grooming.This pre-selection cut on the reconstructed fatjets is imposed to suppress the SM background, especially the large QCD multijet background.As we wish to study the signal which dominantly comes from the production of the heavy SU(2) R gauge bosons which have masses above 5 TeV, the fatjets originating from them are expected to have significantly large P T .This pre-selection helps in removing the background contributions without compromising on the signal events.Two fatjet inclusive searches: In this signal topology, our final state contains at least two high P T fatjets, with each containing one high P T lepton lying inside the fatjet radius.The five different production channels that contribute to this final state in our model are i.Signal 1: pp → W R , W R → E 4 ν 4 , E 4 → ν 4 j j (both the ν 4 then decay into e/µ with additional two jets), ii.Signal 2: pp → W R , W R → E 6 ν 5 , E 6 → ν 5 j j (both ν 5 will then decay into e/µ with additional two jets), iii.Signal 3: pp → W R , W R → E 6 ν 5 , E 6 → µ j j (ν 5 then decays into e/µ with additional two jets), iv.Signal 4: pp → Z R , Z R → (ν 4 ν 4 )/(ν 5 ν 5 ) (ν 4 /ν 5 then decay into e and µ with additional two jets), v. Signal 5: pp → Z R , Z R → E 4 Ē4 , E 4 → ν 4 j j (ν 4 then decays into e/µ with additional two jets).Three-fatjet inclusive searches: Here our final state contains at least three high P T fatjets, where the two leading (P T -ordered) fatjets must contain one high P T lepton in them.The three production channels that dominantly contribute to this three-fatjet signal are given by i. Signal 1: pp → W R , W R → E 4 ν 4 , E 4 → ν 4 j j (all ν 4 then decay into e/µ with additional two jets), a representative Feynman diagram for this channel is shown in Fig. 3, ii.Signal 2: pp → W R , W R → E 6 ν 5 , E 6 → ν 5 j j (both the ν 5 then decay into e/µ with additional two jets),

>
tioned variables for two and three inclusive fatjet searches, where the signal and dominant SM background contributions are shown together to characterize the differences between their distributions.The normalized distributions are shown with events that have passed the following selection cuts: 200 GeV; |η fatjet | < 4. (28)In Fig.4 and 5, we show the relevant distributions of the selected kinematic variables in the two-fatjet inclusive search for BP2.In the top-panel of Fig.4, we have shown the P T distribution of the leading and subleading fatjet where a pre-selection cut of P T > 200 GeV (before grooming for two leading fatjets) is in place.The distribution for the signal events peak at higher

Fig. 4
Fig. 4 Normalized distributions of the signal and different SM backgrounds as a function of the transverse momentum (P T ) of the two leading fatjets (top-panel), the mass (m j ) of two leading fatjets (middle-panel), and invariant mass (m j 1 j 2 ) of the two leading fatjets (bottom-panel) in the two-fatjet inclusive searches for BP2.

Fig. 5
Fig. 5 Normalized distributions of the signal and different SM backgrounds as a function of the transverse momentum of the leptonic sub-jet (P T s j l ) of the fatjets (top-panel), the lepton sub-jet fraction variable involving the leading (LSF 1 ) and sub-leading (LSF 2 ) charged lepton (middle-panel) and the lepton mass drop variable with the leading (LMD 1 ) and sub-leading (LMD 2 ) charged lepton subtracted out (bottom-panel) in the two-fatjet inclusive searches for BP2.

Fig. 5
Fig.5highlights the relevance of the substructure variables used in our analysis.The top-panel shows the event distribution of the signal and background as a function of P T of the sub-jets inclusive of the leading and sub-leading charged lepton.As expected, the signal fatjets which carry much larger P T than the fatjets in the SM background events also give larger P T sub-jets within them, leading to the sub-jets also getting a hard P T distribution compared to the background.The middleand bottom-panels of Fig.5display the event distributions as a function of LSF and LMD variables, respectively.The LSF variables are found to sharply peak around 1 for the signal and the t t (leptonic and semileptonic) background.This happens when the lepton is boosted enough and becomes one of the hardest constituents in the fatjet.As the signal and t t (leptonic and semileptonic) channels are most likely to give a fatjet with a non-isolated charged lepton within, these modes are the only ones that show the peak behavior around 1. Note that the peak has vanished in the distribution of LSF 2 for the t t (semileptonic), as only one top decays in the leptonic channel.It is therefore less likely to give a second charged lepton inclusive fatjet.In addition, the charge leptons do not suffer from soft radiations into the sub-jets, and we can reconstruct the energy

Three-fatjet inclusive
searches: We have already discussed the major signal contribution to three-fatjet inclusive searches for our scenario in Section 4.1.1.The distributions of various variables are also shown in Fig. 7 and Fig. 8.To analyze the threefatjet final state events, we consider three different sets of kinematic selections (Cut-A, Cut-B and Cut-C).A summary of all cuts mentioned below is also summarized in Tables 8.
originating from boosted particles that come from the decay of heavy resonances in a left-right symmetric extension of the SM.The model contains heavy leptons, heavy quarks, and heavy neutrinos, which play a fundamental role in achieving a universal seesaw mechanism for the generation of all the SM fermion masses.The fatjets, which include non-isolated leptons, originate from heavy neutrinos (few 100 GeV) that decay from a heavier (few TeV) right-handed gauge boson.We employ substructure-based variables LSF and LMD, together with hard kinematic cuts, to search for fatjets with associated charged leptons, to reduce SM background while keeping enough statistics for the signals.

• |∆ η| < 4 of two ungroomed fatjet • P T s j l 1 > 1 Cut
200 GeV.Table 11 Signal significance of the four-fatjet signals at the LHC with a center-of-mass energy √ s = 14 TeV with L = 3000 fb −1 and L = 600 fb −1 .After applying Cut L = 3000 fb −1 L = 600 fb −We analyze a multi-fatjet signal topology that would be typical in the model when the heavy fermionic states are produced and finally decay to SM particles.We have focused mainly on two-fatjet and three-fatjet final state searches while commenting on the possibility of a four-fatjet scenario with very limited signal sensitivity for our parameter choices.A critical aspect of the analysis included the presence of a charged lepton in the fatjets, helping us identify the signal over the SM background by employing the variables LSF and LMD.In the two-fatjet final state, both jets include a charged lepton, while for three-fatjet, two of them include a charged lepton.We find that the two-fatjet signals, which also include leptons, can be discovered at the LHC with a moderate luminosity of L ≃ 600 fb −1 for m E ≃ 1 TeV, m W R = 4.5 TeV (BP1) and m E ≃ 0.72 TeV, m W R = 5 TeV (BP3), while the signal can be discovered with L ≃ 3000 fb −1 for m E ≃ 1 TeV, m W R ≃ 5 TeV (BP2), as shown in Tables 7. On the other hand, a three-fatjet signal can be discovered with L ≃ 600 fb −1 (see Table 9 ) for m W R = 4.5 TeV (BP1).In the case of m W R = 5 TeV (BP2 and BP3), we require L ≃ 3000 fb −1 for a 5σ discovery.We finally conclude by noting that although these signals of fatjet can have origins from an altogether different underlying theoretical framework compared to left-right symmetry or a seesaw framework, their observation will be still crucial in the search of new physics at LHC and will provide hints on the presence of new heavy states beyond the SM.

Fig. 8
Fig. 8 Normalized distributions of the signal and different SM backgrounds as a function of LSF and LMD in the three-fatjet inclusive searches for BP2.

Table 3
(25)heavy quarks, charged and neutral lepton masses for the benchmark points with the input parameters shown in Eqs.(24)and(25).

Table 7
Signal significance for the two-fatjet signals for different cuts at the LHC with a center-of-mass energy √ s = 14 TeV with L = 3000 fb −1 , L = 600 fb −1 and L = 300 fb −1 .

Table 8
Cut summary for three-fatjet inclusive searches.

Table 9
Signal significance for the three-fatjet signals for different cuts at the LHC with a center-of-mass energy √ s = 14 TeV with L = 3000 fb −1 , L = 600 fb −1 and L = 300 fb −1 .

Table 10
Cut summary for four-fatjet inclusive searches.

Table 12
Decay branching ratios for W R and Z R , heavy charged leptons E i , and the lightest three heavy neutrinos ν i , (i = 4, 5, 6) for BP2 and BP3.

Table 14
Signal and background events surviving after applying the pre-selection criteria as well as selection cuts at √ s = 14 TeV and L = 3000 fb −1 for two-fatjet searches in BP3.The individual contributions of SM subprocesses to the background are as shown in Table6.

Table 15
Signal and background events surviving after applying the pre-selection criteria as well as selection cuts at √ s = 14 TeV and L = 3000 fb −1 for three-fatjet searches in BP1.

Table 16
Signal and background events surviving after applying the pre-selection criteria as well as selection cuts at √ s = 14 TeV and L = 3000 fb −1 for three-fatjet searches in BP2.The individual contributions of SM subprocesses to the background are as shown in Table15.

Table 17
Signal and background events surviving after applying the pre-selection criteria as well as selection cuts at √ s = 14 TeV and L = 3000 fb −1 for three-fatjet searches in BP3.The individual contributions of SM subprocesses to the background are as shown in Table15.