Next-to-leading order QCD factorization for semi-inclusive deep inelastic scattering at twist 4.