In these proceedings, we apply the recently developed S-ACOT-MPS factorization scheme at the next-to-leading order to prompt charm production at hadron colliders. It provides a good agreement with experimental data on charm meson production measured by LHCb at 7 and 13 TeV. The low-$p_T$ data are on the margins of the theoretical error bands, emphasizing the importance of including contributions beyond the next-to-leading order.