Simulation of ATLAS Collaboration’s Measurement Results of Jet T p Correlations in Pb + Pb Collisions at NN s = 2.76 TeV

In this article, a model based on Glauber Monte Carlo is built to simulate the procedure of jet quenching in Quark-Gluon Plasma (QGP). In this model, energy loss of jets in QGP is parametrized by two quantities: path length of jets in QGP, L, and the initial transverse momentum of the jet T p . The path length of each branch of the jet are labeled L 1 and L 2 . As input data, original jet energy data of p + p collisions were obtained from CMS measurement. After being processed by our model, simulated Pb + Pb jet energy data could be given and were compared to the data of ATLAS’s measurement in Ref. [1]. Distributions of ( J where , also noted as “frequency”, are presented as a function of T p 1 and collision centrality. As the final result, two different forms of energy loss formula were found, both of which have good adaptation to certain centrality and T p ranges.


INTRODUCTION
Relativistic heavy-ion collisions can produce a medium made up of free quarks and gluons, known as the Quark-Gluon Plasma (QGP). Hard scattering processes occurring in these collisions produce high transverse momentum ( T p ) partons called jets which can be one or multiple. Due to the reason that jets carry color charges, they strongly couple with the QGP, thus lose energy as they propagate through the medium, resulting in the phenomenon of "jet quenching", which is widely concerned and studied by various scholars [2][3][4].
A number of research groups have already published their measurement results related to energy of jets produced in these heavy-ion collisions at various collision energy levels. Some experimental data were Open Access Natural Science compared to theoretical results [5, 6] while others remained unexplained by theory. In this article, a model based on Glauber Monte Carlo is built to simulate the measurement results of dijet events T p correlations in Ref. [1] by the ATLAS collaboration, successful theoretical simulation of which has not been published before. As the final result of this model, the energy loss formula reveals the specific way that energies of jets are modified by QGP. The completion of this model has filled the blank of theoretical prediction of jet energy distributions in this specific energy level of collision. Following ATLAS's way of illustrating the data, the distribution of jet energy is plotted against a variable defined as In the Glauber Monte Carlo model section of this article, frequently used terms include impact parameter, b defined as the distance between the center of two colliding heavy-ions and number of nucleons that are hit by at least one nucleon in heavy-ion collisions, part N . These nucleons are also called participants.

Probability Density as a Function of Impact Parameter
In the real circumstance of a collider, two beams of heavy nuclei move in opposite direction. They randomly distribute in a cross-section area (of mm in dimension) far greater than the cross-section area of one nucleon (of fm in dimension).
To determine the probability density of a collision taking part at impact parameter range , first switch to a reference frame where one nucleon in one of the two beams is located at the origin. The positive direction of z axis is selected as the same as the velocity of this nucleon. In this reference frame, incident nuclei are evenly distributed in the whole xy plane.
The number of incident nuclei that fall on the thin ring corresponding to impact parameter range where σ is the area density of incident nucleon. Then the probability density function

Applying Monte Carlo Glauber Model
The specific procedure of generating a Monte Carlo Glauber Model is well explained in Ref.
[7]. Fermi distribution is used in our model, with the expression of nuclear charge density In the computer program, parameter r should be randomly generated according to the probability density function where the specific ration index does not make a difference. For Pb + Pb collision, the Fermi parameters are listed in Table 1.

Results of Glauber Monte Carlo Model
Implementing the Monte Carlo Glauber, we can plot a 2-D graph showing the distribution of all Natural Science nucleons of two Pb nuclei on the transverse plane, as shown in Figure 1.
After that, we can get the number of participants-impact parameter and number of collisionsnumber of participants relationship, as shown in Figure 2.

Relationship of Impact Parameter to Centrality Ranges
In the procedure of Monte Carlo Glauber model, impact parameter b, the distance between the two centers of Pb nucleon, is an independent variable. In the ATLAS measurement Ref. [1], however, the variable centrality is used to represent the degree of overlapping between two nucleons. To utilize the MC Glauber model for simulating the results of ATLAS's experiments, quantitative relationship of impact parameter to centrality needs to be determined.
The definition of centrality cannot be written explicitly in formula. It is based on the division of the histogram of part N as shown in Figure 3. The curve of part N histogram and x-axis form a certain shape which corresponds to a certain area. In this graph, we can find a series of critical values 0 1 2 9 10 , , , , , x x x x x such that is a discrete variable. The integrals above degenerate into summation. The specific values of these critical points are listed in Table 2.

Calculation of Path Length
To simulate the energy loss of jets produced in QGP, the path length of each branch of the jet first needs to be calculated. Note that the QGP forms as a result of a number of nucleons colliding together, the shape of QGP is likely to be irregular and its boundary does not have a clear definition. It is thus more practical to define the path length based on position of participants, the origin of the jet and the equation of the jet line.
However, there still exists various ways to define jet length. The definition we adopt is the number of nucleon that each branch of the jet passes through. Under such definition, the path length is a discrete function. Another possible definition of the path length is the geometric distance from the origin of the jet and the point where the jet leaves the QGP. This definition was finally abandoned because of the following reasons. For one thing, as mentioned above, the boundary of the QGP is not clear. More importantly, it Natural Science has already been proved that the density in each part of QGP is not homogeneous [8]. The density is higher where there are originally more participants and vice versa. So, the energy loss of the jet per unit length is not constant.
The specific criterion to judge whether a participant contribute to the jet and if it does, which branch of the jet it contribute to can be explained by Figure 4.
In this graph, all participants in the orange shadow area make one-unit contribution to upper branch of the jet and all participants in the blue shadow area make one-unit contribution to lower branch of the jet. The width of the shadow area is D.
In each centrality range, the histogram of path length of each branch of the jet, noted as 1 L and 2 L , is shown in Figure 5.

Possible Forms of Energy Loss Formula
As the final step of our computer model, the form of the energy loss formula needs to be determined. As a jet propagates through QGP, the amount of energy it loses may be correlated to the path length and the initial energy of the jet. It is worth mentioning that the jets studied in our investigation have momenta of 100 GeV order of magnitude and thus are all ultra-relativistic. Their energy is approximately proportional to their momenta: . E pc = Thus energy loss is directly proportional to momentum loss. Before the energy loss formula is finally determined, we can only come up with possible forms of it.
where C is a constant parameter to adjust the overall magnitude of energy loss, L is path length and R is a randomly generated parameter to simulate fluctuations. Each term of the formula may have the following forms:   where the sign p . The minima of this 3-dimensional parameter space is relatively easy to look for, and we search for it by traversing the parameter space, evenly selecting points first in larger intervals and then in smaller intervals.

Implementing the Energy Loss Formula to pp Data
As input data to our model, original data of jet events in pp collisions at 2.7 V 6 Te

Effects of Some Parameters on the Shape of the Curve
Even though the energy loss formula consists of several independent parameters, and they act on energy of unquenched jets as a whole. It is still discovered that different parameters have different modulation effect on the shape of simulated curve.
Firstly, it is found that the greater the variance of the random term R, the less obvious the peak at centrality range 0 -10% and T p range 100 -126 GeV is. Let R follow Gaussian distribution with different variances and 1 R = as average value, the result of the curve is shown in Figure 6.
Secondly, we found that ( )

Adjustment of Parameters
In the procedure of adjusting the energy loss formula, ( ) f L is usually decided at first. After that, two parameters C and 0 p are adjusted to match the data in Ref. [1]. To describe the degree of resemblance of our predicted Pb + Pb data to that in ATLAS's measurement, peaks in each plots of our predicted data and that in ATLAS's measurement were found and compared. Also, we define a new variable "deviation" as:   The simulated data is considered to better resemble ATLAS's experimental results when average deviation D is smaller.
In our investigation, there are basically two ways to select the parameters. One is to minimize average deviation of centrality ranges 0 -10% and 10% -20% at T p range 100 -126 GeV, recorded as 2 D . The conspicuous peak of ATLAS's measurement is located in this region.
The other is to minimize average deviation of all 9 different combinations of centrality ranges and T p ranges, recorded as 9 D . The reason why these two ways of parameter selection were developed is that when 9 D reaches its minumum, tendency of the simulated curve of centrality ranges 0 -10% and 10% -20% deviates significantly from the measured result, failing to match the most significant feature of ATLAS's measurement. As shown in the Result section, these two ways of parameter selection will lead to two different forms of energy loss formula with better fitness at different centrality and T p ranges.

RESULTS
By the end of our investigation, the two ways of parameter selection did not give similar forms energy loss formula, leading to two two different forms of energy loss formula with better fitness at different centrality and T p ranges instead. And the formula that has the best fitness under two criteria, minimum 2 D and 9 D are shown separately. Predicted data is shown in Figure 9 and Figure  10.

Optimization Result Considering All Centrality and pT Ranges
For the second way of parameter selection, the relationship of minimum 9 D to α and the dependence of 9 D on C and 0 p at 0.8 α = is shown in Figure 11. The energy loss formula with best fitness is: corresponding to deviation 9 1.84 D = and 2 3.56 D = . Predicted data is shown in Figure 12 and Figure  13.

Discussion of the Results
For the first way of optimization, the simulated curves deviate from experimental data at centrality ranges 20% -30% and 30% -40% most significantly and for the second way of optimization, simulated curves fail to resemble the conspicuous peak at centrality range 0 -10% and when 1 100 1 6 GeV 2 T p < < . It is not satisfying enough that an energy loss formula that can simulate all centrality and T p ranges very well has not been found. So possibility remains that 1) the procedure of jet quenching in QGP cannot be described by a simple form of energy loss formula in Equation (4), and 2) the unfolding procedure done by the ATLAS collaboration cannot fully represent physical reality as the peculiar peak at low centrality ranges was not reported in measurements done by other experimental groups [11].

CONCLUSION
In conclusion, our investigation has simulated the procedure of jet quenching in Quark-Gluon Plasma using Glauber Monte Carlo Model, using energy data of jets produced in pp collision obtained from the CMS collaboration. The research has characterised the energy loss of jets in QGP with path length of jets, L, and their initial transverse momentum, T p . Simulated data were compared with ATLAS's measurement in Ref. [1] to determine best combinations of parameters in our model. Two different forms of energy loss formula were found to adapt for certain centrality ranges and T p ranges well. Our result