In this contribution we present new analytic and numerical results for two-loop QCD amplitudes
mediated by massive-quark loops, relevant for dijet and 𝛾+jet production at hadron colliders. Building on recent developments in elliptic Feynman integrals, we construct a generalized canonical basis for the full set of master integrals. This basis allows us to write the amplitudes in terms of iterated integrals, providing a convenient framework to study their analytic properties and revealing several simplifications for the amplitudes up to the finite part. Moreover, we discuss how fast and reliable numerical evaluations in the physical region can be obtained through series expansion representations of the amplitudes.

