Supersymmetry Without Prejudice at the 7 TeV LHC
SLACPUB14382 ANLHEPPR1113 NUHEPTH/11014 BONNTH201105
Supersymmetry Without Prejudice at the 7 TeV LHC
John A. Conley
Physikalisches Institut, Universität Bonn, Nußallee 12, 53113 Bonn, Germany
James S. Gainer
High Energy Physics Division, Argonne National Laboratory, Argonne, IL 60439 USA and
Department of Physics and Astronomy, Northwestern University, Evanston, IL 60208 USA
JoAnne L. Hewett, My Phuong Le, Thomas G. Rizzo
SLAC National Accelerator Laboratory, 2575 Sand Hill Rd., Menlo Park, CA 94025, USA
We investigate the model independent nature of the Supersymmetry search strategies at the 7 TeV LHC. To this end, we study the missingtransverseenergybased searches developed by the ATLAS Collaboration that were essentially designed for mSUGRA. We simulate the signals for k models in the 19dimensional parameter space of the pMSSM. These models have been found to satisfy existing experimental and theoretical constraints and provide insight into general features of the MSSM without reference to a particular SUSY breaking scenario or any other assumptions at the GUT scale. Using backgrounds generated by ATLAS, we find that imprecise knowledge of these estimated backgrounds is a limiting factor in the potential discovery of these models and that some channels become systematicslimited at larger luminosities. As this systematic error is varied between 20100%, roughly half to 90% of this model sample is observable with significance for 1 fb of integrated luminosity. We then examine the model characteristics for the cases which cannot be discovered and find several contributing factors. We find that a blanket statement that squarks and gluinos are excluded with masses below a specific value cannot be made. We next explore possible modifications to the kinematic cuts in these analyses that may improve the pMSSM model coverage. Lastly, we examine the implications of a null search at the 7 TeV LHC in terms of the degree of finetuning that would be present in this model set and for sparticle production at the 500 GeV and 1 TeV Linear Collider.
1 Introduction and Background
The LHC has had an initial run at 7 TeV with both the ATLAS and CMS experiments collecting pb of useful data. Even with this low integrated luminosity these experiments have been able to extend searches far beyond the reach of the Tevatron for many new physics scenarios with, so far, null results [1, 2]. This clearly demonstrates the power of increasing the center of mass energy in the search for new physics at hadron colliders. Starting soon, the LHC is to begin a longer run at 7 TeV and is expected to collect of order fb of data over the next 2 years. Such a data set will allow for a first exploration of the TeV mass scale, and if new strongly interacting particles exist in this kinematic regime, they should be observed.
A wellmotivated, and perhaps most popular, possibility for new physics that may be discovered during this coming LHC run is Supersymmetry (SUSY) [3]. Both ATLAS [4] and CMS [5] have designed detailed searches for many of the SUSY partners of the Standard Model (SM) particles; these are mostly (but not exclusively) based on the assumption of mSUGRA/CMSSMlike soft breaking within the Minimal Supersymmetric Standard Model (MSSM) framework. This assumption greatly simplifies the exploration of the vast Supersymmetric parameter space. While these searches are designed to well cover the parameter space of these SUSYbreaking scenarios it is important to ascertain their discovery potential in a more general MSSM context. This is particularly worrisome in light of results from the Tevatron, where it has been realized [6, 7] that relatively light gluinos may have escaped undetected. The question then arises whether these specific scenarios and associated searches adequately describe the true breadth of the MSSM and its possible collider signatures, and whether the LHC searches as presently designed could fail to observe sparticle production. This has prompted several studies of more model independent search strategies as well as the development of simplified models [8].
Recently, we have addressed [9] this question by investigating the capability of the 14 TeV LHC to explore a more general MSSM model parameter space, i.e., that of the pMSSM (phenomenological MSSM) [7], to be described below, from the point of view of the ATLAS detector. In particular, we examined the performance of the planned ATLAS SUSY searches in exploring this more general MSSM scenario. This analysis provides insight into general features of the MSSM without reference to a particular SUSY breaking scenario or any other assumptions at the GUT scale. We found that the ATLAS mSUGRAinspired searches, based on missing , did surprisingly well at covering the kinematically accessible portions of this model space; we also found that some interesting exceptions can arise in these more general models. Given the lowerenergy run of the LHC over the next 2 years it behooves us to determine how well the corresponding mSUGRA motivated searches designed by ATLAS would perform under these conditions, since this is the situation presently before us. This is the goal of the present paper. We note that recently there have been several mSUGRAbased studies evaluating the capability of the 7 TeV LHC run to probe that parameter space [10]
As is well known, soft SUSY breaking within the MSSM in all generality leads to a theory with over a priori free parameters
which prohibits a detailed study of this theory. A number of theoretically possible scenarios exist which describe the breaking of Supersymmetry; maybe
even multiple mechanisms are simultaneously responsible. Practically speaking, there are two ways to approach reducing this large number of a priori unknown
parameters to something more manageable. One approach is to consider only specific, wellmotivated SUSY breaking scenarios, such as mSUGRA or others.
This leads to a drastic reduction in the number of free parameters (to only ) so that detailed analyses of the
resulting parameter space can be easily achieved. A problem with performing such studies is that they may result in a bias as to the nature
of SUSY signals when searching for collider or other SUSY signatures. An alternative approach is to be less prejudicial and to instead follow a bottomup
analysis which we have employed in a number of recent works [9, 11]
and will make use of here. By imposing a set of theoretically
and experimentally wellmotivated constraints on the general MSSM (to be described below), without making any reference to the specific mechanism of SUSY breaking,
we arrive at a theory with 19 TeVscale parameters. This is known as the pMSSM, which is
significantly more
manageable than the full Supersymmetric parameter space, and yet allows for more breadth than is present in , e.g., mSUGRA.
To this end, we examine the based SUSY searches developed by the ATLAS collaboration for the 7 TeV LHC [12]. We simulate the pMSSM signal for roughly 71k pMSSM models (hereafter ‘model’ refers to a point in the 19dimensional pMSSM parameter space) that we generated in our previous work [7]. We employ SM backgrounds provided by the ATLAS Collaboration. In Section 2, we describe our SUSY model generation and LHC analysis procedure. It is important to note that we strictly adhere to the analyses as designed by ATLAS. While numerous, and perhaps improved, SUSY collider search techniques have been discussed in the literature [13], it is not our present purpose here to discuss or employ them. Section 3 contains our main results. We find the systematic error in determining the SM background to SUSY production is a limiting factor in the potential discovery of these models; in fact, some channels become systematics limited at larger luminosities. In this section we determine the fraction of our pMSSM model set that is discoverable at the 7 TeV LHC. We then examine the model characteristics in some detail that render some of the models undetectable. We find that the observability of models depends on the precise details of the sparticle spectrum and that a blanket statement of constraints on the mass of, say the gluino or squarks, cannot be made. In this Section, we also explore potential modifications to the kinematic cuts in these analyses that may improve model observability. In Section 4, we examine the implications of a null SUSY search during this run with respect to the degree of finetuning present in these models, as well as the implications for sparticle production at a high energy Linear Collider. A summary and our conclusions can be found in Section 5.
2 Analysis Procedure for Inclusive SUSY Production at the LHC
The purpose of this work is to explore how well mSUGRAinspired inclusive SUSY searches (in particular the set proposed by the ATLAS collaboration [12]) apply to the larger and much more general pMSSM parameter space for the 7 TeV run of the LHC. This is similar in spirit to to Ref. [9], which explored this question for the more powerful 14 TeV LHC. The pMSSM model sample that we study was generated in Ref. [7]; here, we briefly review the procedure employed to generate this sample. We then describe our procedure for generating the signal events, comparing to background, and determining the statistical criteria for discovery. We will show that we faithfully reproduce the ATLAS results in each analysis channel for their benchmark SUSY model.
As stated above, we study the 19dimensional parameter space of the pMSSM. This set of parameters was arrived at [14] by imposing the following set of requirements onto the general RParity conserving MSSM: () the soft parameters are taken to be real so that there are no new CPviolating sources beyond those in the usual CKM matrix; () Minimal Flavor Violation(MFV) [15] is assumed to be valid at the TeV scale; () the first two generations of sfermions having the same quantum numbers are taken to be degenerate and to have negligible Yukawa couplings and () the lightest neutralino is assumed to be the Lightest Supersymmetric Particle(LSP) and is a stable thermal WIMP. Most of these assumptions are applied in order to avoid issues associated with flavor physics constraints. With these conditions, the remaining 19 free softbreaking parameters are given by the three gaugino masses, , the ten sfermion masses , the three terms for the third generation fermions (), and the usual Higgs sector parameters , and .
To generate the specific pMSSM parameter points that we study below (hereafter referred to as our set of models), we performed numerical scans over the space formed by these 19 parameters. This required both a selection of the parameter range intervals as well as an assumption about the nature of the scan prior for how points are chosen within these intervals. These issues are both described in detail in our previous works [7, 9, 11] and the interested reader should refer to them directly. Here, we simply note that two scans were performed: one employing a flat prior beginning with points and one with a log prior employing points. The main distinctions between these two scans directly relevant to our analysis here are that () all SUSY mass parameters were restricted to be TeV for the flat prior case, while for the log case the upper limit on mass parameters was raised to 3 TeV, and (b) the choice of the log prior generally leads to a more compressed sparticle spectrum than does the flat prior case. Note that the restriction on the upper limit for the mass parameters ensures relatively large production cross sections at the LHC for the case of the flat prior model sample.
Once these points were generated, we demanded that they be consistent with a large
number of both theoretical and experimental constraints in order to ensure that the model
sets are valid to study. We mention the most important of these restrictions
here
We now turn our attention to the analysis procedure that we followed in generating and analyzing the signal events from sparticle production at the 7 TeV LHC. Throughout our analysis, we adhere to the search strategies developed by ATLAS [12] as closely as possible. In this reference, ATLAS considers 10 MET search channels, including selections where the minimum number of jets is 2, 3, or 4 and the number of leptons is 0, 1, or 2. In the dilepton case, oppositesign (OSDL) and samesign (SSDL) pairs are considered separately, with SSDL only being considered in association with 2 jets. Note that flavor tagged final states are not considered here. We consider 85 SUSY production processes that contribute to these 10 signatures.
Accurate estimates of the SM backgrounds for the various channels are crucial to the validity of this study. We obtained details of the background distributions presented in Ref. [12] directly from the ATLAS SUSY Group [18]. These backgrounds were produced with stateoftheart Monte Carlo event generators and the full ATLAS detector simulation. Employing these ATLAS computed backgrounds in our analysis allows us to concentrate on generating and analyzing signal events for each of the 71k parameter space points in our pMSSM model sample.
2.1 Generation of the Signal Events
The steps involved in the generation of the signal events are very similar to those detailed in Ref. [9]. Here we will briefly summarize the procedure and point out any differences in the present analysis. For the generation and analysis of events for a single model, the workflow is:

Using PYTHIA 6.418, events were generated, fragmented, showered, and hadronized for each of 85 SUSY production processes, with each process being weighted by its Kfactor.

Detector effects were simulated using an ATLAStuned version of the fast detector simulation PGS4 [22] with the default isolation cuts removed.

The simulated events were then analyzed using the analysis cuts for the 10 ATLAS analyses listed above, as well as the isolation cuts described in Ref. [12].
We note that as in our previous work, a subset (about 1%) of the models suffered serious enough errors that the Pythia event generation halted. These “PYSTOP” models are excluded from our results.
2.2 Analysis Cuts
For the reader’s convenience, we provide here the full set of kinematic cuts for each analysis channel, summarizing the information given in Section 4 of Ref. [12]. All channels have a missing energy cut of GeV, and all analyses except the SSDL channel have a transverse sphericity cut of . Table 1 summarizes the cuts for the searches with jets which are independent of the choice of lepton channel except for SSDL, which will be described below.
Number of jets  

Leading jet (GeV)  
Other jets (GeV)  
min.  
We complete the description of the kinematic cuts by specifying the additional cuts that are specific to the various lepton channels. For the case with zero leptons, events are rejected that have at least one lepton with GeV. For the onelepton channels, one lepton with GeV is required, no additional leptons with GeV are allowed, and the transverse mass of the selected lepton and the missing energy vector must satisfy GeV. (The definition of can be found in Ref. [4].) For the OSDL channel, exactly two leptons with GeV are required, and they must have opposite charge.
As mentioned above, the SSDL channel has distinct cuts; unlike all the other analyses, the three different jet selection options specified in Table 1 are not employed. Instead, two jets with GeV are required. In addition, two leptons with GeV, same charge, and invariant mass GeV must be present, and there is a veto on additional leptons with GeV. The transverse mass of the leading lepton and the missing energy vector must satisfy GeV.
Lastly, when performing our statistical analysis, the cut is optimized for each channel, and for each pMSSM model, in steps of 400 GeV [12].
2.3 Statistical Procedure
To compute the significance of the signal for each search channel, we follow the statistical procedure described in detail in Refs. [9, 4], which is that employed by the ATLAS collaboration. The probability that the expected background fluctuates to the number of observed events is computed assuming that the systematic error on the background is Gaussian and the statistical error is Poissonian. A significance of is required for the observation of a signal. As mentioned above, the cut is optimized for each channel, and for each pMSSM model, in steps of 400 GeV [12]. As was discussed in Ref. [9], and will be further demonstrated here, the accuracy of the background estimation has a profound impact on the signal significance and the resulting search reach. In order to quantify this, we will present results assuming a 20, 50, and 100% systematic error on the background. We will also consider integrated luminosities of 0.1, 1, and 10 fb.
2.4 Comparison with ATLAS Benchmark Models
The ATLAS SUSY group has published signal rates for a single mSUGRA benchmark point (SU4) in their study of Supersymmetry at the 7 TeV LHC [12]. It is imperative for us to check the results of our analysis against these published results for this benchmark model before proceeding to apply our analysis to the pMSSM model set.
Here it is important to remind the reader that our SUSY signal generation, as described in detail in Ref. [9], differs slightly from the procedure employed by ATLAS. In particular, the numerical programs used to compute the SUSY spectrum and decay tables, as well as event generation, are different. Furthermore, we use, by necessity, a fast detector simulation as opposed to the ATLAS full GEANTbased simulation. Therefore, a small degree of discrepancy can be expected. The comparisons shown in Figures 15, however, indicate that we are indeed able to faithfully reproduce the results obtained by ATLAS for this benchmark model for all of the various inclusive analyses.
3 Results of the 7 TeV Analysis
In this Section, we relate the results of our study on the effectiveness of the ATLAS 7 TeV analyses in detecting our pMSSM model sample. We first discuss the impact of the size of the background systematic errors on SUSY searches, and then turn to the discovery coverage of the pMSSM. We examine the characteristics which cause some models to be undetectable as well as study the effects of modifying the ATLAS SUSY analysis cuts. We remind the reader that our sample of k models is not intended to be a full description of the 19dimensional pMSSM parameter space. However, the sample does contain numerous models which exhibit properties that are quite different than those expected in mSUGRA and thus provides insight into general features of the full MSSM.
3.1 Influence of Background Systematic Errors
As mentioned above, the size of the SM background systematic errors plays an important role in the ability of the ATLAS searches to discover Supersymmetry, including the pMSSM. This is not surprising as the number of signal events necessary to reach critically depends upon both the size of the estimated background itself as well as the background uncertainty. For a fixed systematic uncertainty, search channels with large backgrounds clearly require a large number of signal events in order to claim a discovery. To get a feel for this in the case of the ATLAS analyses studied here, we determine the necessary number of signal events to reach the level in each analysis as function of the fractional background uncertainty. We remind the reader that the SM backgrounds for each channel were supplied to us by the ATLAS SUSY working group [18]. In performing these calculations we exactly follow the discussion as given by ATLAS in Ref. [4]. Our results are displayed in Figs. 6 and 7 for the ten ATLAS channels assuming 1 fb of integrated luminosity. Here, we see the number of signal events that are required to obtain the discovery criterion of for various values of the final cut. In the case of the nj0l channel, which has the largest SM background, we note that the required number of signal events is particularly large and is quite sensitive to the value of the cut. Note that as the systematic error increases, the number of required signal events can rise drastically, in some cases by an order of magnitude or more. In particular, the difference between a reasonably low 20% systematic error and taking a 0% error (i.e., ignoring this effect) is substantial and theoretical analyses that do not include this error are thus wildly optimistic.
We will use these numerical results in our subsequent analyses of the pMSSM model coverage in these based searches in the next subsection. They indicate the importance of reducing background systematic errors in order to increase the coverage of new physics parameter spaces.
3.2 pMSSM Model Coverage
We now run each of our pMSSM models through the analysis chain described above. The first question we address is how well do the various search analyses cover the pMSSM model sample, or, more precisely, what fraction of these models can be discovered (or not) by these searches. Further, we also determine which of the analyses provide the best model discovery capabilities. Clearly the answers to these questions will be highly sensitive to the assumed values of both the integrated luminosity and the estimated SM background uncertainty. Figure 8 shows the fraction of the pMSSM models that can be discovered with in each of the ATLAS channels for the flat prior model set as a function of the integrated luminosity assuming three different choices for the background systematic error. The corresponding results obtained in the case of the log prior model set can be found in Fig. 9. Here, we again emphasize that our pMSSM sample is not meant to provide full coverage of the 19dimensional parameter space (such coverage would be computationally prohibited). However both the very large number of models in our pMSSM sample, and the distinct characteristics they possess, make this sample an ideal testbed for this set of mSUGRA designed search strategies. Our fractional results based on our pMSSM model set are thus indicative of the behavior of the MSSM under these search routines.
These figures reveal a number of interesting results: () The size of the background systematic errors makes a significant impact on model coverage for all search channels and integrated luminosities. For the searches with significant SM backgrounds, i.e., the nj0l and nj1l channels, variation in the background uncertainty leads to substantial changes in the capability to observe the model sample. The search least affected by systematics is 2jSSDL since the backgrounds in this case are quite small. This behavior confirms the results of the previous subsection. () The model coverage in almost all cases is significantly better for the flat prior model set than for the log prior sample. This, too, is not very surprising since the masses of the sparticles in the log prior case extend out to much larger values and the sparticle spectrum is generally more compressed in this set [7]. The latter leads to softer jets and leptons in the corresponding cascade decays which have a more difficult time passing the analysis cuts. () For all values of the background systematic error, the nj0l channels yield the best model space coverage with 4j0l affording the best discovery opportunity. In fact, we see that the channels which require more jets to be present have a better chance of being observed for the nj0l and nj1l searches. For the flat prior sample, the 4j0l analysis with low background systematics is observed to cover a very large fraction of the model set by itself once significant integrated luminosities are obtained. () As the number of leptons required to be present in the final state increases, the model coverage is found to decrease significantly, especially for smaller values of the integrated luminosity. This is due to the fact that the branching fractions for leptons to appear in squark and gluino induced cascade decays are generally not very large in our model sample, as we have seen in our earlier work [9]. () Independently of the specific search, as the background systematic errors become large, the pMSSM model coverage is seen to increase more slowly with the integrated luminosity. Some of the search channels nearly saturate at high luminosity due to the large background uncertainties and thus become systematics dominated.
As discussed above, the final step in the ATLAS analyses is to apply a cut on , where the particular value of that is chosen (in units of 400 GeV) is the one that maximizes the signal significance given the SM background and its corresponding uncertainty. This choice not only depends upon the particular channel but also on the amount of integrated luminosity. It is important to note that if this cut is taken to be too large when maximizing the signal, then the analysis will be very sensitive to the detailed shape in the tails of both the signal and expected background distributions, especially with higher luminosities. This happens when there are very few events with large values of , e.g., 2 TeV. In this situation, small fluctuations in the SM background and/or SUSY signal expectations due to limited Monte Carlo statistics can lead to inconsistencies in whether a given model is observable in a specific analysis or even whether or not it would be detected overall. This effect only occurs in the case of the search analyses designed for the 7 TeV run, as the cut applied in the planned 14 TeV ATLAS analyses was fixed at relatively low values.
Figure 10 shows the optimized value for the cut for the nj0l analyses, as an example, for the flat prior model sample. The three things we see here are: () for large background systematic errors, a harder cut is required to optimize the search, () as the number of required jets in the final state decreases, the strength of the cut can be reduced. Both of these results are also found to hold for the nj1l and njOSDL searches although the cut itself turns out to be less important as the number of required leptons in the final state increases. () Given the warning about distribution tails in the discussion above, it is a welcome result to see that in the majority of cases only a moderately strong cut is required to optimize the signal significance. Note that there are some cases where the cut does not contribute very much to increase the significance of the SUSY signal; this happens in particular for scenarios where the background is low and the effect of systematic uncertainties is not very significant.
It is interesting to evaluate the fraction of models that can be discovered in multiple analyses. This is important to consider as, first, it is valuable to verify the discovery of new physics in more than one channel, and second, the availability of multiple discovery channels admits for the possibility of further studies that will allow for, e.g., the extraction of sparticle masses. To be specific, Tables 2, 3 and 4 show the fraction of pMSSM models which are observed in (exactly) channels assuming a SM background systematic error of 20, 50, and 100%, respectively. Note that the distribution shifts towards more models being observed in multiple channels as the integrated luminosity increases and the background systematic error decreases, as expected..
# anl.  Flat  Flat  Flat  Log  Log  Log 

0  38.172  7.5501  0.9965  63.64  43.988  22.92 
1  9.2928  4.1988  0.90862  5.376  4.8674  5.8482 
2  8.7432  4.6665  1.6102  3.6687  5.6665  6.0298 
3  41.836  59.878  39.573  26.008  34.907  35.38 
4  0.65686  4.9257  7.9422  0.25427  2.2158  6.4657 
5  0.53472  4.2629  6.7163  0.47221  2.0341  4.8311 
6  0.54366  8.5391  13.494  0.32692  3.0875  6.5383 
7  0.067026  2.5217  8.9044  0.21794  1.453  4.1773 
8  0.062558  1.2288  5.6364  0.036324  0.72648  2.2884 
9  0.077452  1.2958  6.548  0  0.58118  2.9422 
10  0.013405  0.93241  7.6711  0  0.47221  2.579 
# anl.  Flat  Flat  Flat  Log  Log  Log 

0  54.756  21.772  4.8782  71.558  55.903  32.546 
1  14.143  10.547  4.847  8.1729  7.3011  9.8801 
2  7.8435  11.453  9.959  5.0854  7.1195  12.532 
3  22.552  42.949  40.705  14.857  24.228  28.478 
4  0.29938  4.1407  8.3533  0.18162  1.7436  4.5768 
5  0.15788  3.1562  7.619  0  1.3803  3.4871 
6  0.1415  3.3036  9.1487  0.072648  1.0534  3.4871 
7  0.061068  1.4075  6.049  0.036324  0.79913  1.9615 
8  0.031279  0.58536  3.6166  0.036324  0.32692  1.4166 
9  0.013405  0.43493  2.9716  0  0.036324  1.235 
10  0.0014895  0.25172  1.853  0  0.10897  0.39956 
# anl.  Flat  Flat  Flat  Log  Log  Log 

0  74.112  47.23  17.635  81.911  69.016  47.875 
1  13.894  16.834  15.996  9.8438  10.825  14.094 
2  4.4759  13.331  21.917  3.1602  7.7007  15.91 
3  7.3282  18.166  26.186  4.9401  10.607  15.365 
4  0.10575  1.8827  6.478  0.036324  0.79913  2.6153 
5  0.037237  1.0322  4.7174  0.036324  0.32692  1.7436 
6  0.023832  0.7075  3.3008  0  0.32692  1.0171 
7  0.019363  0.57345  2.1703  0.072648  0.36324  0.87178 
8  0.0029789  0.16831  1.0025  0  0.036324  0.36324 
9  0.0014895  0.064047  0.46474  0  0  0.10897 
10  0  0.011916  0.13257  0  0  0.036324 
By combining our results for these searches, we can also determine the fraction of pMSSM models that are undetected in all of the 7 TeV search analyses designed by ATLAS; this corresponds to the case of in these Tables. Figure 11 presents the fraction of pMSSM models which are undetected in all of the search channels as a function of integrated luminosity for both the flat and log prior sets. In the flat prior case we see that as the integrated luminosity increases, the model coverage substantially improves, and approaches (or exceeds) for 10 fb with SM background uncertainties of 50% or less. In the log prior case, the improvement in pMSSM model coverage as the luminosity increases is much more gradual as we expected. However, even in this case, at high integrated luminosities substantial model coverage is seen to be obtainable at 7 TeV.
This figure also shows the important playoff between increasing the integrated luminosity and decreasing the SM background systematic error in terms of pMSSM model coverage. (Of course, increased luminosity often results in decreased systematic errors, up to a point.) For example, it is interesting to compare the effectiveness of the analyses for the flat prior set assuming and fb with other values. Figure 11 shows that taking and fb or with fb produces essentially identical model coverage. This demonstrates that small reductions in the SM background uncertainty can be worth a significant amount of increased integrated luminosity in terms of pMSSM model coverage.
Another very important message to take home from this figure is that for the integrated luminosity collected during the 2010 LHC run at 7 TeV ( pb for ATLAS), some significant fraction of these pMSSM models should already have been observed at the level. Explicitly, for a background systematic error of 20(50,100)% we find that 46(28,15)% of the flat prior model set should have been discovered; in the log prior case the corresponding results are found to be 30(24,14)%. Since these are discovery results, an even greater portion of the pMSSM model sample would be expected to be excluded by these analyses. This shows the incredible power of going from Tevatron energies up to the 7 TeV LHC in performing searches for massive objects, such as SUSY sparticles, even when only small amount of integrated luminosity is available.
3.3 Why are Models not Detected by the ATLAS Searches?
Here, we investigate the main reasons why some pMSSM models are not discoverable in the ATLAS searches at 7 TeV. We addressed this question in some detail in our earlier work for the case of the analyses designed for the 14 TeV LHC [9], so our discussion here will not be as extensive. Clearly, many of our previous results will carry over qualitatively into the present 7 TeV analysis.
There are multiple explanations as to why some pMSSM models are undetected by the searches, the most obvious one being small production cross sections for the colored sparticles that initiate the familiar SUSY cascades. As an example, we note that for our pMSSM models in the flat prior set, the cross sections for the production of gluino and squark pairs are found to cover an enormous range of several orders of magnitude as can be seen in Fig. 12 (recall that the upper limit on sparticle masses in our flat model set is TeV). Here we see that the large (or small) values of the gluino pair cross section is completely uncorrelated with the corresponding values for first generation squarks within a particular model. Furthermore, by summing over all of the QCD production channels involving gluinos and/or first generation squarks (i.e., , , and ) we obtain an approximate handle on the total overall rate for SUSY production which we see ranges over four orders of magnitude. Note that for any particular value of the squark or gluino mass, the corresponding production cross section itself can vary by up to an order of magnitude or more depending upon the remainder of the pMSSM model spectrum.
While it is certainly clear from this figure that some models have too small a cross section to be discovered, the bottomright panel indicates that this cannot be the entire explanation. Here we show the search significance, , of the 4j0l channel (as it is the most powerful channel in terms of discovery capability) as a function of the total NLO QCD production cross section assuming fb and . Here we observe that () there are models with cross sections pb which are missed by this analysis, while () there are models with cross sections fb which are discovered. () For any given value of the cross section, the range of the significance is large and can be up to two orders of magnitude or more. This validates the claim that there are reasons other than small production cross sections that render models unobservable by these analyses. () For any given gluino mass there is a strong correlation of the signal significance with the mass of the LSP. Clearly when these two masses are close the average of the jets will be softer and this will make it more difficult to pass analysis cuts. Also if the LSP mass is large then that implies even larger squark and gluino masses that will result on average in smaller production cross sections. In the case of the nj0l analyses, a larger number of signal events is required for discovery due to the sizeable SM backgrounds and hence such models will be missed by these analyses. Visibility will then require the production of leptons with significant in cascade decays in order to pass the lower background nj1l, njOSDL and 2jSSDL searches. Unfortunately, lepton branching fractions are low in these cascades in our model sample (since, e.g., sleptons are heavy) and thus some models will be missed entirely. Of course, models that are leptonrich will automatically fail all the nj0l analyses since they veto events with high leptons, but will be picked up by the searches containing leptons.
Though the production cross section for SUSY particles is reasonably well correlated with their masses, we can ask more directly if larger sparticle masses lead to their nonobservation in these searches. Figure 13 shows that, indeed, models with lighter squarks or gluinos tend to lead to signals with greater significance in the 4j0l channel. However, as we can also see from this Figure, this is not true universally, e.g., there are many models with gluino (lightest squark) masses below 300 (200) GeV that have , while conversely there are models with 1 TeV gluinos that have . We see that for any given squark or gluino mass the value of can vary significantly. The top panel of Fig. 14 displays this property even more strongly where we see that these results hold even when all of the searches are combined. This Figure shows the set of flat prior models that are unobservable in all the search channels in average light squark mass–gluino mass plane. Note that there are a significant number of these models which contain light squarks and gluinos. Thus while the masses of the colored sparticles do play an important role in model observability clearly there are additional important factors.
In the lower panels of Fig. 14 we see that the mass splitting between squarks and/or gluinos and the LSP can play an important role in determining
model observability as was first noted in Ref. [6] and was seen explicitly in our earlier work on the generation of
the pMSSM models [7]
and the 14 TeV ATLAS SUSY analyses [9]. The obvious reasoning here is that as the degeneracy in the spectrum
increases and mass splittings become smaller, the values of, e.g., the of the jets, will be
reduced so that it will be more difficult to satisfy any of the analysis
cuts. These figures show this result explicitly. Note, however, that, e.g., in the case of light gluinos with small gluinoLSP splittings there are many models which
are still discoverable in the 4j0l channel. The reason for this is that while the
efficiency
Another cause for models being undetected is the occurrence of detectorstable sparticles at the end of gluino or squark induced decay chains instead of the LSP [9]. This happens with reasonable frequency in both the 7 and 14 TeV analyses. In such cases, the amount of that is produced is substantially decreased which reduces the capability of the relevant models to pass any analysis requirements. Most commonly, these sparticles are actually longlived charginos that are reasonably degenerate with the LSP in wino or Higgsinolike LSP scenarios. In these cases, searches for longlived sparticles, as discussed below, will be a important supplement to the conventional searches. Of course, a loss of the signature can happen in other ways. For example, if the initial squark or gluino produces a very long decay chain then the particles produced at the end of such a chain will be somewhat soft. In some cases this may lead to the inability to pass the necessary and/or requirements for the various based searches and the model will not be observed. Such long decay chains were shown to occur with a reasonable frequency in our earlier work [9] and can be a contributor to models failing to pass the analyses requirements.
There are, of course, other reasons that prevent models from being discovered. As noted above, subtleties in any sector of the sparticle spectrum can make a significant difference as to whether a given model is observed by various analyses. Here, we will discuss a couple of examples where this occurs. As in Ref. [9], the approach we follow is to compare a model which fails to be observed in all search channels to one with a similar spectrum (dubbed a ‘sister’ model) that is detected in at least one channel and then examine the difference between them. For this study, we concentrate on the more difficult cases by taking the flat prior model sample and assume fb and ; this leaves only 670 models that are not observed. In order to avoid the statistical issues associated with the tails of the distribution discussed above, and to further reduce this model set to a more manageable size, we will only consider models whose optimized value is GeV. We note that a large fraction of models in this set have relatively heavy LSPs with masses in excess of 400 GeV. We now discuss two brief examples of these comparisons.
Figure 16 shows a comparison of two similar models, 8944 (observed in the 3,4jOSDL channels) and 21089 (missed by all analyses), that have comparable total colored sparticle production rates (3.4 and 4.6 pb, respectively). Both of these models are not observed in the nj0l searches since the lighter squarks are too close in mass to the LSP to produce hard jets. The gaugino sectors of these two models are quite similar (with the LSP and being Higgsinolike and being binolike), while their colored sparticle spectra are somewhat different. In either model the decay of allows for the OSDL production through an intermediate slepton which has sufficient to pass the analysis requirements. However, while 8944 has a with mass below that of the gluino (which is light enough to give a reasonable cross section), allowing for the decay into , only is (sufficiently) heavier than the in model 21089. In this case, is also much more massive than the gluino, through which it will dominantly decay, and so it will not have a large enough branching fraction into to produce the OSDL signature.
Figure 17 compares models 9781 (discovered in the 2jSSDL channel) and 20875 (completely missed) that have total colored sparticle production cross sections of 1.3 and 1.1 pb, respectively, but yet produce too few hard jet plus events to be found in the nj0l channels due to spectrum compression. Model 9781 is quite interesting as the charginos and lightest three neutralinos are highly mixed combinations of winos and Higgsinos. In this model, (which is relatively light) decays to with a branching fraction. Since has a large bino content, it decays of the time through sleptons which subsequently decay directly to the LSP with a branching fraction of . Thus model 9781 can easily populate the leptonic final state and since the neutralinos are Majorana fermions, the 2jSSDL final state becomes accessible, a channel with an extremely small background. On the other hand, model 20875 does not allow for the generation of a typical leptonic signal. This is because the are quite heavy so that the neutralinos only allow for decay to leptons via an intermediate onshell . This is correlated with the relative lightness of the as well as the Higgsinolike LSP. In addition, the lightest squarks decay primarily directly to the LSP which will not produce any high leptons.
These two examples demonstrate that the full sparticle spectrum may conspire to render a model undiscoverable in the based analysis suite, even if the colored sparticle cross section is large. Discovery is not based on the value of the squark and/or gluino masses alone and blanket limits that claim are ruled out below some value cannot be set.
3.4 Detector Stable Sparticles in Cascades
As we mentioned above, one of the reasons that some pMSSM models may not be observed in the ATLAS analyses is that squark and gluino cascade
decays can sometimes lead to a final state with low . In many cases this is due to the
existence of longlived sparticles, usually charginos, which appear with sizable
branching fractions in such cascade chains and are essentially detector stable. Particles that decay outside the detector when produced at the LHC would
provide a dramatic signal of new physics (see [23] and references therein).
In fact, data from the LHC are already extending the mass limits on such detectorstable
sparticles [24]. Since the inclusive analyses discussed above do not consider such sparticles, we will briefly sketch the TeV
discovery prospects for detectorstable sparticles in our model sets. As discussed in our earlier work [9], the existence of such longlived
states is relatively common in our pMSSM model sample as can be seen in
Fig 18.
Sparticle  LHC Reach pb GeV  LHC Reach fb  LHC Reach fb 

(Winolike)  GeV  GeV  GeV 
(Higgsinolike)  GeV  GeV  GeV 
GeV  GeV  GeV  
GeV  GeV  GeV  
GeV  GeV  GeV 
Our estimation of the TeV LHC mass reach for each longlived sparticle is shown in Table 5 for the specified integrated luminosities, assuming direct pair production. These results are deduced from Figure 1 of [25] by taking the geometric mean of the search reach for a and TeV LHC at the specified luminosities and are interpolated to the luminosities considered here where necessary. These results are somewhat conservative as only detectorstable sparticle production in the hard process is considered; additional detectorstable sparticles could be produced through cascade decays as discussed below.
The number of detectorstable sparticles of various species in our model sample is shown in Table 6. This Table also shows the number of detectorstable sparticles which will not be discovered at LHC with pb and fb of integrated luminosity, using the approximate mass reaches presented in Table 5. We assume here that the mass reach is roughly generationindependent and that it is the same for, e.g., stops and sbottoms. This assumption is reasonable except where there could be significant channel production for the first or (to a lesser extent) second generation, for instance in the case of up or down squarks.
Sparticle  In Model Set  LHC Reach pb  LHC Reach fb  LHC Reach fb 

We now focus on the specific case of longlived charginos, which are by far the most common longlived sparticles in our model sets. If the production cross section for colored sparticles (times relevant branching fractions into charginos) are sufficiently large, these stable charginos should be found in searches for (effectively) stable charged particles occurring at the end of a cascade decay chain. In Fig. 18 we display the estimated value for for the production of detectorstable charginos in cascade decays in our flat prior set. (Note that this does not include the direct contribution arising from direct chargino pair production.) Here, we see that roughly of models with detectorstable charginos lead to values in excess of 10 fb at 7 TeV and so we expect them to be observable in the upcoming run of the LHC. In this estimation, we assumed the largest contribution to the production cross section arises from the production of gluino and light squark species. Using the decay tables generated for each model with detectorstable charginos, we calculated the branching fraction for the gluino and light squarks to produce a stable chargino at the end of each possible decay chain and then weighted them by their corresponding production crosssections. Note that mass information for neither the mother particles or the daughter chargino is used to indicate how likely it is that the chargino will pass the trigger criteria for detection. As we used PYTHIA to compute the LO crosssection for the light squarks, we do not separately generate the production for the various light squark species. Thus in order to obtain our estimate, we make the assumption that the overall cross section is 100% dominated by that arising from the lightest squark. This assumption will break down when the squarks are nearly degenerate, especially when their branching fractions to charginos are vastly different due to the complexities in the gaugino sector. Many of our models have a large production crosssection for colored sparticles, but the relevant branching fractions to charginos can be simultaneously quite small. This can result in very small overall production rates for stable charginos well below fb.
3.5 SUSY Mass Scale From
In our earlier work on the 14 TeV ATLAS analyses, we demonstrated that the relationship between and the mass of the lightest colored sparticle found in mSUGRA, i.e., (where LCP stands for Lightest Colored Particle), proposed long ago [26] does not necessarily hold in the pMSSM. This possible relationship is important as it might be used to get the first handle on the overall mass scale of the sparticle spectrum. Here, we briefly note that this result remains valid for the 7 TeV ATLAS analyses as can be seen in Fig. 19. For both the 4j0l and 2j0l channels we see explicitly that the values of lie mostly above the expected value of , especially in the low sparticle mass region. However, for lightest colored sparticle masses in excess of GeV we see that, indeed, the relationship provides a fairly good estimate in both of these searches.
3.6 Modifying ATLAS SUSY Analysis Cuts
Given the properties of the various sparticles in our model sets, we can try to determine whether the canonical cuts employed in the ATLAS based search analyses can be strengthened to reduce SM backgrounds without any significant loss in the coverage of our pMSSM model space. This is certainly a nontrivial issue and the structure of our analysis, being based on the fixed ATLAS analyses cuts, is not directly set up to obtain completely definitive answers. However, it is possible to make some reasonable estimates based on the information that we do have available. We will concentrate on the three nj0l analyses as they generally provide the greatest pMSSM model coverage and have large statistics. The most important kinematic quantities for these searches are the requirements on the leading jet and the required amount of . Here, we make use of the average values of the distributions in these quantities for our pMSSM model sample, as well as the corresponding fitted width of the part of the distribution below this average value for preselected events. This information then provides us with an estimate of where these two kinematic distributions ‘turn on,’ which we take to be the average value minus this width, on the low energy/momentum side below their peak average values.
First consider the cut on the leading jet, , for the 2(3,4)j0l analyses; ATLAS chooses the value for this cut to be 180(100,100) GeV, respectively. Fig. 20 shows the distribution of the ‘turnon’ values for these three ATLAS analyses obtained from analyzing the flat prior model set. For the 2j0l analysis, we see that the lower edge of the ‘turnon’ values lies somewhat below the cut value of 180 GeV from which we can conclude that this cut is already reasonably hard and cannot be increased without a loss of model coverage. However, for both the 3j0l and 4j0l analyses we instead observe that the ‘turnon’ values lie above those of the ATLAS cuts by GeV suggesting that that the cut in these two channels may be raised without impacting coverage rates for the pMSSM.
Similarly, the nominal cut imposed by ATLAS for the nj0l analyses is 80 GeV. However, there is an additional subsequent cut imposed by ATLAS based on the value of , i.e., where for the analyses. Fig. 21 shows the distributions for the ‘turnon’ values of for these three channels employing the flat prior model set. Here we see that the lower edge of this distribution occurs at GeV for , respectively. This is suggestive that the nominal cut made by ATLAS may be increased for these three analyses without losing significant pMSSM model coverage.
4 Implications of the 7 TeV Run
In this section, we explore some implications of a null search for Supersymmetry at the 7 TeV LHC. We examine the degree of finetuning that would be placed on our pMSSM model sample and we discuss the resulting expectations for sparticle production at a 500 GeV and 1 TeV Linear Collider.
4.1 Finetuning in the Undiscovered pMSSM Models
As has recently been discussed in the mSUGRA/CMSSM context [27, 28], it is apparent that if SUSY signatures are not discovered at the
7 TeV LHC as the integrated luminosity accumulates it is likely that the SUSY parameter space must become more finetuned, and hence more problematic as a solution
to the hierarchy problem. Since we know which models in our sample are discoverable (or not) by the ATLAS search analyses,
we can ask whether this same result also holds in the
case of our pMSSM model sets.
Figure 22 display the results of this analysis assuming a background systematic error of 50% for both the flat and log prior model samples. Clearly, in the flat prior case, one sees that as the integrated luminosity increases and more models can be discovered by ATLAS, those remaining yet undiscovered tend to be more finetuned as expected. In other words, the fractional loss of models from the full distribution occurs more rapidly with increasing luminosity for models with smaller amounts of finetuning. This is not too surprising as, overall, models with less tuning tend to have lighter SUSY sparticle spectra and are thus more easily discovered at the LHC. On the other hand, the results from the log prior model set appear to be affected somewhat differently in that the overall shape of the finetuning distribution does not appear to change very much by removing models that should have been already been discovered by ATLAS as the luminosity increases. In this case, we see that there is not much of an increase in the amount of finetuning as the set of undiscovered models shrinks. This represents one of the few apparent differences between these two different model sets. This can be explained by the fact that while the log prior model set tend to have light sparticle spectra (though they extend out to larger mass values than do those for the flat prior models) and are thus less finetuned to begin with, these same spectra are generally compressed making these models more difficult to discover at the LHC as was discussed above. This would imply that the models missed by the ATLAS analyses in the log prior case are generally not much more finetuned than those appearing in the originally generated model set. Thus we find that the amount of finetuning that remains in the LHCundiscovered pMSSM model sample can depend upon the prior used to generate the original model set.
4.2 Implications of pMSSM Searches for the Linear Collider
If there are no clear SUSY signals as the 7 TeV LHC integrates more luminosity, the question arises whether the production (and study) of charged sparticles
remains viable at the proposed 500 GeV Linear Collider (LC). Based on mSUGRA/CMSSM model coverage projections from both the ATLAS [12] and
CMS [29] Collaborations at 7 TeV with a 1 fb integrated luminosity (as well as their results from the 2010 SUSY searches), it would seem very
unlikely that either light sleptons or gauginos (other
than perhaps the LSP itself) will remain kinematically accessible at a 500 GeV LC if nothing is observed. Here, we address the question whether
this expectation also remains true for
our pMSSM model sets. To this end we examine the set of flat and log prior models which are not detected by any of the ATLAS search analyses for
assumed values of both the integrated luminosity and background systematic error and then determine the part of the sparticle spectrum within these models which is
kinematically accessible at a 500 GeV LC.
We present the results of this analysis in various ways. Tables 7 and 8 show the number of sparticles of various species that are kinematically accessible at a 500 GeV or 1 TeV LC within the subset of ATLASundetected models assuming fb with , respectively. These two cases represent both a possibly conservative and a more optimistic performance for the LHC over the 201112 running period. Here we see several things: () In addition to the gauginos and sleptons, sparticles such as stops, sbottoms and other squarks are potentially almost as likely to also be kinematically accessible at the LC [32]. () The number of models with kinematically accessible sparticles and their variety is significantly greater in the log prior sample as these models are more likely to have a lighter and more compressed sparticle spectrum () The difference between the two cases presented in these Tables is quite significant; in particular, we see that for the flat prior model set there is a huge depletion in the number of unobserved models at higher luminosity and with lower background systematics. There are extremely few flat prior models remaining at high luminosity with any accessible sparticles at a 500 GeV LC. () Going from a 500 GeV to a 1 TeV LC substantially increases the number of models with kinematically accessible sparticles, especially in the flat prior case. It is clear that, at least for the flat prior model sample with luminosities in excess of 1 fb at the LHC, that the 500 GeV LC does not seem to be a good place to study our pMSSM models if no signal for SUSY is found at the LHC in 201112.
GeV  TeV  

Sparticle  Flat  Log  Flat  Log 
107  101  3052  347  
260  209  3938  565  
730  381  7431  869  
30  36  1288  207  
151  117  3168  356  
386  236  4366  553  
5487  1312  14,510  1539  
2738  1035  10,714  1395  
429  352  5667  903  
10  18  1267  202  
4856  1208  13,561  1495  
94  54  3412  456  
0  0  1088  65  
35  11  2459  117  
220  96  3630  526  
52  16  2545  123  
124  64  3581  273  
289  75  5553  590  
1  0  409  21  
93  9  3727  217  
0  0  2  0 
GeV  TeV  

Sparticle  Flat  Log  Flat  Log 
0  37  63  142  
0  72  53  223  
2  142  165  338  
0  11  9  69  
0  42  64  146  
0  85  81  236  
26  507  587  626  
4  397  352  557  
0  136  57  357  
0  5  5  66  
25  467  505  608  
0  17  16  170  
0  0  27  5  
0  3  73  24  
1  18  63  157  
0  5  81  24  
0  14  86  79  
0  20  103  189  
0  0  3  4  
1  2  94  58  
0  0  0  0 
In order to study these LC results in more detail we examine their dependence on the LHC integrated luminosity; this is shown for the 500 GeV LC in Fig. 23 and for the 1 TeV LC in Fig. 24. Here we display the fraction of the unobserved set of models that have a kinematically accessible sparticle of a particular variety. At a 500 GeV LC, this fraction for charginos, stops, sbottoms, selectrons (or smuons) and staus in the flat prior model set is seen to decrease significantly as the LHC accumulates integrated luminosity without observation of a signal for SUSY. However, note that for the corresponding log prior model set, while the number of surviving models decreases with any corresponding increase in the LHC integrated luminosity (or with any decrease in the SM background uncertainty) as one would expect, the fraction of the surviving models with a kinematically accessible sparticle changes very little, if at all! For charginos at a 500 GeV LC this fraction is found to be quite large, , but it is found to be somewhat smaller for the other sparticles, e.g., for and for .
One possible explanation of this unexpected behavior in the log prior sample is as follows: As we saw in the previous discussion of finetuning, in the log prior case, the removal of pMSSM models from the log prior set (as they are ‘discovered’ by ATLAS) must affect the various sparticle mass distributions in a roughly uniform manner. Otherwise the observed amount of finetuning would necessarily increase. However, in the flat prior case, models with lighter sparticles are preferentially ‘discovered’ by ATLAS searches. This hypothesis can qualitatively explain why there is no significant reduction in the fraction of the log prior models with kinematically accessible sparticles at a 500 GeV (or 1 TeV) LC (especially in the case of noncolored sparticles as we have seen above). It also simultaneously explains why the apparent amount of finetuning does not change appreciably as the LHC covers more of the log scan parameter space. Of course, as can be seen from these figures, at a 1 TeV LC a substantial portion of both the log and flat prior model sets which remain undiscovered at the LHC have sparticles which are kinematically accessible even at high LHC integrated luminosities. We further note that the fraction of LHCunobserved models where no SUSY sparticles whatsoever are accessible at a LC (not even the LSP) is quite large for the flat model set at a 500 GeV LC, but is only in the 1220% range for the corresponding log prior model set. On the other hand, at a 1 TeV LC, this value is seen to lie below for the log prior models, while for the flat prior set it remains below . Clearly, a 1 TeV LC will be far better at accessing the sparticles in our pMSSM model sets.
In order to test this hypothesis, we show the mass distributions for the and in both the flat and log prior model sets in
Fig. 25.
5 Summary and Conclusions
In this paper we have analyzed the capability of the ATLAS based SUSY analyses to discover Supersymmetry in a model independent fashion at the 7 TeV LHC. To this end, we tested these search channels on a large set of model points, k, in the 19dimensional parameter space of the pMSSM. This model sample contains a wide variety of properties and characteristics and provides a framework to explore the breadth of possible SUSY signatures at colliders and elsewhere. These models were generated in a previous work and comply with a set of minimal theoretical assumptions as well as the global precision electroweak, heavy flavor, collider searches, and astrophysical data sets. We simulated ten ATLAS search channels, which were designed in the context of mSUGRAbased SUSY, and employed the SM backgrounds as provided directly by the ATLAS SUSY working group. We first checked that our analyses were in agreement with ATLAS results for the mSUGRA benchmark point that the collaboration had previously simulated.
We passed our model sample through the ATLAS analysis chain and computed the significance of the signal for each model in each search channel. A significance of was used as the criteria for discovery in each channel; we employed the same numerical technique that ATLAS does for calculating this value. We found that the systematic error due to uncertainties in the size of the expected background made a substantial impact on model discovery for the range of expected integrated luminosities. In fact, some channels become systematics dominated at luminosities of order fb. Overall, for 1(10) fb of integrated luminosity roughly 80(95)% of the flat prior model sample is discoverable, assuming a 50% background systematic error. Larger (or smaller) systematic errors greatly reduce (or increase) this model coverage. We found that the 4j0l channel is the most powerful in terms of observing a signal, whereas the leptonic channels had a much reduced model coverage. This is due to the suppression of leptonic cascade decays appearing in our model sample compared to expectations from e.g., mSUGRA. The model coverage was worse for the log prior sample due to kinematic reasons.
We explored the characteristics that caused a model to not be observed in these search channels. While production cross section values as related to the sparticle mass obviously plays a role, it does not tell the whole story. There are cases with low mass gluinos and/or squarks with large cross sections that are missed by these search analyses, while models with heavy masses and small cross sections are sometimes observed. We found that the mass splitting between the gluino/squarks and the LSP plays an important role in detecting models, and that this can sometimes be compensated by very large production rates or ISR. We also saw that subtleties in the sparticle spectrum can conspire to render a model to not be detected. A fraction of our model set contains detectorstable sparticles which appear at the end of their cascade decay chains and hence are not detected by the based searches. We studied the effectiveness of the planned stable charged particle searches in these cases and found that some, but not all, of these models will be discovered.
We briefly considered potential modifications to the ATLAS kinematic cuts in these analyses that would improve their discovery potential. We studied the optimal cut on as well as for the of the leading jet and overall . Our results indicate that the cuts for both the leading jet transverse momentum and the could be increased from their nominal value without seriously impacting model coverage.
Lastly, we studied the implications of a null result from the 7 TeV LHC run. We found that the degree of finetuning that would be imposed on the pMSSM depended on the choice of priors which generated the model sample, but overall would not be as large as in the case of mSUGRA. However, the expectations for sparticle production at a high energy Linear Collider would be greatly impacted if Supersymmetry is not discovered during this LHC run. Basic kinematics would essentially exclude sparticle production at a 500 GeV Linear Collider, and would point towards the need for a higher energy machine in order to study Supersymmetry.
In summary, we find that the mSUGRA motivated based searches for Supersymmetry perform well over a larger and more complicated SUSY parameter space such as the pMSSM. However, there are some exceptions and coverage is not perfect. The details of the full sparticle spectrum play a very important role in the observability of a model. There are no blanket statements regarding the potential for discovery, or in setting a mass limit, that that can honestly be made.
We anxiously await the discovery of Supersymmetry in the near future.
Acknowledgments
This project would not have been possible without the assistance and input from many people.
The authors would like to thank the members of the ATLAS SUSY group, in particular, S. Caron, P. de Jong, G. Polesello and G. Redlinger, for discussions and, most importantly, for providing us with additional information detailing the ATLASgenerated SM background distributions for their SUSY studies and the corresponding results for their mSUGRA benchmark analyses.
We would also like to thank the SLAC ATLAS group for assistance with computing resources.
We would like to thank T. Plehn for his help overcoming the special problems we had with the implementation of PROSPINO for our pMSSM model set. We would also like to thank A. Djoudai and J. Conway for similar help with the implementations of SDECAY and PGS, respectively.
We would like to thank L. Dixon for discussions about the theoretical assessment of the size of the systematic errors for the SM backgrounds to SUSY signals.
We would like to thank R. Cotta for important discussions related to our analysis.
We would like to thank J. Cogan for his assistance with setting up the ATLAS analysis code during the earlier part of this work.
Work supported by the Department of Energy, Division of High Energy Physics, Contracts DEAC0276SF00515 DEAC0206CH11357, and DEFG0291ER40684, and by the BMBF “Verbundprojekt HEPTheorie” under contract 05H09PDE.
Footnotes
 Even in such a case a full exploration of this large parameter space is at best difficult if not impossible with present computing power.
 For full details, see Ref. [7]
 Here, efficiency is defined as the fraction of generated signal events that pass the analysis cuts.
 Recall that we will define a particle to be detectorstable if its unboosted decay length in at least 20m. Note that typical values of for longlived particles resulting from cascade decays are in the 23 range. Further note that the dependence of the number of detectorstable sparticles of various species in this model set on the value of the decay width that is assumed is discussed in detail in [9].
 Note that the amount of finetuning in both our model sets was examined in some detail in our earlier work [7]
 Note that we have not performed any analysis here to ascertain whether or not a given kinematically accessible sparticle is actually observable at such a LC. See, however, the work in Refs. [30, 31].
 Note that there is nothing special about the choice of these two particular noncolored sparticles and the features that we will now describe are also found in the mass distributions of other sparticles.
References
 For recent overviews of LHC searches for new physics, see E. Gross (ATLAS) and H. Flaecher (CMS), talks given at LaThuile 2011, Les Rencont res de Physique de la Vallée d’Aoste, 27 Feb5 Mar, 2011.
 V. Khachatryan et al. [ CMS Collaboration ], [arXiv:1101.1628 [hepex]]; CMS Collaboration, [arXiv:1103.0953 [hepex]]; CMS Collaboration, [arXiv:1103.1348 [hepex]]; G. Aad et al. [ Atlas Collaboration ], [arXiv:1102.2357 [hepex]]; ATLAS Collaboration, [arXiv:1102.5290 [hepex]].
 H. E. Haber and G. L. Kane, Phys. Rept. 117, 75 (1985); S. P. Martin, arXiv:hepph/9709356; D. J. H. Chung, L. L. Everett, G. L. Kane, S. F. King, J. D. Lykken and L. T. Wang, Phys. Rept. 407, 1 (2005) [arXiv:hepph/0312378]; L. Pape and D. Treille, Rept. Prog. Phys. 69, 2843 (2006); H. K. Dreiner, H. E. Haber, S. P. Martin, Phys. Rept. 494, 1196 (2010); M. Drees, R. Godbole and P. Roy, Hackensack, USA: World Scientific (2004) 555 p; H. Baer and X. Tata, Cambridge, UK: Univ. Pr. (2006) 537 p.
 G. Aad et al. [ The ATLAS Collaboration ], [arXiv:0901.0512 [hepex]].
 G. L. Bayatian et al. [ CMS Collaboration ], J. Phys. G G34, 9951579 (2007).
 J. Alwall, M. P. Le, M. Lisanti, J. G. Wacker, Phys. Lett. B666, 3437 (2008). [arXiv:0803.0019 [hepph]].
 C. F. Berger, J. S. Gainer, J. L. Hewett et al., JHEP 0902, 023 (2009). [arXiv:0812.0980 [hepph]].
 J. Alwall, P. Schuster, N. Toro, Phys. Rev. D79, 075020 (2009). [arXiv:0810.3921 [hepph]]; N. ArkaniHamed, P. Schuster, N. Toro, J. Thaler, L. T. Wang, B. Knuteson, S. Mrenna, [hepph/0703088 [HEPPH]].
 J. A. Conley, J. S. Gainer, J. L. Hewett et al., [arXiv:1009.2539 [hepph]].
 S. Akula, N. Chen, D. Feldman, M. Liu, Z. Liu, P. Nath, G. Peim, [arXiv:1103.1197 [hepph]]; S. Heinemeyer, [arXiv:1103.0952 [hepph]]; B. C. Allanach, T. J. Khoo, C. G. Lester, S. L. Williams, [arXiv:1103.0969 [hepph]]; P. Athron, S. F. King, D. J. Miller, S. Moretti, R. Nevzorov, [arXiv:1102.4363 [hepph]]; O. Buchmueller, R. Cavanaugh, D. Colling, A. De Roeck, M. J. Dolan, J. R. Ellis, H. Flacher, S. Heinemeyer et al., [arXiv:1102.4585 [hepph]]; S. Scopel, S. Choi, N. Fornengo, A. Bottino, [arXiv:1102.4033 [hepph]]; B. C. Allanach, [arXiv:1102.3149 [hepph]]; D. Feldman, K. Freese, P. Nath, B. D. Nelson, G. Peim, [arXiv:1102.2548 [hepph]]; O. Buchmueller, R. Cavanaugh, D. Colling, A. De Roeck, M. J. Dolan, J. R. Ellis, H. Flacher, S. Heinemeyer et al., Eur. Phys. J. C71, 1583 (2011). [arXiv:1011.6118 [hepph]].
 R. C. Cotta, J. A. Conley, J. S. Gainer, J. L. Hewett, T. G. Rizzo, JHEP 1101, 064 (2011). [arXiv:1007.5520 [hepph]]; R. C. Cotta, J. S. Gainer, J. L. Hewett, T. G. Rizzo, Nucl. Phys. Proc. Suppl. 194, 133144 (2009). [arXiv:0909.4088 [hepph]]; R. C. Cotta, J. S. Gainer, J. L. Hewett, T. G. Rizzo, New J. Phys. 11, 105026 (2009). [arXiv:0903.4409 [hepph]].
 ATLAS Note, Prospects for Supersymmetry discovery based on inclusive searches at a 7 TeV centreofmass energy with the ATLAS detector, Tech. Rep. ATLPHYSPUB2010010, CERN, Geneva (2010).
 L. Pape, D. Treille, Rept. Prog. Phys. 69, 28433067 (2006); G. Weiglein et al. [ LHC/LC Study Group Collaboration ], Phys. Rept. 426, 47358 (2006). [hepph/0410364].
 A. Djouadi, J. L. Kneur, G. Moultaka, Comput. Phys. Commun. 176, 426455 (2007). [hepph/0211331].
 G. D’Ambrosio, G. F. Giudice, G. Isidori, A. Strumia, Nucl. Phys. B645, 155187 (2002). [hepph/0207036].
 F. Abe et al. [ CDF Collaboration ], Phys. Rev. D46, 18891894 (1992); V. M. Abazov et al. [ D0 Collaboration ], Phys. Rev. Lett. 102, 161802 (2009). [arXiv:0809.4472 [hepex]].
 V. M. Abazov et al. [ D0 Collaboration ], Phys. Lett. B660, 449457 (2008). [arXiv:0712.3805 [hepex]].
 We would like to thank the members of the ATLAS SUSY Physics Group for providing us with this information.
 A. Djouadi, M. M. Muhlleitner, M. Spira, Acta Phys. Polon. B38, 635644 (2007). [hepph/0609292].
 W. Beenakker, R. Hopker, M. Spira, P. M. Zerwas, Nucl. Phys. B492, 51103 (1997). [hepph/9610490]; W. Beenakker, M. Kramer, T. Plehn, M. Spira, P. M. Zerwas, Nucl. Phys. B515, 314 (1998). [hepph/9710451]; W. Beenakker, M. Klasen, M. Kramer, T. Plehn, M. Spira, P. M. Zerwas, Phys. Rev. Lett. 83, 37803783 (1999). [hepph/9906298]; M. Spira, [hepph/0211145]; T. Plehn, Czech. J. Phys. 55, B213B220 (2005). [hepph/0410063].
 P. M. Nadolsky, H. L. Lai, Q. H. Cao, J. Huston, J. Pumplin, D. Stump, W. K. Tung, C. P. Yuan, Phys. Rev. D78, 013004 (2008). [arXiv:0802.0007 [hepph]].
 J. Conway, PGS4, http://www.physics.ucdavis.edu/ conway/research/software/pgs/pgs.html
 M. Fairbairn, A. C. Kraan, D. A. Milstead, T. Sjostrand, P. Z. Skands, T. Sloan, Phys. Rept. 438, 163 (2007). [hepph/0611040].
 V. Khachatryan et al. [ CMS Collaboration ], [arXiv:1101.1645 [hepex]].
 A. R. Raklev, Mod. Phys. Lett. A24, 19551969 (2009). [arXiv:0908.0315 [hepph]].
 I. Hinchliffe, F. E. Paige, M. D. Shapiro, J. Soderqvist and W. Yao, Phys. Rev. D 55, 5520 (1997) [arXiv:hepph/9610544].
 A. Strumia, [arXiv:1101.2195 [hepph]].
 S. Cassel, D. M. Ghilencea, S. Kraml et al., [arXiv:1101.4664 [hepph]].
 CMS Note 2010/008.
 C. F. Berger, J. S. Gainer, J. L. Hewett et al., Phys. Rev. D80, 095018 (2009). [arXiv:0712.2965 [hepph]].
 C. F. Berger, J. S. Gainer, J. L. Hewett et al., Phys. Lett. B677, 4853 (2009). [arXiv:0711.1374 [hepph]].
 T. G. Rizzo, [arXiv:1005.4626 [hepph]].