  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.
