Charged Higgs induced 5 and 6 lepton signatures from heavy neutrinos at the LHC

We propose an anomaly free gauged U(1) extension of the SM where three right handed heavy neutrinos, being charged under the general U(1) gauge group, are introduced to explain the origin of the tiny neutrino mass through the seesaw mechanism after the general U(1) symmetry is broken. Due to the breaking of the general U(1) symmetry a neutral beyond the standard model gauge boson Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime $$\end{document} acquires mass. There are two Higgs doublets in this model where one interacts with the SM fermions and the other one interacts with the right handed heavy neutrinos and charged leptons. The charged multiplet of the second Higgs can completely decay into the heavy neutrinos and charged lepton in the neutrinophilic limit of the model parameters. The charged Higgs pair production can be influenced due to presence of the Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime $$\end{document} boson at the High Luminosity LHC (HL-LHC) in addition to the neutral SM gauge bosons. The pair produced charged Higgs bosons decay into SM charged leptons and heavy neutrinos. Following the leading decay modes of the heavy neutrinos into charged leptons and W boson we study the 5 and 6 lepton final states after the leptonic and hadronic decay of the W bosons considering solely muons and electrons in the final state. Combining the electron and muon final states we estimate the significance of the 5 and 6 charged lepton processes in the mH±-mN\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{H^\pm }-m_N$$\end{document} plane for different benchmark points of mZ′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{Z^\prime }$$\end{document}. It is found that the 5 (6) charged lepton processes could be probed at the High Luminosity LHC (HL-LHC) with at least 5 (3) σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} significance and there are parameter regions where the significance could be larger.


I. INTRODUCTION
A strong indication of the existence of beyond the Standard Model (BSM) physics has been indicated by the observation of the neutrino oscillation data [1] which lead us to extend the SM.There are many proposed ways and amongst them seesaw scenario which could be the simplest proposal where SM is extended by a SM-singlet heavy Majorana Right Handed Neutrinos (RHNs) [2][3][4][5][6][7].Being SM gauge singlet RHNs mix with the SM light neutrinos to interact with the SM gauge sector.Hence the RHNs being can only be produced at the high energy colliders through the light-heavy mixings which can be naturally small O(10 −6 ) for the TeV scale Majorana type RHNs to reproduce light neutrino mass around 0.1 eV.Attempts have been made to study such TeV scale or lighter Majorana neutrinos at the Large Hadron Collider (LHC) where same sign dilpeton signature could be distinctively an interesting channel from a variety of production modes [8][9][10].With a general parameterization of the Dirac mass matrix [11] this light-heavy mixing can be larger.However, it becomes less than 0.01 in order to simultaneously satisfy different experimental constraints including the neutrino oscillation data, the electroweak precision measurements and lepton flavor violation [12].As a result production of RHNs at the LHC is extremely suppressed.In addition to that RHN can be produced at the electron positron collider in association with neutrinos and in electron proton collider in association with jet.However, those production modes are also suppressed by the square of the mixing.The study of RHN production in these colliders from a variety of final states and the limits on the light heavy mixing has been studied in [13][14][15][16][17][18][19].In these studies light-heavy mixing has a strong implication because production cross section is proportional to the mixing squared.Therefore smallness of the light-heavy mixing makes RHN production from the SM gauge bosons become challenging at the high energy colliders.On the other hand RHN neutrinos can be produced in pair from the neutral current interaction which is suppressed by fourth power of light-heavy mixing.As a result the production cross section further suppressed [20] making heavy neutrino search more challenging.
Apart from the seesaw scenario there is another interesting possibility where a gauged B−L (Baryon minus Lepton) extension of the SM could be a more compelling scenario [21][22][23][24][25][26][27][28][29][30][31] to explain the origin of tiny neutrino mass.Here the RHNs play crucial role to cancel the gauge and mixed gauge-gravity anomalies.After the spontaneous breaking of the B−L gauge symmetry, the RHNs acquire Majorana masses which automatically implement the seesaw mechanism after the electroweak symmetry breaking.Due to the U(1) B−L gauge extension there exists a neutral, BSM gauge boson Z ′ having direct interaction with the RHNs.Due to the nonzero B−L charges the Z ′ boson couples with the SM fermions.As contribution of the seesaw mechanism, the RHNs decay into the SM particles through the light-heavy mixing.There is an alternative approach where SM can be extended by a general U(1) X gauge group.There are three generations of the RHN present in this scenario in order to solve the mixed gauge and gauge gravity anomalies.As a result the general U(1) X charges of the left and right handed SM fermions becomes different which is manifested by their interactions with Z ′ boson showing a chiral nature.A variety of phenomenological studies have been performed including pair production of heavy neutrino at the high energy colliders from Z ′ using prompt [32,33], displaced [34][35][36] and boosted [37] decays of heavy neutrinos, pair production of charged leptons at the electron positron colliders [38] and light Z ′ searches at beam dump experiments [39] respectively.Due to the presence of the general U(1) X charges the production of the heavy neutrino pair from Z ′ can be enhanced which has been proposed in [32,33].
In this paper we consider a general U(1) X extension of the SM based on SU(3) c ⊗ SU(2) L ⊗ U(1) Y ⊗ U(1) X where U(1) X is a generalization of U(1) B−L such that the U(1) X charges of the particles are realized as a linear combination of the SM U(1) Y and U(1) B−L charges so called non-exotic U(1) X model [32,33,40].In the U(1) X model flavor universal charges are assigned for the three generations of the RHNs [32], however, there is another charge assignment for the RHNs to make model anomaly free.In this alternative assignment two RHNs have U(1) X charge −4 and the remaining one is assigned with a +5 U(1) X charge.With a minimal extension of the Higgs sector two RHNs with U(1) X charge −4 are involved in a minimal seesaw mechanism while the other RHN with U(1) X charge +5 can not couple with SM particles because of the U(1) X charges of the SM particles [33].The connection between additional Higgs doublet and heavy neutrinos have been explored in [41][42][43][44][45][46][47].Due to the fact of seesaw mechanism light neutrino masses are generated.This scenario also involves a neutral BSM gauge boson Z ′ which interacts with the RHNs and SM particles of this model and these aspects have been studied in [32][33][34][35][36][37].
Apart from Z ′ , there is another interesting fact that the RHNs with U(1) X charge −4, involved in the seesaw mechanism, can interact with the charged multiplet of a second Higgs doublet involved belonged to the extended Higgs sector.We found that dilepton searches give stringent constraints on the U(1) X gauge coupling.Consider High Luminosity LHC (HL-LHC) at √ s = 14 TeV with 3000 fb −1 luminosity we propose a parameter space on charged Higgs and RHN masses at different signal significance.Due to the fact that all the Higgs fields are charged under the U(1) X gauge group, they can interact with Z ′ .As a result if kinematically allowed, the charged multiplets of the Higgs bosons can be produced from Z ′ in pair and each charged Higgs can further decay in to an RHN and charged lepton.The RHNs can decay into SM gauge bosons and leptons.Depending on the hadronic and leptonic decays of the SM gauge bosons originated after the RHN decay, we obtain multilpeton signature containing five or six leptons in the final state/s in association with missing momentum and jets coming from the SM gauge bosons.Such rare multilepton mode has not been explored in the context of U(1) X scenario before.Additionally there is another interesting fact in this model due to the presence of general U(1) X symmetry.The left and right handed fermions differently couple with the Z ′ due to the presence U(1) X gauge symmetry.As a result here the Z ′ has chiral nature and it can be manifested when Z ′ interacts with quarks and leptons.As a result the charged Higgs pair production from the left handed and right handed quarks will be different.In addition to that the charged Higgs pair production cross section from Z ′ will depend on the U(1) X charges of the quarks and the charged Higgs.Due to the dependence on U(1) X charges, such a chiral behavior can be manifested in charged Higgs production which has not been explored in any previous literature.Depending on the U(1) X charges, the branching ratio of the Z ′ into charged Higgs can be maximized comparing with that of Z ′ into SM charged leptons.The general U(1) X gauge coupling can be constrained by the dilepton, dijet searches from the LHC [48][49][50][51], ATLAS technical design report (TDR) [52] and LEP-II [53][54][55].
The article is arranged as follows.We discuss the model in Sec.II.Signal and backgrounds of the charged Higgs pair production and its decay into multilepton modes through the RHNs at the LHC have been studied in Sec.III.Finally we conclude the article in Sec.IV.

II. MODEL
The U(1) X extension investigated in this article consists of three generations of SM-singlet RHNs with non-universal charge assignments helping to cancel gauge mixed gauge gravity anomalies.First two generations of the RHNs have charge as −4 and third generation has U(1) X charge as 5. Due to the non-universal charge assignment the scalar sector of the SM needs to be extended with one SU(2) L doublet (H 2 ) and three singlet SM-singlet scalers (φ A,B,C ).Depending on the U(1) X charge assignments of the fermions in the model we find that the RHNs interact only with the extended Higgs sector.The particle content and the corresponding U (1) X charges of the particles are given as (x ′ f , f = {q, u, d, ℓ, e}) in Tab.I.The U(1) X charges of the SM fermions are the same for three generations.The general charges can be related to each other from the following gauge and mixed gauge-gravity anomaly cancellation conditions respectively.Using SM ⊗ U (1) X gauge symmetry we write the Yukawa interactions as where we find that N R1,2 interact with the new Higgs doublet (H 2 ) to generate the neutrino Dirac mass term while the SM singlet scalars (φ A ) generate the Majorana mass terms for N R1,2 after the U(1) X symmetry breaking respectively.The Majorana mass term of N R3 is generated from the VEV of φ B .Due to the U (1) X charges no Dirac mass term is generated for N R3 as a result N R3 does not participate in the neutrino mass generation mechanism at the tree level as a result it may be considered as a potential DM candidate.These Yukawa interaction terms impose to obtain the U(1) X charges using the anomaly cancellation given in Eq. 1.We assume a basis in which Y A N of Eq. 2 is diagonal without the loss of generality.
Apart from the Z ′ interactions with the fermions of this model, the scalar sector in this model is also interesting.The scalar potential for two Higgs doublets and three SM singlet scalars under the SU (3) × SU (2) × U (1) Y × U (1) X gauge group can be given by We parameterize the scalar fields as with j = A, B, C, where v 1,2 and v (A,B,C) are the VEVs of the scalar fields in the model satisfying the condition For simplicity we consider the mixed quartic couplings between the scalars to be very small so that the Higgs doublet sector can be effectively separated from the singlet Higgs sector allowing the higher order mixing effect between the RHNs after the U (1) X symmetry breaking to be highly suppressed.Hence the singlet and doublet scalars communicate through the µ Using this term we further impose on the U(1) X charges of Tab.I to obtain the final form.The charge assignment of the particles after the anomaly cancellation is given in Tab.I where the charge assignment for the SM fermions are linear combination of the U(1) Y and U(1) B−L .We find that for changing x H = −2, − 1, − 0.5, 0.5, 1, 2 the left and right handed charges under the U(1) X differ which is observed in their interactions with Z ′ manifesting its chiral nature.In case of x H = −2 the left handed fermions do not interact with Z ′ which is an alternative U(1) R scenario.For x H = −0.5 and 1 right handed up and down type quarks do not interact with Z ′ .Finally for x H = −1, the right handed electron do not interact with Z ′ .For x H = 0 we notice that the fermion sector and H 1 in the model manifest the B−L charge assignment.Therefore we call it an alternative B−L scenario.We consider the collider constraints the triple coupling has no significant effect on determining v 1 , v 2 , v c when we arrange the parameters in the scalar potential to have v A ∼ v B ∼ v C and µ < v A .When φ C generates VEV, then mixing mass term between H 1 and H 2 are generated.Hence the potential of the Higgs doublet sector effectively becomes the potential for the Higgs potential of the two Higgs doublet model.Since there is no mixing mixing mass term among φ A,B,C , there exists two physical Nambu-Goldstone (NG) modes in our model which are not physically dangerous.We fix the SM singlet scalars to be heavier than Z ′ boson so that Z ′ boson can not decay into the NG modes.

A. Scalar Sector
The vacuum expectation values (VEVs) can be obtained by minimizing the scalar potential given by ∂V /∂v 1 = ∂V /∂v 2 = ∂V /∂v j = 0 where j = A, B, C. Minimizing the scalar potential given in Eq. 4 using the scalar field from Eq. 5 we obtain the following stationary conditions Using the minimization conditions, the mass matrix for the charged scalars is given by which when diagonalized to obtain the mass of the charged Higgs (H ± ) and massless charged Goldstone boson G ± such that and we transform the charged counterparts in the following way where The mass matrix of the CP-odd scalar sector is a 3×3 symmetric matrix which can be diagonalized to obtain two neutral Goldstone bosons G 1,2 which are eaten up by Z, Z ′ gauge bosons and A as the CP-odd neutral scalar.The CP-even scalar sector has five physical degrees of freedom (φ 3,7,A R ,B R ,C R ) with a 5 × 5 symmetric matrix which can be diagonalized numerically to obtain five CP even scalars χ 1,2,3,4 and h.We identify h as the lightest observed Higgs boson of mass 125 GeV [56,57].The mass splitting between charged Higgs and other neutral scalars can lead to interesting scenarios for example the bosonic decay modes of charged like H ± → χ 1,2,3,4 (h, A)W + can overcome the neutrinophilic behavior of charged Higgs.Several studies of bosonic decays of charged Higgs are done in two Higgs doublet model [58][59][60][61][62][63][64].However, we chose the simplest scenario where these bosonic decays are suppressed either kinematically or by suitable mixing angles so that the neutrinophilic nature of charged Higgs gets prominence.We refrain further study on the scalar sector assuming only the neutrinophilic scenario of charged Higgs.In this context we write down the constraints from the stability of the scalar potential as [65,66] and the conditions for perturbativity and unitarity as [66][67][68] where a 1,2,3 are solutions of the cubic equation

B. Gauge Sector
The kinetic terms of the scalar fields are where the covariant derivatives of the scalar fields can be written as and X H1 , X H2 , X φ A , X φ B , X φ C are the U (1) X charges of the the fields H 1 , H 2 , φ A , φ B , φ C respectively as given in Table I.Here g, g ′ and g ′′ are the gauge couplings associated with the SU (2) L , U (1) Y and U (1) X gauge symmetries respectively.The mass of W ± gauge bosons are given by m 2 W ± = 1 4 g 2 v 2 and the mass matrix for the neutral gauge boson is given by where we defined the parameters a µ , B µ ) T by Weinberg angle to obtain the massless photon A µ such that where cos θ W = g/ g 2 + g ′2 and sin θ W = g ′ / g 2 + g ′2 as in SM.We obtain the mass term of the neutral gauge bosons in the basis of ( Zµ , B ′ µ ) basis as The mass matrix can be diagonalized by rotation matrix to obtain the neutral gauge bosons (Z µ , Z ′ µ ) Zµ The physical masses of the neutral gauge bosons are In the limit ∆ 2 → 0 we approximate m Z ≈ M Z,SM and m Z ′ ≈ M Z ′ .Here we have ignored the kinetic mixing of U (1) Y and U (1) X gauge bosons assuming that the mixing effect is small [69].

C. Neutrino Mass
The mass term for the neutral fermions coming from the Yukawa Lagrangian in Eq. 2 after the U (1) X symmetry is spontaneously broken.Hence the neutrino mass term can be written as where where the m ij D is the Dirac mass term generated after the electroweak symmetry breaking and m N is the Majorana mass term generated after the U (1) X symmetry breaking.Finally we write the neutrino mass matrix as where the first two generations of RHNs, N R1,2 are involved whereas the third generation, which can be treated as a where U is the PMNS mixing matrix to leading order (if we ignore the non-unitarity effects for simplicity).Hence the charged-current (CC) interaction in the lepton sector is then given by where g is the SU (2) L gauge coupling and P L = (1 − γ 5 )/2 is the left-chiral projection operator.Similarly, the neutral-current (NC) interaction is given by where θ w is the weak mixing angle.In this analysis we consider a neutrinophilic scenario.The small neutrino mass originated from the seesaw mechanism requires the Dirac Yukawa coupling to be small which further requires tan β ≤ 10 −3 to reproduce the neutrino oscillation data.When we consider the RHNs are heavier than the SM bosons so that they can decay into ℓW , ν ℓ Z, and ν ℓ h on-shell modes.The corresponding partial decay widths are given by respectively where R αi is the corresponding light-heavy mixing element.The corresponding branching ratios are shown in Fig. 1.

III. SIGNAL AND BACKGROUND ANALYSIS
In this model the charged Higgs (H ± ) couples with the SM gauge bosons, Z ′ through the scalar covariant derivatives given in Eq. 15 which involve the corresponding U (1) Y and U (1) X charges.At the LHC we consider pair production of H ± from the neutral gauge bosons.Due to the neutrinophilic nature of H 2 , the charged Higgs H ± decays into the ℓ ± N mode.Considering the decays of the N we obtain 5 and 6 charged lepton final states in association with jets and/or missing momentum.

A. Constraints on the model parameters
In the context of our model scenario the charged Higgs interaction with the quark sector can be given by The interaction is similar to Type-I two Higgs doublet model with cot β replaced with tan β.The low values of tan β, i.e., v 1 < v 2 are more favorable to satisfy the seesaw hierarchy by explaining the small neutrino mass.As a result the charged Higgs comes dominantly from the second Higgs doublet and the quark sector Yukawa couplings with charged Higgs are highly suppressed.In LHC the searches of charged Higgs are motivated by QCD induced H ± tb associated production.Thereafter H ± decays into fermionic or bosonic channels.The suppressed coupling/s of H ± to the quark sector reduces the strong production of charged Higgs.Hence there are no exclusion limits on H ± parameter space.With similar arguments we can evade the flavor physics constraints like b → sγ [70].
For charged Higgs heavier than top quark, the fermionic decay mode of the charged Higgs manifests a dominant decay into H ± → tb and H ± → ℓN i channels.Other than the fermionic channel, charged Higgs can decay into the bosonic channels W ± h, χ 1,2,3,4 , A. The mass degeneracy between the H ± , χ 1,2,3,4 and A due to the T -parameter constraint restricts the decay H ± → W ± χ 1,2,3,4 , A. By suitable choice of the mixing angles α 1,2,3 we completely turn off H ± → W ± h decay mode.In our analysis we consider charged Higgs decay only into the fermionic sector.From Eq. 2, we get the charged Higgs interaction with leptons and right handed neutrinos as Depending on the choice of the tan β and m H ± charged Higgs decays into tb mode or τ ν mode along with the ℓ ± N signature.For the low tan β the second term of Eq. 29 dominates over the other decay modes since it is proportional to cos β.Hence charged Higgs leptonic interaction is dominant into leptons and heavy neutrinos (N 1,2 ).The charged Higgs-quark coupling is proportional to tan β hence H ± becomes neutrinophilic for tan β << 1.The partial decay width of the charged Higgs into the tb channel is given as where N C (= 3) is the color factor, m t is the top quark mass, λ(x, y) 1/2 = (1 + x 2 + y 2 − 2x − 2y − 2xy) 1/2 is the kinematic factor, G F is the Fermi constant and V tb is the CKM matrix element respectively.For the neutrinophilic region the decay width of charged Higgs into charged lepton and heavy neutrino is where i = 1, 2, 3 and j = 1, 2. In Fig. 2 we show the branching ratio comparison between these two channels for different benchmark points.We find that when m H ± is less than top quark mass (m t ) (left panel), then m H ± can decay into ℓ ± N and τ ± ν modes where ℓ ± N mode is dominant when tan β < 0.1.Similar behavior is found in the case where m H ± is greater than m t (right panel) where dominant production of the ℓ ± N mode can be obtained for tan β < 0.007.Hence the charged Higgs becomes neutrinophilic below this limit.Hence in the further analysis we consider tan β = 10 −4 .In the neutrinophilic region H ± will completely decay into the ℓ ± N mode.The partial decay width of Z ′ → H + H − in the limit m Z ′ >> m Z and sin θ → 0 is given by The partial decay widths of the Z ′ into a pair of charged fermions and heavy neutrinos can be given by and respectively.The corresponding branching ratios of the Z ′ are shown in the left panel of Fig. 3 considering negligible Z − Z ′ mixing taking [69] into account.We find that Z ′ can dominantly decay into a pair of heavy neutrinos and the next dominant decay mode can be considered as Z ′ → H + H − respectively.In the right panel we find that BR(Z ′ → H + H − ) is dominant over the dilepton decay mode from Z ′ at x H = −1.24which is considered to study the

FIG. 3: Branching ratio of Z ′ into different modes (left). The ratio of BR(Z
charged Higgs pair production.As a result if we consider x H = −1.24, the Z ′ couplings with the left and right handed SM fermions will be different showing the chiral nature of Z ′ .In addition to that branching ratio of Z ′ into charged Higgs pair will be 30 times more than that of the dilepton mode.As a result the dominant constraints from the dilepton search at the LHC will be less constrained for x H = −1.24.Fixing x H = −1.24we find the U(1) X charge of the quark doublet is 0.54, that for the right handed up and down type quarks are 1.16 and -0.08 respectively whereas the U(1) X charges for the left handed and right handed leptons are −1.16 and −2.24 respectively.The U(1) X charge of H 1 and H 2 are −0.62 and 2.38 respectively.Hence the chiral nature of the model persists for the charge assignment when charged Higgs pair production dominates over the dilepton mode.
To estimate these bounds from the LHC we have considered the case when all the RHNs are heavier than M Z ′ 2 which are shown by the dashed lines by solid Red, solid Green and dashed Blue to represented the LHC dilepton searches from ATLAS-dilepton, CMS-dilepton and ATLAS-TDR (2-electron) respectively.Due to this fact the Z ′ can not decay into the RHNs which will enhance the other decay modes of Z ′ allowing the strongest bound on g ′′ from the dilepton.We use narrow width approximation at 139(140) fb −1 luminosity at the LHC where the Z ′ production cross section is proportional to g ′′ 2 .Hence we estimate the bounds on g ′′ using where g ′′ Model is the coupling considered to estimate the dilepton production cross section from our model.At the High Luminosity LHC (HL-LHC) of 3000 fb −1 we can scale the limits from ATLAS (CMS) considering M N > M Z ′ 2 case and obtain that the limit may get uniformly 0.215 (0.216) times stronger than the current luminosity following at a future luminosity of L future .The limits are shown in Fig. 4. In the similar way we estimate the projected limits on the g ′′ − m Z ′ plane from the ATLAS technical design report (TDR) [52] and the LEP-II limits for x H = −1.24are estimated following [38].
We show the production cross sections of the charged Higgs pair production at the 14 TeV LHC in our model from the s− channel photon, Z and Z ′ mediated processes including the interference effect.In Fig. 5 we show the charged Higgs pair production at leading order as a function of m H ± for m Z ′ = 3 TeV and 5 TeV considering g ′′ = 0.114 and 0.615 respectively from Fig. 4

B. Signal Background Analyses
In this paper we consider the charged Higgs pair production and the charged Higgs dominantly decays into ℓ ± N mode.We consider each heavy neutrino decays into dominant mode N → ℓW as depicted in Fig. 1 followed by the hadronic decay of one of the W bosons and leptonic decay of the other.Hence we find a final state of 5 charged lepton with two jets in association with missing momentum.On the other hand there is another possibility where both the W bosons decay into leptons giving rise to a completely leptonic final state with 6 charged leptons in association with missing momentum.Hence we concentrate only on the leading decay mode of the heavy neutrinos.Implementing the model in FeynRules [71,72] we generate the events using MadGraph [73,74] applying CTEQ6L [75] parton distribution function fixing the factorization scale µ F as the default MadGraph option followed by the showering, fragmentation and hadronization of the signal and SM backgrounds by the PYTHIA8 [76].Finally the detector simulation of the showered events were performed by Delphes [77].

Five charged lepton final state
First we discuss the 5 charged lepton final state from the RHN pair production induced by H ± .We generate irreducible backgrounds for 5 charged leptons which could be fully muons or fully electrons.In this case we have the events with 5 charged leptons, two jets in association with missing momentum.The irreducible backgrounds are: (i) ZZW ± where Z decays leptonically and W ± decays leptonically, (ii) t tZZ where Z decays leptonically, t t produce W ± with b-jets and one W decays leptonically and the remaining one decays hadronically, (iii) t tW Z: the W and Z decay leptonically and t t produce leptonic decay from the W ± boson with b-jets, (iv) ZZW + W − : the Z bosons decay leptonically followed by the leptonic decay of one W boson whereas the other W boson decays hadronically and (v) ZW ± W ± W ∓ : the Z and W bosons decay leptonically to produce a five lepton final state in association with missing momentum.The jets in this events come from the initial state radiations.In addition to the irreducible backgrounds we mention about a reducible background in the form of ZZ+jets.In this case four charged leptons are coming from the Z bosons and other lepton can be misidentified from the jets.In this case a Z-veto could be applied to reduce this background.In our analysis we did not include N → νZ mode therefore we safely ignore this background.To study the signal and backgrounds at the 14 TeV LHC we impose the following selection cuts: 1. Basic Cuts: We select events with 5 leptons (fully muons/electrons) in the final state with minimum transverse momentum p T (ℓ) > 10 GeV and |η ℓ | < 2.4.Minimum transverse momentum for the accompanying radiated jets to be, p T (j) > 20 GeV and |η j | < 2.5.

Momentum cuts:
To reduce the backgrounds further we impose a hard cut on the transverse momentum of the leptons.We impose p T (ℓ 1,2 ) > 30 GeV and p T (ℓ 3,4,5 ) > 20 GeV respectively where leptons are ordered according to the transverse momentum.Along with that we impose a missing momentum cut of p miss T > 30 GeV to reduce events with jets as fake leptons further.
Additional cuts like b-veto and Z-veto reduce the top quark and Z boson mediated backgrounds significantly, however since the backgrounds are already negligible due to the low cross sections, the above mentioned cuts are sufficient enough.The background cross sections based on the above cuts are mentioned in Tab.II for the case of 5 fully muons/electrons final states.For the signal we choose the benchmark parameters as m Z ′ = 3 TeV and 5 TeV, m H ± = 150 GeV and 650 GeV with m N = 100 GeV and 300 GeV respectively.We recognize Benchmark Points (BPs): BP1 as m H ± = 150 GeV, m N = 100 GeV and BP2 as m H ± = 650 GeV, m N = 300 GeV respectively for each choice of m Z ′ .The cut flow for the SM backgrounds and the corresponding signals are given in Tab.III where we estimate the combined significance for the signal of the corresponding BPs using where s and b stand for signal and background events at 3000 fb −1 luminosity.Studying the signal and background processes over grids of 150 GeV ≤ m H ± ≤ 1.1 TeV and 100 GeV ≤ m N ≤ 1 TeV we estimate the significance contours for the combined 5 charged leptons in Fig. 6 at 3000 fb −1 luminosity.We show the contours for m Z ′ = 3 (5) TeV in the left (right) panel.The boundary of the light blue region shows 5σ significance and the light blue region represents m H ± and m N which could be probed with a significance more than 5σ but below 7σ where the light cyan region represents m H ± and m N which could be probed with a significance more than 7σ but below 10σ.The boundary of the blue shaded region represents m H ± and m N which could be probed with 10σ, however, the blue region represents m H ± and m N with more than 10σ.From Tab.III we find that for m Z ′ = 3 TeV BP1 and BP2 reside in the blue and cyan regions respectively.Similar results are obtained for m Z ′ = 5 TeV.Above the black line H ± → ℓ ± N is kinematically suppressed.

Six charged lepton final state
We discuss the 6 charged lepton final state from the heavy neutrino pair production induced by H ± and followed by the dominant N → ℓW mode where W bosons decay leptonically.In addition to that we consider the N → νZ mode from both the heavy neutrinos followed by leptonic decay of the Z boson.Due to the lepton multiplicity, the 6 charged lepton final state will have small cross section.Therefore N → Zν mode is added to boost the signal because BR(N → Zν) is roughly 25% for m N ≥ 141 GeV which can be found in Fig. 1.
We generate irreducible backgrounds for 6 charged leptons which could be fully muons or fully electrons.Here we have the events with 6 charged leptons in association with missing momentum where jets may appear from the initial state radiations.The irreducible backgrounds are: (i) ZZZ where three Z bosons decay leptonically, (ii) t tZZ where Z decays leptonically, t t produce W ± with b-jets and the W bosons decay leptonically, (iii) t tW W Z: the W and Z decay leptonically and t t produce W ± in association with b−jets and finally the W ± bosons decay leptonically, (iv) ZZW + W − : the Z and W bosons decay leptonically.In this analysis additionally we consider a reducible background coming from ZZ+jets.In this case 4 charged leptons are coming from the Z bosons and the remaining ones come from the misidentified jets.Identification of two mistagged leptons from jets is negligibly small.Hence without Z-veto, the ZZ+jets background will be very small for 6 charged lepton final state.Hence we ignore this background.To study the signal and backgrounds at the 14 TeV LHC we impose the following selection cuts: 1. Basic Cuts: We select events with 6 leptons (fully muons/electrons) in the final state with minimum transverse momentum p T (ℓ) > 10 GeV and |η ℓ | < 2.4.Minimum transverse momentum for the accompanying radiated jets to be, p T (j) > 20 GeV and |η j | < 2.5.

Momentum cuts:
To reduce the backgrounds further we impose a hard cut on the transverse momentum of the leptons.We impose p T (ℓ 1,2 ) > 30 GeV for the leading leptons and p T (ℓ 3,4,5,6 ) > 20 GeV for the trailing leptons respectively where leptons are ordered according to the transverse momentum.Along with that we impose a missing momentum cut of p miss T > 30 GeV to reduce events with jets as fake leptons further.V where we estimate the significance of the signal using Eq.37 at 3000 fb −1 luminosity.
Studying the signal and background processes over grids of 150 GeV ≤ m H ± ≤ 1.1 TeV and 100 GeV ≤ m N ≤ 1 TeV we estimate the significance contours for the combined 6 charged leptons in Fig. 7 at 3000 fb −1 luminosity.We show the contours for m Z ′ = 3 (5) TeV in the left (right) panel.The boundary of the light blue region shows 3σ significance and the light blue region represents m H ± and m N which could be probed with a significance more than 3σ but below 5σ where the light cyan region represents m H ± and m N which could be probed with a significance more than 5σ but below 7σ.The boundary of the blue shaded region represents m H ± and m N which could be probed with 7σ, however, the blue region represents m H ± and m N with more than 7σ.From Tab.V we find that for m Z ′ = 3 TeV BP1 and BP2 reside in the blue and cyan regions respectively.Similar results are obtained for m Z ′ = 5 TeV.Above the black line H ± → ℓ ± N is kinematically suppressed.

IV. CONCLUSION
In this paper we consider a scenario where SM is extended by a general U(1) X gauge group.To cancel the gauge and mixed gauge gravity anomalies we introduce three generations of RHNs with nonuniversal charge assignments under U(1) X gauge group.Due to the U(1) X gauge symmetry the RHNs interact with the extended Higgs sector of the SM giving rise to the seesaw mechanism to generate the light neutrino mass after the U(1) X symmetry is broken.In this scenario a neutral gauge boson Z ′ is evolved where the left and right handed fermions, being charged differently under general U(1) X , interact with the Z ′ manifesting the chiral nature of the model.Thanks to the general   V: Signal cross sections and significance for the 6 charged lepton events at √ s = 14 TeV LHC after each selection cuts for the benchmark points BP1 and BP2 at a luminosity of 3000 fb −1 .U(1) X charge assignment which is responsible for the dominant decay mode of Z ′ into a pair of charged Higgs over its dilepton decay mode.Hence for that U(1) X charge assignment, the experimental bounds on the gauge coupling for different M Z ′ obtained from the dilepton searches at the LHC and LEP-II will be comparatively loose.It will help to enhance charged Higgs pair production cross section.Meanwhile, in the limit of low tan β we find that the second Higgs becomes neutrinophilic and its charged multiplet completely decays into a charged lepton and RHN.Producing the charged Higgs in pair and studying its decay into RHNs we obtain rare multilepton final states.Here we consider the decay of the RHNs into the dominant mode ℓW .Following the hadronic and/or leptonic decay of the W boson we obtain five or six charged leptons in the final state.This leptonic multiplicity helps us to probe such rare multilepton final states from the RHNs.Comparing the signals and irreducible backgrounds for different regions in m H ± − m N plane we obtain significantly high discovery potential for these final states at the 14 TeV LHC with 3000 fb −1 luminosity.The six charged lepton final state is weaker than the five charged lepton final state beacsue leptonic decay of the W boson reduces cross section.We comment that in this article we have studied purely electron and muon signals.However, it is possible to study electron and muon both in the final state which may come from the charged Higgs, heavy neutrino and the leptonic decay of the W boson.This will definitely increase the significance of the signal over the backgrounds, however, we already have a clean final state due to lepton multiplicity.Hence due to simplicity we did not consider the mixed signal.In future several interesting aspect of including collimated objects, displaced scenarios could be tested from the charged Higgs induced RHN productions in the context of hadron colliders.The reconstruction of the Z ′ boson will be another interesting aspect at the high energy colliders where all the final states are visible involving leptons and jets.
Particle content of the alternative U(1) X model with general U(1) X charges before and after anomaly cancellation, and α = 1, 2, 3 stands for three generations of the fermions.The charges of the RHNs are non-universal.Considering different benchmark values of the x H we obtain different U(1) X charges of left and right handed fermions of the model manifesting the chiral nature of the model.Here, x H = 0 is an alternative B−L case, which is a vector-like scenario.
potential dark matter candidate, gets the mass m B .Assuming the seesaw hierarchy m N >> m D and diagonalizing Eq. 23 to the light neutrino mass as m ν ≃ −m D m −1 N m T D where we define R = m D m −1 N as the mixing between the light and heavy neutrinos.Due to the light-heavy neutrino mixing, a light neutrino flavor eigenstate (ν α ) is a linear combination of the light (ν m ) and heavy (N m ) neutrino mass eigenstates

FIG. 2 :
FIG. 2: Branching ratios of the charged Higgs into different modes as a function of tan β with m N1,2 = 100 (300)GeV in the left (right) panel.

24 FIG. 4 :
FIG.4: Limits on the U(1) X coupling as a function of m Z ′ where the shaded region is ruled by the existing experimental data.

FIG. 5 :
FIG. 5: Leading order H ± production cross section at the 14 TeV LHC as a function of m H ± for m Z ′ = 3 TeV and 5 TeV respectively considering x H = −1.24,tan β = 10 −4 and m N = 1 TeV.The cross section without the influence of Z ′ falls rapidly at high m H ± .

FIG. 6 :
FIG.6: Different significance contours for the 5 charged lepton final state with 3000 fb −1 luminosity.The black straight line represents m H ± = m N , below which the on-shell decay of H ± occurs.

3 .
Lepton separation cuts: We impose the lepton-lepton separation in the η − ϕ plane as ∆R ℓℓ > 0.4 and if there are any accompanied jets, then ∆R ℓj > 0.4 and ∆R jj > 0.4 have been considered where ∆R = ∆η 2 + ∆ϕ 2 .Additional cuts like b-veto and Z-veto reduce the top quark and Z boson mediated backgrounds significantly, however since the backgrounds are already negligible due to the low cross sections, the above mentioned cuts are sufficient enough.The background cross sections based on the above cuts are mentioned in Tab.IV for the case of 5 fully muons/electrons final states.For the signal we choose the benchmark parameters as m Z ′ = 3 TeV and 5 TeV, m H ± = 150 GeV and 650 GeV with m N = 100 GeV and 300 GeV respectively.We recognize Benchmark Points (BPs): BP1 as m H ± = 150 GeV, m N = 100 GeV and BP2 as m H ± = 650 GeV, m N = 300 GeV respectively for each choice of m Z ′ .The cut flow for the SM backgrounds and the corresponding signals are given in Tab.

3 FIG. 7 :
FIG.7: Different significance contours for the 6 charged lepton final state with 3000 fb −1 luminosity.The black straight line represents m H ± = m N , below which the on-shell decay of H ± occurs.

TABLE II :
Cross sections for the backgrounds of 5 charged leptons at √ s = 14 TeV LHC after each selection cuts.

TABLE III :
Signal cross sections and significance for the 5 charged lepton events at √ s = 14 TeV LHC after each selection cuts for the benchmark points BP1 and BP2 at a luminosity of 3000 fb −1 .

TABLE IV :
Cross sections for the backgrounds of 6 charged leptons at √ s = 14 TeV LHC after each selection cuts.