We consider the O ( α α s ) corrections to single on-shell gauge boson production at hadron colliders. We concentrate on the contribution of all the subprocesses where the gauge boson is accompanied by the emission of two additional real partons and we evaluate the corresponding total cross sections. The latter are divergent quantities, because of soft and collinear emissions, and are expressed as Laurent series in the dimensional regularization parameter. The total cross sections are evaluated by means of reverse unitarity, i.e. expressing the phase-space integrals in terms of two-loop forward box integrals with cuts on the final-state particles. The results are reduced to a combination of master integrals, which eventually are evaluated in terms of generalized polylogarithms. The presence of internal massive lines in the Feynman diagrams, due to the exchange of electroweak gauge bosons, causes the appearance of 14 master integrals which were not previously known in the literature and have been evaluated via differential equations.