Electron–proton collisions are an important mode for, e.g., determination of parton distribution functions, making it necessary to understand the dynamics and kinematics of these collisions. To achieve highest precision in the different observables it is necessary to develop simulations of these events which leverage the Fixed Order calculations and match them to the parton showers.
I will review the importance of photo-production as the complementary regime to DIS and present the progress in the Sherpa event generator to achieve NLO accuracy matched to the parton shower for photo-production and present the corresponding procedure. I will then show validation against data from OPAL and ZEUS and review the quality of parton distributions inside the photon and the interplay with photo-production phenomenology. Last but not least, predictions for the upcoming EIC are presented for event shapes and jet observables.