SIDIS at small $x$ at next-to-leading order: gluon contribution

We calculate the contribution of gluons to single inclusive hadron production at next-to-leading order (NLO) accuracy in Deep Inelastic Scattering (DIS) at small $x$ using the Color Glass Condensate formalism. It is shown that the only divergence present is the standard collinear divergence which is absorbed into evolution of gluon-hadron fragmentation function. Our calculations are performed at finite $N_c$ and we provide general finite $N_c$ expressions for the structure of Wilson lines appearing in inclusive dihadron and single hadron production cross sections. We also comment on how one can obtain rapidity distribution of hadron multiplicities from our results.