Probing Supersymmetry in Rare B Decays

We determine the ability of future experiments to observe supersymmetric contributions to the rare decays $B \to X_s \gamma$ and $B \to X_s l^+l^-$. A global fit to the Wilson coefficients which contribute to these decays is performed from Monte Carlo generated data. This fit is then compared to supersymmetric predictions for several different patterns of the superpartner spectrum.