# Modification of Jet Structure in High-Multiplicity pp Collisions due to Multiple-Parton Interactions and Observing a Multiplicity-Independent Characteristic Jet Size

Authors: Zoltán Varga; Róbert Vértesi; Gergely Gábor Barnaföldi
Journal: Advances in High Energy Physics
 (2019)

---

## Abstract

We study the multiplicity dependence of jet structures in pp collisions using Monte Carlo event generators. We give predictions for multiplicity-differential jet structures and present evidence for a nontrivial jet shape dependence on charged hadron event multiplicity that can be used as a sensitive tool to experimentally differentiate between equally well-preforming simulation tunes. We also propose a way to validate the presence and extent of effects such as multiple-parton interactions (MPI) or color reconnection (CR), based on the detection of nontrivial jet shape modification in high-multiplicity events at highpT. Using multiplicity-dependent jet structure observables in various pT windows might also help understanding the interplay between jet particles and the underlying event (UE). We introduce a multiplicity-independent characteristic jet size measure and use a simplistic model to aid its physical interpretation.


---

## Body

## 1. Introduction

High-multiplicity events of small colliding systems at high center-of-mass energies show similar collective features to those observed in events of heavy-ion collisions with comparable multiplicities, such as long-range near-side correlations andvn (“flow”) coefficients [1, 2]. Whether this behavior may be attributed to the presence of a deconfined state in small systems is an open question. However, possible medium-like effects in pp may question the widely exploited assumption that pp collisions are safe to use as a reference for heavy-ion systems. On the other hand, recent studies showed that flow patterns may emerge from features different than hydrodynamics. For instance, radial flow in pp collisions may be explained by pure QCD mechanisms such as multiple-parton interactions (MPI) [3]. Alternative explanations also exist for the observed collectivity in small colliding systems, such as calculations based on CGC initial states with Lund fragmentation [4] or non-Abelian Bremsstrahlung [5]. These scenarios question signatures previously considered as definite signs of the QGP. Recent analyses of pp and p-Pb collisions also show a universal enhancement of heavy-flavour particles, that is usually attributed to MPI and higher gluon radiation associated with short distance production processes [6]. However, we lack the qualitative understanding of these effects. While we cannot expect to observe direct modification of particle yields by any medium created in collisions of small systems (because of the small volume of such a medium), phenomena that act in the soft-hard transitional regime should in principle pose an effect on hard processes as well. A modification in the shapes of developing jets can in principle be accessible by existing experiments.Jet profile measurements in hadron colliders have long been suggested as sensitive probes of QCD parton splitting and showering calculations [7–9], and even as an indicator of the QGP [9]. A recent study suggests to verify a possible existence of a QGP-droplet by measuring properties of jets in association with a Z-boson in ultra-central pp collisions [10]. Another suggestion is to measure the structures of photon-tagged recoil jets in comparison to inclusive jets. This would allow for the detailed understanding of the color-charge effect in fragmentation as well as the elimination of surface bias in jet quenching [11]. In experiment, jet structure observables with full jet reconstruction have been studied in different collisional systems at HERA, the Tevatron, and the LHC [12–19], among others. It is especially important to gain a detailed understanding on multiplicity dependence of the jet structures up to high momenta with the recent advent of machine learning classification techniques in jet studies [20, 21], since these rely heavily on the modelling of parton shower and fragmentation and their connection to the underlying event, in order to avoid possible selection biases.We use the PYTHIA event generator [22] to extensively study the multiplicity-dependent jet shapes, using different tunes and setups of PYTHIA to examine the possible effects of MPI on jets. We provide predictions for pp collisions at s=7 TeV to motivate similar, multiplicity-dependent jet structure measurements at the LHC. In models with string hadronization, the recombination of overlapping color strings (color reconnection or CR) influence fragmentation and are also known to produce collective-like patterns such as radial flow [23]. We investigate the effects caused by the choice of the CR scheme within PYTHIA on the simulated jet structures. As a reference point in our investigations, we decided to use a set of jet structure measurements by the CMS experiment at s=7TeV, carried out in a wide jet momentum range from 15 GeV/c up to 1000 GeV/c [18]. A previous CMS study has investigated multiplicity-differential jet structures, albeit momentum-inclusively with a pTjet>5 GeV/c jet transverse momentum threshold, to understand the influence of underlying events (UE) on jets [19].Our paper is organized as follows. In Section2, we describe our analysis method in detail, and show its validation on CMS data. In Section 3, we present and discuss our results, complemented by simplistic model calculations that aid the understanding of the physics implications. Finally we summarize our results in Section 4.


## 2. Simulation and Analysis

We used the PYTHIA 8.226 [22] event generator to generate random pp collisions at a center-of-mass energy of s=7TeV. We allowed any hard pQCD process, but in order to decrease simulation time we limited the phase space by requiring a certain minimum invariant transverse momentum pT^ of the hardest 2→2 process in an event. We chose pT^>5 GeV/c, pT^>20 GeV/c, pT^>40 GeV/c, and pT^>80 GeV/c for the evaluation of jets with pTjet>15 GeV/c, pTjet>50 GeV/c, pTjet>80 GeV/c, and pTjet>125 GeV/c, respectively. These cutoffs were determined so that they do not have influence on the shape of the reconstructed pTjet spectrum. We simulated 5 million events with each of the settings.Since many physical details cannot be derived from first principles due to our limited understanding of Nature, the MC event generators, including PYTHIA, require extra input parameters. Determining these parameters are far from trivial, and a given set of the parameters are generally sufficient only for reproducing certain experimental data. A given configuration of these parameters, optimized for reproducing experimental, results in certain physical aspects, which are called tunes. Besides the default tune Monash 2013 (Monash) we also investigated two others, theMonash∗ and 4C tunes. The Monash tune, which uses the NNPDF2.3LO PDF set [24], is specifically configured to both e+e- and pp/pp¯ data [25]. Monash∗ (or CUETP8M1-NNPDF2.3LO) is an underlying-event tune based on the Monash tune and was configured to CMS data [26]. The 4C tune is a newer one introduced with PYTHIA version 8.145 [27]. It is based on the tune 2C, but it uses the CTEQ6L1 PDF set [28] and has further changes including a reduced cross section for diffraction and modified multiparton interaction parameters to produce a higher and more rapidly increasing charged pseudorapidity plateau for better agreement with some early key LHC numbers [29]. Using the Monash tune as a starting point we also did investigations where we changed some settings in PYTHIA to directly study their effect on the jet structure. There are continuously developed models of multiple-parton interactions implemented in PYTHIA [30, 31]. To understand the multiplicity-dependent jet modification by MPI we used data samples where we switched this effect on and off.We also investigated different color reconnection schemes provided by PYTHIA, including turning off this feature. Color reconnection is a built-in mechanism in PYTHIA that allows interactions between partons originating in MPI and initial/final state radiations, by minimizing color string length. Since this procedure is quite ambiguous, several models are implemented. The original MPI-based scheme used in PYTHIA 8.226 (that we denote CR0 in the followings) relies on the parton shower-like configuration of the beam remnant. In an additional step, it merges the gluons of a lower-pT MPI system with gluons of a higher-pT MPI system. A newer QCD-based scheme [32] (CR1) relies, however, on the full QCD color configuration in the beam remnant. Then the color reconnection is made by minimizing the potential string energy. The QCD color rules are incorporated in the CR to determine the probability that a reconnection is allowed. This model also allows the creation of junction structures. Besides the above-mentioned CR schemes, a so-called gluon move scheme [33] (CR2) has been implemented to PYTHIA recently, in which gluons can be moved from one location to another so as to reduce the total string length.We carried out a full jet reconstruction including both charged and neutral particles, using three popular algorithms, the anti-kT [34], kT [35, 36], and Cambridge-Aachen [37, 38] algorithms, provided by the FASTJET [39] software package. All of them are sequential clustering algorithms, meaning that the closest particle tracks in momentum space are sequentially merged one-by-one according to the minimum of a distance measure between the particle four-momenta. While all three algorithms are infrared and collinear-safe, in high-multiplicity environments the clusterization outcomes will be rather different. Anti-kT is popular because it is only slightly susceptible to pile-up and underlying events, and it clusterizes hard jets into nearly perfect cones with a resolution parameter R even in high-multiplicity events, in accordance with the general image of how a jet should look like. The other two algorithms are more suitable for jet substructure studies but provide jets of irregular shape that are not uniform in area, especially the kT algorithm, where the area of the jets fluctuates considerably [39]. Similarly to the CMS analysis [18, 19], we selected inclusive jets, with a resolution parameter R=0.7. We considered constituent particles, with a transverse momentum threshold pT,track>0.15GeV/c, at the generator level. Our experience matches earlier findings that the detector effects, after corrections, do not change the simulated jet observables significantly [18]. We examined jets in the pseudorapidity window ηjet<1. We restricted our investigations to the 15GeV/c<pTjet<400GeV/c jet momentum range, where multiplicity-differential studies are feasible in the near future.For the investigation of a possible jet shape modification we analyze the transverse momentum profile of the jets. Two widely used observables are the differential jet shape (ρ) and the integral jet shape (ψ). As our main focus we carry out a detailed study on the multiplicity dependence of ρ, but we also use ψ evaluated at fixed radii as it allows for easier experimental comparison. The differential jet shape describes the radial transverse momentum distribution inside the jet cone and is defined as follows:(1)ρr=1δr1pTjet∑ra<ri<rbpTi, where pTi is the transverse momentum of a particle inside a δr wide annulus with inner radius ra=r-δr/2 and outer radius rb=r+δr/2 around the jet axis and pTjet is the transverse momentum of the whole jet. The distance of a given particle from the jet axis is given by ri=(ϕi-ϕjet)2+(ηi-ηjet)2, where ϕ is the azimuthal angle and η is the pseudorapidity. The integral jet shape gives the average fraction of the jet transverse momentum contained inside a cone of radius r around the jet axis and is calculated as(2)ψr=1pTjet∑ri<rpTi, where the symbols denote the same quantities as for the differential jet shape.As a first step we showed that our simulations reproduce CMS data [18] within uncertainty throughout this range. We show examples in three different pTjet windows in Figure 1. For harder jets, the calculated momentum density distribution gets steeper in the central (small-r) region of the jets, in qualitative accordance with the calculations of [9].Figure 1
Differential jet structureρ(r) measured by the CMS experiment in pp collisions at s=7 TeV [18], compared to different PYTHIA tunes, for 20GeV/c<pTjet<25 GeV/c(left), 50GeV/c<pTjet<60GeV/c(center) and 110GeV/c<pTjet<125GeV/c(right).We investigate the jet structure for different charged hadron multiplicity (Nch) classes. Generally, PYTHIA is known to reproduce multiplicities in LHC data with little differences over a broad pT range. Charged hadron multiplicities at midrapidity are well reproduced by the 4C tune except for a slight discrepancy at very high Nch values [19]. The CUETP8M1 and Z2∗ tunes reproduce pion and kaon average pT versus track multiplicities within errors [40]. The D6T and Z2 tunes show a marginal agreement with the CMS jet-multiplicity data, with about 5% higher predictions than the mean values, flat in pT [18]. We use charged hadron multiplicity at mid-rapidity (referred to as multiplicity in the followings for the sake of simplicity), defined as the number of the charged final state particles with η<1 in a given event. We show the multiplicity distributions in Figure 2 for the jet momentum window 110GeV/c<pTjet<125GeV/c as an example. As shown in the left panel, distributions of the multiplicity are very similar for the different tunes. However, when considering the multiplicity distribution from different settings of the Monash tune, shown on the right panel, a substantial difference can be seen between the settings with and without MPI or CR. Disabling MPI (and CR, which assumes MPI) causes the distribution to shift towards lower values, while keeping a similar shape. Disabling CR only, on the other hand, causes the multiplicity distribution to extend towards higher values. This means that care should be taken when one compares distributions with MPI or CR settings on and off, as it may be biased when the chosen multiplicity class is too wide. We note that multiplicity distributions from different color reconnection schemes do not differ significantly. The pTjet dependence of the mean and RMS values of the multiplicity distribution is compared in Figure 3 for different tunes, as well as for different settings in the case of the Monash tune. The three tunes predict very similar mean and RMS values throughout the pTjet range. While the means of the 4C and Monash tunes overlap, Monash∗ predicts slightly lower multiplicities. The pTjet dependence of multiplicity distributions is a key observable for validating the strength of multiple-parton interaction and color reconnection effects in Monte Carlo models, as both MPI and CR have a grave effect on the distributions. Switching off MPI causes a downward shift of about 15 to 25 in mean Nch at any pTjet, or almost a factor of three at low pTjet values, while switching off CR alone causes a somewhat less drastic increase of about 10 to 20 in mean Nch counts. The three examined color reconnection models provide very similar distributions except for at low pTjet, where the gluon-moving CR scheme predicts slightly narrower multiplicity distributions and a somewhat lower mean value.Figure 2
Distributions of event multiplicity for jets in the110GeV/c<pTjet<125GeV/c window, compared for the Monash, Monash∗, and 4C tunes(left) and for the Monash tune with the CR0, CR1, CR2, and settings as well as CR turned off and MPI turned off(right).Figure 3
Mean and RMS values of the event multiplicity distributions for jets depending onpTjet, compared for different tunes(left) and settings(right). The uncertainties of the mean and RMS values are smaller than the symbol size.


## 3. Results

In this section we present our results and consider the possible physical implications. As a first step we computeρ(r) similarly to Figure 1, but this time while dividing up the data into two multiplicity classes, Nch≤50 and Nch>50, respectively. We see a multiplicity dependence in the jet shapes in Figure 4. Namely, the jets contain a higher fraction of their transverse momentum closer to their axis and a lower fraction further away from their axis in the case of low multiplicity. For high multiplicity the jet shape behaves in the opposite way. This is a trivial, expected multiplicity dependence arising from two reasons. The first one is that event multiplicity is correlated with jet multiplicity, resulting in a higher fraction of narrow jets in low-Nch events. The second reason is the UE background, which affects the jet structure more at higher r values, and its effect is stronger in the case of high-Nch events.Figure 4
Simulated differential jet structureρ(r) compared in multiplicity-integrated events (black), low-multiplicity (Nch≤50,red), and high-multiplicity events (Nch>50,green), for 20GeV/c<pTjet<25GeV/c(left), 50GeV/c<pTjet<60GeV/c(center), and 110GeV/c<pTjet<125GeV/c(right).Measurements by the CMS experiment [19] that compare five multiplicity classes within the range 10<Nch≤140 and reconstruct jets at momenta pTjet>5GeV/c saw a remarkable difference between low and high-multiplicity ρ(r) at low r values. We can make the same observation at relatively low pTjet values (Figure 4 left panel). Dividing ρ(r) for both the high- and low-multiplicity classes with the multiplicity-integrated ρMI(r) (no condition on Nch), shown in Figure 5, highlights this trend. The curves are much more apart at small pTjet for low r values, while there is relatively little difference between different pTjet windows at high r. This suggests that jets in high-multiplicity events contain much more contribution from the soft regime, and soft physics is selected by a lower choice of momentum range.Figure 5
Ratioρ/ρMI of differential jet structure in low-multiplicity (Nch≤50,red) and high-multiplicity events (Nch>50,green) over multiplicity-integrated events, for 20GeV/c<pTjet<25GeV/c(left), 50GeV/c<pTjet<60GeV/c(center), and 110GeV/c<pTjet<125GeV/c(right).Jets in low-multiplicity events are on average narrower than in high-multiplicity events, hence the correspondingρ(r)/ρMI(r) ratio is above unity, while for high-multiplicity events this ratio is below unity. At high r values, where UE tracks give a nonnegligible contribution especially in the high-multiplicity events, the situation is just the opposite. In between there is a point at a given r value where the two curves intersect each other at unity, meaning that at that radius the jets are just average. In Figure 5 we see three examples in different pTjet windows and we can observe that the intersection point is dependent on the jet momentum. This is not unexpected since harder jets are narrower and UE is significant already at smaller radii. To have a closer look at this behavior we evaluate ρ(r)/ρMI(r) in a more refined division of data with seven multiplicity classes in the range 1≤Nch≤250. We find that the curves intersect unity at virtually the same location for a given pTjet value. This statement holds even if we compare different PYTHIA tunes and MPI or CR settings, as shown on the examples in Figure 6 for the Monash and 4C tunes as well as the Monash tune without color reconnection.Figure 6
Ratioρ/ρMI of differential jet structure in several multiplicity classes (see legend) over multiplicity-integrated events, for jets within the 110GeV/c<pTjet<125GeV/c window. In theleft hand side,central andright panels show events generated using the Monash tune, the 4C tune, and the Monash tune without CR, respectively.In the left and center panels of Figure7 we plot the r dependence of the intersection radius with respect to the jet transverse momentum for different tunes, as well as the different settings of the Monash tune. The intersection radius is computed using a linear interpolation between the two nearest points of ρ(r)/ρMI(r), and its uncertainty is estimated by taking both the high and the low-multiplicity classes, moving the points to the upper and lower edge of their error bars in both cases, and determining the maximal and minimal values of the intersection radius from these cases. We observe that for all tunes and settings that we tested, the intersection radii are consistent within uncertainties for any chosen pTjet value. There is additional uncertainty on the obtained intersection radius stemming from the linear interpolation between finite, δr=0.1 wide bins. In order to estimate this, we repeated the analysis with the three tunes in δr=0.05 wide bins. While the statistical fluctuations increase, the points move a maximum of 4% upwards or 28% downwards in a strongly correlated manner (see Figure 7). Nevertheless, the overall shape of the curves remains very similar and statistically consistent between different tunes point-by-point. Therefore, we suggest that the intersection radius Rfix=rρ=ρMI be considered as a characteristic jet size measure specific for a given jet transverse momentum. We note that the value of Rfix should not be compared to the resolution parameter R, that is typically chosen so that most of the jet momentum is contained within the radius R. In contrast, Rfix is defined as a radius where the momentum density of the jet from events of any multiplicity is just like in the average jet, and substantial fraction of jet momentum falls towards smaller as well as towards larger radii.Figure 7
Evolution of the intersection radiusRfix with the jet momentum pTjet, compared for several PYTHIA tunes(left), different settings(center), and different jet reconstruction algorithms in the case of the Monash tune(right). The range indicated with the shaded band is the absolute uncertainty arising from the choice of bin width (see text for details).Jet shapes depend on the jet reconstruction algorithm, so we investigated whether the observed stability of the intersection radius can be an artifact of the jet reconstruction algorithm itself. Besides the anti-kT algorithm we first used, we have reprocessed all the data with using the kT and the Cambridge-Aachen algorithms. We do not find a significant difference beyond the statistics-driven fluctuations between data reconstructed by different clusterization algorithms in any of the tunes or MPI/CR settings. In the right panel of Figure 7 we show a comparison of Rfix(pTjet) for the Monash tune with the three different algorithms.Jets are more collimated with increasing transverse momentum. In a simple picture this can be linked to Lorentz-boost, i.e., the momentum of the initiating parton in the laboratory system. ThepTjet-dependent evolution of Rfix may also be explained by the Lorentz-boost that high-pT jets undergo (see the illustration in the left panel of Figure 8). In order to gain an effect-level understanding, we use a simplistic model. We consider particles radiating from a point in a plane with momenta of equal absolute value p0. We boost these particles along the axis perpendicular to their plane, with a certain momentum pboost. The resulting particles will form a cone around the boost axis in the lab system, representing our “jet”. In the right panel of Figure 8 we see that the resulting size of the “jet”, Rcone, depends on pboost in a qualitatively similar manner to how the intersection radius Rfix depends on pTjet. This attests to the assumption that the universal behavior can, at least partially, be understood by the narrowing by Lorentz-boost of high-pT jets. It is to be noted however, that our model is a simplistic treatment of jet narrowing, and therefore an exact agreement is not expected. A proper description requires a perturbative QCD-based approach where the narrowing is a consequence of the running of the strong coupling constant αs with kT [7–9]. Since there is no angular cut-off that would limit the jet sizes at low momenta, Rcone is allowed to blow up at low pboost values. Also, one cannot expect real jets to go below a certain size because after certain point the clustering algorithms will be driven by the presence of UE. This may explain the apparent convergence of the Rfix curves to a finite value at high pT. As mentioned before, Rfix at high-pT is also influenced by the choice of δr. A particularly interesting question is whether Rfix can be generalized to the larger and more complex systems produced in heavy-ion collisions. To see that, one would need to do simulations in heavy-ion collisions and verify the outcome with data. In case Rfix is representative of the jet size in heavy-ion collisions, it would provide a handy observable for the exploration of medium modification of jets.Figure 8
Left:cartoon illustrating the narrowing of a jet cone by Lorentz-boost.Right:cone radii according to the simplistic boosted-cone model for different p0 assumptions (see text for details).We plotted in Figure9 the differential jet structure for various PYTHIA tunes in a particular pTjet window to compare them in the low and high-multiplicity regions. In the right panel of Figure 9 we take the differential jet shapes for the above-mentioned low- and high-multiplicity classes and divide them with each other to highlight the differences for the different tunes. As expected, jets from low-multiplicity events have a more steeply falling momentum density distribution than the ones from high-multiplicity events, which is also reflected in a falling ratio. However, there are also certain significant differences between the selected tunes that are beyond this trivial effect.Figure 9
Differential jet structure compared for several PYTHIA tunes, for110GeV/c<pTjet<125GeV/c. Jet ρ(r) from events with a lower multiplicity of Nch≤25(left) is compared to jets from events with a higher multiplicity of 100<Nch≤250(center). The ratio of the two is shown in theright panel.To highlight the differences between the jet structures from different tunes we compute the double ratio,(3)DRr=ρlow/ρhighρlow/ρhighref.tune, where we divide the former ratio of the high- and low-multiplicity classes with the very same ratio calculated for the Monash tune. After the trivial effect is gone, a rather sizable effect in the order of a factor of 2 can be seen for both the 4C and the Monash∗ tunes, with respect to the Monash as the reference tune. The right panel of Figure 10 shows the same calculations for the 4C tune, for several different choices of high and low-multiplicity class pairs. In this selected pTjet range all show similar structures, and generally the effect is larger when the separation in multiplicity is larger. It is very important to note that these curves are derived from statistically independent samples, hence cannot be explained by fluctuations. Since on Figures 9 and 10 we calculate ratios of binned data without a bin center correction, we tested its possible effect by decreasing the bin size from δr=0.1 to δr=0.05. We did not find any difference beyond statistical uncertainties.Figure 10
Left: the double ratio with Monash as the reference tune for the same multiplicity classes.Right: the double ratio shown for different selections of high- and low-multiplicity classes (see text in figure).In order to understand the dependence of the effect onpTjet, one might wish to describe the deviations for each pTjet value with a single number. Therefore we compute the squared sum of the bin-by-bin deviations of the double ratio from the Monash tune, i.e.,(4)RSD=∑0<ri<RDRri-12 at a given pTjet. In Figure 11 we show the results for different tunes as well as for different selections of multiplicity class pairs. Again we see a rather parallel behavior of the 4C and Monash∗ tunes (or, in other words, the Monash tune is the one that deviates from these two). The behavior versus pTjet is nontrivial with several minima and maxima, and is not easily explained without taking into account peculiar details of each tune. However, one sees again a very strong correlation between curves of different multiplicity selections calculated independently from each other, and that the amplitude strongly depends on the separation between the low- and high-multiplicity classes. Thus we can conclude that the multiplicity-dependent analysis of jet structures in a wide pTjet range has the potential of evaluating the goodness of tunes that otherwise preform equally well in several tests.Figure 11
Left: square sum of the bins in the DR(r)-1 diagram for the 4C (green) and Monash∗ (red) tunes with respect to the Monash reference tune, depending on the pTjet, for 0<Nch≤25 as low-multiplicity, and 80<Nch≤100 as high-multiplicity selections.Right: square sum of the bins in the DR(r)-1 diagram for the 4C tunes with respect to the Monash reference tune, depending on the pTjet, for various low- and high-multiplicity selections, as listed in the legend. (The pTjet range is restricted to omit parts with large fluctuations.)Lastly, in Figure12, the integral jet shape is plotted in terms of the multiplicity at high momenta, 180GeV/c<pTjet<200GeV/c. In theleft panel, where tunes are compared, there is no observable effect in the integral structure between the tunes Monash, Monash∗, and 4C. We present the effects of different MPI and CR settings on the integral jet structure in theright panel of Figure 12. Different color reconnection schemes do not lead to significant differences, but there is a slight deviation at high Nch values when color reconnection is turned off. However, the lack of MPI causes a significant difference within the same multiplicity class, that grows approximately linearly with Nch, which suggests that the MPI has a strong influence on the jet structure, especially at high Nch values. It is to be noted that the effect is less significant in case of lower pTjet windows and in case of larger r values. At lower multiplicties, MPI and CR cause little difference in the integrated jet shape. That the ψ(r) values at high Nch are lower in the case the MPI is turned off, means the jets are more concentrated in a narrow cone. This can be understood by a higher relative fraction of soft tracks coming from the UE in case when there is no MPI, compared to the MPI case with the same multiplicity where there is a more relevant contribution from tracks that come from the jet itself. Note that the points in Figure 12 are not at the bin centers, but they are placed to represent the weight of the Nch distribution in a given bin, to eliminate the possible bias stemming from the different Nch distributions within multiplicity classes.Figure 12
Evolution of the integral jet structureψ(r=0.2) with event multiplicity Nch, at 180GeV/c<pTjet<200GeV/c, with a jet resolution parameter R=0.7. Several PYTHIA tunes(left) and settings(right) are compared, as explained in the text. The points are placed according to the weight of the distribution in each multiplicity class.Understanding the observed dependence of the integrated jet structure on the multiplicity needs further analysis supported by experimental data. The above observation, if compared to real data, may provide a control over the extent of MPI effects. Further studies are needed to identify MPI/CR effects and separate them from the UE, also using other observables that are less sensitive to the UE.


## 4. Summary

We performed a novel jet shape analysis ins=7 TeV pp collisions to explore the multiplicity and pTjet-dependence of differential and integrated jet structure observables. We used several models implemented in the PYTHIA 8.226 event generator. We demonstrated that the simulations describe CMS data, and we gave predictions for the jet structure observables in several multiplicity classes, over a wide momentum range. We found that there is a given radius Rfix where jet momentum density is independent of multiplicity. This radius is insensitive to the choice of simulation settings (choice of tune, presence, and modelling of MPI and CR) within the investigated model class and even of jet clustering algorithms, and its pT-dependence qualitatively follows a Lorentz-boost curve. These observations suggest that Rfix is an inherent property of jets that is characteristic to the spatial development of the parton shower at a given momentum.We compared the multiplicity dependence of jet structure variables for three popular PYTHIA tunes as well as different MPI and CR models in severalpT bins. We found that the evolution of the differential jet structure ρ(r) with multiplicity significantly differs in several pTjet ranges for the Monash, Monash∗, and 4C tunes. The shape of the difference is nontrivial in pTjet, but persistent through all tested choices of multiplicity selections. With this we demonstrated that the multiplicity-dependent analysis of jet momentum profiles can differentiate among otherwise well-established models. This lack of understanding may have grave consequences on studies based on classification by jet properties. Our observation highlights the need of extending multiplicity-dependent jet structure measurements such as in [19] to higher pTjet regimes.We also see that the integrated jet structure variableψ(r=0.2) shows a rather different Nch-dependence when MPI are turned off. This attests to the important role of multiple-parton interactions in higher multiplicity events and the need for their detailed understanding in order to develop accurate models in jet physics.



---

*Source: 6731362-2019-03-18.xml*