In high energy hadronic collisions, a general CP -violating Higgs boson Φ with accompanying jets can be efficiently produced via gluon fusion, which is mediated by heavy quark loops. In this article, we study the dominant sub-channel gg → gggΦ of the gluon fusion production process with triple real emission corrections at order α 5 s . We go beyond the heavy top-quark approximation and include the full mass dependence of the top-and bottom-quark contributions. Furthermore, in a specific model we demonstrate the features of our program and show the impact of bottom-quark loop contributions in combination with large values of tan β on differential distributions sensitive to CP -measurements of the Higgs boson.