  The generator BCVEGPY is upgraded by adding the hadroproduction of the
$P$-wave excited $B_c$ states (denoted by $B_{cJ,L=1}^*$ or by $h_{B_c}$ and
$\chi_{B_c}$) and by improving some features of the original one as well. We
denote it as BCVEGPY2.0. The $P$-wave production is also calculated by taking
only the dominant gluon-gluon fusion mechanism (with the subprocess $gg\to
B_{cJ,L=1}^*+\bar{c}+b$ being dominated) into account as that for $S$-wave. In
order to make the addendum piece of the upgraded generator as compact as
possible so as to increase its efficiency, we manipulate the amplitude as
compact as possible with FDC (a software for generating Feynman diagrams and
the algebra amplitudes, and for manipulating algebra formulae analytically etc)
and certain simplification techniques. The correctness of the program is tested
by checking the gauge invariance of the amplitude and by comparing the
numerical results with the existent ones in the literature carefully.
