  The most general model-independent analysis of the rare $B$ decay, $\Bsll$,
is presented. There are ten independent local four-Fermi interactions which may
contribute to this process. The branching ratio, the forward-backward
asymmetry, and the double differential rate are written as functions of the
Wilson coefficients of the ten operators. We also study the correlation between
the branching ratio and the forward-backward asymmetry by changing each
coefficient. This procedure tells us which types of operator contribute to the
process, and it will be very useful to pin down new physics systematically,
once we have the experimental data with high statistics and the deviation from
the Standard Model is found.
