In Fig. REF  the mean and standard deviation of the \(p_{\mathrm {T}}\)  spectra are compared for pp, p==Pb and Pb==Pb collisions at the same centre-of-mass energy per nucleon pair of \(\sqrt{s_\mathrm {NN}\,} = 5.02~\mathrm {TeV} \) .
All three collision systems have similar values at \(N_{\mathrm {ch}} = 1\)  and then coincide until Pb==Pb deviates at \(N_{\mathrm {ch}} \approx 5\)  and p==Pb deviates at \(N_{\mathrm {ch}} \approx 25\)  from the trend observed in pp.
This observation is consistent with an earlier comparison of the \(\langle p_\mathrm {T}\rangle \) –\(N_{\mathrm {ch}}\)  correlation for the three systems at different centre-of-mass energies [1]}.
Figure REF  shows the mean (left) and standard deviation (right) of the transverse momentum spectra as a function of the charged-particle multiplicity \(N_{\mathrm {ch}}\)  for pp (top), p==Pb (middle), and AA (bottom) collisions at different centre-of-mass energies per nucleon pair.
For all collision systems, a clear ordering of \(\langle p_\mathrm {T}\rangle \)  as well as \(\sigma (p_\mathrm {T})\)  with collision energy is observed, which can be attributed to the larger momentum transfers involved at higher \(\sqrt{s_\mathrm {NN}\,}\) .
For pp collisions at all centre-of-mass energies, the average transverse momentum increases monotonically with an almost linear trend up to \(N_{\mathrm {ch}}\)  \(\approx 16\)  and beyond that continues with an again almost linear dependence on \(N_{\mathrm {ch}}\)  but reduced slope.
In p==Pb collisions, a similar multiplicity dependence is observed up to \(N_{\mathrm {ch}} \approx 25\) .
At higher multiplicities, the increase in \(\langle p_\mathrm {T}\rangle \)  is slower than in pp collisions.
In both pp and p==Pb, \(\sigma (p_\mathrm {T})\)  follows a similar trend as \(\langle p_\mathrm {T}\rangle \) .
On the other hand, for AA collisions one observes an increase of \(\langle p_\mathrm {T}\rangle \)  with multiplicity up to about one third of the measured range, followed by a constant trend for the rest of the \(N_{\mathrm {ch}}\)  range.
The \(\sigma (p_\mathrm {T})\)  increases for \(N_{\mathrm {ch}} \lesssim 100\)  to a maximum and decreases afterwards.
This is unique to large collision systems and is presumably a consequence of flow and jet quenching [2]}.
The high \(N_{\mathrm {ch}}\)  resolution of this measurement makes it possible to spot differences between the spectral evolution with multiplicity in Xe==Xe and Pb==Pb collisions at \(\sqrt{s_\mathrm {NN}\,} = 5.44~\mathrm {TeV} \)  and \(\sqrt{s_\mathrm {NN}\,} = 5.02~\mathrm {TeV} \) , respectively.
The observed difference in the trends might be a result of the slightly deformed Xe nuclei [3]}.
In Fig. REF  both the mean (left) and standard deviation (right) of the \(p_{\mathrm {T}}\)  spectra as a function of \(N_{\mathrm {ch}}\)  are summarised for all data sets (top panels) and then shown as a function of relative multiplicity \(N_{\mathrm {ch}}/ \langle N_{\mathrm {ch}}\rangle  \)  (middle panels) as well as divided by their respective multiplicity-integrated values (bottom panels).
In the latter scaling, the overall energy dependent increase of average kinematic energy and number of produced particles are accounted for.
As a result, the values for each collision system align almost perfectly for the \(\langle p_\mathrm {T}\rangle / \langle p_\mathrm {T}\rangle  _{\mathrm {incl}}\) .
For \(\sigma (p_\mathrm {T})/ \sigma (p_\mathrm {T}) _{\mathrm {incl}}\) , for a given system the scaled value increases slightly with centre-of-mass energy.
This small remaining ordering may be due to the larger available phase space at higher centre-of-mass energies.
<FIGURE><FIGURE><FIGURE>