A method is suggested to combine tree level QCD matrix elements 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 \cite{Catani:2001cc}.