Ad Lucem: QED parton distribution functions in the MMHT framework