Master Integrals for the mixed EW-QCD corrections to the Drell-Yan production of a massive lepton pair

We showcase the calculation of the master integrals needed for the two loop mixed QCD-Electroweak virtual corrections to the neutral current Drell-Yan process $(q\bar{q}\rightarrow l^+ l^-)$. After establishing a basis of 51 master integrals, we cast the latter into canonical form by using the Magnus algorithm. The dependence on the lepton mass is then expanded such that potentially large logarithmic contributions are kept. After determining all boundary constants, we give the coefficients of the Taylor series around four space-time dimensions in terms of generalized polylogarithms up to weight four.