We perform a phenomenological study of the invariant mass distribution of hadronic jets produced in pp collisions, in conjunction with a groomer (modified MassDrop Tagger and Soft Drop). Our calculation resums large logarithms of the jet mass and includes the full dependence on the groomer’s energy threshold parameter, and it is matched to fixed-order QCD matrix elements at next-to-leading order. We accounted for non-perturbative contributions by including a correction factor derived from Monte Carlo parton-shower simulations