The main deliverable of this paper are \(\text{NLO} +\text{NLL}  ^\prime \)  predictions for jet angularities in \(Z\) +jets events in
proton–proton collisions that can be immediately exploited for data–theory comparisons. However, this
work has already shown us the path towards higher precision that we need to pursue. In this context, we
can identify four complementary lines of research. First of all, although we do not expect large corrections,
one should improve the accuracy of the resummation, in order to reduce the theoretical uncertainty. The
framework for NNLL calculations has been worked out using SCET, although
explicit results only exist, to the best of our knowledge, for the case \(\alpha =2\)  [1]}, [2]}.
Second, it would be desirable to match the resummation to NNLO distributions. However, we recognise that this endeavour
is much more formidable as it requires calculating \(Z\)  plus two
partons at NNLO accuracy.
Some first steps have however already been taken in this direction [3]}, [4]}, [5]}, [6]}.
Next, given that the size of non-perturbative corrections can remain
large even for groomed angularities, e.g. for \(\alpha < 1\) , it
would be interesting to study additional grooming strategies, such as
complementing the SoftDrop procedure with a filtering
step [7]} or using the Recursive
SoftDrop algorithm [8]}, although these tools would most
likely further complicate the structure of the analytic resummation.
Finally, a SCET framework
to deal with non-perturbative effects has been recently developed [9]}, [10]}
and it would be interesting to quantitatively compare its predictions to the phenomenological model employed here.
