Total: 1
We develop CGC-py, a Monte Carlo event generator for deep-inelastic scattering. It couples the full Color Glass Condensate (CGC) cross section for $γ^*p(A)\to q\bar q+X$ to a Parton-Branching transverse-momentum-dependent backward initial-state shower, while \textsc{Pythia}~8 handles final-state radiation and hadronization. CGC-py retains the complete target-elastic and target-inelastic contributions without taking the back-to-back correlation limit, allowing single- and di-hadron observables to be generated consistently from the same event sample. We validate the generator through an analytic closure test of the single-inclusive quark spectrum and a comparison of charged-hadron spectra in $ep$ collisions with H1 data, finding excellent agreement. The predicted nuclear modification factor $R_{e\mathrm{Au}}^h$ shows the expected saturation pattern: suppression at low $p_T^*$ followed by a rise toward unity at higher $p_T^*$. A comparison with a \textsc{Pythia}~6 baseline, together with an $x_g$-rescaling study, indicates that small-$x$ CGC evolution and collinear DGLAP dynamics contribute comparably to the growth of the dihadron away-side width with energy. Genuine saturation-driven broadening emerges only at the highest energies considered. Within CGC-py, $e\mathrm{Au}$ collisions exhibit an enhanced away-side width and a suppressed back-to-back yield relative to $ep$ collisions. These nuclear effects remain modest over EIC kinematics, motivating measurements at the most forward accessible kinematics and the use of complementary observables to maximize sensitivity to gluon saturation.