III. PARAMETERS 4 in the speed of the code. For the processes involving two quark-antiquark pairs, the matrix elements squared given in Ref. [16] are implemented in MCFM, rather than the amplitudes of Ref. [23], because they lead to faster code. The values of the amplitudes calculated by the new analytic code and the previous semi-numerical code [12] are in full numerical agreement for all amplitudes. The improvement in the performance of our numerical code means that it is appropriate to revisit the phenomenology of Higgs + 2 jet production and to extend it in a number of ways. The improvement in the speed of the code means that it is possible to include the decays of the Higgs boson, specifically for the processes: where h$_{1}$ , h$_{2}$ represent partons inside the incident hadron beams. All four of these processes are included in MCFM v5.7. Throughout this paper we make use of the MSTW2008 parton distribution functions [25], using the LO fit ( α$_{s}$ ( M$_{Z}$ ) = 0 . 13939 and 1-loop running) for the lowest order calculation and the NLO fit ( α$_{s}$ ( M$_{Z}$ ) = 0 . 12018 and 2-loop running) at NLO. The W mass and width are chosen to be, The mass is taken from Ref. [26]. The total width given in Eq. (9) is derived from the measured branching ratio for W → ℓ ¯ ν , 10 . 80 ± 0 . 09% [26] by using a lowest order calculation m$_{W}$ = 80 . 398 GeV , Γ$_{W}$ = 2 . 1054 GeV . (9) h$_{1}$ + h$_{2}$ → H + j$_{1}$ + j$_{2}$ → τ + + τ − + j$_{1}$ + j$_{2}$ (5) h$_{1}$ + h$_{2}$ → H + j$_{1}$ + j$_{2}$ → b + ¯ b + j$_{1}$ + j$_{2}$ (6) h$_{1}$ + h$_{2}$ → H + j$_{1}$ + j$_{2}$ → W − + W + + j$_{1}$ + j$_{2}$ | | → ν + e + | → e − + ¯ ν (7) h$_{1}$ + h$_{2}$ → H + j$_{1}$ + j$_{2}$ → Z + Z + j$_{1}$ + j$_{2}$ | | → e − + e + | → µ − + µ + (8)