# What Can We Learn from (Pseudo)Rapidity Distribution in High Energy Collisions?

Authors: Fu-Hu Liu; Tian Tian; Jian-Xin Sun; Bao-Chun Li
Journal: Advances in High Energy Physics
 (2014)

---

## Abstract

Based on the (pseudo)rapidity distribution of final-state particles produced in proton-proton (pp) collisions at high energy, the probability distributions of momenta, longitudinal momenta, transverse momenta (transverse masses), energies, velocities, longitudinal velocities, transverse velocities, and emission angles of the considered particles are obtained in the framework of a multisource thermal model. The number density distributions of particles in coordinate and momentum spaces and related transverse planes, the particle dispersion plots in longitudinal and transverse coordinate spaces, and the particle dispersion plots in transverse momentum plane at the stage of freeze out in high energy pp collisions are also obtained.


---

## Body

## 1. Introduction

Due to the complexity of high energy collisions, it is impossible to measure all quantities and characteristics. Instead, some quantities and characteristics can be measured by experiment [1], while others can be obtained from the theoretical explanations based on the experimental results.In experiments, high energy proton-proton (p
p), proton-nucleus (p
A), deuteron-nucleus (d
A), and nucleus-nucleus (A
A) collisions have been performed at available accelerators such as the Super Proton Synchrotron (SPS) [2–4] and colliders such as the Relativistic Heavy Ion Collider (RHIC) [5–8] and the Large Hadron Collider (LHC) [9–12]. Some quantities and characteristics such as the pseudorapidity (or rapidity) distributions, transverse momentum distributions, azimuthal correlations, flow effects, and some other interesting information have been obtained [13].A lot of models were introduced in the past several decades. Different phenomenological mechanisms including initial interactions, intermediate processes, and final-state statistical laws have been proposed [14–19]. In a workshop held a few years ago at the CERN Theory Institute [20], many models reported their predictions for the heavy ion program at the LHC based on the explanations of the experimental results at the RHIC. Recently, some of the models are tested by p
p, p
A, and A
A collisions at the LHC. Meanwhile, these models are further tested in p
p, p
A, d
A, and A
A collisions at the RHIC and SPS.In the past years, we suggested a multisource thermal model [21–23] to describe the (pseudo)rapidity distributions, multiplicity distributions, transverse momentum distributions, azimuthal distributions, flow effects, and so forth. Some interesting quantities such as the temperature, speed of sound, and number of sources are obtained. It is noticed that some quantities and characteristics are related to others. Usually, based on a few distributions, other distributions can be obtained due to some modelling assumptions.Very recently, the NA61/SHINE Collaboration reported the negatively charged pion productions in inelasticp
p collisions at 20–158 GeV/c at the SPS [1]. The rapidity spectrums are parameterized by the sum of two Gaussian functions symmetrically displaced with respect to midrapidity. It is interesting for us to give a further test for the multisource thermal model by using the rapidity spectrums in p
p collisions at SPS momenta (energies). According to the rapidity spectrums, we hope some other distributions can be obtained by the model.In this paper, in the framework of the multisource thermal model [21–23], we analyze the rapidity spectrums in p
p collisions measured by the NA61/SHINE Collaboration at the SPS [1]. Then, a series of other distributions are obtained due to the description of rapidity spectrums.


## 2. The Model and Calculation Method

In the multisource thermal model [21–23], we assume that many emission sources are formed in high energy collisions. In rapidity space in the laboratory or center-of-mass reference frame, these sources distribute at different rapidity y
x and form a target cylinder in rapidity interval [
y
T
min
⁡
,
y
T
max
⁡
] and a projectile cylinder in rapidity interval [
y
P
min
⁡
,
y
P
max
⁡
], respectively. Because of the symmetry in p
p collision, we have [
y
T
min
⁡
=
-
y
P
max
⁡
] and [
y
T
max
⁡
=
-
y
P
min
⁡
].In the rest frame of a given source, we can use different formalizations to describe the production of particles. For example, we can use the relativistic ideal gas model to give a description for the source. Then, we have the momentump
′ distribution of particles to be
(1)
f
p
′
(
p
′
)
=
1
n
d
n
d
p
′
=
C
p
′
2
exp
⁡
(
-
p
′
2
+
m
0
2
T
)
,
where C
=
(
1
/
m
0
2
T
)
(
1
/
K
2
(
m
0
/
T
)
) is the normalization constant, K
2
(
m
0
/
T
) is the modified Bessel function of order 2, m
0 is the rest mass of the considered particle, n is the particle number, and T is the temperature parameter of the source. The particles are assumed to emit isotropically in the rest frame of the source. Then, the distributions of emission angle θ
′ and the azimuth ϕ
′ can be given by f
θ
′
(
θ
′
)
=
(
1
/
2
)
sin
θ
′ and f
ϕ
′
(
ϕ
′
)
=
1
/
(
2
π
), respectively.Letp
z
′, p
T
′, p
x
′, p
y
′, m
T
′, E
′, y
′, and η
′ denote, respectively, the longitudinal momentum, transverse momentum, x-component of momentum, y-component of momentum, transverse mass, energy, rapidity, and pseudorapidity of a considered particle in the rest frame of the considered source. In the Monte Carlo calculation, let r
1,2
,
3 denote, respectively, random numbers distributed evenly in [
0,1
]. The value of p
′ satisfies
(2)
∫
0
p
′
‍
f
p
′
(
p
′
)
d
p
′
<
r
1
<
∫
0
p
′
+
d
p
′
‍
f
p
′
(
p
′
)
d
p
′
.For other quantities, we haveθ
′
=
arctan
[
2
r
2
(
1
-
r
2
)
/
(
1
-
2
r
2
)
] if the right side of the equation is greater than 0; otherwise θ
′
=
arctan
[
2
r
2
(
1
-
r
2
)
/
(
1
-
2
r
2
)
]
+
π. Further, we have some quantities such as ϕ
′
=
2
π
r
3, p
z
′
=
p
′
cos
⁡
θ
′, p
T
′
=
p
′
sin
θ
′, p
x
′
=
p
T
′
cos
⁡
ϕ
′, p
y
′
=
p
T
′
sin
ϕ
′, m
T
′
=
p
T
′
2
+
m
0
2, E
′
=
p
′
2
+
m
0
2
=
p
z
′
2
+
p
T
′
2
+
m
0
2
=
p
z
′
2
+
m
T
′
2, and y
′
≡
(
1
/
2
)
ln
⁡
[
(
E
′
+
p
z
′
)
/
(
E
′
-
p
z
′
)
].In the laboratory or center-of-mass reference frame, letr
5,6 denote, respectively, random numbers distributed evenly in [
0,1
]; we have the related quantities y
x
=
(
y
T
max
⁡
-
y
T
min
⁡
)
r
5
+
y
T
min
⁡ if r
6
<
0.5, or y
x
=
(
y
P
max
⁡
-
y
P
min
⁡
)
r
5
+
y
P
min
⁡ if r
6
>
0.5. Further, we have y
=
y
x
+
y
′, p
T
=
p
T
′, p
x
=
p
x
′, p
y
=
p
y
′, m
T
=
m
T
′, p
z
=
m
T
sinh
⁡
y, p
=
p
T
2
+
p
z
2, E
=
m
T
cosh
⁡
y, and θ
=
arctan
(
p
T
/
p
z
) if p
z
>
0, or θ
=
arctan
(
p
T
/
p
z
)
+
π if p
z
<
0. The pseudorapidity η
≡
-
ln
⁡
tan
(
θ
/
2
). Then, the probability distributions of the considered quantities and their correlations can be obtained. The space density distributions of p and p
T are ρ
(
p
)
=
[
1
/
(
4
π
p
2
n
)
]
(
d
n
/
d
p
) and ρ
(
p
T
)
=
[
1
/
(
2
π
p
T
n
)
]
(
d
n
/
d
p
T
), respectively. Besides, some exchangeable quantities in the above expressions are p
z
=
p
z
′
cosh
⁡
y
x
+
E
′
sinh
⁡
y
x, E
=
p
z
′
sinh
⁡
y
x
+
E
′
cosh
⁡
⁡
y
x, and y
≡
(
1
/
2
)
ln
⁡
[
(
E
+
p
z
)
/
(
E
-
p
z
)
].Further, the velocity and corresponding longitudinal and transverse components of the particle can be given byβ
=
p
/
E, β
z
=
p
z
/
E, and β
T
=
p
T
/
E, respectively. If we assume that the time interval from initial collision to the stage of freeze-out is t
0, then we have the space coordinates x
=
t
0
β
T
cos
⁡
ϕ, y
=
t
0
β
T
sin
ϕ, z
=
r
z
=
t
0
β
z, r
T
=
x
2
+
y
2
=
t
0
β
T, and r
=
t
0
β. Thus, the probability distributions of the considered quantities and their correlations can be obtained. The space density distributions of r and r
T are ρ
(
r
)
=
[
1
/
(
4
π
r
2
n
)
]
(
d
n
/
d
r
) and ρ
(
r
T
)
=
[
1
/
(
2
π
r
T
n
)
]
(
d
n
/
d
r
T
), respectively.


## 3. Comparison and Extraction

The NA61/SHINE Collaboration has measured the rapidity distributions and other features of negatively charged pions produced in 20–158 GeV/cp
p collisions at the SPS [1]. It is shown that the rapidity spectrums are parameterized by the sum of two Gaussian functions symmetrically displaced with respect to midrapidity:
(3)
d
n
d
y
=
n
2
σ
0
2
π
·
[
exp
⁡
(
-
(
y
-
y
0
)
2
2
σ
0
2
)
+
exp
⁡
(
-
(
y
+
y
0
)
2
2
σ
0
2
)
]
,
where n, σ
0, and y
0 are the normalization constant, distribution width for one Gaussian, and peak position parameter, respectively. The values with errors of the three parameters at different momenta are obtained by the NA61/SHINE Collaboration [1].We can compare directly our modelling results with the experimental distribution function of the NA61/SHINE Collaboration [1]. In Figure 1, the rapidity and pseudorapidity distributions of negatively charged pions produced in (a) 158, (b) 80, (c) 40, and (d) 20 GeV/c p
p collisions are displayed. The two-dotted curves represent the rapidity distribution range of the NA61/SHINE experiment [1]. The solid curves which fall into the experimental distribution range are our modelling results obtained by using the multisource thermal model. The modelling results in fact give almost whole superpositions to the two-Gaussian functions which are not presented in the panels. Correspondingly, the modelling results on the pseudorapidity distributions are given by the dashed curves. In the calculation, we have used the Monte Carlo technique and χ
2-testing method. The values of parameters T, y
P
max
⁡
(
=
-
y
T
min
⁡
), and y
P
min
⁡
(
=
-
y
T
max
⁡
) for Figures 1(a), 1(b), 1(c), and 1(d) are listed in Table 1. One can see that the model describes the experimental rapidity distributions of negatively charged pions produced in p
p collisions at the SPS. The pseudorapidity distribution has a lower peak and wider range than those of the rapidity distribution. The temperature increases and the rapidity shift y
P
max
⁡ increases with increase of the incident momentum, and the rapidity shift y
P
min
⁡ does not change with the momentum.Table 1
Values of parameters corresponding to the solid curves in Figure1.
Figure
T (GeV)
y
P
max
⁡
y
P
min
⁡
Figure1(a)
0.150
±
0.012
2.00
±
0.02
-
0.21
±
0.01
Figure1(b)
0.145
±
0.012
1.71
±
0.02
-
0.21
±
0.01
Figure1(c)
0.141
±
0.011
1.43
±
0.01
-
0.21
±
0.01
Figure1(d)
0.138
±
0.011
1.20
±
0.01
-
0.21
±
0.01Rapidity and pseudorapidity distributions of negatively charged pions produced in (a) 158, (b) 80, (c) 40, and (d) 20 GeV/cp
p collisions. The two-dotted curves represent the rapidity distribution range of the NA61/SHINE experiment [1]. The solid and dashed curves are our modelling results on the rapidity and pseudorapidity distributions.
(a)
158 GeV/cp
p, π
-
(b)
80 GeV/cp
p, π
-
(c)
40 GeV/cp
p, π
-
(d)
20 GeV/cp
p, π
-By using the above parameter values, we can obtain some interesting features for other quantities. Figure2 presents the distributions of (a) p, (b) p
z, (c) p
T, and (d) m
T
-
m
0 of negatively charged pions produced in p
p collisions at the SPS. The solid, dotted, dashed, and dotted-dashed curves correspond to the results at 158, 80, 40, and 20 GeV/c, respectively. The circles in Figure 2(b) represent the contribution of the target cylinder to p
z distribution at 158 GeV/c, while the contribution of the projectile cylinder is symmetrical. For the p, p
T, and m
T
-
m
0 distributions, both the contributions of the target and projectile cylinders are the same and equal to a half of the distributions. Similarly, Figure 3 presents the distributions of (a) E, (b) β, (c) β
z, and (d) β
T of the negatively charged pions. The circles in Figure 3(c) represent the contribution of the target cylinder to β
z distribution at 158 GeV/c, while the contribution of the projectile cylinder is symmetrical. For the E, β, and β
T distributions, both the contributions of the target and projectile cylinders are the same and equal to a half of the distributions.Distributions of (a)p, (b) p
z, (c) p
T, and (d) m
T
-
m
0 of negatively charged pions produced in p
p collisions at the SPS. The solid, dotted, dashed, and dotted-dashed curves correspond to the results at 158, 80, 40, and 20 GeV/c, respectively. The circles in Figure 2(b) represent the contribution of the target cylinder to p
z distribution at 158 GeV/c.
(a)
(b)
(c)
(d)The same as Figure2, but showing the distributions of (a) E, (b) β, (c) β
z, and (d) β
T. The circles in Figure 3(c) represent the contribution of the target cylinder to β
z distribution at 158 GeV/c.
(a)
(b)
(c)
(d)The space angular distributions of negatively charged pions produced inp
p collisions at the SPS are displayed in Figure 4(a). The solid, dotted, dashed, and dotted-dashed curves represent the results at 158, 80, 40, and 20 GeV/c, respectively. The open circles are the contribution of the target cylinder at 158 GeV/c, while the contribution of the projectile cylinder is symmetrical. For comparison, the result of an isotropic emission is given by the closed circles. To compare the pseudorapidity distributions at different incident momenta, the dashed curves in Figures 1(a)–1(d) are plotted together in Figure 4(b) by the solid, dotted, dashed, and dotted-dashed curves, respectively. The contribution of the target cylinder at 158 GeV/c is shown by the open circles, while the contribution of the projectile cylinder is symmetrical. For comparison, the result of an isotropic emission is given by the closed circles. One can see that the final-state particles are far from the isotropic emission.Distributions of (a)θ and (b) η of negatively charged pions produced in p
p collisions at the SPS. The solid, dotted, dashed, and dotted-dashed curves correspond to the results at 158, 80, 40, and 20 GeV/c, respectively. The closed circles show the results of the isotropic emission. The open circles represent the contribution of the target cylinder at 158 GeV/c.
(a)
(b)The space density distributions of (a)r, (b) r
T, (c) p, and (d) p
T are presented in Figure 5, where the value of t
0 is taken to be 2 fm/c for convenience. The solid, dotted, dashed, and dotted-dashed curves correspond to the results at 158, 80, 40, and 20 GeV/c, respectively. One can see that the maximum ρ
(
r
) appears at the maximum r and the minimum ρ
(
r
) appears in the region close to 0. The maximum values of ρ
(
r
T
), ρ
(
p
), and ρ
(
p
T
) appear at 0 and the minimum values appear at the maximum coordinate or momentum (transverse momentum) values.Space density distributions of (a)r, (b) r
T, (c) p, and (d) p
T of negatively charged pions produced in p
p collisions at the SPS. The solid, dotted, dashed, and dotted-dashed curves correspond to the results at 158, 80, 40, and 20 GeV/c, respectively.
(a)
(b)
(c)
(d)The dispersion plots inx
(
y
)
-
z, y
-
x, and p
x
(
p
y
)
-
p
z planes for negatively charged pions produced in (a) 158, (b) 80, (c), 40, and (d) 20 GeV/c p
p collisions are presented in Figures 6, 7, and 8, respectively. The circles and squares correspond, respectively, to the contributions of the target and projectile cylinders in 500 events. We can see that the density in large |
z
| and small |
x
(
y
)
| region is greater than those in other regions. The density in small |
p
x
(
p
y
)
| and |
p
z
| region is greater than those in other regions. The contribution points of the two cylinders are mixed in small longitudinal momentum region. The particle number densities in coordinate space and momentum space increase obviously with increase of the incident momentum.Dispersion plots inx
(
y
)
-
z planes for negatively charged pions produced in (a) 158, (b) 80, (c), 40, and (d) 20 GeV/c p
p collisions. The circles and squares correspond, respectively, to the contributions of the target and projectile cylinders in 500 events.
(a)
(b)
(c)
(d)The same as Figure6, but showing the dispersion plots in y
-
x planes.
(a)
(b)
(c)
(d)The same as Figure6, but showing the dispersion plots in p
x
(
p
y
)
-
p
z planes.
(a)
(b)
(c)
(d)


## 4. Conclusions

From the above discussions, we obtain the following conclusions.(a) The multisource thermal model is used to describe the rapidity distributions of negatively charged pions produced inp
p collisions at SPS energies (momenta). A target cylinder and a projectile cylinder are assumed to form in the collisions. For each source, a relativistic ideal gas model is used. The calculated results are in agreement with the rapidity distributions measured by the NA61/SHINE Collaboration [1]. In our calculation, the main parameters are the temperature, the maximum rapidity shift, and the minimum rapidity shift. The temperature and the maximum rapidity shift increase with increase of the incident momentum, and the minimum rapidity shift does not change with the momentum.(b) According to the parameter values obtained from the rapidity distributions, the other distributions are obtained in the framework of the multisource thermal model. The other distributions include the probability distributions of momenta, longitudinal momenta, transverse momenta (transverse masses), energies, velocities, longitudinal velocities, transverse velocities, space angles, and pseudorapidities, as well as the number density distributions of particles in coordinate and momentum spaces. The final-state particles are far from the isotropic emission. The maximumρ
(
r
) appears at the maximum r and the minimum ρ
(
r
) appears in the region close to 0. The maximum values of ρ
(
r
T
), ρ
(
p
), and ρ
(
p
T
) appear at 0 and the minimum values appear at the maximum coordinate or momentum (transverse momentum) values.(c) According to the parameter values obtained from the rapidity distributions, the particle dispersion plots inx
(
y
)
-
z, y
-
x, and p
x
(
p
y
)
-
p
z planes are also obtained in the framework of the model. The density in large |
z
| and small |
x
(
y
)
| region is greater than those in other regions, and the density in small |
p
x
(
p
y
)
| and |
p
z
| region is greater than those in other regions. The contributions of the target cylinder and the projectile cylinder are mixed in small longitudinal momentum region.



---

*Source: 863863-2014-04-08.xml*