Determining the parameters of the Minimal Supergravity Model from 2l + E_T^{miss} + (jets) final states at LHC

We analyse the events with two same-flavour, opposite-sign leptons + $E_T^{miss}$ + (jets) as expected in pp collisions at LHC within the framework of the minimal Supergravity Model. The objective is the determination of the parameters m_0 and m_{1/2} of this model (for a given value of $\tan\beta$). The signature $l^+ l^- + E_T^{miss}$ + (jets) selects the leptonic decays of $\tilde{\chi}^0_2$, $\tilde{\chi}^0_2 \to \tilde{\chi}^0_1 l^+ l^- $, $\tilde{\chi}^0_2 \to \tilde{l}_{L,R}^{\pm} l^{\mp} \to \tilde{\chi}^0_1 l^+ l^-$. We exploit the fact that the invariant dilepton mass distribution has a pronounced structure with a sharp edge at the kinematical endpoint even in such an inclusive final state over a significant part of parameter space. We determine the domain of parameter space where the edge is expected to be visible. We show that a measurement of this edge already constrains the model parameters essentially to three lines in the ($m_0, m_{1/2}$) parameter plane. We work out a strategy to discriminate between the three-body leptonic decays of $\tilde{\chi}^0_2$ and the decays into sleptons $\tilde{l}_{L,R}$. This procedure may make it possible to get information on SUSY particle masses already with low luminosity, L_{int}=10^3 pb^{-1}.


Introduction
If 'low-energy' supersymmetry (SUSY) is realised in Nature it should show up at the Large Hadron Collider (LHC). Strongly interacting particles as gluinos and squarks will be most likely the first SUSY particles to be seen at LHC. Gluinos of mass less than ∼ 2 TeV and squarks of mass less than ∼ 1.5 TeV [1][2][3] can be detected, covering in such a way the whole theoretically motivated parameter space. LHC is also a good laboratory for the search of electro-weakly interacting particles, e.g. sleptons [4,5]. In a recent paper [6] it was shown within the minimal supergravity (mSUGRA) [7] model that sleptons in the mass range of ∼ 100 to 400 GeV can be detected at LHC by investigating the signature two leptons + E miss T + no jets. However, this final state where direct production (Drell-Yan) of sleptons predominates requires high luminosity, L int = 10 5 pb −1 . Sleptons can be also produced indirectly in the decays of charginos and neutralinos, especially inχ 0 2 →l L,R l decays. The charginos and neutralinos, can in turn be produced directly or come from gluinos and/or squarks. This leads to final states with ≥ 2 leptons + E miss T + (jets). Actually, this indirect slepton production throughg,q decays has the largest cross-section in a sizable region of the parameter space accessible at LHC and could allow sleptons to be already revealed at L int = 10 3 pb −1 , i.e. simultaneously with strongly interacting sparticles. Thus, indirect production of sleptons can be more important for observing a slepton signal than direct one [8]. Moreover, in such a way the mass reach for sleptons search can be extended up to ml L ∼ 740 GeV. Having evidence for SUSY at LHC, one of the next tasks will be to find out the underlying model and to determine the model parameters. In this paper we work out a method to determine SUSY parameters and suggest a strategy for getting information on masses of SUSY particles by means of the signature l + l − + E miss T + (jets). Our study is made within the framework of the minimal supergravity model (mSUGRA) [7]. In this model all scalar particles (sfermions and Higgs bosons) have a common mass m 0 at M GUT ≈ 10 16 GeV. The gaugino masses M 1 , M 2 , M 3 (corresponding to U(1), SU (2), and SU (3), respectively) unify to a common gaugino mass m 1/2 , and all trilinear coupling parameters A ijk have the same value A 0 at M GUT . One also has unification of the electroweak and strong coupling parameters α i , i = 1, 2, 3 [9]. A further reduction of the parameters is given by invoking 'radiative symmetry breaking'. As a consequence, one has only the following input parameters: m 0 , m 1/2 , A 0 , tan β, sign(µ). Here tan β = v 2 v 1 , the ratio of the two vacuum expectation values of the two Higgs doublets, and µ is the Higgsino mass parameter. The whole SUSY particle spectrum can then be calculated by making use of renormalization group equations (RGE). This model is also incorporated in the Monte-Carlo generator ISASUSY [10] which is used in our analysis. This paper is aimed at determining the parameters m 0 and m 1/2 with fixed tanβ. Knowing these parameters, we can calculate the masses of the superpartners using RGE. For this purpose, we study in detail the leptonic decays ofχ 0 2 which have some useful features. Within the mSUGRA model,χ 0 2 has two-body decays,χ 0 2 →l ± L,R l ∓ , in the region m 0 < ∼ 0.5 · m 1/2 of the parameter space, whereas in the region m 0 > ∼ 0.5 · m 1/2 , m 1/2 < ∼ 200 GeV theχ 0 2 has three-body decays,χ 0 2 → l + l −χ0 1 . In both regions the invariant dilepton mass spectrum (M l + l − ) has a maximum M max l + l − , and therefore a pronounced structure with a sharp edge can be seen at the kinematical endpoint. This property was discussed first in ref. [11] in the case of three-body decays ofχ 0 2 and then in ref. [12] in the case of two-body decays. The generality of this feature, i.e. the observability of an edge in the M l + l − spectrum even in inclusive l ± l ∓ l +− and l + l − + E miss T final states in a large part of the parameter space was shown in [13]. We will show how much the parameters are constrained by a measurement of the M max l + l − value of the dilepton mass spectra. Moreover, we will discuss a method, based on the analysis of the M l + l − spectrum, to find out whether the observed edge is due to the two-body or three-body decays ofχ 0 2 .

Sparticle masses in mSUGRA
Within the Minimal Supersymmetric Standard Model (MSSM) the masses of the neutralinos are determined by the parameters M = m 1/2 (M Z ) , µ, and tanβ using M 1 ≃ 5 3 tan 2 θ W ≃ 0.5M (M 1 being the U(1) gaugino mass). In the following, we fix tanβ = 2 (we assume that tanβ could be known from previous experiments) and take A 0 = 0. In mSUGRA | µ | quite generally turns out to be | µ |> M, so that in this case In supergravity the slepton masses are given by [14]: Analogous equations exist for squarks [14]. Therefore, when the parameters m 0 , m 1/2 and tan β are known we can calculate all sparticle masses. A special case is the third generation of squarks and sleptons, where L−R mixing plays a crucial rôle.
3 Production and leptonic decay ofχ 0 2 Neutralinosχ 0 2 can be produced at the LHC through a Drell-Yan mechanism (direct production), in association with strongly interacting sparticles, or in the decay chain of gluinos and squarks (indirect production). Gluino and squark pair production processes are the dominant source ofχ 0 2 's because of large strong interaction cross-sections. The branching ratios of gluino and squark decays intoχ 0 2 are also sizable and are shown in fig.1 as a function of the model parameters m 0 and m 1/2 . One can see that the regions in (m 0 ,m 1/2 ) plane where the decaysg →χ 0 2 + X andq →χ 0 2 + X are open, are complementary. In the region m 0 > ∼ 1.47 · m 1/2 gluinos are lighter than squarks and can decay intoχ 0 2 , while squarks prefer to decay into gluinos. In the region m 0 < ∼ 1.47 · m 1/2 squarks are lighter than gluinos and can decay intoχ 0 2 , but gluinos decay into squarks. Hence the decaysg →χ 0 2 + X (q →χ 0 2 + X) andq →g →χ 0 2 + X (g →q →χ 0 2 + X) can coexist (see also fig.1). Fig.2(a,b) shows σ × Br for indirectχ 0 2 production from gluinos and squarks as a function of m 0 and m 1/2 .
2 Branching ratios ofχ 0 2 decays into leptons directly and via sleptons are shown in fig.3 as a function of m 0 and m 1/2 . The regions where these decays are kinematically allowed are complementary in the parameter plane, depending on whetherχ 0 2 is lighter or heavier thanl L,R . Thus, one can distinguish three domains in the (m 0 ,m 1/2 ) plane, which are (also see fig.4): In domain III, the decayχ 0 2 →l R l would also be kinematically allowed, but since the B − ino component ofχ 0 2 is very small, the coupling tol R l is also small. Therefore, the decayχ 0 2 →l R l is very much suppressed in the whole domain. In fig.5 we show the regions for σ × Br(χ 0 .002 pb in the (m 0 , m 1/2 ) plane from indirect and associatedχ 0 2 production followed by decays toχ 0 1 l + l − final states directly or via sleptons. One can see that there are regions in domains II and III where the mentioned decays coexist. Finally, in fig.6 we show σ × Br for indirect and associated production ofχ 0 2 decaying into leptons directly or via sleptons.

Determination of m 0 , m 1/2 and sparticle masses
In order to determine the parameters m 0 and m 1/2 we will exploit in the following the kinematical features of the two-and three-body leptonic decays ofχ 0 2 . As pointed out in [11], the decayχ 0 2 → l + l −χ0 1 (domain I) has the useful kinematical property that the invariant mass of the two leptons M l + l − has a maximum at whereas for the decaysχ 0 2 →l ± L,R l ∓ → l + l −χ0 1 (domain II and III) the maximum of M l + l − is given by [12]: Thus, the M l + l − distribution has a very characteristic shape with a sharp edge at the kinematical endpoint M max l + l − . As the main source ofχ 0 2 's is their indirect production in gluino and squark decays, the most suitable signature for selecting theχ 0 2 decays is provided by the topology with two same-flavour opposite-sign leptons accompanied by large missing transverse energy and usually accompanied by a high multiplicity of jets. In this paper we thus concentrate on the two same-f lavour, opposite-sign leptons + E miss T + (jets) channel, where the final state leptons are electrons and muons. 3

CMS detector simulation
The simulations are done at the particle level, with parametrised detector responses based on detailed detector simulations. These parametrisations are adequate for the level of detector properties we want to investigate, and are the only practical ones in view of the multiplicity and complexity of the final state signal and background channels investigated. The essential ingredients for the investigation of SUSY channels are the response to jets, E miss T , the lepton identification and isolation capabilities of the detector, and the capability to tag b-jets.
The CMS detector simulation program CMSJET 3.2 [15] is used. It incorporates the full electro-magnetic (ECAL) and hadronic (HCAL) calorimeter granularity, and includes main calorimeter system cracks in rapidity and azimuth. The energy resolutions for muons, electrons (photons), hadrons and jets are parametrised. Transverse and longitudinal shower profiles are also included through appropriate parametrisations. The main detector features incorporated in the Monte-Carlo description are: • Hadronic tracks, muons and electrons are measured up to | η |=2.4 • Deflection of charged particles due to the 4 T magnetic field is included.
• The resolution for the muon system is parametrised according to [16].
• The calorimetric coverage goes up to at | η |=5 for the HCAL and | η |=2.6 for the ECAL.
• ECAL energy resolution parametrized as: • HCAL energy resolution is parametrised according to [17] as a function of η; a typical hadron resolution is: • Energy resolution for very forward calorimeter (VFCAL), in the parallel plate chambers option: • Granularity of calorimeters: • E threshold on cells: • a modified UA1-jet finding algorithm with a cone size of ∆R=0.9 (for description see CMSJET 3.5 [15]) is used for jet reconstruction.

Observability of edges in invariant dilepton mass distributions
In this chapter we determine the regions in the (m 0 , m 1/2 ) parameter plane, where the characteristic edge in the M l + l − distribution can be observed in inclusive final states with two same-f lavour, opposite-sign leptons + E miss T + (jets) with different luminosities at LHC.
Standard Model background processes are generated with PYTHIA 5.7 [18]. We use CTEQ2L structure functions. The largest background is due to tt production, with both W 's decaying into leptons, or one of the leptons from a W decay and the other from the b-decay of the same t-quark. We also considered other SM backgrounds: W + jets, WW, WZ, bb and τ τ -pair production, with decays into electrons and muons. Chargino pair productionχ ± 1χ ∓ 1 is the largest SUSY background but gives a small contribution compared to the signal.
To observe an edge in the M l + l − distributions with the statistics provided by an integrated luminosity L int = 10 3 pb −1 in a significant part of the (m 0 , m 1/2 ) parameter plane, it is enough to require two hard isolated leptons (p  Fig.7 shows the invariant mass spectra of the two leptons at various (m 0 , m 1/2 ) points from domains I, II and III, respectively. The observability of the "edge" varies from 77σ and signal to background ratio 31 at point (200,160) to 27σ and a signal to background ratio 2.3 at point (60,230). The appearance of the edges in the distributions is sufficiently pronounced already with L int = 10 3 pb −1 in a significant part of (m 0 , m 1/2 ) parameter plane, see fig.9. The edge position can be measured with a precision of ∼ 0.5 GeV.
With increasing m 0 and m 1/2 cross-sections are decreasing, therefore higher luminosity and harder cuts are needed. To achieve maximal reach in m 1/2 with L int = 10 4 pb −1 for points from domain III, a cut up to E miss T > 300 GeV is necessary to suppress the background sufficiently. For points with large m 0 (domain I) the transverse momentum p T of the leptons and E miss T are not very large, but there are more hard jets due to gluino and squark decays. Thus for these points we keep the same cuts for leptons and missing energy as before (p l 1,2 T > 15 GeV, E miss T > 100 GeV) and require in addition a jet multiplicity N jet ≥ 3, with energy E jet T > 100 GeV, in the rapidity range | η jet |< 3.5. To optimise the edge visibility we also apply an azimuthal angle cut, ∆φ(l + l − ) < 120 0 . For points from domain II, the jet multiplicity requirement is also helpful. Right sleptons are too light to provide large lepton p T and E miss T , and to use cuts on p l T and E miss T alone is not very advantageous. With L int = 10 5 pb −1 , to suppress the background at larger accessible m 0 , m 1/2 values, we have to require at least 2 or 3 jets, depending on the m 0 , m 1/2 region to be explored. Fig.8 shows invariant dilepton mass distributions at some (m 0 , m 1/2 ) points close to maximum reach with L int = 10 4 pb −1 and L int = 10 5 pb −1 respectively.
The regions of the (m 0 , m 1/2 ) parameter plane where an edge in the M l + l − spectra can be observed at different luminosities are shown in fig.9. In fig.10 we show separately the three domains where an edge due toχ 0 2 → llχ 0 1 ,l R l andl L l decays can be observed at L int = 10 3 pb −1 . One can notice a small overlapping region, where we expect to observe two edges, due toχ 0 2 → l + l −χ0 1 and toχ 0 2 →l ± R l ∓ → l + l −χ0 1 decays (case 1). With increasing luminosity and correspondingly higher statistics, this overlapping region increases, see figs.11 and 12. These plots show the same as fig.10, but for L int = 10 4 pb −1 and L int = 10 5 pb −1 , respectively. An additional region appears where two edges can be observed simultaneously, due toχ 0 2 →l ± R l ∓ → l + l −χ0 1 andχ 0 2 →l ± L l ∓ → l + l −χ0 1 decays (case 2). These regions (case 1 and 2) are due to the coexistence of differentχ 0 2 decay modes has been seen in fig.5. An example of a M l + l − distribution for case 1 is shown in fig.13. Therefore, to a given integrated luminosity at LHC (L int = 10 3 pb −1 to 10 5 pb −1 ) there corresponds a definite parameter region where the characteristic structure in the M l + l − distribution can be seen. This fact already gives a preliminary information about the parameters m 0 and m 1/2 . The observation of two edges would give even stronger constraints.

M max
l + l − analysis of (m 0 , m 1/2 ) parameter plane More specifically, at low luminosity, L int = 10 3 pb −1 , at the beginning of the LHC operation, the accessible values of M max l + l − lie in the following ranges (see figs. 10 and 14): It follows from the discussion above that a measurement of M max l + l − in the dilepton mass distribution, with a single edge, constrains the parameters in general to three lines in the 6 (m 0 , m 1/2 ) parameter plane. In case of M max l + l − > ∼ 90 GeV the constraint is stronger, there are just two possible lines. The most favourable case is when the measured M max l + l − value is large, M max l + l − > ∼ 180 GeV. Then one is left with a single line in the (m 0 , m 1/2 ) parameter plane. For the present study we have chosen, as an example of the general situation, the case of M max l + l − = 74 ± 1 GeV, with three lines corresponding to the domains I,II and III, respectively. The next step is to find out which line in the (m 0 , m 1/2 ) plane is the right one. To this purpose we have analysed points along these lines, given in tables 1-3.
The study is made for the low luminosity case, L int = 10 3 pb −1 .  One should first notice that the observation of two edges at L int = 10 3 pb −1 would determine the (m 0 , m 1/2 ) point uniquely. This is due to the fact that the set of the edge position values in the M l + l − spectrum is different at each point of the parameter region (case 1), where two edges are expected to be observed at 10 3 pb −1 , see figs.10 and 14. With a luminosity L int = 10 4 pb −1 the positions of the two edges will fix two (m 0 , m 1/2 ) points, belonging to domains II and III and corresponding to case 1 and case 2, respectively. At high luminosity, L int = 10 5 pb −1 , the observation of two edges can give up to three possible (m 0 , m 1/2 ) points. One of them is from domain II, corresponding to case 1. The lines M max l + l − = const corresponding toχ 0 2 →l R l decays have the form of an ellipse and can cross the M max l + l − = const lines corresponding toχ 0 2 →l L l decays twice. Hence, two points with the same set of edge positions in the M l + l − spectrum can be found in domain III corresponding to case 2. A discrimination between these points is possible on basis of the event kinematics, and/or by an analysis of the total event rate and the relative number of events corresponding to the two peaks, (see figs.5,10-12).

Discrimination between differentχ 0 2 leptonic decays
For a definite value of the edge position M max l + l − one expects a different shape of the M l + l − distributions in two-and three-body decays (see fig.13, where the first peak is due to a two-body decay ofχ 0 2 and the second one due to a three-body decay). As we have seen from figs.7 and 8 the signal events contribute in the interval 0 < ∼ M l + l − < ∼ M max l + l − . In the following we only consider events in this mass region. The average value <M l + l − > of signal and background events of this mass region is shown in fig.15 as a function of m 0 . The errors are calculated by taking into account the statistical error, a systematic error in the measurement of the edge position, and a systematic error of 30 % for background uncertainty (the main background is tt). One clearly sees that <M l + l − > is significantly smaller in the case of a direct three-body decayχ 0 2 → llχ 0 1 . Thus the shape of the dilepton mass spectrum already allows one to decide whetherχ 0 2 decays into a slepton or not. In order to distinguish between domains II and III, we suggest to use the fact that in general the contour lines with the same M max l + l − for right and left sleptons have no overlap in m 1/2 in regions of parameter space which are accessible at a given luminosity, see figs.9 and 14 as an example. It means that the masses ofχ 0 1 's are different for these two lines, and hence E miss T is expected to be different. In fig.16 we show the <E miss T > values for events withχ 0 2 →l ± L,R l ∓ → l + l −χ0 1 decays after the cuts p l 1,2 T > 15 GeV and E miss T > 100 GeV, M l + l − < M max l + l − . The errors are calculated by taking into account the statistical error and a systematic error in the measurement of the edge position. As can be seen from fig.16, <E miss T > is larger in the case ofχ 0 2 →l ± L l ∓ → l + l −χ0 1 than in the case of χ 0 2 →l ± R l ∓ → l + l −χ0 1 as expected.

Event rate analysis
When the correct M max l + l − line is chosen, the last step is to find the point (m 0 , m 1/2 ) on this line. In general the cross section falls with increasing m 1/2 and m 0 . Thus, we study the event rate along the corresponding M max l + l − line. We first discuss the domain III, where the situation is simpler. For the event rate analysis at L int = 10 3 pb −1 , to reduce the uncertainties due to background, we use a harder cut on E miss T , E miss T > 130 GeV. The dependence of the expected event rate on m 0 is shown in fig.17. The errors are calculated by taking into account the statistical error and a systematic error of 30 % for background uncertainty. A systematic error due to the precision of the edge position measurement is also taken into account. From the observed event rate we can then determine m 0 with a good accuracy, δm 0 ≃ 4 GeV. The parameter m 1/2 is then given by the M max l + l − −line in the (m 0 , m 1/2 ) plane. The precision obtained in such a way is δm 1/2 ≃ 4 GeV.
In domain II, the event rate along a line of definite M max l + l − is first increasing and then decreasing with m 0 , see fig.18. This is mainly due to the change in the branching ratios (see fig.4). The dependence of the event rate is, however, such that m 0 cannot be determined uniquely. To a given event rate there correspond in general two m 0 values. The ambiguity can, however, be solved at high luminosity L int = 10 5 pb −1 , when two edges in the M l + l − distribution can be observed.
For domain I, the m 0 dependence of the event rate is shown in fig.19a, again for M max l + l − ≃ 74 ± 1 GeV. Notice the steep increase of the rate at m 0 ≃ 120 − 130 GeV. This is due to the fact that the decay channelχ 0 2 → llχ 0 1 is just opening in this region. As can be seen from the curve in fig.19a, there is an ambiguity in the determination of m 0 if the event rate is in the region 3700 < ∼ N EV < ∼ 5600 or 120 GeV < ∼ m 0 < ∼ 240 GeV. Here it helps if we look at the average number of jets <N jet > in the events under study. Fig.19b shows <N jet > as a function of m 0 . <N jet > is rising with m 0 as more jets are produced as the squarks become heavier. With the measured <N jet > we can resolve the ambiguity in the mentioned region 120 GeV < ∼ m 0 < ∼ 240 GeV and thus determine m 0 with δm 0 ≃ 7 − 3 GeV.

Conclusions
In this paper we have performed a detailed analysis of events with the signature l + l − + E miss T + (jets) to be expected in pp collisions at LHC. Our aim has been to determine the parameters m 0 and m 1/2 of the Minimal Supergravity Model and to get information on the mass spectrum of SUSY particles, assuming knowledge of tanβ from previous experiments. We have exploited the property of theχ 0 2 leptonic decaysχ 0 2 → l + l −χ0 1 , χ 0 2 →l L,R l → l + l −χ0 1 that the invariant mass of the two final leptons has a maximum, M max l + l − , clearly visible even in inclusive production. We have determined for different luminosities the regions in the (m 0 , m 1/2 ) parameter plane where one or two edges can be observed in the invariant dilepton mass distributions. These regions already give preliminary information about the model parameters. The appearance of the edges in the M l + l − distributions can be already seen with a luminosity L int = 10 3 pb −1 . Therefore we have concentrated on a low luminosity study. On the other hand, in case no such observation will be made at this luminosity, the corresponding parameter region can be excluded, and the same analysis can be done at higher luminosity.
We have shown that a measurement of the M max l + l − value constrains the parameters mainly to three lines in the (m 0 , m 1/2 ) parameter plane. The lines correspond to the decay modesχ 0 2 → l + l −χ0 1 ,l L l → l + l −χ0 1 ,l R l → l + l −χ0 1 respectively. We have worked out a method to discriminate the three-body from the two-bodyχ 0 2 decays. In the case of three-bodyχ 0 2 decays the parameter m 1/2 can be determined by the measured value of M max l + l − with a precision of ∼ 0.5 GeV. The parameter m 0 can then be determined from the observed event rate with a precision of 7-3 GeV. In the case of two-bodyχ 0 2 decays, a measurement of the missing transverse energy can allow one to distinguish between the two possible decaysχ 0 2 →l L l andχ 0 2 →l R l, but a more detailed study is needed. By an event rate analysis along the corresponding line in the (m 0 , m 1/2 ) plane we can determine m 0 and m 1/2 , δm 0 ∼ δm 1/2 ∼ 4 GeV.
Knowing m 0 , m 1/2 and tanβ, the masses of all SUSY particles (except for the 3rd generation of squarks and sleptons) are calculated by RGE. The precisions which can be achieved are ∼ 1 − 6 GeV.
In such a way it is possible to obtain information about SUSY particle masses already with low luminosity (L = 10 3 pb −1 ) even without having direct experimental evidence for their existence. This is especially important for sleptons in a parameter region where high luminosity would be necessary to detect them through direct production.
This study has been performed for tanβ = 2, but it is also possible for high values of tanβ. Most likely, for large tanβ > ∼ 30 a higher luminosity will be needed because of smaller branching ratios of theχ 0 2 leptonic decays. Let us mention some further interesting aspects of this work. Selecting the twobodyχ 0 2 leptonic decays by our method represents an indirect evidence for sleptons in the framework of mSUGRA. In this way it is possible to probe slepton masses up to ∼ 740 GeV well beyond what is possible in direct [5,6] searches. As it has been shown in this study, the edge in the invariant dilepton mass distributions is expected to appear at M max l + l − > ∼ 10 GeV, being quite generally a signal for a two-or three-body decay of some abundantly produced heavy object. Hence such an observation may serve as a first evidence for physics Beyond the Standard Model, and if observed with significant E miss T it would be a clear evidence for SUSY, more specifically forχ 0 2 production. Figure 1: Decay branching ratios as a function of m 0 and m 1/2 (in GeV) for: a)g →χ 0 2 +X, b)ũ L →χ 0 2 + X and c)g →ũ L + X, d)ũ L →g + X, for tan β = 2, A 0 = 0, µ < 0. Figure 2: Sigma times branching ratios as a function of m 0 and m 1/2 (in GeV) for indirect χ 0 2 production from gluinos (a) and squarks (b), for tan β = 2, A 0 = 0, µ < 0. 13 Figure 3: Branching ratios ofχ 0 2 decays: a)χ 0 2 →χ 0 1 l + l − , b)χ 0 2 →l ± L l ∓ and c)χ 0 2 →l ± R l ∓ as a function of m 0 and m 1/2 , for tan β = 2, A 0 = 0, µ < 0. Figure 4: Domains of the decaysχ 0 2 →χ 0 1 l + l − (dashed line),χ 0 2 →l ± L l ∓ (solid line) and χ 0 2 →l ± R l ∓ (dashed-dotted line) in the (m 0 , m 1/2 ) plane, corresponding to decay branching ratios in excess of 1% and 10% respectively, tan β = 2, A 0 = 0, µ < 0.   . Also shown are the explorable domain in sparticle searches at LEP2 (300 pb −1 ) and the Tevatron (1 fb −1 ), theoretically and experimentally excluded regions [19]. 20 Figure 10: Domains where the observed edge in the M l + l − distribution is due to the decaysχ 0 2 →l ± L l ∓ →χ 0 1 l + l − (solid line),χ 0 2 →l ± R l ∓ →χ 0 1 l + l − (dashed-dotted line), χ 0 2 →χ 0 1 l + l − (dashed line), L int = 10 3 pb −1 .