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}$, $\mathrm{tan}\ensuremath{\beta}$, ${m}_{t}$, ${m}_{b}$ and ${\ensuremath{\alpha}}_{s}({M}_{Z})$ are applied using a Markov chain Monte Carlo algorithm. We use measurements of $b\ensuremath{\rightarrow}s\ensuremath{\gamma}$, $(g\ensuremath{-}2{)}_{\ensuremath{\mu}}$ and ${\ensuremath{\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}\ensuremath{\rightarrow}{\ensuremath{\mu}}^{+}{\ensuremath{\mu}}^{\ensuremath{-}}$ and for ${m}_{\stackrel{\texttildelow{}}{\ensuremath{\tau}}}\ensuremath{-}{m}_{{\ensuremath{\chi}}_{1}^{0}}$. An upper limit of $2\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}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 ${\stackrel{\texttildelow{}}{q}}_{L}\ensuremath{\rightarrow}{\ensuremath{\chi}}_{2}^{0}\ensuremath{\rightarrow}{\stackrel{\texttildelow{}}{l}}_{R}\ensuremath{\rightarrow}{\ensuremath{\chi}}_{1}^{0}$ with a likelihood of $24\ifmmode\pm\else\textpm\fi{}4%$. The stop-coannihilation region is highly disfavored, whereas the light Higgs region is marginally disfavored.