Abstract: Higgs production in association with three jets via gluon fusion is an important channel for the measurement of the CP properties of the Higgs particle at the LHC. In this paper, we go beyond the heavy top effective theory approximation and include at LO the full mass dependence of the top- and bottom-quark contributions. We consider the dominant subchannel gg -> Hggg which involves the manipulation of massive rank-5 hexagon integrals. Furthermore, we present results for several differential distributions and show deviations from the effective theory as large as 100% at high p(T) for light Higgs masses.