Abstract

With the framework of the multisource thermal model, we analyze the experimental transverse momentum spectra of various jets produced in different collisions at high energies. Two energy sources, a projectile participant quark and a target participant quark, are considered. Each energy source (each participant quark) is assumed to contribute to the transverse momentum distribution to be the TP-like function, i.e., a revised Tsallis–Pareto-type function. The contribution of the two participant quarks to the transverse momentum distribution is then the convolution of two TP-like functions. The model distribution can be used to fit the experimental spectra measured by different collaborations. The related parameters such as the entropy index-related, effective temperature, and revised index are then obtained. The trends of these parameters are useful to understand the characteristic of high energy collisions.

1. Introduction

In central heavy ion (nucleus-nucleus, A-A) collisions at high energy, quark-gluon plasma (QGP) is believed to create possibly [13], because the environment of high temperature and density is formed. After the formation, QGP experiences the process of hadronization and then produces lots of final-state particles. Meanwhile, at the early stage of collisions, some products such as various jets are produced and interact subsequently with QGP. Because of the interactions between jets and QGP, jets lost their energies when they go through QGP region. Not only lots of identified particles but also various jets can be measured in experiments at high energies. Indeed, in the abundant data on high energy collisions, the data on various jets are one of the most important constituents. We are interested in analyzing the experimental transverse momentum () spectra of various jets, because they can reflect some information of early collisions of participant quarks or partons.

Generally, the spectra of various jets are wider than those of identified particles. In fact, both the spectra of various jets and identified particles cover a wide range. Even if for the later, one may divide the range into low- and high- regions. It is expected that the spectra in low- region are contributed by the soft excitation process, while the spectra in high- region are contributed by the hard scattering process. In some cases, the spectra in low- and high- regions are still complex. One may divide further the low- and high- regions into very low- and low- regions as well as high- and very high- regions, respectively. It is expected that the spectra in different regions can be analyzed by different functions. This means that one needs two-component or even four-component function to fit the wide spectra.

It is known that perturbative Quantum Chromodynamics (pQCD) successfully describes the processes which involve large momentum transfers. In particular for proton-proton (p-p) collisions, multijet production at high is well described if initial- and final-state radiations are considered (see, e.g., ref. [4]). In addition, multiparton interactions may play a role at low [5]. However, pQCD is very complex, which limits its wider applications in high energy proton-nucleus (p-A) and A-A collisions. We hope to use an alternative and thermal like or statistical method to describe uniformly the spectra of various jets and identified particles in both low- and high- regions in p-p, p-A, and A-A collisions at high energies. As the first step and as an example, we consider the two-component function.

There are two methods to superpose the two components in a function [68]. The first method uses a weighted sum for the two components, and there are correlations between the parameters of the two components, though the point of linkage is smooth. The second method uses a step function to link the two components [8], and there is a nonsmooth linkage between the two components, though the parameters are uncorrelative. It is imaginable that more issues will appear if we consider four components in a function. Although the two-component function is widely used in literature, it is not an ideal treatment method, not to mention the four-component function. We hope to use a method to treat the two or four components uniformly. Even a single component function is used to fit the spectra in wide range.

Fortunately, to search for the single component function for the spectra in wide range is possible, because the similarity, universality, or common law is existent in high energy collisions [915]. To search for the single component function, we have tested many potential functions. Finally, we have found that the convolution of two or more revised Tsallis–Pareto-type functions [16, 17] is a suitable choice. For the purpose of doing a convenient description, we call the revised Tsallis–Pareto-type function the TP-like function in our recent work [17] and this paper. In fact, the Tsallis–Pareto-type function is more proper to restrict only to a Tsallis distribution, once one essentially uses the nonextensive statistical mechanics [18]. Because the Tsallis distribution has more than one form of expression, the term of Tsallis–Pareto-type function is used in our work to mention the concrete form.

The application of the convolution of two or more functions is a general treatment method with the framework of the multisource thermal model [7], where the considered distributions are assumed from the contributions of two or more energy sources. The considered distributions include at least the multiplicity, transverse energy, and transverse momentum (transverse mass) distributions. At least two energy sources are considered in the collisions. Three or more energy sources are not excluded if the two energy sources are not enough to fit the spectra. The concrete number of energy sources is determined by the quality of fits and the scenario of physics. This fit methodology is suitable and unified for the analyses of spectra in different rapidity intervals, centrality classes, and collision systems.

In this paper, in the framework of the multisource thermal model, we assume that a projectile participant quark and a target participant quark take part in the production of various jets, and they contribute to the distribution to be the TP-like function [16, 17]. Then, we may use the convolution of two TP-like functions to fit the experimental spectra of various jets. The related data quoted in this paper are from proton(anti)proton (p-p), deuteron-gold (d-Au), gold-gold (Au-Au), proton-lead (p-Pb), and lead-lead (Pb-Pb) collisions, with different selection conditions, over a center-of-mass energy ( or simplified as for p-p collisions) range from 0.2 to 13 TeV.

The remainder of this paper is structured in the following. The formalism and method are described in Section 2. The results and discussion are given in Section 3. Finally, we give the summary and conclusions in Section 4.

2. The Formalism and Method

According to ref. [16], the Tsallis–Pareto-type function which describes empirically the spectra of particles with rest mass can be given by which is a probability density function and C is the normalization constant because . In Eq. (1), as an entropy index-related parameter, is related to the entropy index because . Generally, or means an equilibrium state. If is close to 1 or is large enough, the system is close to an equilibrium state. The free parameter in Eq. (1) is an effective temperature that describes the excitation and expansion degree of the emission source for particles. We call the effective temperature because both the contributions of random thermal motion and flow effect are included.

Equation (1) is not flexible enough in the description of spectra of particles, in particular for the spectra in low- region. Empirically, Eq. (1) can be revised artificially by adding a revised index that is nondimensional as the power index of . Then, we have the TP-like function to be [17] where is the normalization constant which is different from that in Eq. (1). For the purpose of convenience, two normalization constants in Eqs. (1) and (2) are represented by the same symbol , though they may be different. Although one more parameter is introduced, Eq. (2) is more accurate than Eq. (1). In particular, we can obtain Eq. (1) from Eq. (2) if we use .

With the framework of the multisource thermal model [7], we assume that many quarks or partons take part in the collisions. Each quark or parton is regarded as an energy source. For a given particle or jet, two quarks, i.e., a projectile participant (the first) quark and a target participant (the second) quark, play the main role in the production process. Other quarks that take part in the interactions with weak contributions can be neglected. For the two main quarks, the contribution amount or portion () of each quark to is assumed to obey the TP-like function, where and 2 are for the first and second quarks, respectively. The TP-like function obeyed by is [17] where is empirically the constituent mass of the -th participant quark.

The total amount contributed by the two quarks is the convolution of two TP-like functions. That is [17]. or where the functions and are given by Eq. (3) for various jets which are produced by different collisions which are listed in the table in the next section. In most cases, the convolution of two TP-like functions is suitable for the spectra of various jets. Correspondingly, two heavy flavor quarks such as , , or should be considered due to more effective energy being needed. For two light flavor quarks such as , , or , we do not need to consider them due to too less effective energy for the production of various jets.

The method of the convolution of three TP-like functions is similar to that of two TP-like functions. Firstly, we may obtain the convolution of the first two TP-like functions and . Secondly, we may obtain the convolution of and . Alternatively, we may obtain firstly the convolution of the last two TP-like functions and , and then, we may obtain the convolution of and . The same idea can be used for the convolution of more than three TP-like functions. At present, the convolution of a projectile participant quark and a target participant quark is enough to fit the spectra of of various jets. Temporarily, we do not need to consider the convolution of three or more participant quarks.

Because of the introduction of , Eq. (2) is more accurate and flexible than Eq. (1). By using , the spectra in very low- region can be described reasonable. With the framework of the multisource thermal model, the method of the convolution of two or more probability density functions is applicable for not only the spectra of but also the spectra of multiplicity and transverse energy. In our analysis, the free parameters are , , and . The normalization constant is a parameter but not a free parameter. The convolution does not introduce new free parameters but the source number from the collision picture. In the method, to search for the probability density function contributed by a single participant or contributor or source is a key issue. This participant or contributor or source can be quark if we study the spectra of particles or nucleon if we study the spectra of nuclear fragments.

3. Results and Discussion

3.1. Comparison with Data

Figure 1 shows the transverse momentum spectra of different jets produced in (a) p-p, (b) d-Au, and (c) Au-Au collisions with midpseudorapidity (mid-, for Figures 1(a) and 1(c) and for Figure 1(b)), as well as in (d) p-p collisions with non-mid- () at , where and denote the numbers of jets and events, respectively. The symbols are cited from the experimental data measured by the STAR Collaboration [1921]. In Figures 1(a) and 1(c), the high tower (HT) trigger jets were selected. In Figures 1(a)1(c), the jet events were selected using a cone radius () and anti- algorithm [22], where denotes the transverse momentum. In Figures 1(b) and 1(c), the data of d-Au and Au-Au collisions were measured in 0–20% centrality class. In the figure, the curves are our fitted results with Eq. (4). In the fit process, two participant top quarks with constituent mass of 174 GeV/c2 for each one are considered. The values of free parameters (, , and ), normalization constant (), , and number of degree of freedom (ndof) are listed in Table 1 in which the parameter trend will be analyzed and discussed in Section 3.2. One can see that the spectra of different jets are shown to obey approximately the convolution of two TP-like functions. The values of mean excitation and expansion degree (defined by the effective temperature parameter ) seem to not be related to the size of collision system in the error range.

Figure 2(a) presents the spectra of fast jets produced in p-Pb collisions with different centralities at , where on the vertical axis denotes the cross-section. The spectra of charged jets produced in Pb-Pb collisions with different centralities at are presented in Figure 2(b), where on the vertical axis denotes the number of binary nucleon-nucleon collisions. The symbols are cited from the experimental data measured by the ALICE Collaboration [23, 24]. The jet events were selected with a cone radius () and mid- (). The curves are our fitted results with Eq. (4), in which two participant top quarks are considered. The values of , , , , 2, and ndof are listed in Table 1. One can see that the convolution of two TP-like functions describes approximately the experimental data of the mentioned jets. The effective temperature parameter is the same with changing the centrality percentage in p-Pb collisions. And in Pb-Pb collisions, increases slightly with the increase of centrality percentage, i.e., decreases slightly with the increase of centrality itself.

Figures 3(a) and 3(b) present the spectrum of leading jets corresponding to the () channel in p- collisions at . Figures 3(c) and 3(d) present the spectra of jets (-jets) corresponding to the lepton and dilepton channels in p-p collisions at . The symbols are cited from the experimental data measured by the D0 [25, 26], CMS [27], and ATLAS Collaborations [28]. In Figures 3(a) and 3(b), the jet events were selected with a cone radius () and wide range (). In Figures 3(c) and 3(d), the jet events were selected with and , respectively. The curves are our fitted results with Eq. (4), in which two participant top quarks are considered for Figures 3(a)3(c), and two participant bottom quarks with constituent mass of 4.19 GeV/c2 for each one are considered for Figure 3(d). The values of , , , , 2, and ndof are listed in Table 1. One can see that the convolution of two TP-like functions provides an approximate description on the data. The effective temperature parameter obtained from the spectra with the lepton and dilepton channels is almost the same within the error range.

The spectra of (a, b) leading -jets and (c–f) leading jets corresponding to the (a, c, e) channel and (b, d, f) channel in p-p collisions at (a–d) and (e, f) 8 TeV are presented in Figure 4. The symbols are cited from the experimental data measured by the ATLAS Collaboration [2931]. The jet events were selected with . The curves are our fitted results with Eq. (4), in which two participant bottom quarks are considered for Figures 4(a) and 4(b), and two participant top quarks are considered for Figures 4(c)4(f). The values of , , , , 2, and ndof are listed in Table 1. One can see that the convolution of two TP-like functions provides an approximate description on the data. The effective temperature parameter from the channel and the channel is almost the same within the error range.

The reconstructed jet spectra for the leading a 2nd, 3rd, 4th, and 5th order jets in the channels produced in p-p collisions at  TeV are shown in Figure 5. The symbols are cited from the experimental data measured by the ATLAS Collaboration [32]. The jet events were selected with a cone radius () and wide range (). The curves are our fitted results with Eq. (4), in which two participant top quarks are considered. The experimental data are approximately fitted with the convolution of two TP-like functions, and the values of related parameters are given in Table 1. One can see that the effective temperature parameter decreases with the growth of jet order .

Figure 6 presents the spectra of (a) leading -jets, (b) subleading -jets, (c) leading jets, and (d) subleading jets produced in p-p collisions at  TeV. The symbols are cited from the experimental data measured by the CMS Collaboration [27, 33]. The jet events were selected with (a, b) and (c, d) . The curves are our fitted results with Eq. (4), in which two participant bottom quarks are considered for Figures 6(a) and 6(b), and two participant top quarks are considered for Figures 6(c) and 6(d). The experimental data are approximately fitted with the convolution of two TP-like functions, and the values of related parameters are given in Table 1. The values of from the spectra of leading jets are much larger than those from the spectra of subleading jets.

Figure 7 displays the spectra of (a) leading light jets, (b) subleading light jets, (c) leading jets, and (d) subleading jets produced in p-p collisions at  TeV, where denotes the number of observation. The symbols are cited from the experimental data measured by the CMS [34] and ATLAS Collaborations [35]. The jet events were selected with . The curves are our fitted results with Eq. (4), in which two participant charm quarks with constituent mass of 1.27 GeV/c2 for each one are considered for Figures 7(a) and 7(b), and two participant top quarks are considered for Figures 7(c) and 7(d). The experimental data are approximately fitted with the convolution of two TP-like functions, and the values of related parameters are given in Table 1. Once more, the values of from the spectra of leading jets are much larger than those from the spectra of subleading jets.

The spectra of (a, b) small-R selected and (c, d) large-R jets corresponding to the (a, c) and (b, d) channels produced in p-p collisions at  TeV are shown in Figure 8. The symbols are cited from the experimental data measured by the ATLAS Collaboration [36]. The curves are our fitted results with Eq. (4), in which two participant top quarks are considered. The experimental data are approximately fitted with the convolution of two TP-like functions, and the values of related parameters are given in Table 1. One can see that the values of from the spectra of and channels are almost the same within the error range.

The spectra of (a) the leading jets and (b) the subleading jets in Zjj baseline region as well as (c) the leading jets and (d) the forth jets with prefit produced in p-p collisions at  TeV are presented in Figure 9. The symbols are cited from the experimental data measured by the ATLAS Collaboration [37, 38]. The jet events were selected with for Figures 9(a) and 9(b) and for Figures 9(c) and 9(d). The curves are our fitted results with Eq. (4), in which two participant top quarks are considered. The experimental data are approximately fitted with the convolution of two TP-like functions, and the values of related parameters are given in Table 1, in which the normalization is very large due to the fact that it denotes the event number accumulated but not the cross-section or jet number per event. One can see that the values of from the spectra of the leading jets are much larger than those from the spectra of the subleading and forth jets.

From the above comparison with data, one can see that the data are approximately fitted with the convolution of two TP-like functions. In most cases, the quality of the fits is acceptable due to appropriate 2/ndof. In a few cases, the result of the fits is worthy of further improvement due to large 2/ndof. To improve the result, we need more suitable function for the component and/or three or more functions in the convolution. In addition, there are cases in which some parameters do not change, but others suffer significant changes, in that sense is not possible to conclude about the characteristics of high energy collisions. Indeed, as an application of the uniform method, more comparisons with data are needed in the future. Meanwhile, we hope to improve the result of the fits in the future.

3.2. Parameter Trend and Discussion

Although the comparison with data is not perfect, one can also see some trends in most cases. To show the trends of main parameters, Figure 10(a) presents the relation of the effective temperature and the centrality percentage C in Pb-Pb collisions at  TeV. The symbols represent the values of obtained from Figure 2 and are listed in Table 1. The curve is our fit by an exponential function

in which and are in GeV and %, respectively. One can see that increases slightly with the increase of , or is almost the same within the error range when varies. The relation between and renders that QGP formed in central Pb-Pb collisions has less influence on the jet transport. Or in the transport process of jets in QGP in central Pb-Pb collisions, jets lost less energy.

Figure 10(b) presents the relation of the effective temperature and the jet order in p-p collisions at  TeV. The symbols represent the values of obtained from Figure 5 and are listed in Table 1. The curve is our fit by an exponential function in which is in GeV. One can see that decreases with the growth of . This trend is natural due to the fact that the jet with high order corresponds to the source with less excitation degree.

Figure 11 shows the relations of the effective temperature and (a) the size of interacting system, (b) ℓ and di-ℓ channels, (c) μ (μμ) and e (ee) channels, and (d) leading and subleading jets. The symbols represent the values of obtained from the above figures and are listed in Table 1. One can see that seems to not be related to the system size in the error range. This is in agreement with the conclusion from Figure 10(a) in which central collisions correspond to large system, and peripheral collisions correspond to small system. In the error range, different lepton channels show nearly the same effective temperature, which renders nearly the same excitation degree of source. At the same time, the values of from the spectra of leading jets are much larger than those from the spectra of subleading jets, which is the same as the conclusion from Figure 10(b).

As a parameter determining the curvature in middle- region and the extended range in high- region, is related to the entropy index because . In most cases, which is not close to 1 because which is not large. This implies that the source of jets does not stay at the equilibrium state. In a few cases, is large, and is close to 1. This happens coincidentally but not implies that the source of jets stay at the equilibrium state. This situation is different from the source of identified particles. Generally, the source of identified particles stays approximately at the equilibrium or local equilibrium state. These nonequilibrium and (local) equilibrium states do not mean a contradiction between the productions of various jets and identified particles. In fact, the two types of products are produced at different stages of system evolution. Generally, various jets are produced at the initial stage, and the identified particles are produced at the kinetic freeze-out stage. From the initial to kinetic freeze-out stages, the collision system evolves quickly from nonequilibrium to (local) equilibrium states.

As a parameter determining the slope of the curve in low- region, is elastic from negative to positive values. A negative results in a cocked up distribution and a positive results in a falling distribution. In many cases, which means that it is necessary introducing in the Tsallis–Pareto-type function. Due to the introduction of , the revised Tsallis–Paretotype function, i.e., the TP-like function, becomes more flexible. The convolution of two or more TP-like functions is expected to fit more spectra in high energy collisions. The effective temperature from the spectra of leading jets is much larger than that from the spectra of subleading jets because the leading jets undergone more violent scattering. However, parameter does not change usually with the order of jet because determines only the spectra at low . To determine , the spectra at medium and high play the main role. Generally, the parameters related to various jets are not associated with the flow because various jets are produced at the initial stage where the flow is not formed. The flow is associated with identified particles which are finally produced at the kinetic freeze-out stage where the flow is formed.

It is accepted that regarding central A-A collisions at high energies, jets lose considerable energy as they propagate through the hot and dense medium [39]. Due to the loss of energy, the spectra are reduced, and we obtain seemingly a relative small and/or large (small ) comparing with peripheral A-A collisions. However, due to large error, in A-A collisions also shows a nearly invariant trend. In central p-A collisions, jets do not lose considerable energy, and the spectra are not reduced due to similar small participant volume comparing with peripheral p-A collisions. This renders that the main parameters in p-A collisions are independent of centrality, as what appears in A-A collisions. Mostly, the spectra of jets cited from p-A and A-A collisions are not under the same condition; it is hard to compare the parameters directly.

In A-A collisions, the medium effect on the probe of hard process is commonly referred to as jet quenching. The effect has been observed via deviations from a well-calibrated baseline established in the absence of a medium (e.g., in minimum-bias p-p collisions). The generic strategy is illustrated with the nuclear modification factor (), which evaluates the deviation of a single particle inclusive spectra away from the baseline.

is expected to be one in the absence of medium effect. Jet measurements in central A-A collisions at high energies have shown that is dependent, and it is smaller than one. For peripheral A-A collisions, approaches to unity, and it is nearly flat (see, e.g., [40]). An analogous analysis in p-A collisions gives a which is consistent to one (see, e.g., [23]).

It seems that these features of data on are not captured by the parameters presented in Table 1. For instance, the extracted parameters which characterize the spectral shape (, , and ) are the same for central A-A collisions and 0-100% p-A collisions. Probably, this is not a surprise because based on the reduced 2, the proposed function seems not to describe all the data which are analyzed in the paper. This seems to illustrate that the paper also lacks consistent results. In fact, these opinions cannot be obtained because the parameters between A-A and p-A collisions cannot be compared directly. As we know, in Figure 2, the spectra in A-A and p-A collisions are from different kinds of jets. The same parameters for the two collisions are coincidental due to the similar trends of data quoted.

Before Summary and Conclusions, we would like to emphasize the functions of parameters. As what we discussed in our recent work [17], the power index describes flexibly the shapes of spectra at low . From negative to positive , the spectra bend from up to down over a range from 0 to 1 GeV/c. Correspondingly, the spectra at medium and high change higher due to the result of normalization. With the increase of and by fixing and , the spectra become wider. Meanwhile, with the increase of and by fixing and , the spectra become narrower. However, the changes with the increases of and cannot be offset. Anyhow, the introduction of makes the TP-like function more flexible.

We would like to point out that the effective temperature is one of the reported parameters; however, the kinetic freeze-out temperature () should be more important to disentangle the and flow. We have been worked on it for identified particles in our previous work [41]. In fact, that paper reports an increase of with the collision energy, as well as the mass dependence of the . The present work and our previous work are not contradictory. However, we cannot compare them directly due to different stages of considerations. The identified particles are studied at the stage of kinetic freeze-out, while various jets are produced at the stage of initial collisions. In addition, different fit functions which correspond to different “thermometers” are used, which also increases the degree of difficulty for direct comparisons.

4. Summary and Conclusions

We summarize here our main observations and conclusions. (a)The transverse momentum spectra of various jets selected in different conditions and produced in different collisions over an energy range from 0.2 to 13 TeV are fitted by the convolution of two TP-like functions, where the TP-like function is a revised Tsallis–Pareto-type function. The experimental data recorded by various collaborations are approximately fitted by the mentioned convolution(b)From the fit on the spectra of charged jets produced in Pb-Pb collisions at  TeV with different centrality intervals, we know that the effective temperature increases slightly with increasing the centrality percentage, or is almost the same in the error range when the centrality changes. Meanwhile, from the spectra of jets in p-p, d-Au, and Au-Au collisions at 0.2 TeV does not show the size dependence. This is consistent to the nearly independence of on centrality(c)The values of from the spectra of leading jets are much larger than those from the spectra of subleading jets due to the leading jets undergone more violent scattering. As expected, extracted from the reconstructed jets produced in p-p collisions at  TeV decreases with the growth of the jet order. In addition, from the lepton and dilepton channels is almost the same, which means that these jets have common property(d)The parameter determines the curvature in middle- region and the extended range in high- region. Meanwhile, is related to the entropy index because . Generally, is not too large. This means that is not close to 1, and the source of jets does not stay at the equilibrium state. This is different from the source of identified particles which stays approximately at the equilibrium or local equilibrium state(e)The parameter determines the slope of the curve in low- region. A negative results in a cocked up distribution, and a positive results in a falling distribution. Due to the introduction of in the Tsallis–Pareto-type function, the revised function, i.e., the TP-like function, becomes more flexible. The convolution of two or more TP-like functions is expected to have more applications

Data Availability

The data used to support the findings of this study are included within the article and are cited at relevant places within the text as references.

Ethical Approval

The authors declare that they are in compliance with ethical standards regarding the content of this paper.

Disclosure

The funding agencies have no role in the design of the study; in the collection, analysis, or interpretation of the data; in the writing of the manuscript; or in the decision to publish the results.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The work of the first author (Y.M.T.) was supported by Shanxi University. The work of the second author (P.P.Y.) was supported by the China Scholarship Council (Chinese Government Scholarship) under Grant No. 202008140170 and the Shanxi Provincial Innovative Foundation for Graduate Education under Grant No. 2019SY053. The work of the third author (F.H.L.) was supported by the National Natural Science Foundation of China under Grant Nos. 12047571, 11575103 and 11947418, the Scientific and Technological Innovation Programs of Higher Education Institutions in Shanxi (STIP) under Grant No. 201802017, the Shanxi Provincial Natural Science Foundation under Grant No. 201901D111043, and the Fund for Shanxi “1331 Project” Key Subjects Construction.