  Motivated by the need to increase the precision of theoretical predictions to
test saturation physics at both the LHC and the EIC, we compute the
cross-sections for the diffractive single and di-hadron production at the NLO
in the shockwave formalism.
