The CP Asymmetry in $b \to s l^+ l^-$ Decay

Using the experimental upper bound on the neutron EDM and experimental result on $b\rightarrow s\gamma$ branching ratio we have calculated CP asymmetry and $\Gamma^{2HDM}(b\rightarrow s{\it{l^{+}l^{-}}})/\Gamma^{SM}(b\rightarrow s{\it{l^{+}l^{-}}})$. It is shown that in the invariant dilepton mass $q^{2}$ region $(m_{\psi'}^{2}+0.2\; GeV^{2})<\;q^{2}<\;m^{2}_{b}$ the CP asymmetry is maximal and quite detectable.