  We address the issues of clustering and non-global logarithms for jet shapes
in the process of production of a Higgs/vector boson associated with a single
hard jet at hadron colliders. We perform an analytical fixed-order calculation
up to second order in the coupling as well as an all-orders estimation for the
specific invariant mass distribution of the highest-$p_t$ jet, for various jet
algorithms. Our results are derived in the eikonal (soft) limit and are valid
up to next-to-leading logarithmic accuracy. We perform a matching of the
resummed distribution to next-to-leading order results from MCFM and compare
our findings with the output of the Monte Carlo event generator Pythia 8. After
accounting for non-perturbative effects we compare our results with available
experimental data from the CMS collaboration for the Z + jet production. We
find good agreement over a wide range of the observable.
