  We have written a Fortran programme BCVEGPY, which is an event generator for
the hadronic production of the $B_c$ meson through the dominant hard subprocess
$gg\to B_c(B_c^*) +b+\bar{c}$. To achieve a compact programme, we have written
the amplitude of the subprocess with the particle helicity technique and made
it as symmetric as possible, by decomposing the gluon self couplings and then
applying the symmetries. To check the programme, various cross sections of the
subprocess have been computed numerically and compared with those in the
literature. BCVEGPY is written in a PYTHIA-compatible format, thus it is easy
to implement in PYTHIA.
