  Using the Monte Carlo event generator tools Pythia and Herwig, we simulate
the production of bottom/charm meson and antimeson pairs at hadron colliders in
proton-proton/antiproton collisions. With these results, we derive an
order-of-magnitude estimate for the production rates of the bottom analogues
and the spin partner of the $X(3872)$ as hadronic molecules at the LHC and
Tevatron experiments. We find that the cross sections for these processes are
at the nb level, so that the current and future data sets from the Tevatron and
LHC experiments offer a significant discovery potential. We further point out
that the $X_b/X_{b2}$ should be reconstructed in the $\gamma \Upsilon(nS)
(n=1,2,3)$, $\Upsilon(1S)\pi^+\pi^-\pi^0$, or $\chi_{bJ}\pi^+\pi^-$ instead of
the $\Upsilon(nS)\pi^+\pi^-$ final states.
