Multidimensional mSUGRA likelihood maps

We calculate the likelihood map in the full 7-dimensional parameter space of the minimal supersymmetric standard model assuming universal boundary conditions on the supersymmetry breaking terms. Simultaneous variations of $m_0$, $A_0$, $M_{1/2}$, $\tan \beta$, $m_t$, $m_b$ and ${\alpha }_s(M_Z)$ are applied using a Markov chain Monte Carlo algorithm. We use measurements of $b\to s\gamma$, $(g-2{)}_{\mu }$ and ${\Omega }_{DM}{h}^2$ in order to constrain the model. We present likelihood distributions for some of the sparticle masses, for the branching ratio of $B_s^0\to {\mu }^{+}{\mu }^{-}$ and for $m_{\tilde{\tau }}-m_{{\chi }_1^0}$. An upper limit of $2\times {10}^{-8}$ on this branching ratio might be achieved at the Tevatron, and would rule out $29\%$ of the currently allowed likelihood. If one allows for non-thermal-neutralino components of dark matter, this fraction becomes $35\%$. The mass ordering allows the important cascade decay ${\tilde{q}}_L\to {\chi }_2^0\to {\tilde{l}}_R\to {\chi }_1^0$ with a likelihood of $24\pm 4\%$. The stop-coannihilation region is highly disfavored, whereas the light Higgs region is marginally disfavored.