The leptonic decays of the heavy gauge bosons W and/or Z provide a clear experimental signature at hadron colliders. The production of accompanying jets is an excellent signal to probe QCD, while also being the main background to many searches for ne
w physics. Describing the complex final state of W or Z + jets is a theoretical challenge with most existing calculations combining matrix elements for high energy jet production with a parton shower for lower energy jet production. We focus on two models: SHERPA, which uses Leading Order matrix elements for boson and jet production; and POWHEG with HERWIG++, which uses a Next-To-Leading Order Matrix element for Z production. In order to isolate the impact of the matrix elements for jet production, it is first essential to constrain the differences in the rest of the calculation in each case: specifically, the Multiple Parton Interaction models, and the tuning of the parton shower interfaced to the matrix elements. We test all three aspects of these models against data from the Tevatron, and perform a study of some basic kinematic variables at the LHC energy.
