  We calculate the next-to-leading QCD corrections to the effective Hamiltonian
for \Bsee in the NDR and HV schemes. We give for the first time analytic
expressions for the Wilson Coefficient of the operator $Q_9 = (\bar s
b)_{V-A}(\bar e e)_V$ in the NDR and HV schemes. Calculating the relevant
matrix elements of local operators in the spectator model we demonstrate the
scheme independence of the resulting short distance contribution to the
physical amplitude. Keeping consistently only leading and next-to-leading
terms, we find an analytic formula for the differential dilepton invariant mass
distribution in the spectator model. Numerical analysis of the $\mt$, $\Lms$
and $\mu \approx {\cal O}(\mb)$ dependences of this formula is presented. We
compare our results with those given in the literature.
