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 ggZ 0/γ*b, consistently including the heavy quark masses and overlap removal.