Electroweak corrections to $gg\rightarrow \gamma\gamma$

We present the electroweak corrections for the production of a photon pair through gluon fusion, focusing on the contribution from the first two generations of quarks. The two-loop amplitude is calculated using a series of projection operators which define scalar form factors. In order to evaluate the Master Integrals which appear in this process we employ both generalized polylogarithms and Chen-iterated integrals. In order to perform a phenomenological study we develop a semi-numerical evaluation of the Master Integrals employing a fitting procedure to speed up the evaluation of burdensome higher weight contributions. We present results for the LHC, finding corrections of around a couple of percent to the leading order $gg \rightarrow \gamma\gamma$ process. Our results are implemented into the parton-level Monte Carlo code \texttt{MCFM}.