Constraining rare B decays by μ+μ−→tc

Motivated by the recent rare $B$ decays measurements, we study the matching procedure of operators ${O}_{9}$, ${O}_{10}$ in the low energy effective Hamiltonian and operators in the Standard Model effective theory (SMEFT). It is noticed that there are more related operators in the SMEFT whose coefficients cannot be determined only from the low-energy data from $B$ physics. We demonstrate how to determine these coefficients with some new physics models, like the ${Z}^{\ensuremath{'}}$ model and leptoquark models, and then consider how to probe these operators of SMEFT at high energy by using the process ${\ensuremath{\mu}}^{+}{\ensuremath{\mu}}^{\ensuremath{-}}\ensuremath{\rightarrow}tc$ at future muon colliders, which can provide complementary information except for ${\ensuremath{\mu}}^{+}{\ensuremath{\mu}}^{\ensuremath{-}}\ensuremath{\rightarrow}bs$ on the underlying models which lead to rare $B$ decay processes. We perform a Monte Carlo study (a hadron level analysis) to show how to separate the signal events from the SM background events and estimate the sensitivity to the Wilson coefficients for different models.