Using the Optical Glauber Monte Carlo Model to Calculate R AA Related Results Produced by the ATLAS Collaboration

The energy loss during jet quenching due to the existence of Quark Gluon Plasma (QGP) is calculated by Optical Glaube Monte Carlo model with data collected by ATLAS Collaboration using the LHC detector. An energy loss formula for this situation was modeled and took the form . The nuclear modification factor, AA R , for jets in a 208 Pb + 208 Pb nucleus collision with rapidity interval of y = 2.8 and the initial transverse momentum of 50 GeV ≤ p T ≤ 1000 GeV , are compared with various data plots produced by ATLAS Collaboration. AA R results are plotted in different centrality bins, which are de-fined by the distribution of number of participating nucleons N part . The AA R value was found to slowly increase at lower transverse momenta and flatten out at higher transverse momenta. The model’s


INTRODUCTION
High energy ultra relativistic heavy-ion collisions between two atomic nuclei in particle accelerators can result in a condensed substance that exhibits fluid-like properties known as Quark Gluon Plasma (QGP) [1], which is mainly constituted of colour charge carriers (e.g. quarks and gluons). Quantum Chromodynamics (QCD) describes the strong interactions between elementary particles such as quark and gluons, which carry colour charges. Unlike in Quantum Electrodynamics (QED), where field lines originating from charged particles spread infinitely throughout space, field lines between interacting quarks only exist in a thin tube between them; this has been determined experimentally and through simulations. Therefore, independent free Open Access Natural Science quarks can only exist under special circumstances, as the energy required to separate confined quarks found in hadrons would scale indefinitely with the distance of separation, and the energy inputted would make it favourable for the creation of a quark anti-quark pair due to vacuum fluctuations.
To better understand the properties and the compositions of QGP, particle colliders are used both to generate the plasma and to read their signatures through processes similar to scattering experiments, which can be done by observing and measuring properties of hard and soft probes exiting out of the QGP [2]. Jets, like vector-bosons and quarkonia, are a form of hard probes, for their signatures going through QGP are different compared to their signatures when moving through vacuum. Jets produced at nucleon-nucleon collision sites are subjected to energy modifications due to the effects imposed by the QGP. It has been observed that the energies of jets produced in Pb + Pb collisions are lower than those of jets produced in pp collisions; this phenomenon has been termed as jet quenching [3,4], as the jets' passage through QGP lowers their energies. The energy by which a jet loses is dependent on its initial transverse momentum T p and the distance through the QGP that the jet travels, which is modeled by the distance through the participating nucleons that the jet penetrates (L). This energy suppression could be observed by comparing the jet T p distribution in Pb + Pb collisions to the jet yield of pp collisions. The characterization of the distribution of the number of participating nucleons, part N , achieved by grouping these events into different centrality bins that correspond to different ranges of impact parameters (for example, smaller impact parameters can mean higher numbers of collisions and larger impact parameters can mean smaller numbers of collisions), can categorize data for further calculation regarding part N as shown in Table 1. Figure 1 displays the distribution of part N over 100,000 samples, separated into different centrality bins.
The measurement of the jet yield and nuclear modification factor R AA represents the relationship between the jet energy loss and the initial transverse momenta with respect to different centralities in a collision range is 90% to 100%; the second range is 80% to 90%, etc. The lowest range is 0% -10%.) N part has a mode at 2, and the frequency at which each N part occurs gradually decreases as N part increases, until falling to 0 after N part reaches 416, the maximum number of participants for a 208 Pb + 208 Pb collision.
Note the use of a logarithmic scale for the y-axis. event. The definition of R AA presented in the paper by the ATLAS Collaboration [5] is given in Equation In order to incorporate the Optical Glauber Monte Carlo Model, the following interpretation of AA R , given in Equation (2), is used to simplify and optimize the calculations.
where ( ) T p ψ is the frequency distribution of transverse momentum, T p . Using data provided by the ATLAS Collaboration from the LHC detector, this paper aims to analyze the relationship between T p and AA R under different centrality ranges, attempt to find a formula for the jet energy loss, and reproduce some of the plots presented in the paper by the ATLAS Collaboration [5] through the use of computer simulations built under the Optical Monte Carlo Model.

METHOD
The program setup is based on the Optical Glauber Monte Carlo approach [6][7][8], which is a simplistic, geometry based method of approximating part N (the number of participating nucleons) and coll N (the number of binary nucleon-nucleon collisions), given an impact parameter b and a certain nucleus radius r of a nucleus-nucleus collision. Within this model, it is assumed that the velocities of the two colliding lead nuclei are exactly opposite of each other. For the purposes of this project, it was assumed that the protons and neutrons are geometrically identical and that the charge that the protons carry has no effects on the collisions. The initial locations of nucleons are determined by assuming that the density distribution of the nucleons with respect to the distance to the center of the nucleus is rectangular, but disregarding the effects of nuclear charge density and minimum inter-nucleon separation distance. As the experiment using the ATLAS detector at the LHC indicated, the 208 Pb nucleus, containing 208 nucleons each with a radius of 0.546 fm, has a radius of 6.62 fm. Natural Science The event of a Pb-Pb collision is simulated using the computer program ROOT, which utilizes random number generators to generate the positions of the nucleons in two spheres to simulate two Lead nuclei. The position of each nucleon is determined randomly by selecting arbitrary points within a sphere of radius 6.62 fm; this step is repeated for each of the 208 nucleons in each nucleus. Coordinates of these three-dimensional spherical nuclei, separated by the impact parameter (which is determined randomly given the distribution d d 2 , are then projected onto a two-dimensional plane. If the distance between the centers of two arbitrary nucleons is less than the diameter of a nucleon (which is 1.092 fm), they are registered as collided and are participants. If a given nucleon does not collide with any other nucleon in the other nucleus, the nucleon would be regarded to as a spectator. Figure 2 shows an example of a simulated collision with this model. This model treats the nuclei collision as numbers of independent binary nucleon-nucleon collisions. The model also assumes that all the nucleons carry ultra-relativistic momentum and that they travel in a straight trajectory according to their initial positions and can collide with multiple nucleons without any path change or deflection. These properties indicate that the inelastic nucleon-nucleon cross section is not disturbed by other collisions made by the same nucleon. Due to the fact that the elastic cross sections do not contribute to a significant amount of energy loss, they are not considered here in this model.
The Optical Glauber Monte Carlo model allows for the derivation of centrality ranges from the N part distribution, which characterize whether the collision is central or peripheral. Dijets, two jets that start at the same position but facing the opposite directions, are produced, initiating at the centers of the nucleon-nucleon collisions going at a random angle θ. It is important to note that, in this model, these jets are infinitely thin and will not be deflected or curved during their contact within the Quark Gluon Plasma. The amount of QGP that each jet goes through is calculated by taking the sum of the total distance of participating nucleons the jet penetrates through (L 1 and L 2 for the distances that the dijet goes through on each side). Figure 3 shows the distribution of L 1 and L 2 over 4 different centrality bins: 0% -10%, 20% -  30%, 40% -50%, and 60% -70%. However, the model has certain downsides regarding the physics behind it. It assumes that all the nuclei and nucleons are completely spherical, thus disregarding most of the ultra-relativistic and quantum affects. According to experimental data, it is theorized that, during process of jets traveling through the QGP, different kinds of interactions will contribute to an energy loss. At present this process is not completely understood, but many researchers have done numerous works on this energy loss formula [9,10]. We propose an energy loss formula written in the form as given by Equation (3) in order to compensate for all the above causes. The use of the equation will be demonstrated later in the paper.

RESULTS AND DISCUSSIONS
Equation (3) presents the formula for the energy lost for a jet with transverse momentum T p traveling through QGP a distance of L. Using this energy loss formula, we produced different graphs of the amount of energy lost based on the initial transverse momentum at different centrality ranges. Figure 4 presents several graphs of energy loss at centrality ranges of 0% -10%, 20% -30%, 40% -50%, and 60% -Natural Science 70%. Similar to the results obtained by the ATLAS Collaboration, Figure 5 shows the nuclear modification factor R AA as a function of Jet T p for eight centrality ranges. The R AA value is calculated by jets with 2.8 y = and T p between 50 GeV and 1000 GeV. Following the ATLAS results, we have modified our bins using variable binning to compress higher T p jet intervals before evaluating R AA . R AA values increase slowly as T p values increase for all events at different centrality ranges until flattening out at high T p values. R AA predictions yield by the Optical Glauber Monte Carlo model is fairly close to the actual experimental data collected by the LHC and presented by ATLAS. Most of the data points presented by the model can potentially lie directly on top of the ATLAS data if luminosity uncertainty is considered. This is a proof that the energy loss formula is sufficient to compensate most the ultra-relativistic quantum effects that are not considered in the Optical Glauber Monte Carlo model.  The results seem like that both the initial T p and the amount of QGP will influence the energy loss during the collision. When the initial T p become higher, the influences of QGP on jets will tend to be constant. Therefore, we predict that the QGP will exert a stronger effect on jets with lower transverse momentum, which means the jets with lower T p will lose relative more energy than the jets with higher T p due to the QGP. It is also explicit that the jet energy loss is dependent on the amount of QGP (and subsequenty, the number of participants and the impact parameter), and the amount of QGP will become more crucial factor as T p increases. Natural Science

CONCLUSION
Nuclear Modification Factor (R AA ) predictions, calculated using the Optical Glauber Monte Carlo model and random Dijet generation functions with transverse momentum ( T p ) data for Pb-Pb nucleon-nucleon collisions at center-of-mass energy of 5.02 TeV provided by the ATLAS detector at the LHC, are consistent with the data yield from ATLAS results. From the plots, one can see that the effects of Quark Gluon Plasma (QGP) on jets vary with T p and the different centrality ranges, which correspond to the part N . These results represent the success of the Optical Glauber Monte Carlo Model at replicating the actual experimental data. So the Optical Glauber Monte Carlo Model turns out to be a plausible theoretical model for the evaluation of jet quenching.