We present results of the simulation of electroweak Higgs boson production at the Large Hadron Collider using the NLO multi-jet merging framework provided by the general purpose event generator Herwig 7. For the hard processes, we use the HJets library for the computation of the $mathcal{O}(alpha^3 alpha_{s}^{n-2})$ matrix elements for $pp to h+n$ jet production at LO for $n=2,3,4$ and NLO for $n=2,3$.