McMule is a framework for fully
differential higher-order QED calculations of scattering and decay
processes involving leptons. It keeps finite lepton masses, which
regularises collinear singularities. Soft singularities are treated with
dimensional regularisation and using FKS^\ellℓ
subtraction. We describe the implementation of the framework in
Fortran 95, list the processes that are currently implemented, and give
instructions on how to run the code. In addition, we present new
phenomenological results for muon-electron scattering and lepton-proton
scattering, including the dominant NNLO corrections. While the
applications presented focus on MUonE, MUSE, and P2, the code can be
used for a large number of planned and running experiments.