METing SUSY on the Z peak

Recently the ATLAS experiment announced a 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} excess at the Z-peak consisting of 29 pairs of leptons together with two or more jets, ETmiss\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{\mathrm {T}}^{\mathrm {miss}}$$\end{document} >225\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$> 225$$\end{document} GeV and HT>600\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{\mathrm {T}}\,> 600$$\end{document} GeV, to be compared with 10.6±3.2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10.6 \pm 3.2$$\end{document} expected lepton pairs in the Standard Model. No excess outside the Z-peak was observed. By trying to explain this signal with SUSY we find that only relatively light gluinos, mg~≲1.2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{\tilde{g}} \lesssim 1.2$$\end{document} TeV, together with a heavy neutralino NLSP of mχ~≳400\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{\tilde{\chi }} \gtrsim 400$$\end{document} GeV decaying predominantly to Z-boson plus a light gravitino, such that nearly every gluino produces at least one Z-boson in its decay chain, could reproduce the excess. We construct an explicit general gauge mediation model able to reproduce the observed signal overcoming all the experimental limits. Needless to say, more sophisticated models could also reproduce the signal, however, any model would have to exhibit the following features: light gluinos, or heavy particles with a strong production cross section, producing at least one Z-boson in its decay chain. The implications of our findings for the Run II at LHC with the scaling on the Z peak, as well as for the direct search of gluinos and other SUSY particles, are pointed out.


Introduction
The discovery by ATLAS [1][2][3] and CMS [4] experiments at the LHC Collider [5,6] of a new particle with a mass of 125 GeV [7,8] and with the expected properties of a Higgs boson has marked the program of high-energy physics for the next coming years. On one side, it is mandatory to be precise enough in the measurements of the properties of the new scalar particle in order to definitively ascertain its nature as the Higgs boson remnant of the electroweak symmetry breaking mechanism. On the other, its mass is still compatible with the requirements imposed by supersymmetry (SUSY) at the expense of moving the SUSY scale above TeV enera e-mail: oscar.vives@uv.es gies. Combined with the current LHC constraints (although model dependent in most cases) from data analyses in the first run of LHC, the search for SUSY effects is becoming more restrictive. Nevertheless, the discovery of SUSY would be such an extraordinary event, not only by itself, but for solving pending fundamental experimental and theoretical problems in particle physics, that an intense well-motivated experimental program to search for SUSY effects is of the highest interest.
Within this scenario, the ATLAS Collaboration has recently presented an intriguing excess at the 3 sigma level of e + e − and μ + μ − pairs just at the Z peak [9,10], accompanied by hadronic activity and missing transverse energy (MET). With an integrated luminosity of 20.3 fb −1 at a p − p CM energy of 8 TeV, the experiment observes a total of 29 pairs of electrons and muons with an invariant mass compatible with the Z -boson mass, with an expected background of 10.6 ± 3.2 pairs. No excess over the expected background is observed outside the Z peak. 1 The question that immediately arises is whether SUSY, or some other extension of the Standard Model (SM), is able to explain that excess of Z + MET events taking into account the current limits on beyond the SM physics. A study in those terms within a SUSY framework is presented in this paper.
As we will show, the observed signal can only be explained if one has a large production cross section of heavy SUSY particles (gluinos or squarks) whose decay chain contains about one Z -boson per parent particle. If such an explanation is indeed the answer to the observed excess, our study points out the way to confirm it in the Run II of LHC, as well as cosmological implications, in particular the particle content of dark matter in the universe. The resulting scheme of SUSY particle mass hierarchy, including charginos and neutralinos, will be apparent.

ATLAS excess in l + l − on the Z peak
In order to motivate and clarify the assumptions that will be made in the next sections, a simplified summary of the mentioned ATLAS search is given here.
The main focus of the search in [9] is the decays of squarks and gluinos with two leptons (electrons or muons), jets and E miss T in the final state, where the two leptons originate from a Z -boson.
In order to discriminate between SM background events and a possible signal, the following requirements are applied: • At least two same-flavored leptons with opposite electric charge are required in each event. If more than two leptons are present in the event, the two with the largest values of p T are selected. The leading lepton, i.e. the lepton with highest p T , must have p T > 25 GeV , whereas the subleading lepton p T can be as low as 10 GeV . Their invariant mass must fall within the Z -boson mass window, here considered as 81 < m ll < 101 GeV . • All events are further required to contain at least two jets with p T > 35 GeV and |η| < 2.5, to have E miss T > 225 GeV and H T > 600 GeV , where H T is the p T sum over all the jets with p T > 35 GeV and |η| < 2.5 and the two leading leptons: • Furthermore, the azimuthal angle between each of the two leading jets and E miss T is required to be higher than 0.4.
A great effort has been made to accurately estimate the number of SM events that survive the previous selection. An enumeration of these expected SM processes together with some of their characteristics follows: • The main background, namely tt, together with W W , W t, and Z → τ τ , which add up to ∼60 % of the predicted background, have been estimated using a datadriven method that has been thoroughly cross-checked with different techniques. • Diboson backgrounds with real Z -boson production (∼25 %) and "rare top" (tt + W, tt + Z , tt + W W and t +Z ) backgrounds (<5 %) are estimated using MC simulation. These are subject to carefully assessed theoretical and experimental uncertainties. • Processes with "fake leptons", i.e. leptons originating from the mis-identification of a jet, (∼10 %) are estimated using a data-driven method used regularly in most of ATLAS analyses. • Finally, particular care has been taken to suppress the Z /γ * + jets background as much as possible, given that it could mimic a possible signal (the cut in the azimuthal angle between each of the two leading jets and the E miss T has been applied to serve this purpose). Nonetheless, a data-driven technique has been used to estimate this small (<1 %) but important background.
The total number of expected SM model events passing all the requirements is 10.6 ± 3.2 and the number of observed data events is 29. This corresponds to a 3.0 σ significance.

Z-boson production in the MSSM
As explained above, if the observed excess is confirmed, it would clearly point to a new non-standard process producing additional Z -bosons at LHC energies. Z -bosons are regularly produced in the decay chains of most of the SM extensions. Still, this signal would require a significant production of Z -bosons without conflicting with all other experimental searches of beyond the SM particles. In fact, using the central value for the expected background and taking into account the Z -branching ratio to muon and electron, this would imply that we have produced 273 ± 48 additional Z -bosons (with 20.3 fb −1 ). Assuming the Z -bosons are produced in the decay chains of beyond the SM particles, Y , produced in the collision, we need to produce at least 273/N (Y → Z ) Y particles, with N (Y → Z ) the average number of Z -bosons produced in the decay of a Y particle. On the other hand, as we will see next, the experimental cuts used in the experiment, namely n jets ≥ 2, E miss T > 225 GeV and H T > 600 GeV , define further the characteristics of the Y particle and its decays. Now, the question is whether it is possible to produce these extra Z -bosons with the associated event attributes in some SUSY extension of the SM, while at the same time all the constraints imposed by present new-physics searches at LHC and other experiments are satisfied. In this section, we will assume that the acceptance of the applied selection, also taking into account the reconstruction efficiencies is ideally equal to unity. In Sect. 4, we will perform a "realistic" simulation in a SUSY model using Delphes [12] to take into account the signal acceptance and detector efficiencies.

Production cross section of MSSM particles
We need to produce 273 (225 at 1 σ ) Z -bosons if we want to accommodate the observed excess in lepton pairs. Assuming that R-parity is conserved, supersymmetric particles are produced in pairs in processes of the type pp → YȲ . Thus, the required cross section for this process would be, where we take into account that two Y particles are produced in each event. So, if we obtained one Z -boson for each Yparticle produced, we would need a production cross section of (6.7 ± 1.1) fb at the LHC with a CM energy of 8 TeV.
Then the first thing we must do is to determine whether it is still possible to have these production cross sections for some supersymmetric particle taking into account the constraints from direct searches at LHC. Here, we consider the production cross sections of different supersymmetric particles separately to identify the relevant processes. However, in a full simulation, as done in Sect. 4, different sparticles contribute to the final Z plus jets plus MET signal and the different contributions must be added. Naively, the first option to consider in a hadron collider would be strong production of squarks or gluinos (assuming they produce Z bosons in their decays). However, current experimental searches of jets plus missing energy at LHC force the masses of these colored particles to be high [13][14][15]. Nevertheless, as we will see below, in some cases we can still find cross sections of the required size.
Production cross sections of gluinos and squarks depend only on their masses and are basically independent of other MSSM parameters. In the case of gluino and squarks of the first generation, the cross section depends both on the squark and gluino masses due to t-channel contributions, but in the case of stops or sbottoms it depends only on the stop or sbottom mass. In Fig. 1 we present the production cross section of gluino pairs and light-flavor squark pairs 2 calculated at NLL+NLO with NLL-fast [16][17][18][19][20][21][22] as a function of the gluino or squark mass. Figure 2 shows the squark plus gluino cross section as a function of the squark or gluino mass with the second mass fixed at different values. The different bands in these figures correspond to the cross section at 1 σ with fixed squark or gluino masses: the blue (dashed) band corresponds to mq ,g = 1000 GeV, the brown (dashdotted) band to mq ,g = 1500 GeV, the orange (dotted) band has mq ,g = 2500 GeV and the red (solid) band corresponds to decoupled squarks or gluinos. In all these plots, the horizontal gray lines represent the required cross section at 1 σ needed to reproduce the referred ATLAS excess.
As we can see in these figures, the required cross section is reached only for light gluino and squark masses. In the case of gluino production, the needed cross section is obtained only for mg 1200 GeV and favors heavy squark masses. In fact, these gluino masses are in the boundary of the allowed region obtained from jets plus missing E T searches at LHC [13][14][15] and would contribute significantly only if everyg produces at least a Z -boson in its decay. For the production of squark pairs, present limits are mq 1400 GeV for heavy gluinos and mq 1650 GeV for degenerate squarks and gluinos [13][14][15]. Under these conditions, σ ( pp →qq) is always well below the required cross section, even for mq 1400 GeV. Here, we do not show the cross section σ ( pp →qq * ) as it is typically one order of magnitude smaller than σ ( pp →qq), and thus irrelevant.
Another important process is the simultaneous production of squark and gluino shown in Fig. 2. However, we see that we would need both squark and gluino to be light, which is not possible if we take into account the bounds from LHC searches [13][14][15]. In summary, the best option seems to be gluino pair production with mg 1200 GeV with relatively heavy squarks mq 3000 GeV, if we can get at least one Z -boson in every gluino decay.
Yet, we can also consider the strong production of stop pairs, where the current bounds on stop masses are much lower, mt 650 GeV [23,24]. The totaltt * production cross section is shown in Fig. 3, calculated at NLL+NLO with NLL-fast. In this case, we can see that we could reach the required cross section for mt 750 GeV which in principle could be achievable in general SUSY models (always assuming that every stop produces a Z -boson in its decay). However, the cuts H T > 600 GeV and E miss show the H T distribution from the decays of stop pairs with mt 1 = 750 GeV. As can be seen here, the H T distribution peaks at H T 400-500 GeV as a consequence of the relatively small stop mass, and only a small fraction of the events are able to overcome the cut on H T . Therefore, we must conclude that it is not possible to generate the required cross section and fulfill the requirements of the observed excess through stop production.
Apart from the production cross sections of colored sparticles, we could still consider the weak production of charginos/neutralinos which can be large enough for light gauginos. Taking into account that the current bounds on chargino and neutralino masses are not very stringent, electroweak production is worth exploring [25,26]. It is well known that the largest electroweak production cross sections are those corresponding toW 0W ± andW +W − , in terms of gauge eigenstates. This would correspond to χ 0 2 χ ± 1 and χ + 1 χ − 1 in terms of mass eigenstates. However, the observed signal requires the Z -boson to be on-shell, at least two jets, E miss T > 225 GeV, and a minimum H T of 600 GeV. As an example, we present, in Fig. 5, the production cross section ) in the case where the only restriction imposed is a minimum p T ≥ 20 GeV for the jets (in blue/dashed), and the same cross section after applying a cut on the hadronic H jets T ≥ 300 GeV. As can be seen, the cross section with charginos of m χ + 1 350 GeV looks, in principle, able to accommodate the required Z production. However, the situation changes after we impose the experimental cuts used in the analysis. The required E miss T is easily obtained if m χ 0 1 150 GeV, but the requirement on H T is very restrictive. In fact, for chargino masses m χ ± 1 350 GeV, the largest contribution to the production cross section would correspond to √ŝ 700 GeV, which can never produce H T of 600 GeV and E miss T > 225 GeV. Then we can expect the cross section σ ( pp → χ 0 1 χ + ) at larger √ŝ values to be strongly reduced. This is shown by the red (solid) line in Fig. 5: a relatively mild cut on the hadronic H T reduces the cross section by more than one order of magnitude.
Thus we must conclude that, although electroweak production could contribute efficiently to the production of additional Z -bosons, these events cannot overcome the experimental cuts and cannot give rise to the observed signal.
With this channel, we have reviewed all relevant production cross sections of different supersymmetric particles that could potentially explain the signal. The next step would be to calculate the average number of Z -bosons per parent particle Y , which we use in Eq. (1).

Decay of supersymmetric particles to Z -bosons
Z -bosons are produced through the decay chains of most MSSM particles, although the number of Z -bosons obtained per each supersymmetric particle produced depends on the identity of the supersymmetric particle initially produced and on the supersymmetric spectrum below its mass. The main sources of Z -bosons are the decays of neutralinos and charginos and also in some squark decays. The couplings of charginos/neutralinos to Z are given by Therefore, we can obtain Z -bosons in the decays of higgsinolike neutralinos and charginos. For instance, we can define a typical mSugra spectrum with gaugino mass unification at the grand unification scale and all the masses fixed by four parameters (m 0 , m 1/2 , A 0 , tan β) plus the sign of the μ parameter. In this spectrum, the gaugino masses at the electroweak scale are approximately given by M 1 0.4m 1/2 , M 2 0.8m 1/2 , M 3 2.6m 1/2 , and the sfermion masses are m 2 q m 2 u m 2d m 2 0 + 6m 2 1/2 , m 2 l m 2 0 + 0.5m 2 1/2 and m 2 e m 2 0 + 0.15m 2 1/2 with relatively small dependence on A 0 and tan β. The value of |μ| 2 has a stronger dependence on A 0 and, specially, on tan β, for tan β ∼ 5, we would have Then, in such a spectrum, the second neutralino will only produce Z -bosons through its (relatively small) higgsino component while the two heavier neutralinos can be expected to produce a sizable number of Z -bosons. On the other hand, charginos can produce Zbosons both through the wino and from the higgsino component but only in decays of the heavier charginos, as the lightest one will only decay to a W -boson and a neutralino (or lepton-slepton if ml ≤ m χ + ).
In Fig. 6 we can see the values of N (χ 0 2,3 → Z ) as a function of m χ 0 2,3 in an mSugra model. Here, we obtain N (χ 2 → Z ) around 0.1, as expected if the higgsino content is relatively small. Then N (χ 3 → Z ) can reach at most 0.45 while the other 50 % of the decays go to χ ± W ∓ . Similarly the heavy charginos produce Z -bosons in their decay as can be seen in Fig. 7. In this case N (χ + 2 → Z ) can be 0.3, while we have similar branching ratios to χ 0 1 W + and χ + 1 h. Besides chargino and neutralino decays, Z -bosons couple also to sfermions through, As we can see, the Z couplings are chirality diagonal and therefore, in decays, they can only enter through chirality mixing. Although these couplings could be flavor changing, flavor mixing is bounded to be small due to the stringent flavor changing neutral current constraints. Therefore we can expect a sizable amount of Z -bosons produced only through chirality mixing in third generation sfermions in decays likẽ t 2 →t 1 + Z , orb 2 →b 1 + Z in the large tan β regime.
In addition, we obtain Z -bosons in the decay chains of strongly produced sparticles. We can obtain Z -bosons at different steps of the decay chain, either through the couplings Therefore taking into account the corresponding branching ratios, this decay chain would contribute with 2 ×BR(g → χ 0 , the number of Z -bosons produced perg produced.
As we can see in Fig. 8, in a typical MSSM spectrum we obtain at most 0.2 Z -bosons per stop or gluino while other squarks produce far fewer Z -bosons per squark. Although Figs. 6, 7, 8 have been obtained from a mSugra spectrum, the expected number of Z -bosons would be very similar in other MSSM versions, as it depends only on the spectrum below the mass of the originally produced Y particle and the content of charginos-neutralinos.
As shown in the previous section, the only possibility to explain the signal is to produce a stop or a gluino pair as the lightest colored sparticle, being all other squarks much heavier and only neutralinos, charginos, and possibly some sleptons can be below the gluino or stop mass. Moreover, given the size of production cross sections consistent with the present searches, we need to obtain nearly one Z -boson per Y particle produced. Therefore, from the expected numbers of Z -bosons that we have seen in this section, we have to conclude that it is not possible to reproduce the observed signal in a MSSM with a stable (and light) neutralino.
Nevertheless, we can still consider different variations of the MSSM: • A first possibility would be to have a light gluino below 1 TeV that can evade the bounds from jets plus missing E T if it decays to a sufficiently heavy neutralino LSP in a sort of compressed spectrum. Under these conditions the gluino would be abundantly produced at LHC and even a small number of Z -bosons per gluino could fulfil the requirements to explain the observed signal. However, this would require strongly non-universal gaugino masses and very heavy LSP's and we will not follow this path here. • A second option is to consider an MSSM where the lightest neutralino decays to a lighter gravitino plus some Zboson. This is the case in gauge mediated MSSM [27][28][29][30] and it could also be possible in gravity mediated MSSM if the gravitino is lighter than the neutralino which then becomes the NLSP [31]. In this case, the neutralino decays to Z -boson and gravitino if it is allowed by phase space and the branching ratio will depend on the lightest neutralino composition. This is the possibility we will explore in the following.
Thus, we will analyze a situation where the neutralino is the next to lightest supersymmetric particle and the LSP is the gravitino. All supersymmetric particles will decay to the lightest neutralino which then decays to gravitino plus a photon, a Z -boson or a Higgs. The decay width of the lightest neutralino to photon, h or Z plus gravitino [31][32][33] is given by with F(x, y, z) a function of the particle masses, irrelevant for our discussion, that can be obtained from [31]. As we can see, if the lightest neutralino is bino-like, N 11 1, and the mass difference between neutralino and gravitino is larger that the Z -mass, the branching ratios are BR(χ 1 →Gγ ) cos 2 θ W 0.8 and BR(χ 1 →GZ) sin 2 θ W 0.2. Similarly if the lightest neutralino is wino-like, the branching ratios get exchanged. From this equation we can also see that it is possible to get a very large BR to Z -bosons as needed to reproduce the observed signal if (−N 11 sin θ W + N 12 cos θ W ) 1, but this is only possible if the lightest neutralino has a very large wino component. Although in gauge mediation models the LSP is always the gravitino, the gaugino masses in minimal models are proportional to the gauge couplings and therefore the LSP is mostly bino with small wino and higgsino components. Then we will have to consider other extensions of the gauge mediation idea, like the so-called general gauge mediation (GGM) where the gaugino and sfermion masses depend on hidden sector current correlators which can be different for different gauge groups or particle representations. As we will show in the following section, in these models it is possible to obtain general neutralino NLSP as required to explain the excess.

A possible explanation in general gauge mediation
Minimal gauge mediation predicts that all scalar and gaugino masses originate from a single scale and powers of the gauge couplings [30]. Recently a model-independent generalization of gauge mediation was proposed under the name of general gauge mediation [34,35], where all the dependence of soft masses on the hidden sector is encoded in three real and three complex parameters obtained from a small set of currentcurrent correlators. In these models the gaugino and sfermion masses are given by with ρ (y) (with ρ = 0, 1/2, 1, corresponding to scalar, fermion and vector) are associated with the currentcurrent correlators in the hidden sector, ζ is a possible Fayet-Illiopoulos term (ζ = 0 in the following), C 2 ( f |r ) the quadratic Casimirs and M s a characteristic SUSY-breaking scale in the hidden sector.
Having six parameters,B 1/2 r (0) andÃ r , to fix the soft masses in the observable sector, it is clear now that we have much more freedom in GGM [36][37][38][39] and, in particular, we have as required to reproduce the observed signal at ATLAS. In particular, we need the NLSP to decay to gravitino and a Zboson with a branching ratio close to one. Fortunately, this is possible in GGM as shown in Refs. [39,40].
In this GGM scenario we have used SPheno-3.3.3 [41,42] to obtain the full supersymmetric spectrum at LHC energies. We define a first parameter point GGM1 with the following parameters: M s = 400 TeV,B  Table 1. With respect to this spectrum, some comments are in order: 1. The two lightest neutralinos and the lightest chargino are very similar in mass, ∼430 GeV and this allows a large neutrino mixing as required. In fact the neutralino mixing matrix is given by On the other hand the relatively large NLSP mass is needed to overcome the E miss T cut. 2. The gluino is relatively light mg = 1088.0 GeV which allows for a sizable production cross section and taking into account the squark masses of order ∼3 TeV. As we can see in Ref. [15], this mass would be allowed at 1 σ in a simplified MSSM with gluino-squark-neutralino. 3. The lightest Higgs mass must reproduce the observed value at LHC of m h 125 GeV and in this spectrum it reaches only 119.4 GeV. This problem (typical in minimal gauge mediation models) can be solved either by increasing the stop masses taking a largerÃ 3 or assuming extra operators in the Higgs sector, as the dimension-5 operators proposed by Dine et al. [43]. Here, we assume that this problem is solved by one of these mechanisms, given it does not affect the observed phenomenology on the Z -peak.
We simulate the production of supersymmetric particles at LHC at 8 TeV (LHC8) with this spectrum using Pythia 8.1 [44] with Prospino2 [16,17,45,46] K-factors and the response of the ATLAS detector using Delphes [12]. The selection of events for this study is done as close as possible to that performed in ATLAS [9]. The dashed red line in Figs. 9 and 10 shows the H T and E miss T distributions respectively for the GGM1 point after applying all the selection detailed in Sect. 2 except for the H T and E miss T cuts. In the H T distribution, we can distinguish the two peaks corresponding to electroweak production at low H T values and gluino production at higher H T . From here, we can expect that the cut on H T will eliminate most of the electroweak production but not the gluino production. This can be seen in Fig. 10, where the E miss T distribution is presented separately for strongly produced events (solid blue line) and for the electroweak component of the production (dashed red line) for the same selection as in Fig. 9 but after applying the H T cut, i.e. the final selection except for the cut on E miss T . The electroweak component is significantly reduced by the H T cut while mainly only events coming from strong production survive the cut, as expected. In fact, in this simulation of point GGM1, 99 % of the gluino points and only 11 % of the electroweak points have survived the H T cut. We see that the peak in the E miss T distribution is approximately at m χ 0 1 /2, and, for m χ 0 1 = 425 GeV, a reasonable fraction of the events will survive the E miss T cut at 225 GeV. In the simulation, 65 % of the gluino point and 53 % of the electroweak points survive this cut. However, due the relatively small production cross section, the final number of events is small. In this simulation and after applying all relevant experimental cuts, an T distribution corresponding to the GGM1 point from strong production (solid blue line) and electroweak production (dashed red line) with the same selection as in Fig. 9, after applying the H T > 600 GeV cut expected signal of 6.34 ± 1.02 lepton pairs is found, to be compared with the observed excess of 19.4 ± 3.2. This number of surviving events was obtained at NLO with Pythia and Prospino2 but, unfortunately, it is still too low to explain the observations.
Trying to obtain a model able to account for the excess, we consider a second point in our GGM scenario with a lighter gluino. The GGM2 point is obtained with M s = 400 TeV,  Table 2.
The neutralino mixing matrix in point GGM2 is similar to the corresponding mixing matrix in GGM1, and the BR(χ 0 1 →G Z) = 0.94. However, gluino is now much lighter and the gluino-gluino cross section is now (41.6 ± 7.5) with Prospino at NLL+NLO, thus we can expect many more T distribution corresponding to the GGM2 point from strong production (solid blue line) and electroweak production (dashed red line) with the same selection as in Fig. 9, after applying the H T > 600 GeV cut gluino pairs to be produced and a larger contribution in the final selection for this point.
The simulation for this GGM2 point is presented by the solid blue line in Fig. 9. The H T distribution peaks at slightly lower values than in the case of GGM1, due to the slightly lower gluino mass, but it is still enough to overcome the H T cut at 600 GeV. The E miss T distributions are similar for both GGM points, due to the very similar neutralino masses in both cases. However, in the case of GGM2 the strong production cross section is larger and much more important in relation with the electroweak production: for the GGM2 point, the strong production represents ∼55 % of the total cross section. As seen in Fig. 11, after applying the selection we obtain an expected signal of 28.0 ± 4.7 events, slightly over the excess reported by ATLAS, showing that a signal point defined along these characteristics can be able to reproduce the excess. We have to emphasize here that it is not difficult to obtain the observed excess for light gluino masses, and a gluino mass in between the two presented examples, mg ∈ [900, 1100] GeV, could reproduce the observed signal. However, these points may be in conflict with direct searches of jets plus E miss T [13][14][15]. There is a tension between this excess and the bounds from gluino searches in jets plus E miss T . To quantify more accurately this tension, we use the program Checkmate [47] that allows one to compare the different points of the model with various experimental analyses [47,48] determining whether the point is excluded or not at 95 % C.L.. We applied the constraints on jets plus E miss T from ATLAS studies: ATLAS.1308.1841 [49,50], ATLAS.1308.2631 [51], ATLAS.1407.0583 [24] and ATLAS.1405.7875 [15]. In fact, we find that this point is indeed excluded by the analysis ATLAS.1405.7875, in the signal region with six (or more) jets, with an r ≡ S−1.96 S S 95 obs = 2.0, where S is the total number of expected signal events, S is the total 1 σ uncertainty on this number and S 95 obs is the experimentally measured 95 % confidence limit on signal events.
Thus, we have to find another parameter point able to provide the signal, but satisfying all the present constraints on jets plus E miss T . The main difficulty in the GGM scenario is that each gluino decay produces typically four or more jets with large p T and E miss T . Therefore, these points clash with constraints from observables with a large number of jets (plus E miss T ), which have low backgrounds in the Standard Model. On the contrary, constraints with fewer jets are easier to satisfy because of larger backgrounds. Then our strategy will be to compress sufficiently the spectrum to reduce the number of (observable) jets. This can be done by increasing the lightest neutralino mass closer to the gluino mass, so that the jets in the decaysg → χ 1 + jets have smaller p T . With this goal, we construct a third GGM point. The  Table 3.
For this point, the BR(χ 0 1 →G Z) = 0.98 and the gluino production cross section is (22.8 ± 3.3) fb with Prospino at  Fig. 9, before (red dashed) and after (solid blue) applying the H T > 600 GeV cut. Notice that electroweak production is negligible for the neutralino masses in GGM3 NLL+NLO. In this case, the signal is lower than in GGM2 but we still get (13.1 ± 2.2) events, as shown in Fig. 12, well above the experimentally measured 95 % confidence limit on signal events.
Concerning the experimental searches on jets plus MET, according to Checkmate this point is allowed with an r = 1.0 with the best signal region being two jets plus E miss T from ATLAS.1405.7875. Moreover, in the equivalent search by CMS of two leptons, jets, and E miss T , CMS-SUS-14-014 [11], we obtain r = 1.2 in Checkmate. This result shows some tension with the ATLAS positive signal from the bin with large E miss T in the CMS analysis, but still consistent with observations slightly above 95 % C.L. Therefore, we have proved that it is indeed possible to construct a supersymmetric model that accommodates the observed excess of lepton pairs on the Z -peak. The simulations presented here are only a proof of existence and the final model may be very different. Nevertheless, this model will have to share the main features of the examples that we presented here.

Prospects for SUSY at LHC13
As we have shown in the previous sections, the excess observed in ATLAS, if due to SUSY, requires a gluino of a mass ∼1 TeV producing nearly one Z -boson per gluino in its decay. This scenario would also require relatively heavy squarks of the first generation with mq 2.5 TeV. If this is indeed the correct explanation to the observed excess, such light gluinos would be abundantly produced at Run II in LHC together with other SUSY particles. Therefore, the results obtained at pp collisions at 13 TeV will immediately confirm or reject this supersymmetric explanation of the ATLAS excess.
Using NLL-Fast with the spectrum of point GGM1, the gluino pair production cross section at 8 TeV was σ ( pp → gg) NLL+NLO LHC8 = 7.6 ± 1.3 fb at NLL+NLO. Similarly the production cross section at 13 TeV would be σ ( pp → gg) NLL+NLO LHC13 = 150 ± 16 fb, that is, we would expect to produce 20 times more gluinos at LHC13 for the same integrated luminosity. Repeating the same exercise with point GGM3, we have a gluino pair production cross section at 8 TeV of σ ( pp →gg) NLL+NLO LHC8 = 22.8 ± 3.3 fb while the production cross section at 13 TeV would be σ ( pp → gg) NLL+NLO LHC13 = 344 ± 44 fb. In this point, the cross section increases by a factor ∼15 in going from 8 to 13 TeV center of mass energy. In any case, for both points this would have unambiguous signatures, both on the Z -peak with a scaling of the signal found at LHC8 and in direct searches for gluinos using jets plus missing E T in the extension of the analysis in [13][14][15].
To confirm that it is indeed SUSY behind the excess found at the Z -peak, we should search for other sparticles at LHC13. As shown in Sect. 3.1, squarks of the first generation are preferred be heavy to increase the gluino production cross section, namely above 2.5 TeV. Then their production cross section at LHC13 would be around 1 fb which would make direct detection very challenging. Charginos and neutralinos are expected to be below the gluino mass and in some cases could be abundantly produced at LHC. Although our analysis does not fix the masses of χ ± 1 and χ 0 1,2 , they are expected to be rather heavy m χ 0 1 300 GeV and degenerate. In some cases, the production cross section could be large. For instance σ ( pp →χ + 1 χ 0 1 ) NLO LHC8 = 60.2 ± 0.9 fb and σ ( pp →χ + 1 χ 0 1 ) NLO LHC13 = 154 ± 5 fb calculated with Prospino for m χ 0 1 = 313 GeV and m χ + 1 = 314 GeV (with largeW component in χ 0 1 ). Notice that, as explained above, electroweak production is also large for these points at 8 TeV, but it is eliminated by the cuts on H T and number of jets. Thus, a large electroweak production could be expected at LHC13 and dedicated searches should be encouraged, specially taking into account the requirement of a large production of Z -bosons in chargino and neutralino decays. Finally, the signal does not constrain the masses of sleptons or third generation squarks and, a priory, we cannot make any prediction on their production at LHC13.
Before closing, we should comment on the nature of dark matter in our scenario. As we have seen the signal seems to prefer a non-stable neutralino decaying to a very light gravitino and a Z -boson. Under these conditions the neutralino mass has no relation with the dark matter abundance of the universe and its role as dark matter component would be played by the gravitino. The gravitino mass is not bounded by the observed signal but regardless of its exact mass, unfortunately no signal of dark matter is to be expected in direct search experiments.

Conclusions
Recently the ATLAS experiment announced a 3 σ excess at the Z -peak consistent of 29 pairs of leptons observed to be compared with 10.6 ± 3.2 expected lepton pairs. No excess outside the Z -peak was observed. By trying to explain this signal with SUSY we found that only relatively light gluinos, mg 1.2 TeV, together with a heavy neutralino NLSP of mχ 400 GeV decaying predominantly to Z -boson plus a light gravitino, such that nearly every gluino produces at least one Z -boson in its decay chain, could do it. Unfortunately, this is not possible withing minimal SUSY models, as mSugra, minimal gauge mediation or anomaly mediation. The requirement of a neutralino NLSP decaying to Z plus gravitino points to models of general gauge mediation as the simplest possibility. We have shown that a model of this class is able to reproduce the observed signal overcoming all the experimental cuts. Needless to say, more sophisticated models could also reproduce the signal, however, they will always share the above mentioned features, i.e. light gluinos (or heavy particles with a strong production cross section) with an effective N (g → Z ) 1.