Multijet matrix elements and shower evolution in hadronic collisions: $W b \bar{b}$ + $n$ jets as a case study

We study in this paper the production, in hadronic collisions, of final states with W gauge bosons, heavy quark pairs and n extra jets (with n up to 4). The complete partonic tree-level QCD matrix elements are evaluated using the ALPHA algorithm, and the events generated at the parton level are then evolved through the QCD shower and eventually hadronised using the coherent shower evolution provided by the HERWIG Monte Carlo. We discuss the details of our Monte Carlo implementation, and present results of phenomenological interest for the Tevatron Collider and for the LHC. We also comment on the impact of our calculation on the backgrounds to W (H -> b bbar) production, when only one b jet is tagged.