Hunting for heavy composite Majorana neutrinos at the LHC

We investigate the search for heavy Majorana neutrinos stemming from a composite model scenario at the upcoming LHC Run II at a center of mass energy of 13 TeV. While previous studies of the composite Majorana neutrino were focussed on gauge interactions via magnetic type transition coupling between ordinary and heavy fermions (with mass $m^*$) here we complement the composite model with contact interactions at the energy scale $\Lambda$ and we find that the production cross sections are dominated by such contact interactions by roughly two/three orders of magnitude. This mechanism provides therefore very interesting rates at the prospected luminosities. We study the same sign di-lepton and di-jet signature ($pp \to \ell\ell jj$) and perform a fast detector simulation based on Delphes. We compute 3$\sigma$ and 5$\sigma$ contour plots of the statistical significance in the parameter space ($\Lambda,m^*$). We find that the potentially excluded regions at $\sqrt{s} =13$ TeV are quite larger than those excluded so far at Run I considering searches with other signatures.


I. INTRODUCTION
The recent discovery [1,2] of the Higgs boson at the CERN Large Hadron Collider (LHC) has certainly crowned in the most spectacular way an almost halfcentury long history of successes of the standard theory of the electroweak interactions, the so called standard model (SM).
In spite of tremendous efforts put in by the experimental collaborations working in the LHC experiments the hunt for new physics (supersymmetry, compositeness, extra dimensions etc..) has so far been unsuccessful. Impressive and increasingly stringent new bounds on the scale of several beyond the standard model (BSM) scenarios are continually being reported.
In this paper we propose to study the like sign di-lepton and di-jet (eejj) signature from a gauge model [3] with an hypotetical heavy Majorana neutrino within the well known scenario of compositeness of quarks and leptons, complemented here with contact interactions.
In this scenario the heavy excited states (q * , e * , ν * ) couple, through gauge interactions, with the ordinary SM fermions via magnetic type couplings. Current bounds on excited lepton masses (generically indicated by m * ) have been recently strengthened by the LHC Run I analyses [4,5] of the γ signature arising from * production (pp → * ), via four fermion contact interactions with a compositeness scale Λ, followed by the decay * → γ. In particular in [4] the ATLAS Collaboration reporting an analysis at √ s = 8 TeV with an integrated luminosity of 13 fb −1 gives a lower bound on the mass of the excited leptons m * > 2.2 TeV (derived within the hypothesis m * = Λ). In [5] the CMS Collaboration reported the results of data collected with 19.7 fb −1 at √ s = 8 TeV and (always assuming m * = Λ) excluded excited electron (muon) masses up to 2.45 (2.48) TeV. Preliminary studies within the compositeness scenario of the like-sign di-lepton and di-jet signature were performed long ago [3], assuming the excited neutrino ν * = N to be a Majorana particle. Here our aim is to complement the composite Majorana neutrino model of ref. [3] with contact interactions which are again a generic expectation of a composite fermion scenario [6]. Based on previous studies related to the production at LHC of exotic doubly charged leptons [7] we expect these contact interactions to be the dominant mechanism for the resonant production of the heavy Majorana neutral particles N in the process pp → N . This expectation is indeed verified by our numerical simulations performed with a custom version of CalcHEP [8,9] where our model has been implemented. The heavy Majorana neutrino is produced resonantly in association with a lepton (pp → N ) and then given the relatively important branching ratio for the decay N → jj we perform a detailed kinematic study of the like-sign di-lepton and di-jet final state: including the relevant standard model backgrounds.
Our study shows clearly that a full fledged analysis of the upcoming data from the Run II of LHC at √ s = 13 TeV has the potential of observing the signature or alternatively excluding larger portions of the model parameter space compared to those already excluded from analyses of Run I [4,5].
We remark however that the CMS Collaboration has recently reported an excess over the SM background expectations in the eejj and ep T / jj final states where p T / is the missing transverse momentum. The analysis in [10] for a search of right-handed gauge boson, W R , based on 19.7 fb −1 of integrated luminosity collected at a center of mass energy of 8 TeV reports a 2.8σ excess in the eejj invariant mass distribution in the interval 1.8 TeV < M eejj < 2.2 TeV. A CMS search [11,12] for first generation lepto-quarks at a center of mass energy of 8 TeV and 19.6 fb −1 of integrated luminosity reported an excess of, respectively 2.4σ and 2.6σ in the eejj and ep T / jj channels. Several scenarios have been proposed in the literature to explain the above CMS excesses in the context of various standard model extensions. For instance in [13,14] the authors propose an explanation of the excesses in the context of W R decay by embedding the conventional left-right symmetric model (LRSM) (g L = g R ) in the SO(10) gauge group. Studies of the eejj excess in the context of W R and Z gauge boson and heavy neutrinos (N ) -coupling mainly to electrons-production and decay appear in [15][16][17]. Similarly [18] discusses a model with pseudo-Dirac heavy neutrinos providing a fit to all excesses in a generic LRSM with arbitrary g R , W -W R boson mixing, heavy neutrino N and the ν-N mixing. In addition, the authors point out the consequences of the excesses for neutrino-less double beta decay 0νββ decay, and find for example that 0νββ actually provides a pretty severe limit on the ν-N mixing assuming the excesses are real.
Other interpretations have been proposed within the context of models with vector-like leptons as in [19] showing that resonant pair production of such vector-like leptons decaying to an electron and two jets leads to kinematic distributions consistent with the observed CMS data. The eejj excess has been shown to arise as well in R-parity violating models through spleton resonant production [20][21][22]. An alternative scenario based on lepto-quarks is proposed in [23,24], discussing also possible connections to dark matter, which fits the data of the excess seen by CMS. In [25] the observed CMS excesses are explained within superstring inspired E 6 models which can also accommodate for the baryon asymmetry of the universe via lepto-genesis. Other studies have emphasized that the observed differences between the eejj, µµjj, same sign (SS) and opposite sign (OS) channels could be addressed including mixing and CP phases of the heavy neutrinos [26]. On the other hand it is well known that the like-sign di-lepton and di-jet (eejj), ∆L = 2 violating final state (Keung-Senjanovic process), is the golden signature to look for heavy Majorana neutrinos at high energy hadron collisions [3,[27][28][29][30][31][32][33][34][35]. Studies of heavy (pseudo-Dirac) neutrino production at the LHC within the inverse see-saw mechanism have been performed [36], also considering the quark-gluon fusion mechanism [37,38].
We show that our heavy composite majorana neutrino model, in its simplest version can reproduce, at least qualitatively, some features of the observed excess in the eejj invariant mass distribution. We discuss how, with some refinement, it has the potential to address also other aspects of the excess, such as the absence of a peak in the second-leading-electron−jj invariant mass distribution, the charge asymmetry of the excess, and the fact that the same excess is not observed in the µµjj channel [10].
The rest of the paper is organized as follows: In Sec. II we review the theoretical aspects of the composite model; in Sec. III we discuss the heavy neutrino production cross sections and decay rates; in Sec. IV we discuss the samesign di-lepton and di-jet signature and the main associated SM backgrounds; in Sec. V we present the results of the fast simulation obtained through the Delphes [39] software; finally Sec. VI gives the conclusions with outlooks.

II. COMPOSITE MODEL(S) WITH GAUGE AND CONTACT INTERACTIONS
In this section we review the composite model of excited fermions investigated in [3] within the hypothesis of a heavy Majorana neutrino. Compositeness of ordinary fermions is one possible scenario beyond the standard model. In this approach quarks and leptons are assumed to have an internal substructure which should become manifest at some sufficiently high energy scale, the compositeness scale Λ. Ordinary fermions are then thought to be bound states of some as yet unobserved fundamental constituents generically referred to as preons. While building a fully consistent composite theory has proven to be quite difficult some, important and model independent features of the compositeness scenario can be phenomenologically addressed. Quite natural properties of this picture are [40][41][42]: (i) the existence of excited states of such low lying bound states of preons q * , e * , ν * ... with masses m * ≤ Λ; and (ii) contact interactions between ordinary fermions and also between ordinary and excited fermions. Let us consider here the various possible composite models with respect to the idea of introducing lepton number violation (LNV) via a composite Majorana neutrino. (a) Homo-doublet model. The homo-doublet model [43,44] contains a left handed excited doublet along with a right handed excited doublet: Typically the left and right handed doublet are assumed to have the same mass. It is known that two left and right Majorana fields with the same mass combine to give a Dirac field (with a Dirac mass) [45]. The homo-doublet model, as laid out, cannot therefore accommodate Majorana excited neutrinos, and hence lepton number violation (LNV). This becomes possible if one is willing to introduce a mass difference between the left and right doublet (ν * L −ν * R mixing) or, in other words, a breaking of the L-R symmetry. Such a possibility has been discussed for instance in ref. [46] where the ν * is possibly a linear combination (with mixing coefficients) of Majorana mass eigenstates.
On the other end, if we do not want to introduce a mass splitting (or mixing) between the left and right components in the homo-doublet model, we can account for LNV advocating different models within the compositeness scenario which naturally can accommodate a Majorana neutrino [46,47]. These are the following: (b) Sequential type Model. The sequential model contains excited states whose left handed components are accommodated in SU(2) doublets while the right handed components are SU(2) singlets: and the notation [ν * R ] means that ν * R is necessary if the excited neutrino is a Dirac particle while it could be absent for a Majorana excited neutrino. The magnetic type interactions in this case can be constructed by coupling the left-handed excited doublet to the SM fermion singlets via the Higgs doublet [46]. This results in coupling strengths suppressed by a factor v/Λ [47] where v ≈ 246 GeV is √ 2 times the vacuum expectation value of the Higgs field. (c) Mirror type Model. It is assumed to contain a right handed doublet and left handed singlets: where we may assume that there is no left handed excited neutrino (ν * L ) so that we can associate to ν * R a Majorana mass term and ν * is a Majorana particle. This model is described by a magnetic type coupling between the lefthanded SM doublet and the right-handed excited doublet via the SU (2) L × U (1) Y gauge fields [46,47]: where L T = (ν L , L ) is the ordinary SU (2) L lepton doublet, g and g are the SU (2) L and U (1) Y gauge couplings and W µν , B µν are the field strength for the SU (2) L and U (1) Y gauge fields; f and f are dimensionless couplings which are typically assumed to be of order unity.
The relevant charged current (gauge) interaction of the excited Majorana neutrino N = ν * is then: The above mirror type model is therefore the model to which we will refer our detailed simulation of the like sign di-lepton signature at the Run II of the LHC.
Incidentally we note that SM extensions involving mirror fermions have been recently considered [48] with respect to the phenomenology of the production of mirror quarks at the LHC.
At last we may add that one could also consider extended isospin composite models [49] where the excited states are grouped in triplets (I W = 1) or quadruplets (I W = 3/2) instead of doublets (I W = 1/2) as considered above. Such extensions of the composite scenario contain exotic charge states like doubly charged leptons and quarks of charge Q = (5/3)e. Some phenomenology of these extensions involving the doubly charged leptons has been addressed recently [7,50]. Such extended weak isospin composite models could also be considered with the additional hypothesis that the excited neutrino is a Majorana particle.
The model with I W = 1 can only couple [49] the triplet T = (L 0 , L − , L −− ) with the right-handed lepton singlet R . Therefore we may assume a sequential type structure with a left-handed triplet and right-handed singlets. If the L 0 R is missing we may assume for the L 0 L a Majorana mass term and so the excited neutral L 0 of the triplet is a Majorana neutrino (N ). The magnetic type interaction reads: where f 1 is an unknown dimensionless coupling in principle different from f appearing in Eqs. (5&6). The relevant charged current interaction of the neutral component of the triplet L 0 is in this case: which differs form the one in Eq. (6) in the chirality of the projection operator. [49] with the left-handed SM doublet, so that assuming a mirror type scenario and that there is no L 0 L we can assign to L 0 R a Majorana mass term so that the L 0 neutral of the quadruplet can be a Majorana neutrino (N ). The magnetic type interaction is [49]: (9) where f 3/2 is an unknown dimensionless coupling in principle different from f, f 1 and C( 3 2 , M |1, m; 1 2 , m ) are Clebsch-Gordan coefficients. In particular in the case of I W = 3/2 the relevant neutrino charged current interaction turns out to have the same structure as in Eq. (6): Therefore the interaction in Eq. (6) effectively describes the charge current interaction of a heavy Majorana neutrino both in the (I W = 1/2) mirror type model or in a composite model with extended weak isospin (I W = 3/2), always of the mirror type, provided that we make the correspondence √ 2f 3/2 / √ 3 = f . Contact interactions between ordinary fermions may arise by constituent exchange, if the fermions have common constituents, and/or by exchange of the binding quanta of the new unknown interaction whenever such binding quanta couple to the constituents of both particles [6,47]. The dominant effect is expected to be given by the dimension 6 four fermion interactions which scale with the inverse square of the compositeness scale Λ: where g 2 * = 4π and the η factors are usually set equal to unity. In this work the right-handed currents will be neglected for simplicity.
The single production qq → N proceeds through flavour conserving but non-diagonal terms, in particular with currents like the third term in Eq. 11b which couple excited states with ordinary fermions: which were not considered in [3,50] while are now fully implemented in our simulations.
The Feynman rules corresponding to the Lagrangians in Eqs. (6,12) have been derived with FeynRules [51], a Mathematica [52] package which allows to derive the Feynman rules of a given quantum field theory model once the Lagrangian is given. While the gauge interactions in Eq. (6) where introduced in the CalcHEP [8,9] generator already in [50] and the contact interactions in Eq. (12) were implemented in our CalcHEP model in [7], in this study we have explicitly implemented the Majorana nature of the excited heavy neutrino N assumed in our model.
We conclude this section with one final remark regarding the assumption that in this work the dimensionless couplings f, f , f 1 , f 3/2 are O(1). The production cross sections and all simulations presented in the following are obtained assuming f = f = f 1 = f 3/2 = 1. This should be recalled when quoting the resulting bounds on the other parameters of the model, namely (m * , Λ). In this regard we point out that the cross section yield in the eejj final state cannot be easily rescaled if f = f , f 1 , f 3/2 = 1 because, although the production mechanism is dominated by contact interactions (which do not depend on these constants) the decay of the heavy Majorana neutrino is affected by both gauge interactionsand hence by the factors f = f , f 1 , f 3/2 -, and contact interactions (independent of the f = f , f 1 , f 3/2 ), see the next section for details. A direct comparison with other studies [53,54] which derived bounds on the mixing parameters for the electron flavour of the heavy neutrinos is therefore not possible at the moment. We would need to implement in the generator a model with the additional parameters f, f , f 1 , f 3/2 .

III. CROSS SECTION AND DECAY WIDTH OF THE COMPOSITE MAJORANA NEUTRINO
Heavy Majorana neutrinos N can be singly produced in association with a lepton in pp collisions. The process pp → N can occur via both gauge (Fig. 1, first diagram in the right-hand side) and contact interactions (Fig. 1, second diagram in the right-hand side).
We now present here the production cross section for the heavy Majorana neutrino N in pp collisions expected at the CERN LHC collider stemming from the partonic collisions. Owing to the QCD factorization theorem, the hadronic cross section are given in terms of convolution of the partonic cross sectionsσ(τ s, m * ), evaluated at the partons center of mass energy √ŝ = √ τ s, and the universal parton distribution functions f i which depend on the parton longitudinal momentum fractions, x, and on the factorization scaleQ: (13) For the calculation of the production cross section in proton-proton collisions at LHC, we have used CTEQ6m parton distribution functions [55]. The factorization and renormalization scale has been set toQ = m * .
In Fig. 2 (left) we present the cross section against the heavy Neutrino mass for Λ = 10 TeV for the LHC center of mass energy √ s = 13 TeV. It is evident that the contact interaction dominates the production of the heavy composite Majorana neutrino by a factor that ranges between two and three orders of magnitude, varying the heavy neutrino mass between 1 and 5 TeV, and for the given choice of the compositeness scale (Λ = 10 TeV). In Fig. 2 (right) we compare the cross sections of pp → + + jj with the one of pp → − − jj, for the special case = e. The cross section for the production of positive di-lepton is larger than that for the production of negative di-leptons as expected in proton-proton collisions due to the larger luminosity of a ud pair (needed to produce + + ) compared to that of aūd (needed to produce + + ).
The heavy Majorana neutrino N can decay again through both gauge and contact interactions. The decay amplitudes are related, via appropriate crossing symmetry exchanges, to those describing the single production and depicted in Fig. 1. The possible decays are: In the first we can have a positive lepton, a down-type quark and an up-type antiquark or a negative lepton an up-type quark and a down-type antiquark; in the second owing to the Majorana character of N we can have either a neutrino or an antineutrino of the same flavor of the heavy neutrino N and accordingly two opposite sign leptons belonging to a family that can be the same or different from the other one, or alternatively a positive The dark grey blob describes the production of on shell heavy Majorana neutrinos N in proton-proton collisions at LHC. The production is possible both with gauge interactions (first diagram in the right-hand side) and four fermion contact interactions (second diagram in the right-hand side). For the calculation we used CTEQ6m parton distribution functions and we put the factorization (renormalization) scale tô Q = mN = m * (negative) lepton of the same family of the heavy neutrino and a negative (positive) lepton and an antineutrino (neutrino) belonging to a family that can be the same or different from the other one; in the third we can have a neutrino or an anti neutrino and a quark and an antiquark both of up-type or both of down-type. In Fig. 3 we present the width Γ and the branching ratio B for N → + qq , the decay that gives the final signature under examination of two like-sign leptons and dijet, pp → + + jj which is a signature well known to be rather clean (due to the low expected SM background). Relevant yields are ensured by the rather large B.
This peculiar final state being a lepton number violating process (∆L = +2) is only possible if the heavy neutrino is of Majorana type. In this work we chose to focus on the specific signature with two positive leptons due to its larger cross section as shown in Fig. 2(right). It is important to remark that the like-sign di-lepton plus di-jet signature can be realized by two distinct classes of Feynman diagrams which are shown in Fig. 4. In Fig. 4(a) there is a t-channel exchange of a virtual heavy Majorana neutrino while in Fig. 4(b) the heavy Neutrino is resonantly produced (s-channel) and its subsequent decay. In Fig. 4, each dark (grey) blob includes both a gauge or a contact interaction term whose Feynman diagrams are those shown in Fig. 1.
The process in Fig. 4(a) is the collider analogue of the neutrinoless double-β decay (0νββ), the well known lepton number violating (∆L = ±2) nuclear rare decay [56,57] which, if detected, would unambiguosly verify the Majorana nature of neutrinos. The half life of the 0νββ is currently bounded as T 0νββ 1/2 ≥ 1.1 × 10 25 yr [58] at 90% confidence level, from the data of the 136 Xe EXO (Enriched Xenon Observatory)-200 experiment. Previous searches with 76 Ge (the GERDA experiment) [59] and with 136 Xe (the KamLAND-Zen experiment) [60] had established a half-life longer than 10 25 years. In an high energy collider a heavy Majorana neutrino can be produced in resonance, Fig. 4(b), if the mass of the neutrinos is kinematically accessible m N < √ŝ , where √ŝ is the energy in the parton center of mass frame. In this case the cross section for the signature pp → jj is approximated by σ(pp → jj) ≈ σ(pp → N )B(N → jj). The resonant production rate is dominant relative to the virtual neutrino exchange contribution. This was demonstrated in [31] for the gauge-only case and it is still true in the current model including also the contact interactions. This has been explicitly verified and is shown explicitly in Fig. 5.

IV. SIGNAL AND BACKGROUND
As is well known in the standard model the lepton number L is strictly conserved and thus processes like those in Eq. (1) with ∆L ± 2 are not allowed. However, within the SM there are several processes that can produce two same sign leptons in association with jets. The following processes are considered as main backgrounds [32] pp → tt → + + ννjets , We discuss here the main kinematic differences between the signal and the background to choose suitable cuts for optimizing the signal/background ratio.  transverse plane). The corresponding ∆R distributions all have a peak at low values of (∆R). Therefore, in the reconstruction process, it is possible to have merging, i.e. the two jets can be often reconstructed as a single jet. The same is not true for the electrons. There is a twofold reason for this. On one hand in the production process pp → N , and N are produced in two opposite regions of the transverse plane of the detector, and consequently the same will be true for and the second lepton ( ) from the decay of the heavy neutrino N . On the other hand, it should be kept in mind that the detectors for the LHC experiments can reconstruct the tracks of charged leptons with a very good resolution which certainly warrants the reconstruction of both leptons. The previous considerations can be better understood by the fact that, as shown clearly in Fig. 3, at a given compositeness scale Λ one can identify two regions of the heavy composite neutrino mass m * where the decay (of the heavy composite neutrino) is dominated by gauge interactions (GI) and another where instead it is dominated by contact interactions (CI). Since we assume m * M W , when we are in the region where GI dominate we expect the 2 jets from the W decay to be highly boosted and merged. This effect is expected to be more pronounced as the mass of the heavy neutrino increases, as checked explicitly in Fig. 6 (bottom-left) where the peak at low values of ∆R moves closer to ∆R = 0 for higher values of m * . On the other hand the lepton from the decay of the heavy neutrino is unrelated to the W gauge boson and we thus expect it to be isolated from the jets. When we are in the region where the CI interactions dominate all decay products of the heavy neutrino are produced without being constrained to a particular direction (precisely because the CI are not mediated by a gauge particle), and again the lepton is expected to be isolated from the jet(s). In this case the fraction of events with well separated jets will also increase at higher heavy neutrino masses where the contact interaction dominates, as shown in Fig. 6 (bottom-left) by the peak near ∆R = 3. Incidentally, the above considerations are corroborated by a simple numerical check in the ∆R distribution of Fig. 6(lower-left). For m * = 5000 GeV the peak near zero has a fraction of events (≈ 20%) which is compatible with the ratio of the gauge and contact widths from Fig. 3 (≈ 15%) given the fact that there is a small tail due to contact interactions.
Let us also mention that information about the mass of the heavy Majorana neutrino can be obtained from the invariant mass distribution of the second-leading lepton and the two jets. Indeed in Fig. 6 (bottom-right) we show that this distribution has a very sharp peak corresponding to the heavy Majorana neutrino mass. This is indeed expected since in the resonant production the heavy Majorana neutrino N is decaying to + qq and the lepton from N is expected to be the second-leading, while the leading lepton is the one produced in association with N , in pp → + N .
Finally we show explicitly that the di-lepton plus dijet signature from a heavy composite Majorana neutrino can easily explain the excess observed by the CMS collaboration [10] in the eejj invariant mass distribution in the interval 1.8 TeV < M eejj < 2.2 TeV. Fig. 7 shows, at generator level and for a particular point of the parameter space (Λ = 10 TeV, m * = 1000 GeV), that the eejj invariant mass distribution can easily accommodate an excess in the interval where it has been claimed by the CMS collaboration.
Let us comment here on general expectations about the shape of the eejj invariant mass distribution. Owing to the QCD factorisation theorem, as in Eq. 13, and well known recursive reduction properties of the multiparticle phase-space [61,62] the invariant mass distribution of our eejj final state is easily established to be given by the following relation (note that the jj invariant mass coincides with the energy of the parton center of mass frame, M 2 jj =ŝ): where Q is the QCD factorisation scale and Q is the vir-tual momentum of the resonantly produced heavy neu- trino N . We see that such invariant mass distribution is the product of two factors. One such factor, the first in the right hand side of Eq. 16, is a (dimensionless) parton distribution luminosity factor that vanishes at very large invariant masses M jj ≈ √ s, while the second factor in the right hand side of Eq. 16, is an integral over the virtuality of the produced neutrino, Q, and vanishes for small values of the invariant mass or M jj M N . Therefore in general we expect an invariant mass distribution characterised by a peak for M jj M N . Such picture is of course not altered by the relative importance that contact and gauge interactions may have in the decay process Γ N * → jj (Q), and furthermore the production cross section,σ qaq b → N (M jj , Q), is always dominated by contact interactions.  generator level (shown in Fig. 7) and at the reconstructed level.
We conclude this section by commenting briefly on the fact that the excess is observed in the electron channel but not in the muon channel. This could be explained in our model simply by invoking a rather natural mass splitting between the excited electron (e * ) and muon (µ * ) instead of assuming full degeneracy between the families, i.e. that m e * ≈ m µ * ≈ m * .

V. FAST DETECTOR SIMULATION AND RECONSTRUCTED OBJECTS
In order to take into account the detector effects, such as efficiency and resolution in reconstructing kinematic variables, we interface the LHE output of CalcHEP with the software DELPHES that simulates the response of a generic detector according to predefined configurations [39]. We use a CMS-like parametrisation. For the signal we consider a scan of the parameter space (Λ, m * ) within the ranges Λ ∈ [8,40] TeV with step of 1 TeV and m * ∈ [500, 5000] GeV with step of 250 GeV. For each signal point and each background we generate 10 5 events in order to have enough statistics to evaluate the reconstruction efficiencies ( s , b ) of the detector and of the cuts previously fixed (see Eq. 15a, 15b).
The leptonic flavour of our signature is determined by the flavour of the excited heavy Majorana neutrino: be it either ν * e or ν * µ , (in this work we do not consider a final state with τ leptons due to the production of ν * τ ). In our simplified model characterised by the parameters (Λ, m * ) we are assuming mass degeneracy between the various flavours of excited states. So in principle if we can produce ν * e we can also produce (ν * µ ) and we could have a di-muon and di-jet signature as well. In other words we expect the same number of same-sign di-electrons or di-muons. However our fast simulation of the detector reconstruction is performed only for the electron signature.
To keep our discussion general enough to include both possibilities we use in the text the notation instead of simply ee or µµ. However all results shown (distributions etc...) refer to the electron case which is the one that we have explicitly simulated.
In addition, to be more precise with respect to the hadronic nature of our signature, we may specify that our signal region is defined requiring to have two leptons (electrons) and at least one jet, which means that there may be one or two jets. This selection warrants a very high signal efficiency, regardless of whether there are indeed one or two jets in the reconstructed events.
We then select events with two positive electrons and at least one jet. The number of jets may be just one, in case of merging of the generated two jets, or two, if there is no merging of the generated two jets. Despite the possibility of having a single jet in the event, in the text we will stick with the notations of the main text and will show the results referring to the two jets, coherently with what is produced at the generator level (eejj).
Once we have the number of the selected events we evaluate the reconstruction efficiencies, then for a given luminosity L it is possible to estimate the expected number of events for the signal (N s ) and for the background (N b ) and finally the statistical significance (S): In Fig. 8 Top-left, Top-right and Bottom-left we show the contour plots of S = 3 and S = 5 in the parameter space (Λ, m * ) for three different values of integrated luminosity L = 30, 300, 3000 fb −1 . The regions below the curves are excluded. The colored filled bands are an estimate of the statistical error. In Fig. 8(bottom-right) we compare the 5-σ curves at the three integrated luminosity values. Finally in Fig. 9 we compare our 3-σ contour plots (S = 3) for the three different values of integrated luminosity L = 30, 300, 3000 fb −1 of Fig. 8 with the 95% confidence level exclusion bounds from two Run I analyses at √ s = 8 TeV: ATLAS with 13 fb −1 [4] and CMS with 19.7fb −1 [5]. The shaded regions below the solid, dashed and dot-dashed lines are the current CMS exclusion at √ s = 8 TeV with 19.7 fb −1 of integrated luminosity (blue) [5] and the ATLAS exclusion at √ s = 8 TeV with 13 fb −1 (yellow) and the region of the parameters where the model is not applicable (grey) i.e. m * > Λ. Such experimental exclusion regions from Run I are compared with the contour plots expected from Run II, considering the process studied in this work 1 . The solid (magenta), dashed (red) and dot-dashed (green) -without shadingare the projected contour maps for S = 3 (3-σ) in the parameter space (Λ, m * ) of the statistical significance for √ s = 13 TeV and for the follwing three values of the integrated luminosity L = 30, 300, 3000 fb −1 .
Therefore an experimental study of the eejj signature at LHC is sensitive to a heavy composite Majorana neutrino up to masses of ≈ few TeV. In the absence of a discovery, it will be possible to increase the excluded regions of the parameter space; more so at increasing integrated luminosities. In the authors' opinion a dedicated experimental analysis of this channel with the data of RunII of the LHC should be undertaken.

VI. DISCUSSION AND CONCLUSIONS
In this work we take up the well known composite model scenario [47] in which ordinary quarks and leptons may have an internal substructure. The essential features of this scenario are: (i) the existence of massive excitations of the SM fermions, the so called excited quarks (q * ) and leptons ( * ), which interact via (effective) gauge interactions of the magnetic type with the light SM quarks and leptons; (ii) the presence of four fermion contact interactions between ordinary fermions and also between ordinary and excited fermions.
These excited states have been searched for in a number of dedicated direct searches. In this study we reconsider the hypothesis that the excited neutrino is of Majorana type. This ansatz had been considered in [3] where a model based on gauge interactions only was used to describe the production and decay rates of the composite Majorana neutrino. We have included the contribution of contact interactions in the phenomenology of the excited Majorana neutrino, hitherto not considered in the literature. The model is implemented in CalcHEP which allows quite extensive simulations at the generator level. The contact interaction mechanism turns out to be dominant in the resonant production of the heavy Majorana neutrino. We have performed a fast simulation study of the same-signed di-lepton plus di-jet signature (eejj) arising from the resonant production of a heavy composite Majorana neutrino and its subsequent decay, analysing in detail both signal and background in order to optimise the statistical significance.
We have performed such phenomenological study of the production of heavy composite Majorana neutrinos at LHC also in view of possible connections with the recent observations by the CMS Collaboration of: (i) a 2.8-σ excess in the eejj channel in a search for W R gauge bosons [10]; and (ii) a 2.4-σ and a 2.6-σ excesses respectively in the eejj and ep T / jj in a search for leptoquarks [11,12].
We find that the invariant mass distribution of the system made up of the second-leading lepton and the two jets is highly correlated to the heavy Majorana neutrino mass, see Fig. 6 (bottom-right). A fast simulation of the detector effects and efficiencies in the reconstruction process is performed using the Delphes [39] package based on a CMS-like configuration.
We scanned the two dimensional parameter space (m * , Λ) for some benchmark values and computed the statistical significance. We provide the contour plots of the statistical significance S at 3-and 5-σ (see Fig. 8). We find for instance that with Λ = 15 TeV the LHC can reach a 3-σ sensitivity for masses up to m * = 1500, 2500, 3000 GeV, respectively for an integrated luminosity of L = 30, 300, 3000 fb −1 .
Finally, in the parameter space (m * , Λ), we compared (see Fig. 9) such 3-σ significance curves with the 95% C. L. exclusion regions from experimental data of Run I c.f. [4,5], (see also footnote in Sec. V). Such analyses have investigated signatures of excited electrons and muons ( = e * , µ * ) being produced by contact interactions (pp → * ) and decaying via * → +γ. Strictly speaking such analyses access the parameters spaces (m e * , Λ) and (m µ * , Λ) that are in principle different from the one presented here (m N , Λ). However, all excited states masses can be assumed to be approximately degenerate, at least at a first order approximation. Under the hypothesis M N ≈ m e * ≈ m µ * ≈ m * the eejj signature discussed in this work provides contour maps that can be considered on the same parameter space of the other analyses based on pp → γ [4,5]. This comparison shows that the eejj signature from a heavy composite Majorana neutrino has the potential to improve sensibly the current constraints on the composite scenario.
Before concluding we would like to comment briefly about another anomaly reported by the ATLAS Collabortion and on recent results from the ATLAS and the CMS Collaborations from RunII, and how our model could interpret them.
In a search [63] for high-mass di-boson resonances with boson-tagged jets at √ 8 TeV the ATLAS Collaboration has reported an excess at around 2 TeV with a global significance of 2.5 standard deviations (note however, that the same search performed by the CMS Collaboration did not observe a similar excess [64]). Our model contains fermion resonances (excited quarks and leptons) which do not couple directly to a pair of gauge bosons. Indeed a fermion cannot decay to a pair of gauge bosons by angular momentum conservation. On general grounds our fermion resonances could produce final states with a pair of gauge bosons but these would always be accompained by other objects such as leptons and jets (SM fermions). As an example one might think to pair produce the excited neutrinos pp → Z * → ν * ν * with the excited neutrinos decaying leptonically ν * → W ± e ∓ . One obtains a final signature of W + W − e + e − which is different from the one considered in the ATLAS search for high mass di-boson resonances consisting of only gauge boson pairs (W W, W Z or ZZ).
However, one might imagine to pair produce the charged excited fermions, for instance e * and/or q * , almost at threshold (if they are very massive). Such pair of heavy fermions could in principle form a 1S bound state (via the known Coulomb and/or color interaction) which in turn could decay to a pair of intermediate vector boson given the high mass of the hypothetical heavy fermions [65,66]. Therefore our model has in principle the potential to reproduce an excess in the di-boson signal. The estimate of such effects is certainly very interesting and it would surely be worth further investigation (one needs for instance to understand whether such bound states could form at all in the first place). However, a quantitative analysis goes beyond the purview of the present work and we postpone it to a future study.
Very recently the CMS and the ATLAS Collaborations have released the first results of Run II of the LHC at √ 13 TeV, with, respectively 42 pb −1 and 80 pb −1 of integrated luminosity [67,68], reporting about a search for hadronic resonances in the di-jet channel and showing ≈ 1-σ excess(es) at an invariant mass of about 5 TeV in the measured di-jet invariant mass distribution.
Such excess(es), if confirmed by further data and statistics, could in principle signal the first hadronic (excited quark) resonance level in a composite model scenario beyond the 2 TeV eejj anomaly. Indeed the analysis in [67] excludes excited quarks masses from around 3 TeV at 95% C.L (if m * = Λ) while in [5] excited lepton masses are excluded at 95% C.L. from m * = 2.5 TeV (again assuming m * = Λ). These experimental bounds would seem to preclude the possibility of an excited fermion bound state with regard to the explanation of the di-boson anomaly at 2 TeV (see above) since a q * (e * ) bound state would need to have a mass of at least 6 TeV (5 TeV). However, the quoted bounds for the excited fermions are for m * = Λ which is the limit of validity of the effective composite model. For values of Λ higher than m * the actual bounds on m * are lower (see for example Fig. 9). In such regions of the parameter space the 2 TeV di-boson anomaly could still be explained (in principle) by a q * q * (or e * ē * ) bound state with a mass m * ≈ 1 TeV.
As a last remark concerning the eejj anomaly we would like to comment on the fact that (i) the same excess is not observed in the µµjj channel and (ii) the observed charge asymmetry of the like sign di-leptons [10]. The absence of the excess in the µµjj channel could be explained by our model simply assuming that the excited muon state (µ * ) is somewhat heavier than the e * and so it would be observable only at higher energies. The observed eejj excess consists indeed of 14 events of which 13 are opposite sign (OS) and only one is same sign (SS). It must be said that our Mirror type composite model with one Majorana neutrino will produce the same yield of OS and SS events. Such feature could be explained within our composite model assuming the existence of an additional Majorana ν * state with a slightly different mass. Indeed it has been shown, albeit within a different (seesaw) model [16,17], that the interference between the contributions of two different Majorana states could depress the SS yield relative to the OS. The interference effect could also explain the absence of a peak in the observed invariant mass distribution of the second leading electron and the two jets (see our Fig. 6 -bottom right-). In view of this it could be worthwhile either to upgrade the CalcHEP implementation of our Mirror model to include other Majorana states or alternatively to reconsider the homo-doublet model with ν * L -ν * R mixing. In order to address quantitatively this issue we would need to build a new model (with more than one Majorana neutrino state) in the CalcHEP generator. This goes beyond the scope of the present work and will be addressed in a future study.
In summary the results presented in this work are quite encouraging and certainly endorse the interest and feasibility of a full fledged analysis of the experimental data of the upcoming LHC Run II for a search of the heavy composite Majorana neutrino, within a Mirror type model, in proton proton collisions at √ s = 13 TeV.

NOTE ADDED IN PROOF
While completing this work we became aware that the CMS collaboration has completed an experimental analysis [69,70] of a search for heavy composite Majorana neutrinos based on the model discussed here.
Using 2.6 fb −1 data of the 2015 Run II at √ s = 13 TeV, heavy composite neutrino masses are excluded, at 95% CL, up to m N = 4.35 TeV and 4.50 TeV for a value of Λ of 5 TeV, from the eeqq channel and the µµqq channel, respectively.