  Theoretical predictions used in experimental analysis of LHC data have an
inherent theoretical uncertainty associated to the relevant order in the
perturbative expansion that the observable is computed. In this talk we briefly
introduce and motivate the topic of QCD radiative corrections stressing the
impact of including higher-order QCD effects in the theory predictions and in
particular we consider two-loop QCD corrections. We present numerical results
at NNLO for gluonic jet production where jets are reconstructed using the
anti-$k_T$ jet algorithm from gluon-gluon and quark-antiquark scattering and
discuss the inclusion of leading-$N_F$ quark contributions in the final-state.
