R=0.2 hadron+jet $\Delta\phi$ distribution for $20 < p_{\mathrm{T,jet}^{ch}} < 30$ GeV/c in pp collisions at $\sqrt{\it{s}_{NN}} = 5.02$ TeV compared to PYTHIA and pQCD calculation

$\Delta_\mathrm{recoil}$ is calculated as:

$\Delta_\mathrm{recoil} = \frac{1}{N^\mathrm{AA}_\mathrm{trig}} \frac{\mathrm{d^3}N^\mathrm{AA}_\mathrm{{jet}}}{\mathrm{d}p^{\mathrm{ch}}_\mathrm{T,jet} \mathrm{d}\Delta\varphi \mathrm{d}\eta_\mathrm{jet}} \bigg|_{p_\mathrm{T,trig} \in \mathrm{TT_{Sig}}} - c_\mathrm{ref} \cdot \frac{1}{N^\mathrm{AA}_\mathrm{trig}} \frac{\mathrm{d^3}N^\mathrm{AA}_\mathrm{{jet}}}{\mathrm{d}p^{\mathrm{ch}}_\mathrm{T,jet} \mathrm{d}\Delta\varphi \mathrm{d}\eta_\mathrm{jet}} \bigg|_{p_\mathrm{T,trig} \in \mathrm{TT_{Ref}}}$

Here $c_\mathrm{ref}$ is a scaling factor to account for conservation of jet density, and the trigger track intervals are in this analysis $TT_{Ref}: 5 <  p_\mathrm{T,trig} < 7$ GeV/c and $TT_{Sig}: 20 <  p_\mathrm{T,trig} < 50$ GeV/c. Jets are reconstructed with charged particles with $p_{T,track} > 0.15$ GeV/c, using the anti-$k_T$ algorithm with resolution parameter R=0.2.

The PYTHIA8 calculation uses the Monash2013 tune [1].

The pQCD@NLO is that described in [2]. The pp prediction includes perturbative parton production with Sudakov resummation. The predictions for the jet-normalised $\Delta\varphi$ distributions are scaled such that the magnitude of the prediction matches that of $\Delta_\mathrm{recoil}$ at $\Delta\varphi \sim \pi$.

[1] EPJC 74 3024 (2014), arXiv:1404.5630 [hep-ph] 
[2] Phys. Lett. B 773 (2017) 672-676, arXiv:1607.01932 [hep-ph]