SModelS: a tool for interpreting simplified-model results from the LHC and its application to supersymmetry

We present a general procedure to decompose Beyond the Standard Model (BSM) collider signatures presenting a Z2 symmetry into Simplified Model Spectrum (SMS) topologies. Our method provides a way to cast BSM predictions for the LHC in a model independent framework, which can be directly confronted with the relevant experimental constraints. Our concrete implementation currently focusses on supersymmetry searches with missing energy, for which a large variety of SMS results from ATLAS and CMS are available. As show-case examples we apply our procedure to two scans of the minimal supersymmetric standard model. We discuss how the SMS limits constrain various particle masses and which regions of parameter space remain unchallenged by the current SMS interpretations of the LHC results.


Introduction
The rest of this paper is organized as follows. In Section 2 we briefly describe the general procedure for the decomposition and the use of the experimental results. A more detailed description is given in Section 3; readers who are not interested in technicalities may skip this section altogether. Section 4 discusses the validation of the SModelS framework. In Section 5 we then apply SModelS in two scans of the MSSM and discuss how the SMS limits constrain various particle masses, and which regions of parameter space remain untested. Conclusions and an outlook are given in Section 6.

General Procedure
Currently most ATLAS and CMS experimental analyses consider specific simplified models to present their constraints on new physics. The number of signal events expected in a given signal region is then obtained using the signal efficiency times acceptance for the specific model assumed. In general, the signal efficiency is model dependent and must be calculated for each specific model considered. Nonetheless, most current experimental analyses aim for model independent constraints and consider sufficiently inclusive signal regions. The guiding principle behind the procedure discussed here is the assumption that the signal efficiencies for most experimental searches for new physics depend mostly on the event kinematics and are just marginally affected by the specific details of the BSM model. This allows us to map the full model's signals to its SMS-equivalent topologies and use the latter to constrain the full model.
Clearly, this assumption, which from now on we call SMS assumption, is not always valid. For instance, it is expected to be violated for searches which strongly rely on shape distributions. The signal efficiencies can moreover depend on spin correlations, or on properties of off-shell states in production channels (s-channel or t-channel production) or in decays. A quantitative discussion of the sensitivity of signal efficiencies to various model properties is analysis dependent and a complicated matter, being out of the scope of the current work. Dedicated studies will be presented in future publications. Generally, it is the responsibility of the user to apply SModelS only to models and experimental results for which the SMS assumption is approximately valid. (Note that it is possible to use only a subset of experimental results in the database.) Let us now outline the general SModelS procedure.
Under the SMS assumption, to a first approximation, it is possible to reduce all the properties of a BSM model to its mass spectrum, production cross sections (σ) and decay branching ratios (B). With this knowledge we can decompose the full BSM signal in a series of independent signal topologies with their specific weights given by the corresponding production cross section times branching ratio (σ × B). Such a decomposition is extremely helpful to cast the theoretical predictions of a specific BSM model in a model-independent framework, which can then be compared against the experimental limit on this σ × B. A schematic view of the working principle is given in Fig. 1.
The first step in our procedure is to compute all signal topologies appearing in the full model 1 and their respective weights, σ × B. Since here we only consider models with a Z 2 symmetry (such as R-Parity, T-Parity or KK-Parity), the possible signal topologies will always arise from pair production of new Z 2 -odd particles, which decay as P → P + SM particles, where P and P are the parent and daughter BSM particles, respectively. 2 Hence all topologies will be of the form represented in Fig. 2, which shows the production of the initial pair of BSM states (represented by a circle with two outgoing legs) and their subsequent cascade decays. In our notation all particles appearing in the SMS topology, both Z 2 -even and Z 2 -odd, are on-shell. The case of off-shell decays is always included as 3-body decays, with no mention to the off-shell states. Therefore all the relevant information (in the SMS framework) of such a diagram can be reduced to three main objects: • the diagram topology: number of vertices and SM final state particles in each vertex; • the masses (mass vector) of the Z-odd BSM particles appearing in the diagram; • the diagram weight (σ × B).
The reduction of a particular process to its equivalent SMS topology is illustrated in Fig. 3. Details of the decomposition procedure and the labeling scheme used are explained in Section 3.1. Note that once the decomposition is done, the full model is reduced to its signal topologies and there is no longer any reference to the specific details of the model, except for the relevant Z 2 -odd masses and the σ × B associated to each topology. In this way we can cast the theoretical predictions in a model-independent way. 3 2 Throughout this work we ignore BSM particles which are Z 2 even, such as heavy Higgs bosons in the MSSM. The light MSSM Higgs h is treated as a SM particle. 3 One has to keep in mind, however, that the color factor of the initially produced BSM particles influences the QCD activity in the final state and may thus significantly affect the signal efficiency. This is not a worry in the following as we did not come across any example where constraints from a experimental result assuming QCD production are used to exclude an EW produced topology, or vice-versa, but one might encounter such cases in principle. 4 P 1 P 2 P 3 P 4 P 5 P 6 P 7 P 8 Figure 2: The general type of SMS topology considered in this paper. The P i label the SM final state particles. The end of each decay chain is always the lightest Z 2 -odd particle which is stable. The next and more involved step is to confront the theoretical predictions obtained from the decomposition with the experimental constraints. For that it is necessary to map the signal topologies produced in the decomposition to the SMS topologies constrained by data. For some experimental analyses this is a trivial matter, since they provide an upper limit for a single topology cross section as a function of the relevant BSM mass vector. Examples are constraints on squark pair production, withq → q +χ 0 1 , which give an upper limit on σ × B as a function of (mq, mχ0 1 ), or gluino pair production, withg → tt +χ 0 1 , which limit σ × B as a function of (mg, mχ0 1 ). However it is often the case that the experimental analysis does not constrain a single topology but rather a sum of several topologies, assuming a specific relative contribution from each of them. As an example, consider the slepton pair production limits, where the interpretation constrains the sum over final state lepton flavors (e's and µ's) under the assumption that each flavor contributes 50% to the signal and that selectrons and smuons are mass degenerate, (mẽ, mχ0 1 ) = (mμ, mχ0 1 ). In order to apply this experimental constraint to the signal topologies obtained from the decomposition, it is necessary to combine all topologies with a single lepton being emitted in each branch and which have the same mass vector. Moreover, in order for the experimental constraint to be valid, it is necessary to verify the analysis conditions: topologies with e's and µ's contribute equally to the final theoretical prediction (σ × B). A more involved example are the constraints from tri-lepton+MET searches: here the SMS results forχ ± 1χ 0 2 production with decays through sleptons assume the intermediate slepton being a 5 selectron, smuon or stau (including or not including sneutrinos) so that the limits on σ × B only apply for specific flavor-democratic, tau-enriched, tau-dominated, etc., cases. Such constraints need to be carefully taken into account when mapping the signal topologies obtained from the decomposition to the experimental results. A detailed description of how the analyses assumptions and constraints are described in a model independent language is presented in Sec. 3.2. The procedure for matching the decomposition results to the analysis constraints is discussed in detail in Sec. 3.3. Finally when all signal topologies are combined according to the assumptions of each experimental result, the resulting theoretical predictions for the cross sections of the combined topologies can be directly compared to the experimental upper limits. Thus it can be decided whether a particular parameter point (a particular BSM spectrum) is excluded or not by the available SMS results.

Detailed Description
In this section, we discuss the technical details of the main building blocks of SModelS: the decomposition procedure, the analysis database, and the matching of theoretical and experimental results.

Decomposition Procedure
As explained in the previous section, under the SMS assumption, all the complexity of the BSM model can be replaced by the knowledge of the signal topologies and their weights, together with the relevant BSM masses. 4 The topologies can then be classified according to the number of vertices in each branch and the SM states appearing in each vertex. In order to properly classify the signal topologies, we introduce a formal labeling scheme, which is i) model independent, ii) general enough to describe any topology, iii) sufficiently concise and iv) allows us to easily combine topologies according to the assumptions and conditions in the experimental analyses.
We choose to use a labeling system based on nested brackets, which corresponds to a textual representation of the topology. As discussed in Section 2, we assume that all signal topologies respect a Z 2 -symmetry, hence BSM states are produced in pairs and cascade decay to a single BSM state and SM particles. Therefore a signal topology always contains two branches (one for each of the initially pair produced BSM states), which we describe by [B 1 , B 2 ]. In each branch there is a series of vertices, which represent the cascade decays. From each vertex there is one outgoing BSM state and a number of outgoing SM particles. Thus, inside each branch bracket, we add an inner bracket for each vertex, containing the lists of outgoing SM particles. Note that there is no mention to the intermediate BSM particles (χ 1 , l ,...), which makes our method explicitly model independent. The only information kept from the BSM states are their masses. To illustrate the labeling scheme just defined, we show in Fig. 4 a signal topology containing a 2-step cascade decay in one branch and a one step decay in the other. Following the above prescription, this topology is described by [ Figure 4: The labeling scheme adopted here applied to an example diagram. In parenthesis we show the respective mass vector for the diagram. Given a BSM model -described by its spectrum, branching ratios and production cross sections-we need to obtain all the possible signal topologies and compute their weights (σ × B). The procedure of computing these corresponds to the decomposition of the full model in terms of simplified model topologies. There are two ways to actually perform the SMS decomposition. The first requires the generation of parton-level Monte Carlo (MC) events, followed by the mapping of each event into the corresponding SMS topology. The second method is purely based on the SLHA [20] spectrum file and decay table, supplemented by theoretically computed cross sections. Below we describe both approaches in more detail.

Monte Carlo based decomposition
The decomposition based on parton level MC events is the most general one, since it can decompose any type of BSM model, as long as it is possible to simulate MC events for it. In this case an event file in the LHE format is used as input and each event is mapped to a simplified model topology. The mapping is obviously not one-to-one, since more than one event can generate the same exact SMS topology. Then the sum of MC weights for all events contributing to the same topology directly gives σ × B for the corresponding topology. The disadvantage of this method is the introduction of MC uncertainties in the decomposition result; this can however be easily solved by increasing the MC statistics. We also note that the recent advances on NLO MC generators will allow to produce decomposition results at NLO even in the MC based decomposition.

SLHA-based decomposition
In the SUSY case, there are a number of public codes for the computation of the mass spectrum and decay branching ratios with output in the SLHA [20] format. The production cross sections can then be computed at leading order (LO) through a Monte Carlo generator, or at next-to-leading order (NLO) with Prospino [21]. For gluinos and squarks, cross sections at next-to-leading log (NLL) precision can be computed using NLL-Fast [22]. These cross sections can be included in the SLHA file [23], which then holds all the required information for the SMS decomposition.
In the SLHA-based decomposition the cross sections for pair production of BSM states and the corresponding branching ratios are used to generate all possible signal topologies. The only theoretical uncertainties in this case come from the cross-section uncertainty, hence being much smaller than the uncertainties in the MC-based method, which depends on the MC statistics. In order to avoid dealing with a large number of irrelevant signal topologies, only the ones with σ × B above a minimal cut value are kept. For the results presented below we take this cut value to be 0.03 fb. 5

Compression of topologies
Although the decomposition described here (both for the MC and SLHA based methods) is fairly straightforward, the model independent language introduced above can be used to perform non-trivial operations on the signal topologies. One possibility is the compression of topologies containing a series of invisible decays at the end of the cascade decay chain, as illustrated by Fig. 5(left). In this case the effective final state BSM momentum is given by the sum of the neutrinos and the final BSM state momenta. Therefore, for the experimental analyses, this topology is equivalent to a compressed one, where the effective BSM final state includes the neutrino emissions, as shown in Fig. 5 ]. More importantly, such "invisible compression" can be automatically performed in a model-independent way using the nested bracket notation.
Another case of interest corresponds to topologies which contain spectra with very small mass splittings. In this case one can again approximately map the original topology into a smaller one, where decays between quasi-degenerate states are omitted. This "mass compression" approximation is only reliable when the energy carried away by the SM particles emitted in the decay of quasi-degenerate states is negligible for all experimental purposes. If this is the case, the original topology is equivalent to a compressed one, as shown in Fig. 6. For the Figure 5: A signal topology with invisible decays (left) and its [invisibly] compressed equivalent (right). M k labels the masses of the BSM states appearing in the topology. results presented in Section 5 we perform the mass compression for quasi-degenerate states if their mass difference is below 5 GeV. Once again, this compression procedure allows us to constrain the original topology using experimental constraints for shorter cascade decays. We also point out that, whenever (invisible and/or mass) compressions are performed, care must be taken to avoid double counting topologies. If both the original and the compressed topologies are kept after the decomposition, they should not be combined later, since this would result in a double counting of the topology weight.

Anatomy of an SMS result
The interpretation of the BSM search results in the context of simplified models has become the de-facto standard for the experimental collaborations. The ATLAS and CMS collaborations typically produce two types of SMS results: for each simplified model, values for the product of the experimental acceptance and efficiency (A × ) are determined to translate a number of signal events after cuts into a signal cross section. From this information, a 95% confidence level upper limit (UL) on the product of the cross section and branching fraction (σ × B) is derived as a function of the BSM masses appearing in the SMS. Finally, assuming a theoretical "reference" cross section for each mass combination, an exclusion curve in the plane of two masses is produced.
All limits at 95% CL Figure 7: Examples of SMS results from the CMS (left) and ATLAS (right) collaborations, taken from [24] and [25], respectively. Figure 7 gives two examples. The plot on the left is from the CMS analysis of charginoneutralino production withχ ± 1 →τ ± ν andχ 0 2 →L ± L ∓ , withL ± → L ±χ0 1 , where L = e, µ, τ . Shown are the 95% CL UL on σ × B together with the expected and observed mass limit curves. The plot on the right is an ATLAS result for the case of slepton pair production and direct decay to the LSP,l → lχ 0 1 . Our approach builds upon the binned UL on σ × B, which is the information collected in our analysis database. Neither the efficiency plots nor the exclusion lines are used in our procedure.
The cross section upper limits reported in this way are however subject to a number of assumptions made in the analysis, which also have to be described in the database. In the SModelS language described in Section 3.1, the topologies appearing in Fig. 7 where L = e, µ, τ and the relevant mass vectors are (mχ± 1 , mL, mχ0 1 ) and (mχ0 2 , mL, mχ0 1 ). The experimental analysis in our example assumes degenerate sleptons and mχ0 2 = mχ± 1 , so there are only three independent masses. Furthermore this analysis constraints the sum over lepton flavors and charges for the topologies listed above under the assumption that each term contributes equally (flavor democratic decays). Moreover, the SMS result illustrated in Fig. 7 (left) is for the particular mass relation ml = (mχ± 1 + mχ0 1 )/2 and can thus only be applied for (approximately) these ml and mν values. However, in Refs. [24,25], results for distinct slepton mass values are also provided, which allows us to interpolate between them and apply these results to more generic models. All these constraints and conditions are described in our analysis database in the form of "metadata". As an example we show the metadata for the CMS result in Fig. 7 in Table 1. The entry 'constraint' lists the (sum of) topologies being constrained by the analysis. If charges and/or flavors are not explicitly listed, a sum over charges and/or flavors is assumed. The additional analysis assumptions (each lepton flavor contributes equally to the total weight, σ × B) can actually be relaxed because the efficiency for µ's is higher than the efficiency for e's, which is higher than the one for τ 's. Therefore, instead of requiring an equal weight (σ × B) from each flavor, it suffices to impose: These assumptions are included in the entry 'fuzzycondition' and the function Cgtr(x, y). This function uses the theoretical predictions for x and y to map the condition x > y into a number in the interval 0−1. 6 If Cgtr(x, y) = 0, the condition is fully satisfied (x > y) and if Cgtr(x, y) = 1, the condition is fully violated (x y). If the conditions are strongly violated (Cgtr(x, y) > minimal value) by the input model, the corresponding analysis should not be used to constrain the model. In this way, the user can decide how strictly the conditions are enforced and ignore all constraints which have a too large value of Cgtr(x, y). In the following we ignore all constraints which have Cgtr(x, y) > 0.2.
Finally, the entry 'axes' describes the available slices of the (mχ± 1 , mL, mχ0 1 ) parameter space, which in this example corresponds to ml = 0.5mχ± 1 + 0.5mχ0 1 , ml = 0.95mχ± 1 + 0.05mχ0 1 and ml = 0.05mχ± 1 + 0.95mχ0 1 . This information is used to interpolate between the values in the experimental results. Since different experimental analyses adopt distinct slicing methods of the parameter space, we use a general interpolation procedure which works for any slicing choice. It is based on a tesselation of the mass vector space and a linear interpolation on each simplex.
Another illustrative example is the ATLAS dilepton search [25] for slepton pair production and decay shown in Fig. 7 (right):l + +l − → (l +χ0 1 ) + (l −χ0 1 ), wherel =ẽ,μ. In this case selectrons and smuons are assumed to be mass degenerate and the experimental collaboration constrains the sum of lepton flavors: . Once again, each lepton flavor is assumed to contribute equally. The metadata for this example is given in Table 2. As in the previous example, we assume the muon efficiency to be higher than the 6 The actual function used is Cgtr(x, y) = (|x − y| − (x − y)) / (2(x + y)).

Matching Theoretical and Experimental results
Once a BSM spectrum is decomposed according to the procedure described in Section 3.1, all the relevant information for confronting the model with the experimental results is encapsulated in the SMS topologies plus their mass vectors and weights. Any specific model dependent information can be dropped at this point. As discussed in Section 3.2, it is however often the case that the experimental result constrains a sum of topologies instead of a single one. Before a direct comparison with the experimental constraints, it is necessary to combine single SMS topologies (which means adding their weights) according to the experimental analysis' assumptions (described in the metadata of the analysis). Furthermore, it is always implicitly assumed that all summed topologies have a common BSM mass vector (for the slepton pair production example this means (mẽ, mχ0 1 ) = (mμ, mχ0 1 )). Therefore, when combining signal topologies according to the analysis constraints, we must ensure that they have similar mass vectors.
Since different analysis have different sensitivities to the mass vector, we do not use the simple mass distance in GeV in order to verify if two vectors are similar or not. Instead we consider the sensitivity of the analysis in question to the difference between the two vectors. In order to quantify this sensitivity we use the analysis upper limit for each individual mass vector. If both upper limits differ by less than a maximal amount (20%), we render the mass vectors as similar with respect to this particular analysis. However, if the upper limits differ by more than 20%, we consider the two vectors as distinct. Moreover, in order to avoid cases where two upper limits are coincidentally equal, but they correspond to completely different mass configurations, we also require the mass values not differ by more than 100%.
Once again we illustrate this procedure using the slepton pair production constraint as an example. In the MSSM we usually have mẽ L = mμ L and mẽ R = mμ R , but mẽ L = mẽ R . ]], we must first group the sleptons (ẽ L ,ẽ R ,μ L ,μ R ) with similar masses before we can combine the topologies. In order to identify the similar mass vectors we first obtain, for the given analysis, the upper limit for each vector ((mẽ L , mχ0 1 ), (mμ R , mχ0 1 ),...) and cluster together the masses with similar upper limit values. If the analysis is sensitive to the left-handed/right-handed slepton mass splitting, the upper limits will differ significantly and the grouped masses will correspond to (ẽ L ,μ L ) and (ẽ R ,μ R ). On the other hand, if the mass spliting is small and the analysis is not sensitive to it, all upper limits will be similar and all the sleptons will be grouped together.
After ]]) for the slepton analysis. These conditions need to be taken into account when interpreting the experimental results, since each topology may have a different signal efficiency. Therefore we must also verify that these conditions are satisfied, otherwise the experimental upper limit can not be applied. Finally, if the experimental assumptions are satisfied, the resulting theoretical predictions (σ × B) obtained after combining the topologies can be directly compared to the corresponding experimental upper limit.
In Fig. 1 we summarize the main steps required to confront the BSM model predictions with the experimental constraints: the SMS decomposition, the combination of SMS topologies with identical or similar masses into the topology sums assumed by the analyses and finally the comparison with the experimental upper limits obtained from the database of LHC results described in Section 3.2.

Validation
A simple and robust way to validate the SModelS procedure described in Section 3 consists in applying it to a simplified model. In this case, the experimental assumptions are exactly satisfied by the full (simplified) model and we should be able to reproduce the exclusion curve given by the experimental collaboration. For instance, to validate the ATLAS-CONF-2013-035 χ ± 1χ 0 2 → W Zχ 0 1χ 0 1 analysis, we assume the full model only consists of a pair of mass-degeneratẽ χ ± 1 andχ 0 2 (both pure Winos) and the neutralino LSP (pure Bino),χ 0 1 , with B(χ ± 1 → W ( * )χ0 1 ) = B(χ 0 2 → Z ( * )χ0 1 ) = 1. All the other particles in the spectrum are taken to be in the multi-TeV scale and decoupled. Scanning over mχ± 1 and mχ0 1 and verifying which points are excluded by the ATLAS-CONF-2013-035 analysis, we then expect to reproduce the official exclusion curve obtained by the experimental collaboration. 7 We performed such "regression tests" for several SMS results implemented in the SModelS database. Some illustrative examples are shown in Fig. 8. In particular, we show regression tests for ATLAS and CMS electroweak-ino searches, searches for gluinos with three-body decays into light quarks (ATLAS) or top quarks (CMS), squark pair-production with direct decay into jet+LSP (ATLAS and CMS), slepton-pair production (ATLAS), as well as stop-pair production giving tt+MET (CMS). The upper limit obtained by interpolating the experimental results is shown by the color coded background and we assumed a constant binning of 25 GeV for all plots. The official exclusion curve obtained by the experimental collaboration is shown in black (solid), while the one obtained through SModelS is shown in red (dashed). In general we find very good agreement with the official exclusion curves once we take into account differences in binning size and the smoothing procedure adopted by the experimental collaboration. (In Fig. 8(c), the red line and the background do not extend above the diagonal because the results for off-shell chargino decays are not yet included in the database.) The exception is the ATLAS qq,q → qχ 0 1 analysis [29]: here there are regions below the official exclusion line which are actually not excluded by the 95% CL UL on the cross section. This is due to large fluctuations in the upper limit values published by the experimental collaboration, clearly visible in the colored bins in Fig. 8(e).

Parameter scans
To demonstrate how the SMS results currently constrain supersymmetry, and to see which scenarios remain untested by the current searches (based on their SMS interpretations), we perform a scan over the MSSM with parameters defined at the weak scale. The parameters which we consider are the following: • the gaugino mass parameters M 1 , M 2 and M 3 , assuming the approximate GUT relation M 1 : M 2 : M 3 = 1 : 2 : 6; • the higgsino mass parameter µ, the pseudoscalar mass m A and tan β = v 2 /v 1 ; • a common mass parameter Mq ≡ MQ 1,2 = MŨ 1,2 = MD 1,2 for the first two generations of left-and right-handed squarks; • mass parameters and trilinear couplings for stops and sbottoms: MQ 3 , MŨ 3 , MD 3 , A t , A b • for left and right sleptons, we take common mass parameters ML, MR for all three generations, supplemented by the trilinear couplings for staus, A τ ; Since we do not consider heavy Higgses in our analysis, we fix m A = 2 TeV. This leaves us with 12 free parameters. We perform two random scans over this parameter space: one focusing on gauginos and sleptons with all squarks assumed to be heavy (Scan-I), and one where sleptons are assumed to be heavy but squarks, and in particular stops and sbottoms, can be light (Scan-II). The scan ranges and fixed parameters for the two scans are summarized in Table 3. 8 Table 3: Scan ranges and values for fixed parameters used in this study; dimensionful quantities are in TeV units. R ≡ max(MQ 3 , MŨ 3 ). The parameters are scanned over randomly assuming a flat distribution.
We thus effectively have 7 free parameters in Scan-I and 9 free parameters in Scan-II. We use SOFTSUSY3.3.9 [64] for the computation of the masses and mixings, SDECAY1.3 [65] for the sparticle decay tables, SUPERISO3.3 [66] for flavor observables, and the masslimits function of micrOMEGAs3.2 [67][68][69] for LEP limits (we also compute dark matter observables with micrOMEGAs3.2 but do not use them here). We require that the LSP be theχ 0 1 , and that all heavier sparticles decay promptly. Concretely we require mχ± 1 − mχ0 1 > 300 MeV 125.5 ± 0.2 +0.5 −0.6 GeV (ATLAS) [73] 3 GeV [74,75] 125.5 ± 3 GeV 125.7 ± 0.3 ± 0.3 GeV (CMS) [76] sparticle masses LEP -(micrOMEGAs [67,68]) to veto long-lived charginos, for which the available SMS results do not apply. Finally, we require that points be consistent with measurements of B(b → sγ), B(B s → µ + µ − ), the muon (g − 2), the Z → invisible width, and the mass of the SM-like Higgs boson. We also apply SUSY mass limits from LEP. The constraints imposed are listed in Table 4. In each scan, we collected more than 30K points which obey all these constraints; these points are then passed through SModelS to test their compatibility with the SMS constraints. SModelS decomposes the spectrum of each point into its SMS topologies and compares it to the experimental constraints according to the procedure described in the previous sections. It returns detailed information on the occurring topologies, which topologies are actually tested (i.e. SMS topologies which have been constrained by one or more experimental analyses in the database), and how their σ × B compares to the experimental UL. For the results presented below, we only considered the 8 TeV analyses in the database.

Results from Scan-I (EW-ino and slepton focus)
The purpose of Scan-I is to test the sensitivity to EW-ino and slepton searches. Moreover, as we assume a GUT relation between the gaugino masses, it will also be susceptible to limits from three-body gluino decays. To begin with, we show in Fig. 9 scatter plots of the scan points in the mχ± 1 versus mχ0 1 , M 2 versus µ, mμ 1 versus mχ0 1 and mg versus mχ0 1 planes. The red points, which form the top layer, are excluded by at least one analysis, while the blue points are not excluded by any single SMS limit. The grey points do not have any experimental limits either because their topologies do not match any of the existing experimental results or because their masses fall outside the ranges considered by the experimental searches. These points are labeled 'not tested' and mostly appear when both the chargino and gluino masses fall outside the grids of the experimental results.
Some comments are in order regarding the differences between the areas covered by the red points and the naive expectations from the official exclusion curves. One issue concerns the gaugino-higgsino mixing. As seen in Figs. 9(a) and 9(d), excluded points are concentrated along the mχ± 1 ∼ 2mχ0 1 (or mg ∼ 7mχ0 1 ) line, which corresponds to model points with a pure wino chargino (χ ± 1 W ± ). Onceχ 0 2 andχ ± 1 acquire a higgsino component, the constraints become significantly weaker. Notice that this happens even for points with large mχ± 1 − mχ0 1 mass splittings, far from the pure higgsino region. This is explicitly visible in Fig. 9(b), which shows the projection onto the M 2 -µ plane. As we can see, almost all excluded points lie in the wino chargino region (µ > M 2 ). In Fig. 9(c) we show the same results, but now in theμ 1 -χ 0 1 mass plane, where we see that for mχ0 1 125 GeV there are excluded points for any slepton mass, due to constraints on light charginos (decaying to W Z) and gluinos. Once the LSP becomes heavier, the constraints for Figure 10: Breakdown of the most constraining analyses for each excluded parameter point from Scan-I, on the left for ATLAS and on the right for CMS results. The top row shows thẽ χ ± 1 versusχ 0 1 , the bottom row shows theμ 1 versusχ 0 1 mass planes. The labels correspond to the SMS topology constraining the respective point. If not shown explicitly, charginos, neutralinos, sleptons and sneutrinos are assumed to decay asχ ± 1 → W ±χ0 1 ,χ 0 2 → Zχ 0 1 ,l → lχ 0 1 andν → νχ 0 1 , respectively. The exclusion curves are the same as in Figs. 9(a) and 9(c). χ ± 1χ 0 2 decaying through sleptons become relevant and exclude the region mμ 1 2mχ0 1 up to mχ0 1 ∼ 250 GeV (mχ± 1 ∼ 500 GeV). Finally, in Fig. 9(d), the gluino versus neutralino mass plane, we see once again that the points are concentrated along the bino LSP line (mg ∼ 7mχ0 1 ). It is also interesting in this respect to look at a more detailed breakdown of how different analyses constrain the parameter space. In Fig. 10 we show the most constraining topology for each excluded point in the chargino-neutralino and slepton-neutralino mass planes. 9 We also split the constraints by ATLAS (left panel) and CMS (right panel) results. The first interesting point to notice is that the constraints coming from gluino decay topologies (mostly fromgg production withg → bbχ 0 1 ) are stronger then those from EW production only at low chargino masses (mχ± 1 250 GeV), while the high mass region is constrained by either slepton pair production orχ ± 1χ 0 2 production followed by decay through on-shell sleptons. Both ATLAS and CMS results show a similar behavior. 10 The main differences come from the ATLAS off-shell WZ analysis, which can extend the constraints to a significant part of the mixed-chargino region (at low chargino masses), and from the ATLASll search, which obtained stronger constraints than the equivalent CMS analysis. On the other hand, the 2012 CMS EW-ino analysis [52] saw an under-fluctuation in the BG, which resulted in stronger constraints for some regions of parameter space than the ATLAS results. 11 We note furthermore that points with the same (mχ± 1 , mχ0 1 ) combination may or may not be excluded depending on whether the EW-inos are gaugino-or higgsino-like, or whether their decays proceed through W , Z, h or on-shell sleptons. In order to verify which region of parameter space is excluded for all values of the scanned parameters, we show in Fig. 11 the 'allowed' points (i.e. points that are not excluded by any of the SMS results) on top of the excluded ones in the mχ± 1 versus mχ0 1 and (mg versus mχ0 1 ) planes. We can see that only a very small fraction of points along the pure wino-chargino line is strictly excluded by the SMS results. The region with mg < 500 GeV (mχ0 1 < 70 GeV), excluded for any combinations of parameters, is mainly constrained by limits onχ ± 1χ 0 2 → W * Z * χ0 1χ 0 1 , as seen in Fig. 10(c). For basically any point with mg > 500 GeV there is at least one combination of parameters which evades all SMS constraints. However, the official exclusion curves seem to indicate that gluino  Figure 12: Maximum gluino decay branching ratios in the gluino versus neutralino mass plane for all points in Scan-I for theg →χ 0 1 + tt, bb, qq (top left) andg →χ ± 1 + tb (top right), g →χ 0 1 + g (bottom left) andg →χ 0 2,3,4 + g (bottom right) decays.
masses up to ∼ 1.2 TeV are excluded by SMS constraints on gluino signal topologies. In order to understand how points with light gluinos evade the SMS constraints, we note that the current SMS results for gluinos only consider direct 2-body decays to on-shell squarks or direct 3-body decays to the LSP plus tt, bb or qq. 12 Since for light gluinos (mg < mt 2 TeV) the decay to on-shell squarks is kinematically forbidden, only the SMS constraints on direct decays to the LSP are relevant for Scan-I. Hence the excluded region will be drastically reduced whenever theg →χ 0 1 + X branching ratio (BR) is suppressed. In Fig. 12, top left plot, we show the maximum value of BR(g →χ 0 1 +tt,bb, qq) for each parameter point in the gluino versus neutralino mass plane. As we can see, all points have BR(g →χ 0 1 +tt,bb, qq) < 20% and the region with the highest values and light gluinos (mg < 1 TeV) coincide with the excluded region in Fig. 9(d). Since the exclusion curves always assume 100% BRs, the smaller BRs lead to an exclusion region significantly smaller than the ones projected by the curves. We also show in Fig. 12 the maximum of BR(g →χ ± 1 + tb), which clearly dominates over the other 3-body decays. However, this decay leads to gluino signal topologies not constrained by the current SMS results. Last but not least, Fig. 12 also shows that the loop decay into a neutralino plus a gluon jet can be important. The direct decay into the LSP,g →χ 0 1 + g, is constrained by the 2 jets + MET searches. However, its BR is always small, see the bottom left plot. Thẽ g →χ 0 2,3,4 + g decays are much more relevant, see the bottom right plot, but as theχ 0 2,3,4 decays will lead to additional jets or leptons, these are again not constrained by the current SMS results. Although the constraints from gluino topologies are suppressed due to the small BRs, the same is not expected for theχ ± 1χ 0 2 signal topologies, sinceχ ± 1 ,χ 0 2 almost always decay to on/offshell W 's, Z's, h's or sleptons/sneutrinos and the majority of these topologies are constrained by at least one experimental analysis. However, as shown in Fig. 9(b), most of the excluded points are in the wino chargino region (µ > M 2 ), while the mixed higgsino-wino and pure higgsino chargino regions remains largely unchallenged. This is mainly due to the smaller higgsino production cross sections, which suppressχ ± 1χ 0 2 production once µ M 2 . This is explicitly shown in Fig. 13 where we plot the ratio of theχ ± 1χ 0 2 production cross section and the pure wino cross section in the chargino versus neutralino mass plane. As we can see, once we deviate from the pure wino scenario, σ(χ ± 1χ 0 2 ) decreases, reaching values as low as 20% of the pure wino cross section. Although the official exclusion curves suggest that a large fraction of the mixed higgsino-wino region is excluded, they have been obtained under the assumption of pure wino production cross sections and do not take into account the cross-section suppression due to the higgsino mixing. Furthermore, points with light right-handed sleptons or light staus can also easily evade the SMS constraints. In these cases, the decay of neutralinos and charginos to τ 's is enhanced, while the decay to e's and µ's is suppressed, resulting in a reduction of the signal efficiency for most experimental searches. Also, once mχ0 2 − mχ0 1 125 GeV, the neutralino decay to h+LSP becomes the dominant decay mode. Although there are experimental analyses which look forχ ± 1χ 0 2 → W h +χ 0 1χ 0 1 , these are currently too weak to constrain models with gaugino mass unification.
Finally, the constraints from direct slepton production can also be suppressed when ml R ml L , since right-handed sleptons have smaller production cross sections. Hence we conclude that the parameter space excluded by the current SMS results can be drastically reduced when compared to the naive expectations from the official exclusion curves. For the slepton and EWino signal topologies this reduction is mostly due to the suppression of production cross sections and/or small sensitivity to some signal topologies (such as neutralino decays to Higgs). On the other hand, the constraints on gluino production can be potentially enhanced if new analyses include SMS results for the most relevant gluino topology:g →χ ± 1 + tb →χ 0 1 + W ± + tb. Last but not least note that we do not combine limits from distinct SMS topologies, which corresponds to a conservative estimate of the excluded region. In particular our results tend to underestimate the exclusion obtained from gluino production: as explained above the light gluino scenarios marked as "allowed" in our plots typically show a rather complicated mix of 2 ) and the pure wino production cross section value (σ(W ±W )) in the chargino-neutralino mass plane for Scan-I.
gluino three-body decays intoχ 0 1,...,4 +qq andχ ± 1,2 +qq , and of loop-induces decays intoχ 0 1,...4 +g, such that the SMS limits may be evaded despite a production cross section of the order of 1 pb. At least part of these points are effectively excluded by the generic multi-jet and/or multi-jet plus leptons searches. 13 A way to improve this situation may be the use of efficiency maps, as discussed at the recent LPCC workshop [77].

Results from Scan-II (gluino and squark focus)
Let us now turn to the case that squarks of the 1st/2nd and also of the 3rd generation are allowed to be light. Here, we assume decoupled sleptons. The EW-ino spectrum still is constrained by the assumption of gaugino-mass unification and µ is allowed to vary within 0.1-1 TeV. The full range of parameters is listed in Table 3. In Fig. 14 we show again the excluded points (in red) on top of the SMS-allowed points (in blue), but now in the mχ0 1 versus mg, mχ0 1 versus mq, mχ0 1 versus mt 1 , and mχ0 1 versus mb 1 planes. Since the first two generations of squarks are degenerate, the average squark mass is given by mq ≡ (mũ L + mũ R + md L + md R )/4. For comparison we show moreover the following official exclusion curves: for gluino (gg) topologies, we show the CMS exclusion curves forg → qq +χ 0 1 (Fig. 6b of [44]) andg →tt +χ 0 1 (Fig. 13 of [41]) as well as the ATLAS exclusion forg →bb +χ 0 1 (Fig. 12a of [30]); regarding squark production, we show the exclusion curves forqq → qq +χ 0 1χ 0 1 from both ATLAS and CMS, Figs. 6a of [40] and 19c of [29]; for the 3rd generation, we show the exclusion curves for tt witht → t +χ 0 1 from CMS (Fig. 10a of [63]) andt → b +χ + 1 → b + W + +χ 0 1 from ATLAS ( Fig. 6 of [62]) and the curves forbb → bb +χ 0 1χ 0 1 from both ATLAS and CMS (Fig. 4 of [62] and Fig. 10c of [39]).
From Fig. 14(a) we once again notice that the exclusion region is mostly concentrated in the bino LSP region (µ > M 1 ), despite the potential presence of points with light squarks. Nonetheless the excluded region in the mχ0 1 versus mg plane is considerably larger than in 13 We thank Lukas Vanelderen and Matthias Schroeder for explicit checks against the CMS RA2 results. Scan-I and it extends well beyond the gluino 3-body decays exclusion curves. This is expected since the presence of light squarks allows points to be excluded even for very heavy gluino masses. This is explicitly visible in Fig. 14(b), where we see that excluded points with heavy gluinos (mg > 1.8 TeV or mχ0 1 > 250 GeV) correspond to light squark masses, mq < 1 TeV. It is also interesting to note that the excluded points in the squark versus neutralino mass plane extend well beyond the exclusion curves for squark production and direct decay to the LSP, even for heavy gluino masses. The reason is that even a 2-3 TeV gluino is not yet decoupled but gives a t-channel contribution to squark-pair production. Finally, in Figs. 14(c) and 14(d), the excluded points are projected onto the LSP-stop and LSP-sbottom mass planes. While the density of excluded points seems to depend weakly on the stop mass, it is concentrated in the light sbottom region, agreeing well with the expectation of the exclusion curves, except for light LSP masses, where constraints from gluino and EW-ino topologies become relevant.
In order to better understand the excluded regions shown in Fig. 14, we present in Figs. 15 and 16 the breakdown into the most constraining topologies for each excluded point. We again show the constraints from ATLAS (left panels) and CMS (right panels) analyses separately. As already noted above, in the mχ0 1 versus mg plane, the excluded points are concentrated in the bino LSP region (mχ0 1 ∼ 7mg), see the top row of plots in Fig. 15. The few points in the mχ0 1 < 7mg region are excluded by constraints on sbottom and stop topologies, which are independent of gluino masses. We also point out that the constraints from gluino topologies are only relevant for mg 1.2 TeV, while all the points with heavier gluinos are excluded by squark topologies, as expected. Furthermore, we see that all the points excluded by thẽ qq → qq +χ 0 1χ 0 1 topology falls along the bino LSP region. This is expected since squarks only couple to higgsinos through their small Yukawa couplings. Hence once the LSP acquires a sizeable higgsino component, theq → q +χ 0 1 decay becomes suppressed (benefitting decays into heavier neutralinos or the wino-like chargino), which rapidly weakens the constraints.
In Figs. 15(c) and 15(d) we show the same points, but in the squark versus neutralino mass plane. The points excluded by squarks topologies (qq → qqχ 0 1χ 0 1 ) extend up to mq ∼ 1.4 TeV, well beyond the naive expectations from the exclusion curves. As mentioned, this is due to the t-channel gluino contribution to the squark production cross section, which enhances σ(qq) with respect to the fully decoupled gluino case. Since the exclusion curves assume decoupled gluinos, the squark cross sections are significantly reduced, resulting in a smaller reach. As the gluino mass increases (for a fixed squark mass), the squark cross section slowly decreases as well as the corresponding reach in squark mass. This is explicitly seen in the plots of squark versus gluino mass in Fig. 15.
The same enhancement is not present for the 3rd generation squarks since the t-channel gluino contribution is strongly suppressed by the negligible bottom and top PDFs. This is also shown by the plots in Fig. 15, where we see that the reach from sbottom and stop topologies are independent of the gluino mass up to mg ≈ 1.5 − 1.7 TeV. For even higher gluino masses, points with a bino LSP have mχ0 1 > 250 GeV and are beyond the reach of the sbottom/stop SMS constraints, as shown by the exclusion curves in Figs. 14(c) and 14(d). The few points above mg ∼ 1.5 − 1.7 TeV excluded by the sbottom and stop signal topologies have a light higgsino LSP.
Finally, in Fig. 16 we show the analyses breakdown in the stop and sbottom versus neutralino mass planes. We see that most of the points excluded by the stop or sbottom signal topologies agree well with the expectations from the exclusion curves. One exception are the few points at mt 650 GeV and mχ0 1 ≈ 150 − 180 GeV in the top-right plot in Fig. 16. These are excluded by the constraints from the CMS razor analysis [41], which in fact has a higher reach than the exclusion curve at low LSP masses (mχ0 1 180 GeV). In the mχ0 1 versus mb 1 plots in Fig. 16 (bottom row) the points excluded by the sbottom topologies are slightly above the exclusion curves due to uncertainties in the production cross section, but they agree within 1σ. The labels correspond to the SMS topology constraining the respective point. If not shown explicitly, charginos, neutralinos and stops are assumed to decay asχ ± 1 → W ±χ0 1 ,χ 0 2 → Zχ 0 1 andt → tχ 0 1 , respectively. The exclusion curves are the same as in the corresponding plots in Fig. 14. 27 Figure 16: Breakdown of the most constraining analyses for each excluded parameter point from Scan-II as in Fig. 15 but in theχ 0 1 versust 1 (top row) andχ 0 1 versusb 1 (bottom row) mass planes. The labels correspond to the SMS topology constraining the respective point. If not shown explicitly, charginos, neutralinos and stops are assumed to decay asχ ± 1 → W ±χ0 1 , χ 0 2 → Zχ 0 1 andt → tχ 0 1 , respectively. The exclusion curves are the same as in the corresponding plots in Fig. 14. Figure 17: SMS-allowed points from Scan-II in the squark versus gluino mass plane for the bino LSP case (µ > 2M 1 ) and the mixed higgsino-bino scenario (µ < 2M 1 ). The dashed lines correspond to mg = mq.
We now turn our attention to the points which have large cross sections (light squark and/or gluino masses) but none the less evade all current SMS constraints. In Fig. 17 we show the allowed points on top of the excluded ones in the squark versus gluino mass plane. It is instructive to consider separately the bino and mixed LSP scenarios, so in the left panel we plot the points with mostly binoχ 0 1 and winoχ ± 1 (µ > 2M 1 ), while in the right panel we show points with a mixed higgsino-bino LSP (µ < 2M 1 ). In the µ > 2M 1 case we see that gluino masses below 500 GeV are excluded for almost any choice of the parameters, in agreement with the results for Scan-I. As already mentioned, because of the relation between M 1 , M 2 , M 3 , low gluino masses are also constrained by limits onχ ± 1χ 0 2 → W Z +χ 0 1χ 0 1 . However for any higher values of mg there are several allowed points as long as mg < mq and/or µ < 2M 1 . The former case was already discussed in the previous Section, since it is equivalent to the parameter space of Scan-I. As mentioned, if gluino decays to squarks are kinematically forbidden, thẽ g → tb +χ ± 1 decay often dominates, which is not constrained by the current SMS results. Furthermore, if squarks are heavier than gluinos, they predominatly decay tog + q and these signal topologies are also not constrained by the current SMS results. Therefore points with mg < mq and mg < mt ,b can always evade the SMS constraints. On the other hand, if mg > mq, squarks decay to EW-inos and the constraints on squark production can be applied. From the last plot in Fig. 17 we see that points with light squarks (mq 1.3 TeV) and mg > mq are always excluded in the bino LSP scenario (µ > 2M 1 ). Once the LSP acquires a sizeable higgsino component, theq → q +χ 0 1 decay becomes suppressed and the constraints are much weaker, resulting in several allowed points in the light squark/gluino region, shown in the right plot in Fig. 17. 29

Conclusions
We presented a new tool, SModelS, to decompose the LHC signatures expected from BSM spectra presenting a Z 2 symmetry into SMS topologies and test the predicted cross sections (σ × B) for each topology against the existing 95% CL upper limits from ATLAS and CMS. The program consists of three parts, the decomposition procedure, which can be Monte-Carlo based or SLHA based, the database of ATLAS and CMS SMS results, and the interface to match model predictions onto the experimental results. Our concrete implementation currently focusses on SUSY searches with missing energy, for which a large variety of SMS results from ATLAS and CMS are available. SModelS can be used "out of the box" for the MSSM and its extensions, like the next-to-MSSM. The approach is however perfectly general and easily extendible to any BSM model to which the experimental SMS results apply. (Of course, the SMS assumption is subject to several caveats, as explained at the beginning of Section 2 and in footnotes 3 and 4.) As a proof of principle, we applied our procedure in two scans of the weak-scale MSSM and discussed which parameter regions are excluded by the SMS limits and which scenarios remain untested. As we showed, SModelS can be useful for several purposes. First of all it is a convenient tool for relatively fast surveys of complex parameter spaces in order to devise the regions that are definitely excluded by data. Second, it can be used to find relevant topologies for which no SMS results exists-this may be helpful for the experimental collaborations for designing new SMS interpretations. An example that we already identified is "mixed" topologies with tb final states, which can arise e.g. fromg → tbχ ± 1 or fromt 1 pair production with one stop decaying into tχ 0 1 and the other one into bχ ± 1 . Third, one may use the decomposition procedure of SModelS independent of the results database in order to identify the most important signal topologies in regions of parameter space. Finally, in the case of a positive BSM signal, the SModelS framework may be helpful for characterizing which new physics scenarios may explain the observations. While SModelS is already a powerful and useful tool for phenomenological studies, there is of course still much room for improvement. For instance, as we have shown, scenarios with complicated decay patterns-as typical for e.g. light gluinos and heavy squarks-are not well constrained with the current framework, which is based on testing the σ × B upper limits topology-by-topology without the possibility of combining results. Here the use of efficiency maps, once available in a systematic fashion, would allow for significant improvements. This is clearly a development which we will follow. We also foresee to extend SModelS to multiple branches, including sub-branches originating from the main cascade decay chains. This is relevant, e.g., for heavy Higgses (or other non-SM "R-even" particles as present in models with extra dimensions) appearing in the decay chains. We also conceive the inclusion of resonant production of new particles, as well as violation of the Z 2 symmetry; in the context of SUSY this means extension to R-parity violation. Regarding the statistical treatment, for identifying the most sensitive analysis (in order not to alter the effective CL), it would be of great help if the experimental collaborations systematically provided also the expected σ × B upper limits in addition to the observed ones.
Last but not least note that the mapping from the full model signal to a sum over simplified model signal topologies is clearly not an exact procedure. It assumes that, for most experimental searches, the BSM model can be approximated by a sum over effective simplified models.
The validity of this approximation, in particular (i) regarding the question how the type of production channel or the nature of the off-shell states mediating the decays affects the signal efficiencies, and (ii) regarding the question of using SMS limits derived in the SUSY context for non-SUSY models, will be a subject of future work.
The publication of the SModelS code with a dedicated online manual is in preparation [18].

Note added
Shortly before submission of this paper another program package, CheckMATE [78], became available for confronting BSM scenarios with LHC Data. CheckMATE determines whether a model is excluded or not at 95% C.L. by comparing to many recent experimental analyses, based on fast simulation. It thus represents an interesting alternative to the SMS approach followed by SModelS.
Moreover, while the SModelS paper was being refereed, another program based on SMS results, Fastlim [79], was published. Fastlim employs efficiency maps and currently takes into account 11 ATLAS analyses, which are mainly focused on stop and sbottom searches. Restricting the SModelS database to the 11 ATLAS analyses implemented in Fastlim, we have verified that for Scan-II about 70% of the points excluded by SModelS are also excluded by Fastlim and vice-versa. Interestingly, from all the points excluded by SModelS, ≈ 30% are not excluded by Fastlim; a part of these points features a light LSP with mass of 45-65 GeV -this region does not seem to be covered by Fastlim. The reverse is also true: from all the points excluded by Fastlim, ≈ 30% are not excluded by SModelS. A detailed comparison of the two approaches is on the way. We recall, however, that SModelS includes in total more than 50 analyses from both ATLAS and CMS, which allows to exclude an additional 54% of points in Scan-II.