Supersymmetric monojets may be produced at the Large Hadron Collider by the process qg -> squark neutralino_1 -> q neutralino_1 neutralino_1, leading to a jet recoiling against missing transverse momentum. We discuss the feasibility and utility of the supersymmetric monojet signal. In particular, we examine the possible precision with which one can ascertain the neutralino_1-squark-quark coupling via the rate for monojet events. Such a coupling contains information on the composition of the neutralino_1 and helps bound dark matter direct detection cross-sections and the dark matter relic density of the neutralino_1. It also provides a check of the supersymmetric relation between gauge couplings and gaugino-quark-squark couplings.

We study the problem of selecting between different mechanisms of supersymmetry breaking in the Minimal Supersymmetric Standard Model using current data. We evaluate the Bayesian evidence of four supersymmetry breaking scenarios: mSUGRA, mGMSB, mAMSB and moduli mediation. The results show a strong dependence on the dark matter assumption. Using the inferred cosmological relic density as an upper bound, minimal anomaly mediation is at least moderately favoured over the CMSSM. Our fits also indicate that evidence for a positive sign of the mu parameter is moderate at best. We present constraints on the anomaly and gauge mediated parameter spaces and some previously unexplored aspects of the dark matter phenomenology of the moduli mediation scenario. We use sparticle searches, indirect observables and dark matter observables in the global fit and quantify robustness with respect to prior choice. We quantify how much information is contained within each constraint.

In most global fits of the constrained minimal supersymmetric model (CMSSM) to indirect data, the a priori likelihoods of any two points in tan beta are treated as equal, and the more fundamental mu and B Higgs potential parameters are fixed by potential minimization conditions. We find that, if instead a flat ("natural")prior measure on mu and B is placed, a strong preference exists for the focus point region from fits to particle physics and cosmological data. In particular, we find that the lightest neutralino is strongly favored to be a mixed bino-higgsino (~10% higgsino). Such mixed neutralinos have large elastic scattering cross sections with nuclei, leading to extremely promising prospects for both underground direct detection experiments and neutrino telescopes. In particular, the majority of the posterior probability distribution falls within parameter space within an order of magnitude of current direct detection constraints. Furthermore, neutralino annihilations in the sun are predicted to generate thousands of neutrino induced muon events per years at IceCube. Thus, assuming the framework of the CMSSM and using the natural prior measure, modulo caveats regarding astrophysical uncertainties, we are likely to be living in a world with good prospects for the direct and indirect detection of neutralino dark matter.

We provide a bestiary of public codes and other algorithmic tools that can be used for analysing supersymmetric phenomenology. We also describe the organisation of the different tools and communication between them. Tools exist that calculate supersymmetric spectra and decay widths, simulate Monte Carlo events as well as those that make predictions of dark matter relic density or that predict precision electroweak or b-observables. Some global fitting tools for use in SUSY phenomenology are also presented. In each case, a description and a link to the relevant web-site is provided. It is hoped that this review could serve as an "entry-gate" and map for prospective users.

We calculate the likelihood map in the full 7 dimensional parameter space of the minimal supersymmetric standard model (MSSM) 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 -> 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 -> mu^+ mu^- and for m_{stau}-m_{chi_1^0}. An upper limit of 2.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 squark_L -> chi_2^0 -> slepton_R -> chi_1^0 with a likelihood of 24+/-4%. The stop coannihilation region is highly disfavoured, whereas the light Higgs region is marginally disfavoured.

We examine the effect of a prior that favours low values of fine-tuning on Bayesian multi-dimensional fits of the constrained minimal supersymmetric standard model (CMSSM or mSUGRA) to current data. The dark matter relic density, the anomalous magnetic moment of the muon and the branching ratio of b->s gamma are all used to constrain the model via a Markov Chain Monte Carlo sampler. As a result of the naturalness prior, posterior probability distributions skew towards lighter higgs and sparticle masses, the effect being most pronounced in the gaugino sector. Interestingly, slepton masses are an exception and skew towards heavier masses. The lightest CP-even Higgs h^0-pole region becomes allowed at the 2 sigma level for the latest combination of measurements of m_t=172.7+/-2.9 GeV, provided we allow for a theoretical error in the prediction of its mass m_{h^0}. m_{h^0} is constrained to be less than 120 GeV at the 95% C.L. Probing the branching ratio of B_s->mu^+ mu^- to the level of 2 x 10^{-8}, as might be achieved by the Tevatron experiments, would cover 32% of the probability density, irrespective of which of the two priors is used.