  Massive on-shell operator matrix elements and self-energy diagrams with outer
gluon lines are calculated analytically at $O(\alpha_s^2)$, using Mellin-Barnes
integrals and representations through generalized hypergeometric functions.
This method allows for a direct evaluation without decomposing the integrals
using the integration-by-parts method.
