The decay of a neutral Bs meson into a muon pair is a very rare process whose rate in principle could be severely affected by new physics beyond the standard model. We now know it is not: given the rate measured by the LHCb experiment, any new contribution to the decay amplitude has to be smaller than the standard model one. There's a medical discussion going on and on about the interpretation of this result in the context of supersymmetry. Indeed, the statements describing the LHCb result as "a blow to supersymmetry" or "putting SUSY into hospital" are silly (if you think it's the most spectacular change of loyalties since Terminator 2, read on till the end ;-) But what is the true meaning of this result?
To answer this question at a quantitative level it pays to start with a model independent approach (and technical too, to filter the audience ;-) B-meson decays are low-energy processes which are properly described within a low-energy theory with heavy particles, like W/Z bosons or new physics, integrated out. That is to say, one can think of the Bs→μμ decay as caused by effective 4-fermion operators with 1 b-quark, 1 s-quark, and 2 muons:
Naively, integrating out a mediator with mass M generates a 4-fermion operator suppressed by M^2. In the standard model, only the first operator is generated with ML,SM≈17 TeV, dominantly by the diagram with the Z-boson exchange pictured here. That scale is much bigger than the Z mass because the diagram is suppressed by a 1-loop factor, and furthermore it is proportional to the CKM matrix element V_ts whose value is 0.04. The remaining operators do not arise in the SM, in particular there are no scalars that could generate MS or MP (the Higgs boson couples to mass, thus by construction it has no flavor violating couplings to quarks).
In terms of the coefficients of these operators, the Bs→μμ branching fraction relative to the SM one is given by
LHCb says that this ratio should not be larger than 2 or smaller than 1/3. This leads to model-independent constraints on the mass scales suppressing the 4-fermion operators. And so, the lower bound on ML and MR is about 30 TeV, that is similar in size of the standard model contribution. The bound on the scalar and pseudoscalar operators is much stronger: MS,MP≳150,200 TeV. \begin{digression} The reason is that the contribution of the vector operators to the Bs→μμ decay is suppressed by the small ratio of muon and Bs masses, which goes under the name of helicity suppression. Bs is spin zero, and a vector particle mediating the decay always couples to 2 muons of the same chirality. In the limit mμ=0, when chirality=helicity, the muons spins add up, which forbids the decay by spin conservation \end{digression}.
Consequently, the LHCb result can be interpreted as a constraint on new physics capable of generating the 4-fermion operators listed above. For example, a generic pseudoscalar with order 1 couplings and flavor violating couplings to quarks and leptons must be heavier than about 100 TeV. It may sound surprising that the LHC can probe physics above 100 TeV, even if indirectly. But this is in fact typical for B-physics: observables related to CP violation and mixing of B-mesons are sensitive to similar energy scales (see e.g Table I of this paper). Notice however that 100 TeV is not a hard bound on new pseudoscalars. If the new physics has a built-in mechanism suppressing the flavor violating couplings then even weak scale masses may be allowed.
Now, what happens in SUSY? The bitch always comes in package with an extended Higgs sector, and the exchange of the heavier cousins of the Higgs boson can generate the operators MS and MP. However, bounds on the heavy Higgs masses from Bs→μμ will always be much weaker than 100 TeV quoted above. Firstly, the Higgses couple to mass, thus the Yukawa couplings relevant for this decay are much smaller than one. Secondly, the Higgses have flavor conserving couplings at tree-level, and flavor violation is generated only at 1 loop. Finally, models of low-energy SUSY always assume some mechanism to suppress flavor violation (otherwise all hell breaks loose); in typical realizations flavor violating amplitudes will be suppressed by the CKM matrix elements, much as in the standard model. All in all, SUSY appears less interesting in this context than other new physics models, and SUSY contributions to Bs→μμ are typically smaller than the standard model ones.
But then SUSY has many knobs and buttons. The one called tanβ -- the ratio of the vacuum values of the two Higgs fields -- is useful here because the Yukawa couplings of the heavy Higgses to down-type quarks and leptons happen to be proportional to tanβ. Some SUSY contributions to the branching fraction are proportional to the 6th power of tanβ. It is then possible to pump up tanβ such that the SUSY contribution to Bs→μμ exceeds the standard model one and becomes observable. For this reason, Bs→μμ was hailed as a probe of SUSY. But, at the end of the day, the bound from Bs→μμ on the heavy Higgs masses is relevant only in the specific corner of the parameter space (large tanβ), and even then the SUSY contribution crucially depends on other tunable parameters: Higgsino and gaugino masses, mass splittings in the squark sector, the size of the A-terms, etc. This is illustrated by the plot on the right where the bounds (red) change significantly for different assumptions about the μ-term and the sign of the A-term. Thus, the bound may be an issue in some (artificially) constrained SUSY scenarios like mSUGRA, but it can be easily dodged in more the general case.
To conclude, you should interpret the LHCb measurement of the Bs→μμ branching fraction as a strong bound on theories on new physics coupled to leptons and, in a flavor violating way, to quarks. In the context of SUSY, however, there are far better reasons to believe her dead (flavor and CP, little hierarchy problem, direct searches). So one should not view Bs→μμ as the SUSY killer, but as just another handful of earth upon the coffin ;-)
Some pictures borrowed from Mathieu Perrin-Terrin's talk.