Jet Structure Studies in Small Systems

A study investigating a possible jet shape dependence on the charged event multiplicity was performed on collision samples generated by Monte-Carlo (MC) event generators PYTHIA and HIJING++. We calculated the integral jet shape and found a significant modification caused by multiple-parton interactions. By interchanging and enabling different model ingredients in the simulations and analyzing the results in several $p_T$ bins and event multiplicity classes, we found a characteristic jet size measure that was independent of the chosen tunes, settings, and jet reconstruction algorithms.


Introduction
The discovery of collective-like behavior in high-multiplicity proton-proton (pp) and proton-nucleus (pA) collisions was one of the major surprises in early LHCresults [1,2]. The collective-like behavior previously found in large systems, manifested in long-range correlations and a sizeable azimuthal anisotropy, have traditionally been considered as a signature proving the presence of the quark-gluon plasma (QGP). The creation of the QGP in high-energy nucleus-nucleus (AA) collisions can also be investigated by studying the structure of jets and their modification that leads to the well-known jet quenching phenomenon [3,4]. Whether QGP is created in small systems like pp collisions is still an open question [5]. However, the presence of the QGP is not necessary to explaining collectivity: Relatively soft vacuum-QCDeffects such as multiple-parton interactions (MPI) and color reconnection (CR) can also produce a similar behavior [6,7]. These interactions, at least in principle, can modify the jet shapes in even small systems. Although experimental confirmation is not yet available, a recent phenomenology study also suggests the modification of hard processes by soft vacuum-QCD effects in a high-multiplicity environment [8].
MPI is also expected to depend on flavor [9]. Fragmentation of heavy-flavor jets is expected to differ from light-flavor jets because of color charge and mass effects. The internal structures of heavy-flavor jets may therefore provide a deeper insight into the flavor-dependent development of jets and their connection to the underlying event (UE). This paper continues our previous studies [10,11] aimed at the evolution of jet structure patterns and their dependence on simulation components.

Analysis
We used PYTHIA 8.226 and HIJING++ Monte-Carlo (MC) generators to simulate pp collision events at √ s = 7 TeV [12,13]. Three different PYTHIA tunes were investigated, the Monash 2013 tune with the NNPDF2.3LO PDF set [14,15], the Monash* tune with NNPDF2.3LO [16], and tune 4C with the CTEQ6L1 PDF set [17,18]. Collisions with and without multiple-parton interactions and color reconnection were simulated and compared to each other [10,11]. We only considered particles above the transverse momentum threshold p track T > 0.15 GeV/c. We carried out a full jet reconstruction using the anti-k T , k T , and Cambridge-Aachen jet reconstruction algorithms, which are part of the FASTJET software package [19]. These choices are typical in jet shape analyses [20]. The multiplicity-integrated jet shape studies of CMSwere used as a benchmark for our current multiplicity-differential studies [10,20]. Therefore we chose the jet resolution parameter as R = 0.7 and applied a fiducial cut so that the jets were contained in the CMS acceptance |η| < 1. We did not apply underlying event subtraction in the jet cones. We investigated the jets within the 15 GeV/c < p jet T < 400 GeV/c transverse momentum range, where multiplicity-differential studies on real data are feasible in the near future.
We chose the following two jet shape measures to study the multiplicity-dependent behavior of the jet structure: The Ψ integral jet shape (or momentum fraction) and the ρ differential jet shape (or momentum density fraction) [21,22]. The former one gives the average fraction of the jet transverse momentum contained inside a sub-cone of radius r around the jet axis, the latter one is the momentum profile of the jet, i.e., the average transverse momentum of the particles contained inside an annulus with a δr width and boundaries r a and r b . The exact formulae are given by: respectively, where p i T is the transverse momentum of the selected particle and p jet T is the transverse momentum of the jet. The distance r i of the given particle from the jet axis is calculated as where φ is the azimuthal angle and η is the pseudorapidity. For a better understanding, we noted that the aforementioned observables are connected with the equations: where R is the jet resolution parameter. The differential and integral jet shapes were calculated for the above mentioned p jet T range. The three tunes reproduced the CMS results within statistical errors [10,11].

Results
The multiplicity distributions (multiplicity is defined as the number of charged final state particles in a given collision event) were very similar for the three tunes used in our simulations, as shown in the left panel of Figure 1. However, a significant change is observed when we do not consider the effects of the color reconnection and/or the multiple-parton interactions. Without the CR, the multiplicity distribution becomes wider. If we also switch off the MPI, the multiplicity distribution becomes much narrower compared to the setting where both of them are applied. The width of the multiplicity distribution, however, is not sensitive to the choice of the CR model. In the right panel of Figure 1, we plot the mean values of the multiplicity distributions in events where we reconstructed a jet of a given transverse momentum. We did not exclude non-leading jets in our analysis, but we also investigate the effects of selecting only leading jets later in this paper. As expected, the average event multiplicity grew with the transverse momentum of the selected jet. The integral jet shape with r = 0.2 is shown in Figure 2 to compare the effects of the different tunes and settings. As expected [23], we see similar trends for the multiplicity distributions of the tunes. However, there is a substantial difference between the different MPI and CR settings. The most significant difference in the jet shapes was caused by turning off the MPI, which further supports the current view that the MPI contributions need to be included to correctly describe the jet shapes. We note that we did not subtract the underlying event from the jet. Further investigation is necessary to understand whether the interplay between the UE and the observed hard process is significantly modified by the MPI.  In order to gain a more comprehensive picture, we use the currentdifferential jet shape to investigate the multiplicity dependence of the jet shapes. On the left side of Figure 3, we categorize the events into a high and low multiplicity class. The differential jet shapes of both categories are compared to the multiplicity-integrated momentum density (ρ MI ), computed without any selection in multiplicity. A multiplicity dependence is observed. The jets in the higher multiplicity bin appear to be wider while the jets in the lower multiplicity bin appear to be narrower. This is a trivial multiplicity dependence since the event multiplicity strongly correlates with the jet multiplicity. In our previous work, we canceled this trivial multiplicity dependence by applying a double ratio across the different tunes [10].
In the right panel of Figure 3, we plotted the ρ(r)/ρ MI ratio for the low-and high-multiplicty bins divided by the ρ MI curve. In this way, the difference compared to the ρ MI (black) curve is more visible. We also observed that the ρ(r)/ρ MI curves obtained from the low-and high-multiplicty classes intersect each other at a particular radius inside the jet cone. We also investigated whether the choice of multiplicity bins affects the location of the intersection point by categorizing the events into several smaller event multiplicity bins. We found that all the ρ/ρ MI curves corresponding to these different multiplicity classes intersect each other at approximately the same radius, as can be seen on the left panel of Figure 4. Therefore we named this specific r value R fix . The R fix depends, however, on the transverse momentum of the jet, as shown in the right panel of Figure 4. The R fix (p jet T ) curve goes in a fashion expected by a Lorentz boost and converges to a constant value at higher energies [10].
Since we use a linear interpolation between the points of ρ/ρ MI , the value of R fix will have a slight dependence on the bin width in r (δr), especially at higher p T values where the value of R fix is smaller. To make sure that the results are robust enough, we repeated the analysis with narrower bins (δr = 0.05). In the left panel of Figure 5, we show ρ with this finer r binning, while in the right panel of the same figure we plotted the p jet T dependence of R fix . We can see that the effects by the choice of bin width is rather small and does not change the conclusions, so we do not have to sacrifice the current statistics for finer binnings. Using different settings in PYTHIA changes the physics enough to expect different jet structures after jet reconstruction. All of the k T , Cambridge-Aachen and anti-k T jet clustering algorithms that we investigated reconstructed the jet structures differently since they had different susceptibility to the underlying event. We therefore investigated the effects of varying the physics settings and used three different jet reconstruction algorithms ( Figure 6). We conclude that the presence of R fix was very robust for the different physics selections, as its stability was neither an artifact of the particular choice of jet reconstruction algorithms, nor did it depend strongly on the underlying event. We note that R fix is localized to lower r values, therefore we do not expect significant influence from effects that are mostly visible at higher r values, such as the modification of the UE or higher order corrections to the parton shower [24,25]. We also computed the differential jet shapes using a different MC generator. We selected HIJING++ for this purpose because it implements a mechanism of creating the underlying event and QCD effects on the soft-hard boundary that is different from PYTHIA. HIJING++ uses the PYTHIA jet fragmentation, therefore we do not expect any difference during the later stages. Instead of MPI as implemented in PYTHIA, HIJING++ uses minijet production. Differences at lower momenta may arise below the minijet cutoff. In case of the p T and multiplicity distributions, these effects do not exceed the variation caused by applying different tunes in PYTHIA.
In the case of HIJING++, we used two different PDFsets. The results show quantitatively the same R fix dependence on p jet T within systematic errors, as can be seen in the left panel of Figure 7. Jets originating from different flavors undergo different fragmentation due to both the color-charge effect and the dead-cone effect [26]. We compared flavor-inclusive jets to heavy-flavor (beauty and charm) jets in the right panel of Figure 7. We ensured that heavy flavor comes from the initial stages by only enabling leading order processes in PYTHIA. We compared these to leading and subleading flavor-inclusive jets.
It should be noted that the effect of non-leading jets is negligible on the R fix . Although the overall tendency of heavy-flavor R fix is similar to that observed for light flavor, there is also a clear quantitative difference between heavy and light flavors, which points to a different jet structure. The leading bjets differ for higher p jet T and the leading c jets for lower p jet T . This suggests that the interplay between the mass and color-charge effects is non-trivial and needs further investigation. One possibility for that would be a parallel study of the UE and the fragmentation region corresponding to a heavy-flavor trigger, in a similar manner to [27]. From all the above, we can assume that R fix is a property of the jets that is associated with the final state.
We also repeated the same analysis by selecting only the leading and sub-leading jets. We observed no significant difference compared to the case where all the jets were used (right panel of Figure 7).

Conclusions
We conducted a systematic study on jet structures for pp collisions using MC generators PYTHIA and HIJING++. We investigated the effects of CR and MPI on jet shapes and showed that not considering the effects of CR and MPI causes significant jet shape modification. We introduced a characteristic jet size (R fix ) that depends only on the p T of the jets but is independent of the jet reconstruction algorithms, MC generators, parton density functions, and even the choice of simulation parameters such as CR and MPI [10,11]. We have also shown that the choice of δr does not change our conclusions about R fix . These observations suggest that R fix is an inherent property of the jets and is characteristic to the space-time evolution of the parton shower at a given momentum.
However, R fix does depend on the flavor of the jet. Flavor-dependent jet structure studies may be a way to access mass versus color charge effects that is complimentary to hadron-or jet-production cross-section measurements. We believe that our findings motivate further phenomenology studies as well as cross-checks with real data to gain a deeper understanding on flavor-dependent jet fragmentation. Another direction for future research could be to investigate the effects of MPI on jets without the underlying event. This could be done either by choosing an observable that depends very weakly on the underlying event [9], or by the parallel understanding of the underlying event and the fragmentation region [27].