Analysing parameter space correlations of recent 13 TeV gluino and squark searches in the pMSSM

This paper examines unexplored correlations in the parameter spaces probed by recent ATLAS analyses for gluinos and squarks, addressing various shortcomings in the literature. Six 13 TeV ATLAS analyses based on 3.2 fb$^{-1}$ of integrated luminosity are interpreted in the 19-parameter R-parity conserving phenomenological minimal supersymmetric extension to the Standard Model (pMSSM). The distinct regions covered by each search are independent of prior, and we reveal particularly striking complementarity between the 2-6 jets and Multi-b searches. In the leptonic searches, we identify better sensitivity to models than those used for analysis optimisation, notably a squark-slepton-wino scenario for the SS/3L search. Further, we show how collider searches for coloured states probe the structure of the pMSSM dark sector more extensively than the Monojet analysis alone, with sensitivity to parameter spaces that are challenging for direct detection experiments.

This paper extends the ATLAS pMSSM study [53] using fast detector simulation to include combined constraints from six early 13 TeV searches, and we make this exclusion information available online [55]. Instead of preliminary results with more data, we opt to use published results based on 3.2 fb −1 of integrated luminosity, which sets the scope of our study to gluinos and light flavour squarks. Our analysis reveals previously unexamined correlations in the pMSSM space that are already sufficiently elaborate to merit detailed evaluation of search strategies. These results are also of interest to the CMS collaboration and direct dark matter detection experiments. We organise our discussion around the following questions that address various shortcomings in the literature: (a) How distinct are the regions of parameter space being probed by individual analyses? Interpretations using the pMSSM often present combined constraints from multiple searches as fractions of models excluded [39,48,53,54,56]. Overlap matrices were recently used in the literature [53,55] to quantify the complementarity of these searches, namely the fractional exclusion of the same subset of points by two analyses. However, this marginalisation not only obscures which analyses had greatest sensitivity to different pMSSM subspaces, but also depends on the prior distribution from the parameter scan and non-LHC constraints.
To what extent are analyses over-optimising to a set of simplified models, which may preclude sensitivity to a wider class of scenarios? The pMSSM offers a greater variety of decay chains, such as those with suppressed branching fractions and placing those with more intermediate sparticles on-shell, which alter kinematics. How well these simplified model oriented searches are capturing the wider classes of signatures remains relatively unexplored. (c) What neutralino DM scenarios can be probed competitively by 13 TeV collider searches for coloured sparticles? Simplified dark matter models are stimulating the collider frontier for DM searches [57][58][59][60][61][62][63][64][65][66][67]. ATLAS uses these to perform an explicit DM interpretation for the Monojet search [6]. However, its sensitivity is greatly altered in richer dark sectors of the pMSSM, shaped by the composition of the neutralinoχ 0 1 lightest SUSY particle (LSP) [68][69][70][71][72] being dominantly bino, Higgsino or wino (defined in Table 2 of Appendix A.1). Furthermore, studies of neutralino DM often combine exclusions of multiple analyses [48,53] or only focus on the electroweakino sector [52,[73][74][75], omitting potentially important coannihilation roles of other sparticles. Section 2 addresses question (a) by directly correlating the most sensitive of the six analyses considered with the masses of the gluino, LSP and lightest squark. We examine the individual searches that provide sensitivity to distinct regions within these mass plane projections. The regions identified are both prior-independent and reveal significantly richer information than the overlap matrices previously used in literature. Subsection 2.3.2 addresses question (b) by considering the analyses that select leptonic events as a case study. We examine to what extent the simplified models used by the 1-lepton analysis [3] and the same-sign or 3-lepton (SS/3L) search [5] map onto pMSSM points, and identify scenarios beyond those considered for analysis optimisation.
Section 3 addresses question (c) by ascribing DM interpretations to each 13 TeV search considered, allowing for comparisons to the Monojet search. Striking correlations are exhibited, and we discuss how distinct decay cascades are influenced by the bino, Higgsino or wino content of the LSP, while coloured sparticles may act as early universe coannihilators [76]. Further, we identify the collider parameter space affected by recent limits from Xenon-target direct detection experiments [77,78]. Section 4 summarises the conclusions of this work. The pMSSM points we investigate were produced for the ATLAS study [53]. The 19 parameters of the R-parity conserving MSSM, where all flavour and CP violation resides in the CKM matrix, were scanned with flat priors with sparticle mass scales capped at 4 TeV [53]. However, the distinct regions of sensitivity we identify in our results are independent of sampling prior. Constraints were imposed from LEP searches [79], precision electroweak measurements [80][81][82][83][84][85][86][87][88][89], heavy flavour physics [90][91][92][93][94][95][96][97], DM direct detection [98][99][100] and the Planck relic density [101] upper bound 1 . Further details of the theoretical assumptions and experimental constraints used by ATLAS [53] may be found in Appendix A.1. We consider the 181.8k points that survive Run 1 constraints from Ref. [53], after having removed models with long-lived (cτ > 1 mm as defined in Ref. [53]) gluinos, squarks and sleptons as these require dedicated Monte-Carlo simulation. The methodology we employed to interpret the six 13 TeV searches using fast detector simulation is detailed in Appendix A.2.

Complementarity between early 13 TeV searches
For the six 13 TeV (3.2 fb −1 ) ATLAS searches considered in Table 1, we present the regions of sensitivity between these analyses. In Appendix A.3, we present results using existing practices in the literature, based on fractions of models excluded (marginalised distributions) and overlap matrices. This section addresses various shortcomings of these approaches as discussed in the Introduction. Subsection 2.1 defines and quantifies the 'most sensitive analysis' used to exclude the points in our interpretation. Following this, Subsection 2.2 projects this information into 2-dimensional subspaces of the pMSSM involving gluinos, squarks and the LSP, and discusses prior features in these planes. Finally, we examine the complementary sensitivity of each analysis to distinct regions of pMSSM parameter space in Subsection 2.3, partitioning our discussion between analyses that veto events with leptons with those that select on them.
Henceforth in this paper, 'squark'q refers to only the lightest superpartner of the left-or right-handed quark of the first or second generationsq ∈ {ũ,d,c,s} L,R . Similarly, 'slepton'˜ refers to the lightest superpartner of the left-or right-handed lepton of the first or second generatioñ ∈ {ẽ L ,ν eL ,μ L ,ν µL ,ẽ R ,μ R }. (1)

Most sensitive analyses used for exclusion
In this study, a point is deemed excluded at 95% confidence level if at least one analysis returned a CLs value less than 0.05, using the CLs prescription [102]. Of the 181.8k points that survived Run 1 [53], a total of N Excluded 13 TeV 3.2 fb −1 = 28.5k were excluded by our interpretation of the six 13 TeV analyses. We take the analysis with the smallest CLs value as the 'most sensitive analysis' used to exclude the point. In Table 1, we normalise the number of points satisfying this for each analysis N Lowest CLs to the total excluded N Excluded 13 TeV 3.2 fb −1 . For example, the Multi-b search was the analysis with the lowest CLs value for 23% of the 28.5k excluded points. Indeed, almost 95% of the excluded points have either the 2-6 jets or Multi-b searches being the most sensitive. However, care must be taken when interpreting these fractions, as they are prior dependent and correlated with non-LHC constraints. The fractions indicate the relative number of points in forthcoming figures.
Less than 4% of the excluded points have two or more analyses associated with the same smallest CLs value. In such cases, the 'most sensitive analysis' is randomly chosen from this subset of analyses with smallest CLs value to minimise systematic selection bias. In the vast majority of these situations, this is done because the analyses share a CLs value of 0.0.
Also displayed in Table 1 is the total number of excluded points by each analysis N Excluded 13 TeV 3.2 fb −1 out of those that survived Run 1 N Survived ATLAS Run 1 . The relative overlap between analyses is quantified in Table 4 in Appendix A.3. Importantly in Table 1, all analyses retain non-zero percentages in the total fraction of points uniquely excluded by each analysis and none of the other five N Excluded Uniquely 13 TeV 3.2 fb −1 . This emphasises the importance of maintaining a broad programme of searches. While most of the searches optimise for gluino production, and some have dedicated signal regions for squarks, each search maintains unique sensitivity to particular classes of signatures within the MSSM, featuring distinct final states and kinematic regimes.

Features in mass plane projections
The distinct regions of sensitivity for each of the six searches become unambiguous when we project into various twodimensional subspaces of the pMSSM. Figures 1 and 2 display each excluded point styled and coloured according to the analysis that returned the lowest CLs value, projected into the mass planes of gluino vs LSP, gluino vs squark and squark vs LSP respectively. Due to their large numbers, the 2-6 jets and Multi-b analyses are allocated smaller markers to improve clarity of other points.
Before discussing the sensitivity of each analysis in turn (Section 2.3), we note that the two-dimensional distributions (but not the regions to be discussed in Sections 2.3 and 3.1) of excluded points do depend on the flat prior and importance sampling ATLAS used to scan the parameter space [53], in addition to non-LHC constraints. This is especially apparent in the mass planes involving the LSP (Figures 1a and 2), and we comment on several features common to these figures to aid interpretation: -Within around 100 GeV of the grey hatched regions where 100% of models in the mass bins were excluded by Run 1 ATLAS searches, very few points are present. The number of models close to this boundary is often less than ten 2 . Typically, such points marginally survived Run 1 constraints and many have since been excluded by the six 13 TeV searches. -There is a visible break in points at m(χ 0 1 ) ∼ 100 GeV. Below this mass, there are few models with Higgsino-or wino-like LSP as these generally have a near-degenerate chargino, which are excluded by LEP direct searches. For the definition of the three LSP classifications, see Table 2 in Appendix A.1.
-The points below m(χ 0 1 ) ∼ 100 GeV are therefore predominantly bino-like LSP models. These LSPs are close to half the mass of the Z 0 or Higgs boson h 0 so they can undergo resonant annihilation in the early universe in so-called 'funnel regions' to satisfy relic abundance constraints.
-A large density of points is visible along the diagonals where the mass splitting between the LSP and the gluino or squark become small ( 50 GeV Table 1. The exclusion information is from Ref. [55]. These are projected into the mass planes described in the captions. Here m(q) denotes the lightest squark of the first or second generations. The markers are styled according to the analysis with the lowest CLs value: 7-10 jets (magenta triangle), 1-lepton (orange plus), SS/3L (green cross), Monojet (cyan ring), 2-6 jets (blue filled square), and Multi-b (red filled circle). Due to their large numbers, the markers for the 2-6 jets and Multi-b analyses are reduced in size for clarity. The hatched grey regions indicate the mass bins where all model points were completely excluded by Run 1 searches considered in Ref. [53]. Models with long-lived gluinos, squarks or sleptons are removed from these figures.  more apparent by using importance sampling to ensure the number of bino-like LSP models is of the same order as those of Higgsino-and wino-like LSP; see Ref. [53] for details.

Mass plane correlations of searches for squarks and gluinos
This subsection discusses the distinct regions of sensitivity correlated with each analysis considered. Importantly, the identified regions in the two-dimensional planes are independent of the prior distribution of points.

Searches with a lepton veto
First, we investigate the excluded points where each of the 2-6 jets, Multi-b, 7-10 jets, and Monojet searches are most sensitive, in that order. These analyses all veto on leptons, and exhibit the most unambiguous correlations in the mass planes of the gluino, squark and LSP (Figures 1 and 2). Generally, these analyses select events with transverse momentum imbalance of magnitude E miss T , varying jet multiplicity, together with discriminants dependent on mass scale or flavour. For full details, refer to the ATLAS references in Table 1. Figure 1a reveals the extensive distribution in the gluino-LSP plane where the 2-6 jets is the most sensitive analysis (blue points). The 2-6 jets search uses the effective mass discriminant, with varying degrees of minimum jet multiplicity to target models where a squark or gluino directly decays to the LSP. Meanwhile, the larger jet multiplicity regions target scenarios where a chargino mediates the decay of the gluino to the LSP. This analysis has almost exclusive sensitivity to points with gluinog masses below 1 TeV, where its mass splittings with the LSP in the range 25 m(g) − m(χ 0 1 ) 500 GeV. Larger mass splitting scenarios for sub-TeV mass gluinos are excluded by Run 1 searches (hatched grey mass regions). For regions with gluinos above 1 TeV with large gluino-LSP mass splittings, we find the 2-6 jets has reduced sensitivity compared with other analyses, especially for LSP masses m(χ 0 1 ) 500 GeV. For gluino masses m(g) 2 TeV, we expect reduced sensitivity to gluinos, but points in this region are correlated with excluded points involving low mass squarks. This is confirmed by comparing with the gluinosquark plane (Figure 1b). Indeed, the 2-6 jets search is also predominantly the most sensitive analysis for light squarks in regions not far beyond Run 1 sensitivity. High mass gluinos in such scenarios can nevertheless contribute to production cross-sections of the squarks as a mediator via t-channel diagrams. The gluino-squark plane (Figure 1b) also reveals a vertical strip with a lower density of points around gluino mass between about 1 TeV and 1.2 TeV. This corresponds to a region where a lower fraction of models is excluded per mass interval (see also Figures 11c and 13c in Appendix A). This reduced sensitivity, and lack of any other dedicated analysis probing this region, is due to gluino-LSP mass being moderately small (between around 25 and 200 GeV). Such scenarios are challenging for traditional 'missing energy plus jets' searches, but represent the greatest potential for high luminosity where novel techniques are being developed to target such regions [104].

2-6 jets
Further interpretations for the 2-6 jets analysis apply in the squark-LSP plane ( Figure 2). The search has sensitivity to a wide variety of squark mass scenarios, gradually reducing above m(q) ∼ 1 TeV. The strip of blue points for m (q) 1.5 TeV and 600 m χ 0 1 900 GeV where we would expect reduced squark sensitivity are correlated with low mass gluinos. These points largely have sub-TeV gluinos, an interpretation, again confirmed in the gluinosquark plane (Figure 1b).

Multi-b
Where the Multi-b analysis has most sensitivity (red points) is localised to regions complementary to the 2-6 jets analysis. In the gluino-LSP plane (Figure 1a), the region is enclosed by gluino masses of 1.2 m(g) 1.7 TeV and LSP masses of 100 m χ 0 1 500 GeV. This sensitivity also extends to smaller gluino-LSP mass splittings, though to a lesser extent.
The distinctiveness of this region of sensitivity remains in the gluino-squark plane ( Figure 1b). Here, points where the Multi-b analysis is most sensitive are highly correlated with squark masses above 1.5 TeV, largely untouched by other analyses, being strikingly separated from 2-6 jets and to some extent the 7-10 jets analysis. Figure 2 also confirms such correlations of the Multi-b in the squark-LSP plane and again sensitivity to squarks above 1.5 TeV is largely from gluino rather than squark production.
These correlations arise from the Multi-b search targetingg → bbχ 0 1 andg → ttχ 0 1 models with heavy gluino production decaying via off-shell stopst and sbottomsb. This analysis selects events with large E miss T and at least three jets originating from bottom quarks. A subset of the signal regions selects loosely on boosted top quarks decaying from gluinos, including events enriched with lepton presence. Figure 3 takes the points where the Multi-b is most sensitive and illuminates the mass distributions of various pertinent sparticles. The light flavour squarksq are centred around 2 TeV with gluinos around 1.4 TeV, where Figure 1b indicates that squarks are predominantly heavier than gluinos. This fact suppresses gluino decays to light flavour quarks, which proceed via the three bodyg → qqχ 0 1 process. By contrast, the mass distribution of sbottomsb 1 peaks around 800 GeV in Figure 3 and have a preference to be lighter than the gluinos, allowing on-shellg →bb decays to be favoured. Indeed, this demonstrates favourable sensitivity to the on-shell counterpart of the simplified models considered for optimisation by ATLAS. Meanwhile, the distribution of stop masses is relatively uniform for m(t 1 ) 900 GeV. The requirement of three or more jets originating from bottom quarks therefore favours such scenarios. We furthermore note that out of the models excluded where the Multi-b is most sensitive, 56% models have a Higgsinos-like LSP while 22% are wino-like, which have light charginos consistent with the corresponding mass distribution in Figure 3. This preference of the Multi-b analysis for Higgsinolike LSP models can be understood by the higher Yukawa couplings to heavy flavour quarks, which also enhance decays of the gluino to bottom and/or top quarks.
From the region in Figure 1 where the density of red points is greatest, we show a representative point with model number 148229034 (Figure 4a). This contains an LSP with a bino-Higgsino mixture at a mass of 175 GeV, a relatively low mass 1.2 TeV gluino and a 1.3 TeV sbottom enabling g → bbχ 0 1 branching ratios to be preferred.

7-10 jets
The concentration of magenta triangles in Figure 1 illustrates that the 7-10 jets search is the most sensitive analysis for models having light LSPs with large mass splitting from the gluino. Most of these LSPs have masses below 100 GeV, which are bino-like as discussed in Section 2.2, and undergo early universe annihilation through a Z 0 or a h 0 boson in so-called 'funnel' regions. Since this annihilation mechanism requires Higgsino or wino admixtures to pro-  where the corresponding analysis in parentheses had greatest sensitivity. These show the most relevant low-mass sparticles and their decays; for the full mass spectra, see Figure 15 in Appendix C. The grey arrows show the branching ratios between two sparticles and proportional to its brightness. The model displayed for the 7-10 jets analysis is the one with the lightest gluino with an LSP mass below 100 GeV, while for the SS/3L search, we selected the model with the lightest squark. For the Multi-b analysis, displayed is a model selected from the region in figures 1 where the density of red points is greatest. The model for the 1-lepton search is chosen based on one that is representative from the mass distributions in Figure 5a.
ceed, such scenarios typically have light charginosχ ± 1 and next-to-lightest neutralinosχ 0 2 , with relatively small mass splittings between one other.
These types of decays are consistent with the simplified models containing long decay chains that this analysis optimises for. The observed missing transverse energy E miss T therefore tends to be smaller than that required by the 2-6 jets or Multi-b searches. Indeed no explicit requirement on E miss T is made by the search (instead the main discriminant is a ratio E miss T / ∑ p jet T involving the missing energy and scalar sum of transverse jet momentum). The 7-10 jet search also has a looser requirement on jets originating from bottom quarks compared with the Multi-b analysis. Together, this allows the 7-10 jets analysis to maintain a unique coverage of models. In the gluino-squark plane (Figure 1b), the 7-10 jets points occupy a similar space as Multi-b, but the longer cascade chains mean the gluino mass reach is lower. The remainder of the spectrum can be relatively decoupled. Figure 4b displays model number 227558023, which is representative of the models where the 7-10 jets analysis is most sensitive, where we selected the one with lowest gluino mass which had an LSP below 100 GeV.

Monojet
The dedicated Monojet analysis selects events with an energetic jet from initial-state radiation recoiling off a system of large missing transverse momentum and up to three additional jets. In terms of SUSY models, this is optimised for scenarios where the mass of the squark is almost equal to that of the LSP, so-called 'compressed scenarios'.
In Figures 1a and 2, the excluded points (cyan rings) where this analysis is most sensitive involve very small mass splittings between the coloured sparticle and LSP. Though there is significant overlap in regions of sensitivity for the Monojet and 2-6 jets analyses, the former is exclusively the most sensitive analysis involve squark-LSP splittings below 30 GeV. We note that the 2-6 jets search includes a similar signal region but requiring a minimum of two jets called '2jm'. The different jet multiplicity requirements ensure the dedicated Monojet maintains a unique sensitivity to the smallest squark-LSP mass splittings, again demonstrating the complementarity of searches. We find a small number (fewer than 10) of scenarios involving small mass splittings between the lightest 3rd generation squark and the LSP in the models for which the Monojet is most sensitive.

Searches selecting one or more leptons
Both the 1-lepton and SS/3L searches also target production gluinos and squarks, but require one or more leptons in selected events. In Figures 1 and 2 involving the gluino, squark and LSP masses, the correlations between where the 1-lepton (orange plus) or SS/3L (green cross) analyses were most sensitive are less obvious. There points tend to cluster below squark masses of m(q) 1.5 TeV, while few points are present for gluino masses m(g) 1.2 TeV where other analyses dominate. Further investigation reveals other correlations driven by the light flavour slepton and gaugino masses not apparent in these figures, which we discuss in what follows. We also identify scenarios where these searches had most sensitivity, which are beyond what ATLAS optimised for. Figure 5a shows the distribution of masses for various sparticles from the excluded points where the 1-lepton analysis had the lowest CLs value. The 1-lepton analysis requires events with exactly 1 electron or muon with various minimum jet multiplicities and large E miss T . Though the analysis was optimised for gluino production, there is also a prevalence of light squarks, whose distribution peaks around m(q) ∼ 1 TeV. The distributions of the chargino and to a less extent next-to-lightest neutralinoχ 0 2 are skewed towards low masses, peaking around 300 GeV, with a tail that extends above 1 TeV. The tendancy for the next-to-lightest neutralino and chargino to be light is characteristic of Higgsino content in the LSP, or a bino-like LSP with wino-like pair near in mass to the LSP.

1-lepton
The mass distribution also reveals a preponderance of light sleptons, peaking around m(˜ ) ∼ 600 GeV. We find these typically reside between the coloured sparticle and the chargino-LSP states. To demonstrate this interpretation, we define the splitting parameter where m ˜ and m (q) are respectively the masses of the lightest slepton and squark among the first or second generations. The red dashed line in Figure 6 then shows the distribution in x of models where 1-lepton was most sensitive in x, and indeed the majority of models in the histogram have 0 < x < 1. The 1-lepton analysis interprets in a simplified model with a gluino-chargino-LSP decay chain and optimises for various mass splittings between these three sparticles. This is noteworthy given we find the 1-lepton search is most sensitive to richer mass spectra, with squarks and intermediate sleptons feature prominently. Given the mass distributions in Figure 5a, we display model number 13382371 (Figure 4c) as a representative point to illustrate the wider sensitivity of the 1-lepton analysis. This point features a 727 GeV squark that can cascade to a slepton doublet and a 173 GeV Higgsino-like LSP. We also note that all the signal regions were involved in the exclusion of these models, indicating sensitivity to a wide variety of splittings. This addresses the question of over-optimisation raised in the Introduction (Section 1). Despite optimising robustly for the gluino-chargino-LSP model, this analysis was sensitive to wider classes of models.

SS/3L
Though there are similarities with the 1-lepton analysis, many salient differences appear for points where the SS/3L search is the most sensitive analysis. This search selects events with at least two leptons, and if there are exactly two, they are required to have the same electric charge.
One prominent difference is that points where the SS/3L is most sensitive again have almost exclusively wino-like LSPs, and so a nearly mass-degenerate chargino. This is consistent with the mass distributions of these models in Figure 5b: compared with the 1-lepton case (Figure 5a), the SS/3L has a smaller tail of high chargino mass while the next-to-lightest neutralinoχ 0 2 distribution is no longer skewed to lower masses. We note that a future Run 2 version of the 'Disappearing Track' analysis [16], not considered in this work, could be sensitive to these models, since the small chargino-LSP mass splitting in such models ensures that the chargino is typically long-lived on collider time scales.
For these points, the gluinos are all relatively heavy and not strongly correlated with a particular mass scale, being fairly uniformly distributed for m(g) 1.5 TeV, in contrast to the 1-lepton discussion. This implicates that the SS/3L is not the most sensitive analysis to light gluinos, where other analyses such as the 2-6 jets are most sensitive. By contrast, squarks retain a peaked distribution centred around 1 TeV.  )]  1 -˜ -q system, parametrised by x as defined in Eq. 2. These are presented for the excluded points where the 1-lepton (dashed red) and SS/3L (solid blue) analyses are most sensitive. Figure 6 demonstrates, in a similar way to the 1-lepton analysis, that the slepton mass is almost always between the chargino and the lightest squark. We find a negligible number of models (a single entry) has x > 1, corresponding to a slepton mass being above that of the lightest squark. Thus, we find that points where the SS/3L has best sensitivity are strongly correlated with one common feature: a squark-slepton-chargino-LSPq-˜ -χ ± 1 -χ 0 1 ordered mass spectrum. The squark can undergo a three-body decay to a quark, lepton and a sleptonq → q ˜ if the intermediate neutralino is off-shell m(χ 0 2 ) > m(q). This hierarchical structure is displayed in model number 11733067 (Figure 4d). Among the models where SS/3L is most sensitive, this has the lightest squark mass at 436 GeV. Moreover, we find one signal region 'SR0b3j' was used to exclude 98% of these models where the analysis was most sensitive. This suggests other analyses had better sensitivity to models targeted by the other three signal regions, for example the Multi-b analysis is particularly sensitive tog → ttχ 0 1 scenarios. The 'SR0b3j' signal region requires at least three leptons and was optimised to capture ag-χ 0 2 -˜ -χ 0 1 decay chain, distinct from theq-˜ -χ ± 1 -χ 0 1 scenario we just identified.
Taken together, we draw two noteworthy conclusions from our findings for the SS/3L search. First, this analysis is the most sensitive analysis for a different scenario we identified in the pMSSM than all the simplified models those used for analysis optimisaton. Second, the exclusivity of signal region used to exclude this scenario indicates that the SS/3L lacked competitive sensitivity to points in the pMSSM corresponding to the simplified models considered.

Probing the dark sector with strong SUSY searches
In the pMSSM, the relic abundance of dark matter is predominantly shaped by the gaugino and Higgsino content of the LSP. The ATLAS Monojet search makes an explicit DM interpretation in a simplified model, but the complexity of the pMSSM dark sector necessitates other collider searches to provide decisive sensitivity. Complementing electroweak SUSY searches [52], this section sets out to demonstrate that 13 TeV ATLAS searches for coloured sparticles can indirectly probe different dark matter scenarios in the pMSSM. Subsection 3.1 ascribes DM interpretations to each of the 13 TeV searches and manifests their sensitivity to distinct regions of DM parameter space. We also consider the impact of recent direct detection results from the PandaX-II/LUX experiments [77,78] on gluino and squark masses in Subsection 3.2. Present constraints from indirect detection experiments searching for LSP annihilation signals in galactic halos are expected to have little sensitivity on the pMSSM space considered 3 .
3.1 Impact of 13 TeV constraints on dark matter observables 3.1.1 Relic density Figure 7a shows each of the the 28.5k points excluded by the six 13 TeV SUSY searches considered, styled and coloured according to the analysis that had the lowest CLs value, projected into the plane of LSP mass vs relic density. Again, the strong correlations in this plane with the six searches are unambiguous. To facilitate interpretation, we discuss the most salient underlying processes that shape several features in this projection: -The LSP relic abundance Ω χ 0 1 h 2 is set by early universe thermal freeze-out, the hallmark of the weakly interacting massive particle (WIMP) paradigm. We do not require the neutralino to be the sole constituent of dark matter, as other well-motivated candidates such as axions can contribute [69], making the points considered more general than Ref. [56]. Thus the Planck measurement of the cold dark matter (CDM) abundance [53] only serves as an upper bound.
-The composition of the LSP strongly influences the early universe annihilation mechanism of the LSP and the resulting cosmological relic density. Figure 8a illustrates the relic density against the mass of the LSP that survived the constraints from six 13 TeV analyses, coloured by the dominant composition to the LSP, as defined in Table 2. Notably, the Higgsino-and wino-like LSP models are concentrated along a straight diagonal line in the plot. This is because the thermally averaged cross-section σ v ∝ m −2 χ 0 1 [70], and the relic density Ωχ0 1 h 2 is therefore nearly proportional to the LSP mass squared. Mass splittings between the LSP and coannihilating chargino are typically a few GeV (∼ 100s MeV) for Higgsino-(wino-) like LSPs. There are also no Higgsino-or wino-like 3 Such indirect detection constraints were examined in Ref. [48], which finds that current constraints from Fermi-LAT are not expected to constrain the pMSSM space. Nevertheless, it was found that the future Cerenkov Telescope Array (CTA) [106] is projected to have sensitivity to high mass LSPs scenarios not accessible to direct detection and collider searches.
LSPs below around 100 GeV. This is because such models have near-degenerate charginos to the LSP, which are excluded by LEP lower bounds. -By contrast, bino-like LSP models have suppressed earlyuniverse thermally-averaged annihilation cross-sections 4 , leading to larger relic abundance, as generally seen in Figure 8a. Therefore to satisfy the Planck bound, binolike LSPs must either have non-negligible mixing with winos and/or Higgsinos, or there must be a near-degenerate next-to-LSP to act as a coannihilator. In Ref. [53], the natures of the coannihilators were displayed in the planes involving relic density. No distinction was made between light flavour squarks and gluinos, yet they have different phenomenological roles as coannihilators. We therefore elucidate this in Figure 9, differentiating between coannihilators being gluinos (light orange), light flavour squarks (dark blue) and 3rd generation squarks (green). We observe that light gluino coannihilators are more stringently excluded, with few points m(g) 800 GeV than for squarks due to powerful Run 1 constraints. -Focusing now on electroweak particles ('Other', light grey), for LSP masses m χ 0 1 250 GeV, the coannihilation mechanism is predominantly via uncoloured sparticles due to stringent LHC constraints on squarks and gluinos. The two peaks centred around m χ 0 1 ≈ 45 and 63 GeV involve resonant annihilation through a Z 0 or Higgs boson. Meanwhile, for 90 m χ 0 1 250 GeV, the coannihilators are predominantly slepton or gauginos, which are bounded from below by the LEP limit.
Returning to the discussion of points excluded by individual analyses in the relic density plane (Figure 7a), the strong SUSY searches considered are sensitive to models with a wide variety of Ωχ0 1 h 2 . As previously discussed in Section 2.3, the 7-10 jets search is particularly sensitive to models with m χ 0 1 100 GeV, where bino-like LSPs are associated with the Z 0 and h 0 funnel region. Meanwhile, the SS/3L analysis is most sensitive to wino-like LSPs models while the Multi-b analysis had preferential sensitivity to Higgsino-like LSP scenarios, as indicated by the clustering of green crosses and red dots along the wino and Higgsino respective diagonal bands (compare with Figure 8a).
The 1-lepton analysis has sensitivity away from the Higgsino and wino diagonal bands, for LSP masses below about 4 Pure binos do not couple to any gauge or Higgs bosons. This is seen from the couplings g Zχ 0 1χ 0 1 (g hχ 1χ1 ) of the neutralino to the Z 0 and (Higgs) bosons are given at tree-level by [108] g Zχ 0 g hχ 0 1χ 0 1 = g 2 cos α (N 11 − N 12 tan θ W ) (N 13 tan α + N 14 ) .
Here we have the neutralino mixing matrix elements N i j defined in Ta Table 1. A marker is styled according to the analysis with the lowest CLs value, as described in Figure 1. This is projected into the planes of LSP mass m(χ 0 1 ) against relic density ( Figure 7a) and spin-independent LSPnucleon cross-section σ SI N-χ 0 1 (Figure 7b). The grey line overlayed on Figure 7a indicates the observed relic density Ω Planck CDM h 2 measured by Planck [101]. In Figure 7b, the cross-section is scaled by the ratio R = Ω χ 0 1 h 2 /Ω Planck CDM h 2 of the neutralino relic density Ω χ 0 1 h 2 to that observed Ω Planck CDM h 2 . Overlayed on Figure 7b are upper limits at 90% confidence level on the crosssection observed by the LUX 2016 [78] (solid light blue) results, together with the projected limit from 1000 days of data taking by the LZ experiment [107] (dashed green). The central value of the limit from LUX 2016 is twice as strong as that recently reported by PandaX-II [77].   Table 2. 300 GeV, where slepton and gaugino coannihilators are prevalent. The extensive presence of blue dots and cyan rings away from these bands shows respectively that the 2-6 jets and Monojet searches are particularly sensitive to gluinos and squarks that have small mass splitting with the bino-like LSP. Crucially, these impact scenarios where such coloured sparticles are the coannihilators (compare with Figure 9) and therefore indirectly probe LSP masses higher than those currently accessible by direct electroweakino searches at the LHC. Nevertheless, the 2-6 jets analysis also has sensitivity covering Higgsino and wino-like LSP points, with clusters along the Higgsino and wino bands.  The strongly interacting coannihilators of points with a bino-like LSP surviving Run 1 constraints, projected into relic density (left) and spin independent LSP-nucleon crosssection (right) vs LSP mass m(χ 0 1 ). The points are coloured according to dominant coannihilator for the LSP: gluino (light orange), light flavour squarks (dark blue), 3rd generation squarks (green) and 'Other' refers to (co)annihilation mechanisms involving the electroweak sector.

Direct detection
The neutralino LSP contributes to the local dark matter density and can induce nuclei recoils in direct detection experiments. Figure 7b illustrates the points excluded by the most sensitive of the six ATLAS searches, allowing us to compare with direct detection sensitivity for each analysis. This projects into the plane defined by the spin-independent 5 cross-section σ SI N-χ 0 1 of the LSP interacting with a Xenon nucleus, normalised per nucleon, and LSP mass. Again we discuss several features, with reference to other projections where appropriate, to aid interpretation: -Direct detection experiments typically interpret results assuming the LSP fully saturates the cold dark matter (CDM) relic abundance measured by Planck. As the LSP in the pMSSM need not be the sole constituent of dark matter, we rescale the direct detection interaction crosssections σ SI This accounts for the reduction in direct detection sensitivity due to a lower local density of neutralino LSPs. -Recently, the PandaX-II [77] and LUX [78] collaborations presented results that extend sensitivity for WIMP masses m χ 0 1 20 GeV by a factor of 2 and 4 respectively beyond the LUX 2015 result [109]. As a guide to the sensitivity of direct detection, the observed limits from LUX 2016 [78] and the projected sensitivity of LZ [107] are overlayed in Figure 7b.
-In defining the permitted pMSSM points, ATLAS conservatively increased by a factor of four the upper limit on the spin-independent cross-section R·σ SI N-χ 0 1 from LUX 2013 [98] to account for uncertainties in nuclear form factors [53], before rejecting points in the pre-selection. This pre-selection constraint carves out the points at the highest R · σ SI N-χ 0 1 cross-sections around 10 −44 cm 2 . - Figure 8b reveals points that survived constraints from the six 13 TeV searches considered in this study, in the plane of spin-independent cross-section against LSP mass. The points are coloured according to the whether the dominant contribution to the LSP is the bino, Higgsino, or wino. The composition of the LSP has a significant effect on R · σ SI N-χ 0 1 . This is due both to the couplings of the LSP to nucleons and the relic density suppression (Eq. 5 and Figure 8a). Notably, most wino-like LSPs have suppressed direct detection cross-sections given the small coupling to the Higgs boson.
Returning to the discussion of the regions of sensitivity from the ATLAS searches (Figure 7b), points above the blue solid line are independently excluded by both LUX 2016 and the 13 TeV ATLAS searches considered. In this context, collider searches are free from astrophysical uncertainties, and therefore provide powerful cross-checks should either approach report tentative signals. On the other hand, the cosmological lifetime of theχ 0 1 dark matter candidate can only be verified by non-collider means. For LSP masses m χ 0 1 1 TeV, the number of models excluded by the AT-LAS searches decreases rapidly. There remain many points surviving in this regime (Figure 8b), and those with the high R · σ SI N-χ 0 1 are within direct detection sensitivity.
The 2-6 jets analysis (blue points) is sensitive to a large class of models, particularly those with gluino or squark coannihilators. As this difference has important phenomenological consequences, we display the strong sector coannihilators in Figure 9b for bino-like LSPs of points that survived Run 1 constraints. In this projection, it is evident that coannihilation points involving light flavour squarks have enhanced cross-section R · σ SI N-χ 0 1 compared with gluinos. This is due to the s-channel diagram involving the quarks and LSP scattering via an intermediate squark, a point we will elaborate further in Section 3.2. Thus, the 2-6 jets and Monojet analyses share sensitivity to many squark coannihilator scenarios with LUX 2016.
As discussed in previous sections, the 7-10 jets is most sensitive to light mass LSPs with significant bino content. Many of the points excluded by this search are below the current LUX 2016 sensitivity. Meanwhile, the Multi-b analysis tends to favour scenarios Higgsino-like LSP scenarios where squark masses are above 2 TeV. There is particular sensitivity to a region centred around cross-section from 10 −46 cm 2 to 10 −49 cm 2 and LSP mass of 300 to 700 GeV (Figure 7b). Many Higgsino-like LSP models inhabit this region and are beginning to be probed by LUX 2016, but the majority of points where the Multi-b analysis is most sensitive are below the LUX limit.
Notably, ATLAS strong SUSY searches are sensitive to scenarios with direct detection cross-section R · σ SI N-χ 0 1 well below even the projected sensitivity of LZ based on 1000 days of data taking [107]. The SS/3L reach into this regime is especially prominent, being most sensitive to wino-like LSP models with highly suppressed cross-sections R · σ SI N-χ 0 1 . Many gluino coannihilator scenarios occupy this region (Figure 9b) and are dominantly probed by the 2-6 jets analysis. The projected LZ sensitivity is within an order of magnitude of the irreducible neutrino background 'floor', which is a challenging regime for Xenon-target direct detection experiments.
Concluding this subsection, we demonstrated the important complementarity of strong SUSY searches for probing models beyond both Monojet-like collider interpretations and the reach of direct detection experiments. This enables colliders to indirectly constrain bino-like LSPs with coloured coannihilators, in addition to Higgsino LSP scenarios for example before electroweak SUSY searches gain direct sensitivity. This motivates construction of simplified DM models based on such interpretations, but is beyond the scope of this work.

Impact of direct detection constraints on squarks and gluinos
Finally, turning the question around, we explore the impact of LUX 2016 constraint on the parameter space of squarks and gluinos relevant to LHC searches. Considering the set of points that survived ATLAS Run 1 constraints [53], we deem any point with scaled cross-section R · σ SI N-χ 0 1 above the 90% confidence level upper limit observed by LUX 2016 to be excluded. We take the observed central value of the upper limit as reported by LUX 2016 [78], without rescaling to account for the nuclear form factor uncertainties 6 (as was done in Ref. [53]). The upper limit derived by LUX 2016 excludes 30.3% of such pMSSM points constraints (with long-lived gluinos, squarks and sleptons removed). Figure 8b shows that LUX is particularly sensitive to pMSSM points with bino-like LSP, due in part to the relic density suppression of Higgsino and wino-like models. Figure 10a projects the fractional exclusion 7 by LUX 2016 into the gluino vs LSP plane. There is a band of modest exclusion for LSP masses m(χ 0 1 ) 250 GeV, which is relatively uncorrelated with gluino masses. This apparent enhancement of sensitivity is partly an artefact due to the systematic oversampling of coannihilators for bino-like LSPs in Ref. [53]. This region is rich in gaugino and slepton coannihilators due to their weak LHC constraints compared with strongly interacting sparticles. Due to smaller annihilation cross-sections of this electroweak process in the early universe, the relic abundance Ωχ0 1 h 2 for a given bino-like LSP mass is larger. These models are thus not scaled down as far by the R = Ωχ0 1 h 2 /Ω Planck CDM h 2 factor. A prominent feature of the points excluded by LUX appears when projecting into the squark vs LSP plane (Figure 10b). The diagonal region where the mass splitting between the LSP and squark is small m(q) − m(χ 0 1 ) 50 GeV shows a distinctly higher exclusion fraction; this was absent in the gluino vs LSP plane. Furthermore, mass bins with higher exclusion fraction are correlated with lower squark masses. This is particular salient along the diagonal, where points are predominantly squark coannihilators. The band of modest exclusion fraction for LSPs with m(χ 0 1 ) 250 GeV remains, as with the gluino vs LSP plane.
Such contrasting features between gluinos and squarks highlight the importance of the s-channel squark exchange diagram in direct detectioñ 6 Had we weakened this constraint by a factor of four [53], the key qualitative features in the discussion are unaffected; only the numerical fraction of models excluded is reduced to 18.0%. 7 Although N Excluded LUX 2016 /N Survived ATLAS Run 1 is prior-dependent, we will momentarily discuss this, while the forthcoming comparison of gluino and squark coannihilation regions has minimal prior dependence.
Here, a squarkq mediates the LSPχ 0 1 scattering off light flavour quarks q inside the nucleons. Resonant scattering occurs when the LSP and squarks are nearly mass degenerate, enhancing the sensitivity of LUX along the squark-LSP diagonal of Figure 10b These observations are consistent with the LSP-nucleon cross-section vs LSP mass plane with coannihilation mechanism identified (Figure 9b). Squarks of the first or second generation (dark blue) have particularly enhanced direct detection cross-sections. The points with gluino coannihilators (light orange) feature R · σ SI N-χ 0 1 primarily below 10 −46 cm 2 and will only begin to be probed at the direct detection frontier by future experiments such as LZ. We also note that coannihilation points involving third generation squarks (light blue) also tend to have suppressed direct detection cross-sections compared with squarks, due to the negligible 3rd generation content in nucleons.
Taken together, this highlights the important implications of LHC searches for squarks in the context of direct detection experiments, especially when the squark-LSP mass splitting is small. Yet these squark coannhilation scenarios are challenging for colliders, where direct detection experiments can provide a complementary probe.

Conclusion
In this work, we interpreted six published 13 TeV (3.2 fb −1 ) ATLAS SUSY searches for gluinos and light flavour squarks in a 19-parameter pMSSM. The purpose of this study was to analyse previously unexamined correlations between the most sensitive analyses with distinct regions of pMSSM parameter space. Our study addressed various shortcomings in the literature, presented under three questions in the Introduction (Section 1), which we now summarise.
Firstly, we examined these correlations in collider parameter spaces, providing substantially richer information than overlap matrices used in the literature. For the twodimensional projections into gluino, LSP and squark masses, the separation in regions probed by the 2-6 jets and Multib analyses were particularly distinct. The Multi-b was the most sensitive analysis for models with larger gluino-LSP mass splittings, where the 2-6 jets search began to lose sensitivity. The regions identified are independent of the priors in the pMSSM points. Further, while the Monojet and 2-6 jets share substantial overlapping sensitivity, the tighter jet requirements of the former is needed for the scenarios where the coloured sparticle and LSP are near mass-degenerate.  Secondly, we identified classes of models beyond those ATLAS used for optimisation. Arguably the most striking realisation of this was the SS/3L search. Despite optimising to four distinct simplified models, we found one signal region to be most sensitive to a different scenario in the pMSSM not considered by the ATLAS search. It involved a light flavour squark cascading to a slepton and wino-like LSP with a nearly mass-degenerate chargino, which could be used by the experimental collaborations to refine future searches. Meanwhile, though the 1-lepton search optimised for a single simplified model, we showed it was sensitive to scenarios that included squark production, as well as intermediate sleptons and gauginos.
Finally, while ATLAS performed an explicit DM interpretation for their Monojet search, our study manifested the prominent role other searches for coloured sparticles have when the dark sector is beyond non-minimal regimes as in the pMSSM. Bino-like LSPs may rely on coloured coannihilators to be consistent with the observed relic abundance, which are probed by the 2-6 jets and Monojet analyses being sensitive to small squark-LSP and gluino-LSP mass splittings. Light flavour squarks enhance LSP-nucleon scattering cross-sections, and squark lower mass bounds can still be below 500 GeV in the pMSSM. In addition, the SS/3L analysis had preferential sensitivity to wino-like LSP scenarios, which are particularly challenging for direct detection experiments.
Using our findings to design novel search strategies or interpretations was beyond the scope of this work and is deferred to future studies. It would also be of interest to perform similar assessments for the third generation squark and electroweak sectors including long-lived sparticles, interpret searches based on simplified models of DM, as well as develop surveys for non-minimal SUSY scenarios. Previously unprobed regions of the pMSSM will be explored further as LHC luminosity continues to rise.  Table 2: Definition of neutralinoχ 0 1 LSP categories from Ref. [53]. In the neutralino mixing parameter N i j , the first index denotes the neutralino mass eigenstateχ 0 i while the second indicates its dominant composition in the order B ,W ,H 1 ,H 2 .
Appendix A.1: Review of the ATLAS pMSSM19 Starting with the MSSM, minimal flavour violation was imposed and CP violation was restricted to the CKM phase in the quark sector. The neutralino was required to be the lightest supersymmetric particle (LSP) and R-parity is exactly conserved. Table 2 displays the three categories of LSP type by their dominant composition, as defined in Ref. [53]. Due to their different resulting phenomenology, ATLAS employed importance sampling [53] to ensure approximately equal proportions of each LSP type are selected. The 19 resulting parameters were scanned with flat priors, taking an upper limit on mass scales to be 4 TeV. The resulting model points were subjected to the following non-LHC constraints. LEP lower mass bounds [79] were imposed together with precision measurements from the electroweak isospin splitting parameter ∆ ρ [81], g − 2 of the muon [82][83][84][85][86][87][88][89], Z 0 invisible width [80], branching fractions of heavy flavour states [90][91][92][93][94][95][96][97] and Higgs boson mass [111] at the time of pMSSM points generation were applied [53]. Dark matter constraints were subsequently considered. To account for uncertainties in nuclear form factors, ATLAS quadrupled the upper limit on the spin-independent crosssection for LSP-nucleon interaction from LUX 2013 [98] before applying this to the points. Similar constraints for the spin-dependent cross-section of LSP-proton interactions were applied from COUPP [99], and XENON100 [100] for the LSP-neutron cross-section. The LSP was not assumed to be the sole constituent of dark matter so only an upper limit was set for the LSP relic abundance Ωχ0 1 h 2 , taken as Ω CDM h 2 = 0.1208, which is the central value plus twice the reported uncertainty from the Planck measurement [101].
The 310.3k pMSSM points surviving all non-LHC constraints underwent evaluation against 22 relevant Run 1 AT-LAS searches for supersymmetry using 7 and/or 8 TeV data using full ATLAS software and reconstruction; for full details, see Ref. [53]. A total of 40.9% of points were excluded at 95% confidence level. The top row of Table 3 Table 3: Viable model points before Run 2 constraints. These are classified by the dominant contribution to the LSP being bino, wino or Higgsino. Long-lived (LL) gluinos, squarks and sleptons with cτ > 1 mm require dedicated Monte-Carlo simulation and are omitted from this study. Event simulation was performed on non-LL models with total strong sparticle production cross-section σ tot ≥ 5 fb.
Appendix A.2: Interpretation of early 13 TeV searches We now summarise the methodology of Ref. [55] to apply the combined constraints of six published 13 TeV AT-LAS searches considered to the 183.8k pMSSM points that survived Run 1. Long-lived squarks, gluinos and sleptons (cτ > 1 mm as defined in Ref. [53]) make up 1.9k of these points. These require dedicated Monte-Carlo simulation and are removed from consideration of 13 TeV sensitivity. Of the remaining 181.8k points, 71.4% have total inclusive production cross-section of any two coloured sparticles σ tot greater than 5 fb at leading order. This is enumerated in the bottom row of Table 3 and underwent particle event generation to be described. The remaining 28.6% was deemed not to have sensitivity with 3.2 fb −1 of integrated luminosity and therefore not excluded.
We used MADGRAPH5 2.3.3 [112,113] for particle event generation involving any two coloured sparticles from two protons, interfaced with PYTHIA 6.428 [114] for hadronisation and showering, using the CTEQ6L parton distribution functions [115]. Up to one additional parton in the matrix element and the MLM prescription [116] was used to match jets, setting the MADGRAPH minimum parton k T parameter to 100 GeV and PYTHIA jet measure cutoff at 120 GeV, in accord with Ref. [53]. The DELPHES 3.3.2 [117] fast detector simulator was employed with FASTJET 3.1.3 [118,119], using the anti-k T clustering algorithm with cone parameter R = 0.4, to parametrise the performance of the ATLAS detector. Table 1 lists the six Run 2 analyses, whose event selection were implemented in the MADANALYSIS5 1.3 [120,121] recasting framework. We adapted codes from the Public Analysis Database [122][123][124] where available, and write our own otherwise. The RECASTINGTOOLS package was used for limit setting using the CLs prescription [102], where a point is deemed excluded if the CLs value is below 0.05. To validate our code, we ensured at least one benchmark point had cutflows agreeing to better than 30% and the observed simplified model limits were reproduced to within the uncertainties published by ATLAS (see Appendix B for further details). Appendix A.3: Summary of sparticle masses excluded by early 13 TeV searches Overall, out of 181.8k pMSSM points considered for 13 TeV (3.2 fb −1 ) sensitivity, 15.7% were excluded by the combined constraints. Figure 11  For gluinosg, the high fractions (above 60%) of models excluded involve masses below 1 TeV, and to a lesser extent around 1.4 TeV, is unambiguous. Reduced sensitivity around 1.1 TeV is due primarily to mass splittings between the gluino and LSP being less than 200 GeV (Figure 12). This corresponds to a localised region of high exclusion around 1.4 TeV gluino and 1.8 TeV squark masses the gluino-squark plane ( Figure 11). Overall, good corroboration with the simplified models is observed.
The lightest squarkq of the first two generations shows a smaller but noticeable extension of sensitivity beyond Run 1 with 3.2 fb −1 of integrated luminosity. This is due to the less advantageous scaling of cross-sections between 8 and 13 TeV compared with gluino production from LHC parton distribution functions. In the squark-LSP plane (Figure 11), the simplified model limit assumes all eight squarks are of the first two generations are mass-degenerate and therefore over-constrains squark scenarios in the pMSSM. Figure 12 separates the masses of the squarks by the chirality of the partner quark. As the left-handed squarksd L ,ũ L form an SU(2) doublet, their mass splittings are typically negligible compared to those of the right-handed squarks d R ,ũ R . In accord with Ref. [53], right-handed down-type squarks are least constrained of the four light flavour squarks. Analyses targeting third generation squarks were not considered and so there is little impact on the lightest sbottomb 1 and stopt 1 masses. Also shown in Figure 12 are uncoloured sparticles, which can participate as intermediate states in cascade decays of the gluino or squark(s). However, it is seen that the sensitivity of the analyses considered in Table 1 have little to no correlation with these sleptons or electroweakinos. The lower mass bounds on charged electroweak sparticles are due to LEP searches. Here, slepton˜ refers to the lightest superpartner of the left-or right-handed charged lepton of the first or second generation. Table 4 shows the overlap matrix between the six searches we considered. This quantifies the fraction of points excluded by one analysis are also excluded by another. For example, of the points excluded by the 7-10 jets search, 76% of them were also excluded by the 2-6 jets analysis. Overall, the sensitivity is complementary, with no analysis completely excluding the points probed by another search.
Appendix A.4: Distributions prior to 13 TeV interpretation When interpreting the results of this work, especially Figures 11, it is instructive to examine the distribution of the 181.8k points that survived Run 1 constraints, before any 13 TeV constraints are applied. Figure 13 displays these prior distributions projected into the three mass planes involving the gluino, squarks and LSP. Along the diagonal where the gluino or squark are similar in mass to the LSP, there is a larger abundance of models. These are primarily models with bino-like LSP, which require the LSP to coannihilate with a near-degenerate gluino to prevent oversaturation of the observed relic abundance. Their density were enhanced from the importance sampling used by ATLAS in the pre-selection of the pMSSM points [53].
Furthermore, the region below 100 GeV with substantially fewer models are also mainly bino-like LSP models. These correspond to where the LSP relic density is reduced by resonance annihilation through a Z 0 or h 0 boson through the so-called 'funnel' mechanism. Models with Higgsinoor Wino-like LSPs typically have near-degenerate charginos which are constrained by LEP bounds, hence reducing the concentration of models in this region. These projections also show very few (less than 10) models in the mass bins coloured green close to the region where 100% of models were excluded in Run 1.

Appendix B: Details of analysis validation
This appendix details the validation of the analyses we implemented in this work. The purpose of the validation was the check our implementation of event selection in MAD-ANALYSIS. We display the cutflows (signal acceptance after each event selection) for our implementation compared with those published in the supplementary material of each ATLAS analysis. We also present a simplified model limit comparison for each of the six analysis.
To generate signal samples for validation, we used the same simulation setup as all the pMSSM points simulated in this study (outlined in Subection Appendix A.2). The main differences here are that we generate signals with up to 2 extra partons in the matrix element and the jet matching scale was set to one quarter of the produced sparticle mass, in accord with ATLAS. The event yields were normalised to the NLO squark or gluino cross sections from the LHC SUSY Cross Section Working Group [125].   Table 1, out of the points that survived Run 1 constraints. In each mass bin, the colour scale denotes the fraction of points excluded at 95% confidence level normalised to the number of points satisfying 'Without long-lived' in Table 3, where black indicates 100% exclusion. This is projected into the mass planes of the gluino-LSPg-χ 0 1 (upper left) and squark-LSPq-χ 0 1 (upper right). Here, m(q) is the mass of the lightest squark among the first two generations. White regions correspond to no models being produced by ATLAS due to prior non-LHC constraints. Hatched grey regions indicate bins where all points were excluded by Run 1 in Ref. [53]. Overlayed grey solid lines are the simplified model limits from the 13 TeV 2-6 jets search [1] for gluinosg → qqχ 0 1 (upper left) and squarksq → qχ 0 1 (upper right). In the latter case, all eight squarks are of the first two generations and assumed to be mass-degenerate. For the gluino-squark plane (lower), the overlayed grey dashed line is taken from the 'gluino-squark-LSP simplified pMSSM' scenario from the 8 TeV 2-6 jets search [15].    Table 1 for different sparticle masses. The exclusion information is from Ref. [55]. The colour scale indicates the fraction of models excluded, normalised to the number satisfying 'Without long-lived' in Table 3. Black regions denote that all points were excluded in the mass bins of 50 GeV size. Here,q (l) is the lightest squark (slepton) of either the first or second generations. The hatched grey regions indicate all models in the mass bins were excluded by Run 1 ATLAS searches from Ref. [53]. Models with long-lived gluinos, squarks or sleptons are removed from these figures.
In summary, the agreement is sufficient for the scope of this study.

Appendix B.1: Cutflow comparison
For the cutflow comparison, we generated signal samples at a benchmark point presented by ATLAS in the public supplementary material of each analysis. The comparison ensures our implementations of each search in MADANALY-SIS are selecting the correct kinematic spaces in all signal regions. Tables 5 to 10 compare the cumulative percentage signal acceptances found by the collaboration 'ATLAS' after each selection criterion 'Cut' with ours 'MA5 Recast'. In these tables, the model nomenclature and variables labelling each cut follows those of the ATLAS supplementary material for ease of comparison; see the references in Table 1 for full details.
We achieved better than 30% agreement in the final yields for the signal regions 'SR' in each analysis. This is reasonable given the systematic uncertainties introduced by using a different version of PYTHIA for showering compared with ATLAS, together with application of the same object isolation, overlap removal, reconstruction and b-tagging efficiencies in the DELPHES detector simulation for all six analyses (to reduce disk usage to a manageable level). Appendix B.2: Simplified model limits Figure 14 shows comparisons of exclusion limits published by ATLAS (solid lines) with our own excluded points (orange squares) using one of the simplified model scenarios considered in each analysis. We varied the appropriate masses in each model SLHA card and generated signal samples for each point. The RECASTINGTOOLS routine in MADANAL-YSIS was used to set 95% confidence level limits.
Small discrepancies arose between our MADANALYSIS implementation and those of ATLAS due to various sources of systematic uncertainty in our simulation setup, such as different PYTHIA versions and fast detector simulation setup. Nonetheless, the general corroboration in the shape of our MADANALYSIS exclusions across the simplified model parameter spaces, together with overall cutflow agreement, gives confidence that the event selections were implemented correctly.