  We present a next-to-leading order calculation of Higgs boson production plus
one and two jets via gluon fusion interfaced to shower Monte Carlo programs,
implemented according to the POWHEG method. For this implementation we have
used a new interface of the POWHEG BOX with MadGraph4, that generates the codes
for generic Born and real processes automatically. The virtual corrections have
been taken from the MCFM code. We carry out a simple phenomenological study of
our generators, comparing them among each other and with fixed next-to-leading
order results.
