Advances in High Energy Physics

Volume 2014 (2014), Article ID 896037, 7 pages

http://dx.doi.org/10.1155/2014/896037

Research Article

Pion Spectra in Collisions as a Function of and Event Multiplicity

Priyanka Sett1,2 and Prashant Shukla1,2

¹Nuclear Physics Division, Bhabha Atomic Research Center, Trombay, Mumbai 400 085, India

²Homi Bhabha National Institute, Anushaktinagar, Mumbai 400094, India

Received 27 June 2014; Accepted 15 August 2014; Published 2 September 2014

Academic Editor: Edward Sarkisyan-Grinbaum

Copyright © 2014 Priyanka Sett and Prashant Shukla. 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³.

Abstract.

We study the charged pion transverse momentum spectra in collisions as a function of collision energy and event multiplicity using Tsallis distribution. This study gives an insight of the pion production process in collisions. The study covers pion spectra measured in collisions at SPS energies (6.27–17.27 GeV), RHIC energies (62.4 GeV and 200 GeV), and LHC energies (900 GeV, 2.76 TeV, and 7 TeV). The Tsallis parameters have been obtained and parameterized as a function of . The study suggests that as we move to higher energy more and more hard processes contribute to the spectra. We also study the charged pion spectra for different event multiplicities in collisions for LHC energies using Tsallis distribution. The variation of the Tsallis parameters as a function of event multiplicity has been obtained and their behavior is found to be independent of collision energy.

1. Introduction

The particle spectra measured in hadronic collisions are of utmost interest because of their fundamental nature and simplicity, which allow verifying pQCD [1, 2] calculations and also help to make comprehensive phenomenological studies. The ratios of the particle yields obtained from the measured spectra allow getting the chemical freeze-out conditions, whereas the spectra themselves reflect the conditions at the kinetic freeze-out. The particle spectra provide useful information about the collision dynamics. The low region of the spectrum corresponds to the particles originating from low momentum transfer and multiscattering processes (nonperturbative QCD), whereas the high region comes from the hard-parton-scattering (pQCD) among the initial partons. The transition of this nonperturbative to perturbative dynamics has no sharp boundary, though one can have an estimate from the “” [3]. Extensive [4, 5] and nonextensive [6–10] statistical approaches have been used to characterize particle spectra in terms of thermodynamic variables. Extensive statistics assume thermal and chemical equilibrium of the system at hadronic phase which lead to an exponential distribution of the particle spectra. In experiments, the particle spectra show a power-law behavior at high . This behavior is reproduced by the nonextensive approach with an additional parameter. In recent times, the Tsallis [6] statistical approach is widely used to describe the particle spectra obtained in high-energy collisions with only two parameters, the temperature and , known as nonextensivity parameter which is a measure of temperature fluctuations or degree of nonequilibrium in the system.

The Tsallis distribution gives an excellent description of spectra of all identified mesons measured in collisions at  GeV [11]. In a recent work [12, 13], the Tsallis distribution has been used to describe the spectra of identified charged hadrons measured in collisions at RHIC and at LHC energies. Such an approach has also been applied to the inclusive charged hadron data in recent publications [14, 15]. It has been shown in [12, 16] that the functional form of the Tsallis distribution with thermodynamic origin is of the same form as the QCD-inspired Hagedorn formula [17, 18]. This could be the reason of success of Tsallis distribution in collisions which is a power law typical of QCD hard scatterings. The hardness of the spectra is thus related to and the parameter governs the contribution from soft collisions.

Using the Tsallis phenomenological function, we review and study the charged pion spectra in collisions in a large energy regime, spanning from SPS [19] (6.27 GeV–17.27 GeV) and RHIC [20] (62.4 and 200 GeV) to LHC [21] (900 GeV, 2.76 TeV and 7 TeV) energies. The object of the present work is to study the behaviour of the Tsallis parameters as a function of collision energy. We also study the charged pion spectra for different event multiplicities in collisions for LHC energies. Among all hadrons, pions are chosen because of their abundance in collisions, simple quark structure, and availability of the data at different energies.

2. Formalism

The transverse momentum spectra of hadrons, obtained from different fixed and collider experiments, have shown that the high region of the spectra can be described successfully by the power law,

where is the normalization constant and is the power which determines the shape of the spectra at high . However, the low region of the particle spectra shows an exponential shape and can be described by the Boltzmann-Gibbs [22, 23] statistical approach,

where is the normalization constant, is the particle energy, and is the temperature of the system.

In the early 80s, Hagedorn [17] proposed a phenomenological function which describes the particle spectra for both the higher and lower regions:

where , , and are the fit parameters. The above equation describes an exponential behavior for low and a power-law behavior for high . Consider

The parameter in this equation is often related to the “power” in the “QCD-inspired” quark interchange model [18].

In the late 80s, Tsallis [6] introduced the idea of the nonextensive statistics in place of thermal Boltzmann-Gibbs statistics. This approach includes a parameter , known as nonextensive parameter which quantifies the temperature fluctuation [24] in the system as . The nonextensive statistics assume Boltzmann-Gibbs form in the limit . In Tsallis approach, the Boltzmann-Gibbs distribution takes the form

where is the normalization factor. One can use the relation = at mid-rapidity and in (5) to obtain

where is the normalization factor. Equation (6) can be rewritten as

The value of can be obtained by integrating (7) over momentum space:

Here the quantity is the integrated yield. Equation (6) with the normalization constant takes the form [11]

Larger values of correspond to smaller values of which imply dominant hard QCD point-like scattering. Both and have been interchangeably used in Tsallis distribution [7, 11, 20, 25, 26]. The Tsallis interpretation of parameters as temperature and as nonextensivity parameter is more suited for heavy ion collisions while for collisions Hagedorn interpretation in terms of power and a parameter which controls soft physics processes is more meaningful. Phenomenological studies suggest that, for quark-quark point scattering, [27, 28] and when multiple scattering centers are involved, grows larger.

There are many other forms of (9), which are used often to describe particle spectra in literature; see, for example, [7, 20, 25, 26, 29–31].

3. Results and Discussions

All the studies are performed with (9) and the fit parameters , , and are obtained. The different experiments, energies, rapidity ranges, and particles used in the analysis are summarized in Table 1. For SPS energies only available data is for measured by NA61 Collaboration [19]; for RHIC and LHC energies we use . All the data used are measured in mid-rapidity and are given for unit rapidity. The difference in rapidity range is not expected to affect the behaviour of the spectra. CMS experiment presented [21] transverse momentum spectra for different events classified on the basis of number of true tracks here referred to as track multiplicity of event or simply multiplicity. Each multiplicity class is represented by average number of tracks ().

Table 1: The center of mass energy and rapidity of the data used for the study.

Experiments
Center of mass

energy (GeV)
Rapidity
Particles studied

SPS
6.27, 7.74, 8.76, 12.32, 17.27
0.0–0.2

RHIC
62.4, 200

,

LHC
900, 2760, 700

,

3.1. Tsallis Parameters as a Function of in System

In this analysis all the Tsallis parameters are obtained for charged pion spectra as a function of in system for SPS [19], RHIC [20], and LHC [21] energies. Similar study is available in [12] using RHIC and LHC data and in [13] for SPS and LHC data.

The pion spectra measured in collisions at different are shown in Figure 1 along with with Tsallis fits (9) shown by solid lines. The spectra are scaled by arbitrary factors (given in figure) for visual clarity. In case of RHIC data, we restrict the range to 1.7 GeV/ to have similar range at all energies. It can be noticed that the spectra become harder with increase in which is depictive of occurrence of harder scatterings at higher collision energy. The per degree of freedom values for all the fits are given in Table 2. The values are 1, which is indicative of good fit quality.

Table 2: Values of the for Tsallis fits of pion spectra at different .

6.27 GeV
6.31/12

7.74 GeV
5.80/12

8.76 GeV
10.68/12

12.32 GeV
9.25/12

17.27 GeV
2.65/12

62 GeV
0.74/11

200 GeV
0.48/11

900 GeV
24.33/19

2.76 TeV
5.59/19

7.00 TeV
13.11/19

Figure 1: The invariant yield spectra of charged pions as a function of () for SPS [19] energies 6.27 GeV, 7.74 GeV, 8.76 GeV, 12.32 GeV, and 17.27 GeV, RHIC [20] energies 62.4 GeV and 200 GeV, and LHC [21] energies 900 GeV, 2.76 TeV, and 7 TeV. The solid lines are the Tsallis function (9). The negative pion yields are plotted for SPS energies and for all other energies; average yield for positive and negative pion is plotted.

The parameters and obtained from this analysis are shown in Figures 2 and 3, respectively, as a function of . The variation of as a function of is shown in Figure 4. The parameter decreases with increasing and starts saturating at LHC energies. The value of also reduces slowly from SPS energies to LHC energies. The integrated yield increases 10 times when going from SPS to highest LHC energy.

Figure 2: The variation of the Tsallis parameter for charged pions as a function of . The solid curve represents the parameterization .

Figure 3: The variation of the Tsallis parameter for charged pions as a function of . The solid curve represents the parameterization .

Figure 4: The variation of the integrated yield for charged pions as a function of . The solid curve represents the parameterization .

Larger value of (also larger value of ) suggests that the spectra have contribution from processes involving small momentum transfer arising due to the rescattering, recombination of partons, fragmentation from strings, and so forth, whereas smaller values of are indicative of harder processes being involved in particle production. Thus the spectra at SPS energies have large softer contribution and as the collision energy increases more and more contribution from hard processes is added.

All the three parameters can be parametrized by a function of type

Here , , and for , , , and for , and , , and for . Using the parameterizations for by (10) we get in the limit . The extrapolated values for , and for  TeV are ,  MeV, and .

3.2. Tsallis Parameters as a Function of Multiplicity () for LHC Energies

The Tsallis parameters for charged pion spectra are studied as a function of event multiplicity for different LHC energies 900 GeV and 2.76 and 7 TeV. The event multiplicity data was also studied in a recent work [10] but our analysis and interpretations are different.

The invariant yield spectra corresponding to different multiplicities are fitted with (9) and are shown by the solid black lines in Figure 5 for 900 GeV, in Figure 6 for 2.76 TeV, and in Figure 7 for 7 TeV center of mass energy. The spectra are scaled up for distinctness. The Tsallis distribution describes all the spectra well, shown by the values given in Table 3. The values are little higher for some of the lower multiplicities due to the deviation of first data point in spectra with the curve.

Table 3: Values of the for the Tsallis fits in different event multiplicities.

values for

900 GeV
2.76 TeV
7.0 TeV

7
48.89/18
57.66/18
48.10/18

16
25.21/18
26.38/18
27.76/18

28
9.97/18
11.47/18
22.83/18

40
7.20/18
6.91/18
19.25/18

52
8.00/18
6.08/18
22.52/18

63
9.18/18
6.44/18
26.30/18

75
15.01/18
9.05/18
21.30/18

86

8.16/18
19.24/18

98

11.91/18
23.59/18

109

20.82/18

120

16.85/18

131

19.77/18

Figure 5: The invariant yield spectra of [21], as a function of for collisions at = 900 GeV. The yields are shown for 7, 16, 28, 40, 52, 63, and 75. The spectra are scaled up for clarity by a factor of , where , and 6. The solid lines show the Tsallis fits.

Figure 6: The invariant yield spectra of [21], as a function of for collisions at = 2.76 TeV. The yields are shown for 7, 16, 28, 40, 52, 63, 75, 86, and 98. The spectra are scaled up for clarity by a factor of , where , and 8. The solid lines show the Tsallis fits.

Figure 7: The invariant yield spectra of [21], as a function of for collisions at = 7.00 TeV. The yields are shown for 7, 16, 28, 40, 52, 63, 75, 86, 98, 109, 120, and 131. The spectra are scaled up for clarity by a factor of , where , and 11. The solid lines show the Tsallis fits.

The parameters and obtained from the fits are shown in Figures 8 and 9, respectively, as a function of . The circles, squares, and triangles correspond to the parameter values obtained from data at 900 GeV, 2.76 TeV, and 7 TeV, respectively. It is seen that both the parameters and decrease rapidly and then start saturating with the increase of for all three energies. This variation (of and ) is very similar to the variation which we find as a function of . It means that events with higher multiplicity have larger contribution from hard processes. The value of for high multiplicity events at 7 TeV is which is depictive of production from point quark-quark scattering. The variation of and as a function of can be described by the same curve given in the figure for all three energies and are parameterized by

Here , , and for and , , and for .

Figure 8: The variation of the Tsallis parameter for charged pions as a function of . The variation is shown for 900 GeV by black circles, 2.76 TeV by red squares, and 7.00 TeV by green triangles. The dashed curve represents the parameterization .

Figure 9: The variation of the Tsallis parameter for charged pions as a function of . The variation is shown for 900 GeV by black circles, 2.76 TeV by red squares, and 7.00 TeV by green triangles. The dashed curve represents the parameterization .

The integrated pion yield distribution in different multiplicity classes is shown in Figure 10 for the three LHC energies. The total integrated pion yield for each energy can be obtained by integrating the above distributions over all multiplicity classes. It is noticed that as the energy increases more and more high mutliplicity events are added in the sample with mean of the distribution shifting towards higher .

Figure 10: The variation of the integrated yield for charged pions as a function of . The variation is shown for 900 GeV by black circles, 2.76 TeV by red squares, and 7.00 TeV by green triangles.

4. Conclusion

This work presented the study of the transverse momentum spectra of the charged pions for different collisional energies and also for different event multiplicities (at LHC energies) using Tsallis distribution. The Tsallis parameter decreases with increasing and starts saturating at LHC energies. The value of also reduces slowly from SPS energies to LHC energies. It means that the spectra at SPS energies have large softer contribution and as the collision energy increases more and more contribution from hard processes is added. The integrated pion yield increases with increasing and becomes 10 times when going from SPS to highest LHC energy. The Tsallis parameters are also obtained as a function of event multiplicity for all three LHC energies which can be described by the same curve. The variation of and as a function of multiplicity is very similar to the variation which we find as a function of . It means that events with higher multiplicity have larger contribution from hard processes. The value of for high multiplicity events at 7 TeV is which is depictive of production from point quark-quark scattering. The integrated pion yield distribution for the three LHC energies shows that as the energy increases, more and more high mutliplicity events are added in the sample with mean of the distribution shifting towards higher multiplicity.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

G. Sterman, J. Smith, J. C. Collins (CTEQ Collaboration), et al., “Handbook of perturbative QCD,” Reviews of Modern Physics, vol. 67, article 157, 1995.

J. F. Owens, E. Reya, and M. Glück, “Detailed quantum-chromodynamic predictions for high-pT processes,” Physical Review D, vol. 18, p. 1501, 1978.

R. Blankenbecler, S. J. Brodsky, and J. F. Gunion, “Inclusive processes at high transverse momentum,” Physics Letters B, vol. 42, no. 4, pp. 461–465, 1972.

E. Fermi, “High energy nuclear events,” Progress of Theoretical Physics, vol. 5, no. 4, pp. 570–583, 1950.

A. Andronic, P. Braun-Munzinger, and J. Stachel, “Hadron production in central nucleus-nucleus collisions at chemical freeze-out,” Nuclear Physics A, vol. 772, no. 3-4, pp. 167–199, 2006.

C. Tsallis, “Possible generalization of Boltzmann-Gibbs statistics,” Journal of Statistical Physics, vol. 52, no. 1-2, pp. 479–487, 1988.

T. S. Biro, G. Purcsel, and K. Urmossy, “Non-extensive approach to quark matter,” The European Physical Journal A, vol. 40, pp. 325–340, 2009.

W. M. Alberico and A. Lavagno, “Non-extensive statistical effects in high-energy collisions,” European Physical Journal A, vol. 40, no. 3, pp. 313–323, 2009.

T. S. Biro and G. Purcsel, “Non-extensive Boltzmann equation and hadronization,” Physical Review Letters, vol. 95, Article ID 162302, 2005.

K. Urmossy, “Multiplicity dependence of Hadron spectra in proton-proton collisions at LHC energies and super-statistics,” http://arxiv.org/abs/1212.0260

A. Adare, S. Afanasiev, C. Aidala, et al., “Measurement of neutral mesons in collisions at   GeV and scaling properties of hadron production,” Physical Review D, vol. 83, Article ID 052004, 2011.

P. K. Khandai, P. Sett, P. Shukla, and V. Singh, “Hadron spectra in p+p collisions at rhic and LHC energies,” International Journal of Modern Physics A, vol. 28, no. 16, Article ID 1350066, 2013.

M. Rybczynski and Z. Wlodarczyk, “Tsallis statistics approach to the transverse momentum distributions in p-p collisions,” The European Physical Journal C, vol. 74, p. 2785, 2014.

M. D. Azmi and J. Cleymans, “Transverse momentum distributions in proton-proton collisions at LHC energies and Tsallis thermodynamics,” Journal of Physics G, vol. 41, Article ID 065001, 2014.

J. Cleymans, G. I. Lykasov, A. S. Parvan, A. S. Sorin, O. V. Teryaev, and D. Worku, “Systematic properties of the tsallis distribution: energy dependence of parameters in high energy p-p collisions,” Physics Letters B, vol. 723, no. 4-5, pp. 351–354, 2013.

C. Y. Wong and G. Wilk, “Tsallis fits to spectra and multiple hard scattering in pp collisions at the LHC,” Physical Review D, vol. 87, Article ID 114007, 2013.

R. Hagedorn, “Multiplicities, pT distributions and the expected hadron to quark-gluon phase transition,” La Rivista del Nuovo Cimento, vol. 6, no. 10, pp. 1–50, 1984.

R. Blankenbecler and S. J. Brodsky, “Unified description of inclusive and exclusive reactions at all momentum transfers,” Physical Review D, vol. 10, no. 9, pp. 2973–2992, 1974.

N. Abgrall, A. Aduszkiewicz, Y. Ali, et al., “Measurement of negatively charged pion spectra in inelastic p+p interactions at p lab = 20, 31, 40, 80 and 158 GeV/c,” The European Physical Journal C, vol. 74, p. 2794, 2014.

A. Adare, S. Afanasiev, C. Aidala (PHENIX Collaboration), et al., “Identified charged hadron production in = 200 and 62.4 GeV,” Physical Review C, vol. 83, Article ID 064903, 2011.

S. Chatrchyan, V. Khachatryan, A. M. Sirunyan, et al., “Study of the inclusive production of charged pions, kaons, and protons in pp collisions at = 0.9,2.76, 7 TeV,” The European Physical Journal C, vol. 72, article 2164, 2012.

R. Hagedorn, “Statistical thermodynamics of strong interactions at high energies,” Nuovo Cimento, supplement 3, pp. 147–186, 1965.

D. B. Walton and J. Rafelski, “Equilibrium distribution of heavy quarks in Fokker-Planck dynamics,” Physical Review Letters, vol. 84, p. 31, 2000.

G. Wilk and Z. Wlodarczyk, “Interpretation of the nonextensivity parameter q in some applications of Tsallis statistics and Lévy distributions,” Physical Review Letters, vol. 84, p. 2770, 2000.

B. I. Abelev, J. Adams, M. M. Aggarwal, et al., “Strange particle production in collisions at  GeV,” Physical Review C, vol. 75, Article ID 064901, 2007.

J. Cleymans and D. Worku, “Relativistic thermodynamics: transverse momentum distributions in high-energy physics,” The European Physical Journal A, vol. 48, article 160, 2012.

R. Blankenbecler, S. J. Brodsky, and J. Gunion, “Analysis of particle production at large transverse momentum,” Physical Review D, vol. 12, no. 11, pp. 3469–3487, 1975.

S. J. Brodsky, H. J. Pirner, and J. Raufeisen, “Scaling properties of high pT inclusive hadron production,” Physics Letters B: Nuclear, Elementary Particle and High-Energy Physics, vol. 637, no. 1-2, pp. 58–63, 2006.

T. S. Biro, K. Urmossy, and G. G. Barnafoldi, “Pion and kaon spectra from distributed mass quark matter,” Journal of Physics G, vol. 35, no. 4, Article ID 044012, 2008.

K. Aamodt, N. Abel, U. Abeysekara, et al., “Production of pions, kaons and protons in pp collisions at  GeV with ALICE at the LHC,” The European Physical Journal C, vol. 71, p. 1655, 2011.

J. Adams, M. M. Aggarwal, Z. Ahammed, et al., “Identified hadron spectra at large transverse momentum in p+p and d+Au collisions  GeV,” Physics Letters B, vol. 637, no. 3, pp. 161–169, 2006.