Lepton pair production at NNLO in QED with EW effects

We present a fully differential calculation of lepton pair production, taking into account the dominant next-to-next-to-leading order QED corrections as well next-to-leading electroweak and polarisation effects. include all masses, hard photon emission, non-perturbative hadronic corrections. The corresponding matrix elements are implemented in Monte Carlo framework McMule. In to obtain numerically stable implementation, we extend next-to-soft stabilisation, universal technique based on next-to-leading-power expansion, calculations with polarised leptons. As an example, show results tailored Belle II detector current setup potential future configuration that includes beams.