Antenna subtraction at NNLO with identified hadrons

We extend the antenna subtraction method to include hadron fragmentation processes up to next-to-next-to-leading order (NNLO) in QCD in $e^+e^-$ collisions. To handle collinear singularities associated with the fragmentation process, we introduce fragmentation antenna functions in final-final kinematics with associated phase space mappings. These antenna functions are integrated over the relevant phase spaces, retaining their dependence on the momentum fraction of the fragmenting parton. The integrated antenna functions are cross-checked against the known NNLO coefficient functions for identified hadron production from $\gamma^*/Z^* \to q\bar{q}$ and $H \to gg$ processes.