Research Article  Open Access
Comparing Multicomponent Erlang Distribution and Lévy Distribution of Particle Transverse Momentums
Abstract
The transverse momentum spectrums of finalstate products produced in nucleusnucleus and protonproton collisions at different centerofmass energies are analyzed by using a multicomponent Erlang distribution and the Lévy distribution. The results calculated by the two models are found in most cases to be in agreement with experimental data from the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). The multicomponent Erlang distribution that resulted from a multisource thermal model seems to give a better description as compared with the Lévy distribution. The temperature parameters of interacting system corresponding to different types of finalstate products are obtained. Light particles correspond to a low temperature emission, and heavy particles correspond to a high temperature emission. Extracted temperature from central collisions is higher than that from peripheral collisions.
1. Introduction
The Relativistic Heavy Ion Collider (RHIC) in USA and the Large Hadron Collider (LHC) in Switzerland have been built to study properties of matters formed in highenergy collisions. These collisions are helpful in understanding particles’ statistical behavior, production process, interaction mechanism, and related phenomenon in highdensity and hightemperature states. Such highenergy collisions offer us opportunities to carry out investigations not only on the Higgs and dark matter [1–3], but also on particle statistical behavior at ultrahigh energy.
Transverse momentum spectrums of finalstate products are very important in highenergy collisions. Many models have been introduced to describe the transverse momentum spectrums of different finalstate products [4]. From the spectrums, one can extract temperature parameter of interacting system. It is expected that temperature parameters extracted from different particle spectrums are different due to different emission stages and regions in collisions. Although we can compare nuclear temperature with classical temperature, they have different physical meanings.
Temperature parameter in highenergy collisions is very important. Generally speaking, temperatures of interacting system at initial, intermediate, and final states are different [5]. Since these temperatures cannot be measured directly, it may, therefore, be interesting to find out an indirect method for obtaining the temperature of the interesting system. Traditionally, temperature can be extracted from measurements of spectrum slopes or double isotopic ratios at lower energies [5, 6]. In some cases, we cannot obtain absolute values of concerned temperature parameters, but relative values corresponding to different particle spectrums.
Multicomponent Erlang distribution derived from multisource thermal model [7, 8] has been applied to collisions in relatively low energy region comparing to RHIC and LHC energies. Energy spectrum of nuclear fragments, multiplicity distribution of charged particles, neutron number distribution of isotope in nuclear fragments, transverse momentum (mass) spectrum of relativistic particles, and so forth were described by the multicomponent Erlang distribution. The Lévy distribution has been also applied to transverse momentum spectrums in highenergy collisions [9–11]. We can study transverse momentum spectrums by using the multicomponent Erlang distribution [7, 8] or the Lévy distribution [9–11] to extract temperature parameters.
In this paper, the transverse momentum spectrums of different finalstate products produced in nucleusnucleus and protonproton collisions at RHIC and LHC energies are studied with the two distributions mentioned above. Temperature parameters are then obtained from fitting experimental data of the STAR, CMS, and ALICE Collaborations.
2. Formalism
The multicomponent Erlang distribution can be derived from the multisource thermal model [7, 8]. In the model, many emission sources of particles are assumed to form in high energy collisions. According to different interaction mechanisms, geometrical relations, selected conditions, or other factors, the emission sources are divided into groups. Source number in the th group is assumed to be . Each source contributes finalstate distribution to be an exponential function. We have the transverse momentum spectrum contributed by the th source in the th group to be where denotes number of finalstate particles and is mean transverse momentum contributed by the sources in the th group.
The transverse momentum spectrum contributed by the th group is the fold of exponential functions; that is, This is an Erlang distribution. In final state, the spectrum contributed by the groups can be written as where is the relative weight contributed by the th group. It is a multicomponent Erlang distribution.
Considering relative contribution of the th group, we have the mean transverse momentum of finalstate particles to be Generally, reflects the mean excitation degree of the emission sources and can be used to describe the source temperature parameter . As in the ideal gas model in which obeys Rayleigh distribution, we have where denotes rest mass and is mean Lorentz factor of considered particles. Further, where and are mean energy and mean momentum of considered particles, respectively. On other hand, as the inverse slope parameter, can be used to describe excitation degree of the emission sources. We define as a new temperature parameter.
The Lévy distributions appear in many branches of physics, mathematics, biology, economy, computer science, and other areas, where the distribution forms may be different in different branches and the scale of fluctuations may be characterized by long tails and an asymptotic powerlawlike behavior. The Lévy distributions are a generalization of the Gaussian distribution. They are similar to the Gaussian distribution and remain stable under the convolution. In fact the Lévy distributions are quite general distributions which contain Gaussian and Cauchy distributions as special cases [12].
Let be the nonextensive parameter. As a probability distribution, the Lévy distribution is commonly the following powerlike distribution [9]: which is just a oneparameter generation of the BoltzmannGibbs exponential formula with , where is the normalization constant and is in the range from 0 to infinity. For the transverse momentum distributions in highenergy collisions, we use directly the function form of Lévy distribution [10]: where is the slope parameter and represents the scale of possible fluctuation in . The parameter can be regarded as the temperature parameter in the Lévy distribution.
3. Comparisons with Experimental Data
The transverse momentum spectrums of finalstate particles produced in CuCu and AuAu collisions at RHIC energy ( TeV) are shown in Figure 1. The symbols represent experimental data of the STAR Collaboration [11]. The solid and dashed curves represent results calculated by the multicomponent Erlang distribution with or 2 and the Lévy distribution, respectively. The results for different centralities (0–10%, 20–30%, and 40–60% in Cu + Cu, as well as 0–5%, 20–40%, and 60–80% in Au + Au) and also for different particles (, , , and in Cu + Cu, as well as and in Au + Au) in central rapidity range () are displayed in different panels. For the sake of convenience, the spectrums are for various centrality bins, with each being scaled by the amount indicated in the legend. The parameter values used in the calculations are shown in Table 1 along with values of per degree of freedom () and extracted temperatures. One can see that the concerned experimental data are described approximately by the two distributions. Light particles correspond to a lower temperature comparing with the heavy particles. The multicomponent Erlang distribution seems to give a better description than the Lévy distribution. We can use the new distribution, the multicomponent Erlang distribution, to describe the transverse momentum spectrums.

(a)
(b)
(c)
(d)
(e)
(f)
In Figure 2, we give the transverse momentum spectrums of leading and subleading jets produced in PbPb and pp collisions at the LHC energy ( or TeV), where the selections of leading and subleading jets can be found in experimental material [13]. The symbols represent experimental data of the CMS Collaboration [13]. The solid and dashed curves represent results calculated by the multicomponent Erlang distribution and the Lévy distribution, respectively. Figures 2(a), 2(b), and 2(c) correspond to different selected conditions shown in the panels, where , , anti, , and Flow denote the integral luminosity, azimuth, sequential recombination algorithm for high particle, resolution parameter, and particle flow, respectively. The parameter values used in the calculations are shown in Table 2 with values of and extracted temperatures. It is again observed that the two distributions describe approximately the concerned experimental data.

(a)
(b)
(c)
In the Lévy distribution, we need to know the rest mass of finalstate product. However, the rest mass of jet is uncertain. In fact, we regarded as a parameter in Figure 2. To see dependence of jet spectrum on in the Lévy distribution, we redraw the Lévy distribution curves for different values in Figure 3, where the same experimental data [13] as those cited in Figure 2 are used. Different values of correspond to different results shown in the figure by different types of curves. All the parameter values with values of are given in Table 3. We see that the temperature extracted from a given jet spectrum decreases with increase of the jet mass and is greater than that extracted from particle spectrums. It should be noticed that the jet mass is the total mass of particles in the jet. For a jet with a given total transverse momentum, a larger mass corresponds to more particle number. Then, the transverse momentum per particle will be smaller, which renders a lower temperature.

(a)
(b)
(c)
In Figure 4, another data sample on spectrums of leading and subleading jets produced in PbPb collisions at TeV is analyzed. The symbols represent experimental data of the CMS Collaboration [13]. The solid and dashed curves represent results calculated by the multicomponent Erlang distribution and the Lévy distribution, respectively. The values of all the parameters along with the values of are given in Table 2. We see that except for a few points the two distributions describe approximately the experimental data. Different spectrums corresponding to different (dijet imbalance parameter) ranges can be described by the same distribution which reflects a common law in the spectrums.
(a)
(b)
(c)
(d)
The spectrums of charged jets produced in PbPb collisions at TeV is given in Figure 5. The symbols represent experimental data of the ALICE Collaboration [14]. The solid and dashed curves represent results calculated by the multicomponent Erlang distribution and the Lévy distribution, respectively. All the parameter values along with values of and extracted temperatures are given in Table 2. One can see that both the distributions describe approximately the experimental data, and the former one gives a better description than the latter one.
The spectrums of charged particles (which can be approximately regarded as ) produced in TeV PbPb collisions in different centrality bins with different multiplications are shown in Figure 6(a). Meanwhile, the spectrums of , , , and produced in central (0–5%) PbPb collisions at the same energy are shown in Figure 6(b). The symbols represent experimental data of the ALICE Collaboration [14, 15] measured in the pseudorapidity range of . The solid and dashed curves represent results calculated by the multicomponent Erlang distribution and the Lévy distribution, respectively. Corresponding to Figures 6(a) and 6(b), the parameter values with values of and extracted temperatures are given in Tables 4 and 5, respectively. One can see that the multicomponent Erlang distribution describes well the spectrums in all the cases. The Lévy distribution describes well the spectrums in some cases, and in other cases it describes approximately the mean trends of the spectrums.


(a)
(b)
Figures 7(a), 7(b), and 7(c) show, respectively, spectrums of finalstate particles , , and produced in TeV PbPb collisions in different centrality bins with different multiplications. Selected condition for is rapidity being in the range of . For the sake of comparison, the results for and produced in 2.76 TeV pp collisions are also given in Figures 7(a) and 7(b), respectively. The symbols represent experimental data of the ALICE Collaboration [16, 17]. The solid and dashed curves represent results calculated by the multicomponent Erlang distribution and the Lévy distribution, respectively. All the parameter values with values of and extracted temperatures are given in Tables 5 (for Figures 7(b) and 7(c)) and 6 (for Figure 7(a)), respectively. One can see that the multicomponent Erlang distribution describes well the spectrums in all the cases. The Lévy distribution describes well the spectrums in some cases, and in other cases it describes approximately mean trends of the spectrums.

(a)
(b)
(c)
The transverse momentum spectrums of and as well as inclusive electrons produced in inelastic pp collision at 7 TeV are given in Figures 8(a) and 8(b), respectively. Experimental data measured by the ALICE Collaboration [15, 18] are shown by the symbols. Results calculated by using the multicomponent Erlang distributions and the Lévy distributions are shown by the solid and dashed curves, respectively. The parameter values used in the calculation are listed in Table 4. We see that both distributions describe approximately the experimental data.
(a)
(b)
The transverse momentum spectrums of , , and ; , , and ; , , , , and produced in pp collisions at 0.9 TeV are displayed in Figures 9(a), 9(b), and 9(c), respectively. The symbols represent experimental data of the ALICE Collaboration [19, 20]. The solid and dashed curves represent results calculated by using the multicomponent Erlang distribution and the Lévy distribution, respectively. The related parameter values are given in Table 5. One can see that both the two distributions describe approximately the experimental data.
(a)
(b)
(c)
In Figure 10, the transverse momentum spectrum of charged particles (which can be approximately regarded as ) produced in nonsingle diffractive (NSD) pp collisions at 0.9 TeV is presented. The symbols represent experimental data measured in the pseudorapidity range of by the ALICE Collaboration [19]. The solid and dashed curves represent results of the multicomponent Erlang distribution and the Lévy distribution, respectively. The related parameter values are given in Table 4. One can see that both the two distributions describe approximately the experimental data.
To see dependences of temperature ( and ) on centrality and , in Figures 11 and 12, we plot different values of and taken from Tables 1–6. The related impacting types, , centralities, and finalstate products are shown in the figures. Figures 11(a), 11(b), 11(c) and 11(d) as well as 11(e) and 11(f) correspond to dependence on centrality for particle productions at 0.2 and 2.76 TeV and jet production at 2.76 TeV, respectively. Figure 12 corresponds to dependence on for particle productions at RHIC and LHC energies. One can see that the extracted temperature for light particles is less than that for heavy particles. Central collisions or high correspond to a relative high temperature. The multicomponent Erlang distribution extracts a relatively high temperature comparing to the Lévy distribution. Besides, from the parameter tables (Tables 1, 2, and 4–6) and (8), one can easily obtain values of which show similar behaviors as those of .
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
4. Conclusions and Discussions
The transverse momentum spectrums of finalstate products produced in highenergy collisions are analysed by using the multicomponent Erlang distribution and the Lévy distribution. In most cases, both the distributions are approximately in agreement with experimental data at RHIC and LHC energies. The multicomponent Erlang distribution seems to give a better description as compared to the Lévy distribution. Although the Lévy distribution is well known to give the transverse momentum spectrums, the multicomponent Erlang distribution gives a new method to describe the transverse momentum spectrums.
The temperature parameters of interacting system corresponding to different types of finalstate products are extracted from transverse momentum spectrums. Light particles correspond to a low temperature emission, and heavy particles correspond to a high temperature emission. For a jet with a given transverse momentum, larger mass corresponds to larger particle number and lesser transverse momentum per particle, which renders a lower temperature. Central collisions or high correspond to a relative high temperature. The multicomponent Erlang distribution extracts a relatively high temperature comparing with the Lévy distribution.
System size dependence of the hadronic spectrums is well described by the two modeling distributions in the present work. We see some correlations between the parameter values and system size. Particularly, the extracted temperature increases with increase of the system size from pp collision to CuCu and AuAu (PbPb) collisions at the same . This renders that the excitation degree of the interacting system increases with increase of the system size. Comparing with light nuclear collisions, a participant nucleon in heavy nuclear collisions takes part in more binary collisions, and more energy per nucleon deposits in heavy nuclear collisions.
It is well known that most of the hadrons in low transverse momentum region are produced in the process dominated by soft interaction, whereas the hadrons with high transverse momentums are produced in the process dominated by hard partonparton scattering. According to the discussions in the present work, the first group of sources in the multicomponent Erlang distribution corresponds generally to the soft interaction, and the second or third group of sources corresponds to the hard scattering. The Lévy distribution does not distinguish the transverse momentum regions of soft interaction and hard scattering.
Although there are more or less differences in both the modeling distributions for the observed transverse momentum spectrums, the multicomponent Erlang distribution and the Lévy distribution describing approximately the transverse momentum spectrums in different systems render that there are some common laws or universality in multihadron production [21, 22], even in general probability distributions. For example, the multicomponent Erlang distribution is also used to describe the probability distributions of some plant seed masses and sizes [23], and the Lévy distribution has more other applications [12, 24]. We are interested in searching new applications of the two distributions.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was partly finished at the State University of New York at Stony Brook, USA. One of the authors (FuHu Liu) thanks Professor Dr. Roy A. Lacey and the members of the Nuclear Chemistry Group of Stony Brook University for their hospitality. The authors acknowledge the supports of the National Natural Science Foundation of China (under Grant no. 10975095, no. 11247250, and no. 11005071), the China National Fundamental Fund of Personnel Training (under Grant no. J1103210), the Open Research Subject of the Chinese Academy of Sciences LargeScale Scientific Facility (under Grant no. 2060205), the Shanxi Scholarship Council of China, and the Overseas Training Project for Teachers at Shanxi University.
References
 P. W. Higgs, “Broken symmetries, massless particles and gauge fields,” Physics Letters, vol. 12, no. 2, pp. 132–133, 1964. View at: Google Scholar
 P. Huang, N. Kersting, and H. H. Yang, “Extracting MSSM masses from heavy Higgs boson decays to four leptons at the CERN LHC,” Physical Review D, vol. 77, no. 7, Article ID 075011, 14 pages, 2008. View at: Publisher Site  Google Scholar
 L. Maiani, G. Parisi, and R. Petronzio, “Bounds on the number and masses of quarks and leptons,” Nuclear Physics B, vol. 136, no. 1, pp. 115–124, 1978. View at: Google Scholar
 S. Abreu, S. V. Akkelin, J. Alam et al., “Heavy ion collisions at the LHC—last call for predictions,” Journal of Physics G, vol. 35, no. 5, Article ID 054001, 185 pages, 2008. View at: Google Scholar
 P. Zhou, W.D. Tian, Y.G. Ma, X.Z. Cai, D.Q. Fang, and H. W. Wang, “Influence of statistical sequential decay on isoscaling and symmetry energy coefficient in a gemini simulation,” Physical Review C, vol. 84, no. 3, Article ID 037605, 4 pages, 2011. View at: Publisher Site  Google Scholar
 C.W. Ma, J. Pu, Y.G. Ma, R. Wada, and S.S. Wang, “Temperature determined by isobaric yield ratios in heavyion collisions,” Physical Review C, vol. 86, no. 5, Article ID 054611, 6 pages, 2012. View at: Google Scholar
 F.H. Liu, “Unified description of multiplicity distributions of finalstate particles produced in collisions at high energies,” Nuclear Physics A, vol. 810, no. 1–4, pp. 159–172, 2008. View at: Publisher Site  Google Scholar
 F.H. Liu and J.S. Li, “Isotopic production cross section of fragments in ^{56}Fe+p and ^{136}Xe(^{124}Xe)+Pb reactions over an energy range from 300A to 1500A MeV,” Physical Review C, vol. 78, no. 4, Article ID 044602, 13 pages, 2008. View at: Publisher Site  Google Scholar
 G. Wilk and Z. Włodarczyk, “Interpretation of the nonextensivity parameter q in some applications of Tsallis statistics and Lévy distributions,” Physical Review Letters, vol. 84, no. 13, pp. 2770–2773, 2000. View at: Publisher Site  Google Scholar
 J. Adams, M. M. Aggarwal, Z. Ahammed et al., “K(892)* resonance production in Au+Au and $p+p$ collisions at $\sqrt{{s}_{NN}}=200$ GeV,” Physical Review C, vol. 71, no. 6, Article ID 064902, 15 pages, 2005. View at: Google Scholar
 H. Agakishiev, M. M. Aggarwal, Z. Ahammed et al., “Strangeness enhancement in Cu+Cu and Au+Au collisions at $\sqrt{{s}_{NN}}=200$ GeV,” http://arxiv.org/abs/1107.2955. View at: Google Scholar
 W. Ebeling, M. Y. Romanovsky, and I. M. Sokolov, “Velocity distributions and kinetic equations for plasmas including Lévy type power law tails,” Contributions to Plasma Physics, vol. 49, no. 10, pp. 704–712, 2009. View at: Publisher Site  Google Scholar
 Y. Yilmaz, “Jet fragmentation functions measured in Pb–Pb collisions with CMS,” Journal of Physics G, vol. 38, no. 12, Article ID 124157, 4 pages, 2011. View at: Publisher Site  Google Scholar
 M. van Leeuwen, “High${p}_{T}$ results from ALICE,” in Proceedings of the Hadron Collider Physics symposium (HCP '11), Paris, France, November 2011, http://arxiv.org/abs/1201.5205. View at: Google Scholar
 M. Floris, “Identified particles in pp and Pb–Pb collisions at LHC energies with the ALICE detector,” Journal of Physics G, vol. 38, no. 12, Article ID 124025, 8 pages, 2011. View at: Publisher Site  Google Scholar
 H. Appelshäuser, “Particle production at large transverse momentum with ALICE,” Journal of Physics G, vol. 38, no. 12, Article ID 124014, 8 pages, 2011. View at: Google Scholar
 R. Preghenella, “Transverse momentum spectra of identified charged hadrons with ALICE detector in Pb–Pb collisions at $\sqrt{{s}_{NN}}=2.76$ TeV,” in Proceedings of the Europhysics Conference on High Energy Physics (EPSHEP '11), Grenoble, France, July 2011, http://arxiv.org/abs/1111.0763. View at: Google Scholar
 S. Masciocchi, “Inclusive electron spectrum from heavyflavour decays in protonproton collisions at $\sqrt{{s}_{NN}}=7$ TeV measured with ALICE at LHC,” Nuclear Physics A, vol. 855, no. 1, pp. 432–435, 2011. View at: Publisher Site  Google Scholar
 M. Kowalski, “First results on charged particle production in alice experiment at LHC,” Acta Physica Polonica B, vol. 42, no. 34, pp. 859–866, 2011. View at: Publisher Site  Google Scholar
 K. Aamodt, A. Abrahantes Quintana, D. Adamová et al., “Strange particle production in protonproton collisions at $\sqrt{s}=0.9$ TeV with ALICE at the LHC,” The European Physical Journal C, vol. 71, no. 3, Article ID 1594, 24 pages, 2011. View at: Google Scholar
 E. K. G. Sarkisyan and A. S. Sakharov, “Relating multihadron production in hadronic and nuclear collisions,” The European Physical Journal C, vol. 70, no. 3, pp. 533–541, 2010. View at: Publisher Site  Google Scholar
 E. K. G. Sarkisyan and A. S. Sakharov, “Multihadron production features in different reactions,” AIP Conference Proceedings, vol. 828, pp. 35–41, 2006. View at: Google Scholar
 S. H. Fan and H. R. Wei, “Multicomponent Erlang distribution of plant seed masses and sizes,” Journal of the Korean Physical Society, vol. 61, no. 11, pp. 1918–1921, 2012. View at: Google Scholar
 T. J. Kozubowski and K. Podgórski, “Distributional properties of the negative binomial Lévy process,” Probability and Mathematical Statistics, vol. 29, no. 1, pp. 43–71, 2009. View at: Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2014 HuaRong Wei et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP^{3}.