Uncertainties in the simulation of Higgs boson production with up to two jets at next-to leading order accuracy are investigated. Traditional uncertainty estimates based on scale variations are extended employing different functional forms for the central scale, and the impact of details in the implementation of the parton shower is discussed.
