MSSM fits to the ATLAS 1 lepton excess

We use the framework of the p19MSSM to perform a fit to the mild excesses over the Standard Model background recently observed in three bins of the ATLAS 1-lepton + (b-)jets + MET search. We find a few types of spectra that can fit the emerging signal and at the same time are not excluded by other LHC searches. They can be grouped roughly in two categories. The first class is characterized by the presence of one stop or stop and sbottoms with mass in the ballpark of 700-800 GeV and a neutralino LSP of mass around 400 GeV, with or without the additional presence of an intermediate chargino. In the second type of scenarios the stop, lightest chargino, sbottom if present, and the neutralino are about or heavier than ~650 GeV and the signal originates from cascade decays of squarks of the 1st and 2nd generation, which should have a mass of 1.1-1.2 TeV. For the best-fit scenarios, we compare the global chi-squared with respect to several ATLAS and CMS searches with the corresponding chi-squared of the Standard Model expectation, showing that the putative signal is also favored globally with respect to the background only hypothesis. We point out that if the observed excess persists in the next round of data, it should be accompanied by associated significant excesses in all-hadronic final state searches.


Introduction
The LHC is now in the course of its second run, characterized by proton-proton collisions at the center-of-mass energy of 13 TeV. By August 2016, the ATLAS and CMS Collaborations had collected approximately 14 fb −1 of data and presented the results of their searches for new physics beyond the Standard Model (BSM) at the ICHEP conference [1] in Chicago. To some disappointment, in almost all of the hundreds of presented channels the number of observed events turned out to be in excellent agreement with the Standard Model (SM) expectation, so that no new physics seems to be present at the tested energy.
On the other hand, as is bound to happen given the enormous number of kinematical bins and final-state channels analyzed by the two collaborations to cover the parameter space of as many BSM models as possible, a few anomalies have emerged to small significance. In particular, the 1-lepton ATLAS search [2] shows a 3.3 σ excess of events in the so-called DM-low bin, and less significant 2.6 σ and 2.2 σ excesses in the bins called bC2x-diag and SR1, respectively. In this letter we dedicate some attention to these excesses, as they are appearing in one of the classic topologies designed to discover supersymmetry (SUSY), which arguably remains the most comprehensive and appealing scenario for BSM physics.
It is obviously extremely premature to get excited about an excess like this, as its statistical significance is relatively low and experience has shown time and again that background fluctuations of comparable strength do happen with some frequency, especially when many channels are analyzed. On the other hand, given the appeal that SUSY has held on the particle community for many years, we also think it is legitimate to wonder whether the ATLAS 1-lepton excess already points toward specific SUSY spectra, and in particular spectra more involved than the simplified models used by the experimental collaborations to interpret their results.
We address the question by performing a global fit to the bins showing the excess in the generic framework of the phenomenological Minimal Supersymmetric Standard Model (p19MSSM) [3], which is characterized by 19 independent parameters. We use a set of model points generated by the ATLAS Collaboration in Ref. [4], which satisfies a set of experimental constraints from the relic density, direct dark matter (DM) searches, Higgs sector measurements, electroweak and flavor physics, as well as includes the exclusion bounds from the 8 TeV LHC SUSY searches. We also use the state-of-the-art exclusion bounds from 11 other ATLAS SUSY searches based on 3.2 fb −1 and ∼14 fb −1 data, which were derived by one of us in Ref. [5]. We repeat that, while our analysis satisfies a genuine curiosity, the reader will have to keep in mind that it might prove a futile effort were the excess to disappear or lose significance in the next round of data.
It is important to note in this regard that the CMS search [6] equivalent to ATLAS 1-lepton does not show any significant excess over the SM background with ∼13 fb −1 of data. By recasting the CMS analysis we will show, however, that this is not in contradiction with the onset of a possible signal observed at ATLAS, as the kinematical variables used by CMS and the subsequent definitions of the signal bins differ from the ones of ATLAS and the compatibility between the two is not straightforward.
Recently, an analysis of the ATLAS 1-lepton excess was performed in Ref. [7]. The excess was fitted with simplified SUSY model spectra (SMS) characterized by moderately light stops decaying into a bino or, in alternative, higgsino lightest SUSY particle (LSP). The originating signal was confronted with two CMS hadronic searches [8,9]. It was concluded in [7] that the signal is consistent at 2σ with a spectrum characterized by a ∼ 750 − 800 GeV stop, a bino LSP and a higgsino next-to-LSP. We extend the analysis here by considering a broader range of possibilities for the SUSY spectrum, as embodied by our choice of using the p19MSSM points to fit the excess. We show that there exist additional spectra and decay chains that can equally well fit the 1-lepton excess, and for each found scenario we compare its global χ 2 with respect to a sample of LHC searches with the χ 2 of the background only hypothesis, showing that all our scenarios are favored globally over the SM to some significance. We additionally provide possible signals for the next batch of LHC data and give a few comments on possible UV completions.
This paper is organized as follows. In Sec. 2 we outline the fitting procedure and we enumerate the constraints we impose on our model spectra. In Sec. 3 we present and discuss the main results of the fit and favored SUSY spectra. In Sec. 4 we evaluate the significance of the best-fit spectra globally with respect to the background only hypothesis. We summarize our findings and conclude in Sec. 5. Additionally, we present two appendices dedicated, respectively, to the validation of the CMS 1-lepton search, and to the detailed breakdown of the chi-squared contributions of every bin in the most important searches considered in this paper.

Fitting procedure and constraints
We investigate in this paper whether the excess in events shown by the ATLAS 1 lepton + (b-)jets search [2] can be interpreted in the generic framework of the p19MSSM.
To fit the three excesses we use a set of model points provided by the ATLAS Collaboration in Ref. [4] and generated using methods similar to those presented in [3,10,11,12]. The set consists of points for which mass and trilinear parameters have been scanned up to a maximum value of 4 TeV, with the exception of the third generation trilinear soft coupling A t , scanned in the range [−8 TeV, 8 TeV] to allow for the correct Higgs mass even when the stop mass is not very large. Phenomenological constraints from the relic density (upper bound), DM searches, Higgs and flavor physics, and electroweak precision data were taken into account, as was the impact of 22 ATLAS LHC searches at √ s = 7 and 8 TeV with integrated luminosity of 20.3 fb −1 . The SUSY spectra were calculated with SOFTSUSY v.3.4.0 [13], and the decay branching ratios with SUSY-HIT v.1.3 [14]. MicrOMEGAs v.3.5.5 [15,16] and SuperIso [17] were used to evaluate dark matter, precision electroweak and flavor observables, while the lightest Higgs boson mass was calculated with FeynHiggs v.2.10.0 [18,19]. The resulting p19MSSM set contained 183,030 allowed model points. More details on the sampling procedure and implementation of the experimental bounds can be found in Ref. [4].
We proceed as follows. We start with the sample of points that is not at present excluded by the current ATLAS 13 TeV searches for SUSY particles, which was constructed in Ref. [5] out of the 183,030 initial points. 1 All points were there tested on-the-fly by recasting for the p19MSSM 11 direct ATLAS searches:  [32].
Detailed description of the recast procedure, as well as validations of the implemented searches, can be found in Ref. [5].
Additionally, there are a few CMS 13 TeV SUSY searches that should be taken into account as they provide results complementary to those of ATLAS, in particular the 1-lepton [6] and 0-lepton [33,34,35] searches. Reference [6] is similar to the ATLAS search presented in [2] but, in contrast to the ATLAS results, it does not show any significant excess in the kinematical bins. As was mentioned in Sec. 1, this does not automatically reflect an inconsistency between the two results, as the definitions of the bins in both searches are quite different, but one needs to make sure that the models that fit the ATLAS excess are not instead excluded by CMS. To evaluate the impact of the CMS 1-lepton search [6] on our model sample, we implemented it in our recast tool as well. The validation of the recast procedure against the experimental results is given in Appendix A.
The 95% C.L. bounds on the stop/neutralino, sbottom/neutralino SMS provided in the CMS 0-lepton analyses are in general more constraining than those of the equivalent ATLAS searches [23,28] in the part of the parameter space where the ATLAS 1-lepton excess can be fitted. Bounds obtained in 0-lepton searches tend to be fairly stable under variations of the decay channels, so that for most of our points we will take the CMS exclusion lines at face value. On the other hand, we will draw attention to a few spectra that could fit the ATLAS excess if they were not excluded by the CMS 0-lepton searches when the SMS bound is taken at face value. Those spectra are characterized by a significant reduction of the branching ratio to all-hadronic final states, so that there is some doubt on whether they are actually excluded or not, and the answer would require a detailed recasting of the CMS searches as was done in Ref. [5] [2]. The name of the bin, the number of observed and background events, as well as the significance of the excess are provided by the experimental collaboration. The number of signal events that fit the excess best is derived according to the procedure described in the text.
We refrain from doing this here, given the still premature nature of the ATLAS excess, and limit ourselves to pointing out in the text the presence of these points "of doubt." We then proceed to perform a maximum likelihood estimate for the ATLAS 1 lepton + (b-)jets + E miss T search. A moderate excess shows up in 3 of the experimental bins [2], see Table 1. The bins are not statistically independent, and the same signal could be responsible for all 3 observed excesses. We limit ourselves to this hypothesis in this paper.
Indicating the observed number of events, the expected SM background, and the collaboration's estimate of the systematic uncertainty for each bin i as o i , b i , and δb i , respectively, we seek the maximum of the likelihood function where s i is the simulated signal for each bin, as a function of the points in the parameter space. The signals s i are simulated in the usual way, following the procedure given, e.g., in [5]. We use PYTHIA 8 [36] for the hard-scattering event and showering and DELPHES 3 [37] for the detector properties' simulation. We fit to all bins by making use of the ∆χ 2 variable: ∆χ 2 i = −2 log(L i /L 0 ), where L 0 corresponds to the backgroundonly hypothesis for the bins that do not show an excess, and to the signal given in Table 1 otherwise. Since the kinematical bins considered in the ATLAS 1-lepton search are not independent we build a global statistics as The sample of viable points fitting the excesses is narrowed down in a two-step procedure. First we identify the points that fit the excesses approximately at the 95% C.L. In other words we require ∆χ 2 1-lep ≤ 3.84. Given that the most significant excess in Table 1 is at most 3.3 σ, this first step allows a large number of points to survive. To further scale down the sample, we assume that the largest excess will yield in the near future a 5 σ discovery. A simple rescaling of the integrated luminosity gives us an approximate estimate of the luminosity L at which this might occur. Equating L /L 0 = 5 σ/3.3 σ yields, given the present luminosity of L 0 = 13.  search [2]. The mass of the light degrees of freedom in GeV, the branching ratio of the dominant decay yielding a hard lepton, and the goodness of fit according to Eq. (2) are shown.
L ≈ 32 fb −1 . Thus, as a final step, we rescale o i and b i in the 3 bins presenting an excess by the target luminosity L and subsequently fit to the maximum likelihood of Eq. (1) at the 95% C.L. We present in the next section the best fit spectra.

Spectra that fit the signal
Out of 183,030 p19MSSM model points allowed by the 8 TeV LHC data, 138,141 are not excluded after inclusion of the 13 TeV results. 1,293 model points fit the excess observed by ATLAS 1 lepton + (b-)jets + E miss T 13.2 fb −1 search at the 95% C.L., while 194 would produce a 5σ signal at the luminosity of ∼ 32 fb −1 . Out of those, 48 points are not excluded when the CMS 1-and 0-lepton searches are taken into account. These will be considered in what follows as the allowed models.
The resulting spectra can be divided into a few categories, depending on the properties of the light particles present in the spectrum and the mechanism that allows to fit the observed excesses. We describe them below, and summarize their features in Table 2 and in Fig. 1.
• Right-chiral stop scenario. The signal is consistent with the SMS employed by the experimental collaboration or similar spectra. The signal can be generated by the pairproduced stops decaying into a top quark and the LSP with BR(t 1 → tχ 0 1 ) ≈ 100%, as exemplified in Table 2 by the benchmark BP1. The table shows also the branching ratio of the mother sparticle into the channel that produces a hard lepton, and the χ 2 value of this point according to Eq. (2) (at the luminosity of 13.2 fb −1 ). The stop in  Figure 1: Benchmark spectra that can fit the ATLAS 1 lepton + (b-)jets + E miss T excesses [2]. The colors are associated with each sparticle as described in the legend. We indicate the SM particle emitted at every step of the decay chain next to the corresponding arrow.
this scenario is moderately light and mostly right-chiral, with mt 1 ≈ 600 − 800 GeV, and the source of the missing energy is a bino-like neutralino LSP with mass in the range ∼ 360 − 410 GeV. Note that spectra similar to this scenario were found to fit the signal in [7].
The model points of the ATLAS set [4] all yield an acceptable value for the relic abundance (at least within an upper bound), which means that when the LSP is binolike the spectrum presents also additional light particles in mass close to m χ 0 1 . In the case of BP1, these are a light right-chiral selectron and smuon with just about the same mass as the neutralino, which contribute to the early Universe co-annihilation with the LSP. Because their mass is almost degenerate with m χ 0 1 , direct SUSY searches are not yet sensitive to the presence of these particles, and since they do not affect the fitted 1-lepton signal we do not indicate their presence in Table 2 and in Fig. 1.
Alternatively, there might be a light chargino (either wino-or higgsino-like) in the spectrum, which participates in co-annihilating with the neutralino in the early Universe, as exemplified in the spectrum of BP2 in Table 2. Although a light chargino is necessary for the relic density, it does not directly participate to producing a signal in the 1-lepton ATLAS search. It has the effect of mostly reducing the branching ratio into the tχ 0 1 final state (because of a non-negligible BR(t 1 → bχ ± 1 )) and thus allowing for a best fit involving a slightly lighter stop and neutralino than in BP1.
• Stop/sbottom scenario. The signal is consistent with the presence of relatively light, mostly left-chiral stops and sbottoms, with mt 1 ≈ mb 1 ≈ 650 − 850 GeV, neu-tralino LSP (bino, higgsino, or an admixture) of mass ∼ 350 − 480 GeV, and the lightest chargino or second-lightest neutralino not far above the LSP. The ATLAS excess is fitted assuming pair-production of left-handed stops and sbottoms, with the subsequent decay chain depending on the mass splitting betweent 1 ,b 1 , and the light electroweakinos. Several possibilities are allowed:t 1 → tχ 0 1 ,t 1 → tχ 0 2,3 , andb 1 → tχ ± 1 , with branching ratios in all cases varying between 50% and 100%. The charginos and neutralinos just above the LSP decay to χ 0 1 and soft jets via an off-shell gauge boson (BR(χ 0 2 , χ ± 1 → χ 0 1 + jets) ≈ 50 − 60%) so that the acceptance/efficiency for finding an isolated hard lepton in the event does not drop significantly. The benchmark point for this scenario is shown in Table 2 as BP3.
BP3 is at risk of being excluded by the CMS α T search [35], as emerges by the interpretation of this search in the SMS T2bb with a light sbottom and the LSP. On the other hand, for BP3 the efficiency of the α T search to all-hadronic final states is reduced by a factor 2.3 with respect to the T2bb SMS, due to the lepton veto. Only a full simulation of [35] can quantitatively assess the impact of the efficiency reduction on the 95% C.L. exclusion bound. As we explained in Sec. 2, performing a full simulation of every available search far exceeds the purpose of this paper, so that for the moment we limit ourselves to pointing out the existence of this possible tension.
The ATLAS signal can also be fitted when the right-chiral sbottom is light as well. In such a case different mass hierarchies can arise betweent 1 ,b 1 andb 2 , with all three sparticles contributing to the signal. The benchmark point for this scenario is shown in Table 2 as BP4.
• Squark decay scenario. There exists the possibility of fitting the signal with spectra moderately heavier than those considered so far. In the case where the right-chiral stop and the LSP lie above ∼ 600 GeV, the signal can be fitted in the presence of squarks of the first two generations with mq L ≈ 1.1 TeV, which can boost the production cross section. An example of this light squark spectrum is shown in Table 2 and Fig. 1 as BP5. The decay chain proceeds throughq L → χ ± 1 q → bt 1 q (with approximately 60% branching ratio) and the lepton tagged by the ATLAS search originates from the subsequent decay of the stop.
A variation of the above scenario, involving left-chiral stops/sbottoms instead, is presented as BP6. The hard lepton comes now from the decay of the chargino or secondlightest neutralino, as the cascade proceeds throughq L → χ ± 1 (χ 0 2 ) q → tb 1 (t 1 ) q, with a final soft decay ofb 1 ort 1 .
Note that, while in BP5 the spectrum must include additional light sleptons or squarks in mass very close to the LSP if the upper bound on the relic density is to be satisfied, BP6 contains all the ingredients in it, as the correct relic abundance results from co-annihilation with the lightest stop/sbottom. This scenario is not excluded by the ATLAS searches described in Sec. 2. Gluino pairproduction provides in this case a cross section of the right magnitude, and the decay chain followsg → tt 1 , with a subsequent decay of the stop to the LSP and a c-jet (with flavor changing loop-generated coupling). An exemplary spectrum is characterized by mg = 1160 GeV, mt 1 = 858 GeV, and m χ 0 1 = 830 GeV. We do not report this scenario in Table 2 and Fig. 1 as it seems to be excluded by CMS all-hadronic searches [33,34,35] when the bound interpreted in the SMS T2bbbb is taken at face value. However, one must be aware that the acceptance/efficiency to all-hadronic final states is in this case reduced by a factor of 1.6 with respect to the T2bbbb SMS.

Approximate global analysis
The excess appearing in the DM-low bin of the ATLAS 1-lepton search shows a local significance greater than 3σ. We have shown in Sec. 3 that the excess admits an explanation in the context of the MSSM for several decay chains. On the other hand, one may wonder whether the emerging signal will maintain a noticeable significance with respect to the SM when other searches are considered and combined in an approximate global analysis. We attempt to do so in this section.
To begin with, let us notice that many of the kinematical bins constructed by the experimental collaborations show either downward-or upward fluctuations in the number of the observed events. It means that the SM (which corresponds to the background-only hypothesis) does not need to present the best fit to the observed data.
To quantify to what extent the fit improves (or worsens) if a putative SUSY signal is assumed, we proceed as follows. For each search we construct a test statistics as in Eq. (2), where again ∆χ 2 i are defined as −2 log(L i /L 0 ). In every bin i the numerator, L i , will either assume the value corresponding to a benchmark SUSY signal, s i , or to the SM background, s i = 0. The denominator, L 0 , will be instead given by the maximum likelihood value, obtained when s i ≈ o i − b i . Note that this approach may seem to differ from the one employed in the previous sections to fit the 1-lepton signal and to derive the exclusion bounds from the other ATLAS searches, where we always normalized the likelihoood function to the background-only hypothesis. We do so for two reasons. First, we would like to have only positive test statistics. Second, it allows us to quantify how far the SM predictions are from the data. One should also keep in mind that changes in the normalization factor result only in a shift of the χ 2 value, which affect in equal ways the signal and background-only hypotheses.
When combining the test statistics for individual bins, we always follow the experimental prescription. Thus, if for a given search the experimental collaboration explicitly states that the presented bins are independent, our test statistics will be the sum of the ∆χ 2 i ; if, on the contrary, the experimental collaboration presents bins that are not independent, we construct the test statistics as is done in the experimental paper to derive the exclusion bounds. Depending on the search, this is given by Eq. (2) or some modifications of it. The quantity of relevance to assess the goodness of fit will be the relative difference between the global test statistics of the SM and that of the benchmark points.  Table 3: Breakdown of the χ 2 contributions due to different searches for the background only hypothesis (SM) and the benchmark points.
In Table 3 we show the breakdown of the χ 2 values for the implemented LHC searches for the SM and for the postulated SUSY signals. In this presentation we do not consider the results from the 3.2 fb −1 analyses, if a higher luminosity update is available, and we also neglect searches that show a negligible χ 2 contribution. The complete χ 2 breakdown for every bin of every search considered in Table 3 can be found in Appendix B.
As expected, in the ATLAS 1-lepton search there is a large contribution to the SM χ 2 coming from the bin with the largest excess, DM-low. At the same time, the BSM signals fit all bins to good accuracy, as was discussed in Sec. 3.
For the corresponding analysis by CMS, the difference between the SM and the putative signals is not as marked, and it is encouraging to see that the benchmark points are not strongly disfavored with respect to the SM even in the absence of a noticeable excess above the background. Note that contributions to the SM χ 2 in the CMS search originate in some bins from an observed downward fluctuation of the background, in others from slight excesses. On the other hand, χ 2 contributions from the bins with downward fluctuations are smoothed out by large uncertainties in the background determination so that they do not significantly spoil the fit to the benchmark signals significantly. At the same time, the benchmark signals fit the small excesses better than the SM. The overall result is that in the CMS search the goodness of the fit is only slightly worse than in the SM, and for BP5 it is actually better.
To derive the total significance of the ATLAS excess after the CMS data is taken into account we combine both searches by adding the individual χ 2 contributions, as the results of the two collaborations can be safely assumed to be independent. The total significance, calculated from the chi-squared difference with the SM, is for all 6 benchmark points greater than 3σ.
When we add the other ATLAS searches to the global fit, the significance of the BSM signal with respect to the SM predictions can either improve or worsen, depending on the particular benchmark model. For the 0 lepton + b-jets search [28] almost all the models describe the experimental results better than the background-only hypothesis. That is once more due to several small excesses observed in this search that are better fitted by the putative SUSY signals than by the SM background. This could be a hint that the 1-lepton excess is a real phenomenon, but most likely the good fit results from the fact that both searches are optimized for stop production and the background determination is to some degree correlated in one and the other search.
All benchmark models show instead some tension with the 1 lepton + jets analysis [25]. The tension originiates from one bin, optimized for squark pair production, in which no excess is observed. Kinematical cuts employed in this bin are very similiar to those in the bin DM-low of the ATLAS 1-lepton search, so large signals for our benchmark models are to be expected. The tension is the strongest for the scenarios BP2 and BP6, where it reaches approximately the 2σ level. Note that these are the models that predict the largest signals in the DM-low bin. Additionally, models BP5 and BP6 are disfavored at 1.9 and 1.5 σ, respectively, by the all hadronic search [23]. That is also to be expected, as both of them present relatively light left-chiral squarks for which this search is optimized.
In the two bottom rows we present the total χ 2 for all the searches listed in Table 3, and the resulting total significance of the ATLAS 1-lepton excess. When combining the results of different analyses, we use two approaches. First, we assume that all ATLAS searches are statistically independent and we simply sum the corresponding χ 2 contributions. We mark the results thus obtained with the subscript "add." As a result, the significance of the 1-lepton excess decreases for BP2, BP5, and BP6.
As we have mentioned above, however, from the outside it is hard to gauge to what extent ATLAS searches, even the ones with different final states topologies, can be considered independent from one another. The ATLAS Collaboration itself, in their global analysis of 20 searches based on the 8 TeV data [4], decided to employ for overall exclusion limits a so called "best-of" strategy, which only uses the result of the analysis with the best sensitivity. To mimic the same approach, we also provide here significances based on the maximum χ 2 only, dubbed here as "max." Since we are not able to quantify the correlations between different ATLAS searches, we think it is safe to assume that the real significance lies somewhere in between the "add" and "max" scenarios.
One might wonder if the scenarios described above give rise to specific signatures at the LHC (besides a clear 1-lepton signal) that might help distinguish them from one another. Since all scenarios involve the production of sparticles with color charge, 0-lepton searches with different kinematical variables should produce a complementary, if weaker, signal when the integrated luminosity is increased by a factor ∼ 3 or more with respect to the amount of data presented at ICHEP.
In particular, by taking the 0-lepton ATLAS search [23] as an example, and rescaling there the SM background to ∼ 50 fb −1 , our numerical simulation shows bins with clear excesses above the background for BP5, BP6, and possibly BP3, where we have ordered the benchmarks by decreasing significance. Whereas, with the same luminosity, BP1, BP2, and BP4 would still be indistinguishable from the background in a search like [23]. With ∼ 100 fb −1 BP5 and BP6 will probably produce a > 5 σ discovery in the 0-lepton search, and significant excesses would appear for BP3 and BP2, in decreasing order of signal strength. Note, however, that even at this high luminosity, BP1 will remain difficult to observe in searches with a hard lepton veto, because of its large branching ratio to the 1-lepton final state, which allows for a clear fitting of our signal in the first place. The reader should also keep in mind that, at least for the data provided at ∼ 13 fb −1 , the corresponding CMS analyses have shown better sensitivity in the parameter space characterized by a moderately heavy LSP so that the benchmark points should appear earlier and to greater significance in a search like the CMS α T .
Finding appropriate UV completions to the scenarios presented above is beyond the purpose of this paper, although it might become a necessary endeavour were the excesses to be confirmed in the next round of data. Models with light stops in general do not fare well with boundary conditions defined at the GUT scale, because the current bound on the gluino mass (mg ≥ 1.6 − 1.8 TeV for the most common SUSY spectra) translates into a large renormalization of the stop mass at the low scale, independently of the initial choice for GUT-scale boundary conditions. Gauge mediation remains a viable option, as it allows some freedom in the choice of the messenger scale, and in particular models with matter-messenger mixing like the one proposed in, e.g., [38] and several other papers, can easily produce right-chiral stops and binos in agreement with the spectrum of BP1. The remaining benchmark points, however, present more involved and often compressed spectra, which might prove more challenging from the model-building point of view.

Summary and conclusions
In this paper we have investigated the possibility that some mild excesses over the SM background that emerged in the bins of the ATLAS 1-lepton + (b-)jets + E miss T search reported at the ICHEP 2016 conference are due to the first appearance of supersymmetry. Working in the framework of the p19MSSM, we have determined a few different types of SUSY spectra than can explain the observed data and are not excluded by other direct searches for SUSY.
The most straightforward possibility (BP1 in Table 2 and Fig. 1) corresponds to the simplified model spectrum used by the ATLAS Collaboration for the interpretation of the search results, which can fit the putative signal with light right-chiral stops of about 600-800 GeV and a bino-like neutralino LSP with a mass of the order of ∼ 400 GeV. One of the advantages of this scenario, beyond its simplicity, is that it can be relatively easily embedded in known UV completions, for example models of gauge mediation with matter-messenger mixing terms in the superpotential that allow for light righthanded stops.
The signal can also be fitted if charginos of mass in between the stop and the neutralino are present in the spectrum (BP2), where in this case the stop can be slightly lighter than in BP1, as a light chargino has the effect of reducing the branching ratio to the 1-lepton final state topology. Additional good fits can be achieved when the signal is produced by left chiral stop/sbottoms, because of the boost in production cross section that arises by including theb 1 → tχ ± 1 topology (BP3), and when the right-chiral sbottom is also added to the light spectrum (BP4).
For an LSP of mass significantly above 400 − 500 GeV the acceptance for the signal produced by relatively light stops and sbottoms drops drastically. More complex spectra, however, involving associated production of light-generation squarks and electroweakinos with a mass of about 1.1 TeV and cascade decays to stops and sbottoms (BP5, BP6) can still produce the right signal even for a ∼ 700 GeV LSP, because of the increase in production cross section due to the presence of four degenerate "valence" squarks. These light valence squarks are in pole position to be probed by all-hadronic searches in the next round of the LHC data.
For all of the presented benchmark-point scenarios we have compared their global χ 2 with respect to the ATLAS and CMS searches listed in Table 3 with the corresponding χ 2 of the background-only hypothesis. Our computation shows that the signal hypothesis is, at least slightly, favored globally with respect to the SM, and this is true under different choice of test statistics.
On the other hand, we conclude by repeating once again that the analysis presented here is meant to investigate a possibility that might reveal itself as a statistical fluctuation. On the positive side, if the signal is real it will be confirmed with more data and the signatures described in Sec. 3 will be useful to distinguish one or the other scenario.

ACKNOWLEDGMENTS
We thank the anonymous Referee for suggesting the global analysis of Sec. 4. K.K. is supported in part by the DFG Research Unit FOR 1873 "Quark Flavour Physics and Effective Field Theories". The use of the CIS computer cluster at the National Centre for Nuclear Research in Warsaw is gratefully acknowledged.

A Validation of the CMS 1 lepton search
In this appendix we provide the validation of our numerical simulation of the CMS search for direct top squark pair production in the single-lepton final state [6]. Validation of all other searches used in this paper and listed in Sec. 2 can be found in Appendix A of Ref. [5].
We check the consistency of our statistical analysis with the official bounds provided by CMS by comparing the 95% C.L. exclusion line, derived by the experimental collaboration in one of the SMS used for their result interpretation, with the line obtained   with our code in the same SMS. The technical details of the experimental analysis can be found in [6].
The SMS analyzed here consists of a stop decaying into top quark and neutralino LSP. The experimental signature includes one lepton, jets and large amount of missing energy. The analysis is performed with the following pre-selection cuts: • exactly one signal lepton with p T > 20 GeV and |η| < 2.4(1.44) for muons (electrons), • at least 2 signal jets with p T > 30 GeV and |η| < 2.4, • at least 1 signal b-jet, • missing energy E miss T > 250 GeV.
The kinematical variables used to discriminate between the signal and the background are: azimuthal angle ∆φ(E miss T , jet) between two leading jets and missing energy; the transverse mass m T of the signal lepton and missing transverse momentum; the magnitude of the negative vectorial sum of the transverse momentum of jets with p T > 20 in the event, H miss T , and the variable M W T2 defined in [39]. To implement the latter in the recast tool, we used the C++ code provided in [39]. Fifteen exclusive signal regions are defined with various jet multiplicities and binned in E miss T and M W T2 . We present in Fig. 2 the validation of our simulation in terms of the exclusion limits in the parameter space (mt 1 , m χ 0 1 ) for the direct stop production scenario and a subsequent decay into top quark and neutralino LSP. Gray diamonds represent the points excluded by our likelihood function at the 99.7% C.L., cyan circles are excluded at the 95.0% C.L., and blue triangles are excluded at the 68.3% C.L. The solid black line shows the 95% C.L. CMS exclusion limit, which we overlap for comparison.

B Summary of all ∆χ 2 contributions
In Tables 4-10 we present the complete χ 2 breakdown for every bin of every search considered in Table 3. The s i column shows the number of signal events simulated according to the procedure of Sec. 2. The total ∆χ 2 for each search is calculated following the prescription given in the experimental paper and summarized in the corresponding caption.

Bin
Obs.  Table 4: Breakdown of the χ 2 contributions due to different signal bins of the CMS 1-lepton search [6] for the background-only hypothesis (SM) and the benchmark points. The total ∆χ 2 is given by the sum of each individual contribution. The collaboration has not labeled the bins by name. Nevertheless, we do provide an indicative label in column 1, based on the number of jets and the E miss T cut adopted.      Table 7: Breakdown of the χ 2 contributions due to different signal bins of the ATLAS 0 leptons + (b-)jets search [28] for the background-only hypothesis (SM) and the benchmark points. Bins 1, 2, and 3 are orthogonal to each other, and so are bins 4, 5, and 6. The total ∆χ 2 is given by max{ 3 i=1 ∆χ 2 i , 6 i=4 ∆χ 2 i , ∆χ 2 7 , ∆χ 2 8 , ∆χ 2 9 , ∆χ 2 10 , ∆χ 2 11 }. The bins are labeled according to their ATLAS name.  Table 8: Breakdown of the χ 2 contributions due to different signal bins of the ATLAS 3 btagged jets search [27] for the background-only hypothesis (SM) and the benchmark points. For the first two bins the largest ∆χ 2 i is chosen. The chi-squared of bins 3 and 4 are instead confronted in a pairwise manner with the chi-squared of bins 5, 6, and 7, and the combination with the largest contribution is selected. The total ∆χ 2 is given by the sum of the two partial chi-squared. The bins are labeled according to their ATLAS name.  Table 9: Breakdown of the χ 2 contributions due to different signal bins of the ATLAS 1 lepton + jets search [25] for the background-only hypothesis (SM) and the benchmark points. The total ∆χ 2 is given by the maximum of the individual contributions. The bins are labeled according to their ATLAS name.  Table 10: Breakdown of the χ 2 contributions due to different signal bins of the ATLAS 2 btagged jets search [31] for the background-only hypothesis (SM) and the benchmark points. The total ∆χ 2 is given by the maximum of the individual contributions. The bins are labeled according to their ATLAS name.