Parton showers are an essential step in the simulation of realistic collider events because the emission of soft and collinear radiation from charged partons is logarithmically enhanced, which spoils the perturbative convergence of fixed order calculations. Using the universal factorisation of these enhanced limits, parton shower Monte Carlo provides a process-independent, numerical resummation of the large logarithms, allowing an accurate description of multiple emissions. To combine a parton shower with a hard calculation beyond leading order, careful matching is required in order to preserve both fixed order and logarithmic accuracy where both calculations cover the same phase space. Furthermore, the arbitrary recoil assignments inherent in a practical parton shower implementation can introduce higher order, but numerically significant, artefacts in resonant processes due to the rapidly changing cross section near the peak. We will use the MC@NLO algorithm to match a QCD+QED parton shower to NLO QCD+EW and explore how this calculation must be modified for resonant processes.