The DIS 1-Jettiness Event Shape at N$^3$LL+${\cal O}(\alpha_s^2)$

We present results for the $\tau_1$ and $\tau_{1a}$ 1-Jettiness global event shape distributions, for Deep Inelastic Scattering (DIS), at the N$^3$LL + ${\cal O}(\alpha_s^2)$ level of accuracy. These event-shape distributions quantify and characterize the pattern of final state radiation in electron-nucleus collisions. They can be used as a probe of nuclear structure functions, nuclear medium effects in jet production, and for a precision extraction of the QCD strong coupling. The results presented here, along with the corresponding numerical codes, can be used for analyses with HERA data, in EIC simulation studies, and for eventual comparison with real EIC data.