Two photon radiation inWandZboson production at the Fermilab Tevatron collider

We present a calculation of two photon radiation in W and Z boson production in hadronic collisions, based on the complete matrix elements for the processes $q{q}^{\prime }\to l^{\pm }\nu \gamma \gamma$ and $q\overline{q}\to l^{+}{l}^{-}\gamma \gamma ,$ including finite charged lepton masses. In order to achieve stable numerical results over the full phase space, multiconfiguration Monte Carlo techniques are used to map the peaks in the differential cross section. Numerical results are presented for the Fermilab Tevatron.