  This paper presents the method of combining parton shower formalism with
perturbative calculations (matrix elements) in form of a Monte-Carlo algorithm
for the process g g -> Z b anti-b,consistenlty including the heavy quark masses
and overlap removal.
