  A method is suggested to combine tree level QCD matrix for the production of
multi jet final states and the parton shower in hadronic interactions. The
method follows closely an algorithm developed recently for the case of $e^+e^-$
annihilations.
