  We present results of the simulation of electroweak Higgs boson production at
the Large Hadron Collider using both the NLO multijet merging and NLO matching
frameworks provided by the general purpose event generator Herwig 7. For the
simulation of the hard scattering processes, we use the HJets library for the
full calculation and VBFNLO for the approximate calculation to compute the
$2\to h+n$ amplitudes at tree-level with $n=2,3,4$ and at one-loop with
$n=2,3$.
