Exclusive J/psi production off a dilute proton within a refined hotspot description

We revisit the calculation of exclusive $\mathrm{J}/\psi$ production in electron-proton scattering within the QCD dipole model, employing a refined description of the proton in terms of gluonic hot spots in the dilute regime. In contrast to earlier studies, we incorporate key missing elements: the event-by-event fluctuations in both the number of hot spots and their local saturation scales, as well as the first relativistic corrections to the $\mathrm{J}/\psi$ wave function. We perform a Bayesian analysis to constrain the model parameters using HERA data. These enhancements significantly improve the agreement with HERA data. Saturation scale fluctuations are found to be more important than fluctuations in the number of hot spots. Including the first relativistic correction is found to be important, especially in order to describe the observed power-law behavior in the incoherent cross section at large $t$.