Transverse momentum spectra of mesons produced in p+p collisions are studied in the framework of a thermalized cylinder model. In the region of high transverse momentum, the considered distributions have a tail part at the maximum energy of RHIC. A two-component distribution based upon the improved cylinder model is used to fit the experimental data of the PHENIX Collaboration. It is found that the improved approach can describe the meson production in the wider range of transverse momenta.
1. Introduction
Relativistic heavy ion collisions are performed to study nuclear matter under extreme temperature and density conditions [1]. Proton-proton collisions are conventionally used as a reference to compare with nuclear collisions and to understand the observed collective effects. And the Large Hadron Collider (LHC) [2] was originally designed to rise center-of-mass (cm) energy to 14 TeV, which is almost 8 times the record of 1.8 TeV achieved by the Fermilab in the United States. After the LHC data obtained in p+p interactions at sNN=900 GeV and 7 TeV [3, 4], the new interest in general features of p+p collisions at ultrarelativistic energies appeared. In p+p collisions, the meson spectra provide insight into the particle production mechanism and interaction in the hadronic and quark gluon plasma (QGP) phases. Furthermore, the detailed study of meson spectra is important because it acts as an ingredient for estimating the hadronic decay backgrounds in the photon, single lepton, and dilepton spectra, which are the penetrating probes of QGP. To isolate phenomena related to the dense and hot medium created in such collisions, it is also important to measure particle production in smaller collision systems like p-p and d-Au. Measurements of transverse momentum spectra for particles produced in p+p collisions are used as a baseline, to which similar measurements for heavy ion collisions are compared. In addition, the nuclear modification factor RAA of several identified hadrons with high transfer momentum is used to probe jet quenching.
The information about the production process is retained by the final-state particle distributions in the collisions [5–9]. Single-hadron production at large transverse momenta in high-energy hadronic and nuclear collisions results from the fragmentation of quarks and gluons issuing from parton-parton scatterings with large momentum transfer. To explain the abundant experimental data, different phenomenological mechanisms of initial coherent multiple interactions and particle transports were proposed and extended in recent years [10–15]. The comparison to model calculations can provide valuable information of the collision evolution and help better understand properties of the QGP. In particular, several theory models of high energy collisions were reported in a workshop held at the CERN Theory Institute [16]. Recently, systematic studies on the production of final-state particles performed by the NA49 collaboration were discussed [17]. Hadronic transport models fail to describe the production of final-state particles, while the results of statistical models are generally in good agreement with the measured particle yields at available energies.
Based on the one-dimensional string model [19] and the fireball model [20], we have developed a thermalized cylinder model, which successfully describes the particle production in heavy ion collisions over an energy range from the Alternating Gradient Synchrotron (AGS) to the relativistic heavy ion collider (RHIC) [21–25]. In our previous work [25], the transverse momentum distributions of strange hadrons produced in Cu+Cu and Au+Au collisions at RHIC energies were explained by a single component distribution. The excitation degree of the emission source is allowed to be determined by studying the transverse momentum spectra, anisotropic flow effects, and their correlations. We found that the single component distribution only describes a narrow transverse momentum rang. Recently, the invariant differential cross-section for production of neutral mesons in p+p collisions at sNN=200 GeV were published by the PHENIX Collaboration [18]. The high pT reach of the transverse momentum is helpful to characterize the mechanisms of truly perturbative parton-parton scatterings and parton fragmentation in different QCD environments. Particle yields observed in these experiments inspired our work. It is interesting for us to analyze the results of pp collisions at the RHIC energies. In order to verify the thermalized cylinder model and describe the broad distribution range, in this paper, we will develop the single component spectra into a two-component distribution.
2. The Formula
According to the cylinder model [25], we assume the projectile and target cylinders to be formed in A+A, d+A and p+p collisions at high energies. The cylinders are wounded sources of particles. The idea comes from the observation that the process of particle production is not instantaneous, which was noted in [26, 27]. In the reference frame where the longitudinal momentum of a produced particle vanishes, the minimal time necessary for its creation is t0⩾1/MT, where MT=pT2+m2. In the laboratory frame, the particle in question acquires some longitudinal momentum and we have
(1)t⩾γt0=Emt2=coshymt,
where γ and E are the Lorentz factor and the energy of the particle, respectively. So, the resolving power of the longitudinal distance (the uncertainty of the distance from the collision point to those at which the particle is created) is
(2)l=vt=sinhymt.
When the rapidity of the produced particle is large enough so that l>Z(b), where Z(b) is the size of the nucleus at a given impact parameter, the particle cannot resolve individual collisions. It is natural to suggest that its creation may be insensitive to the number of collisions in the source. This is the origin of the idea of wounded sources. The concept of the wounded source may practically be applied in the whole rapidity region.
In AA collisions the cylinder is thick, and in p+p collisions the cylinder is thin. Final-state particles are randomly produced from emission sources in the cylinder(s). The excitation degree of the side-surface region is naturally lower than that of the central axis region in the cylinders. From central axis region to side-surface region of the concerned emission source, the excitation degree is assumed to decrease linearly along the transverse axis direction [25]. When the excitation degree increases, the value of the distribution width σ becomes bigger and bigger. So, the excitation degree can be characterized by the momentum distribution width σ. The emission points with the same excitation degree form an emission circle in the transverse momentum space. Therefore, a given σ corresponds to wounded sources which emits a fixed density of particles. Let σs and σc denote the widths of the transverse momentum distribution of particles produced in the side-surface and central axis regions, respectively. The distribution of σ can be given by
(3)f(σ)=C(σc-σ),
with the normalization constant
(4)C=2(σc-σs)2.
According to the single component distribution in the Monte Carlo calculation, the transverse momentum distribution of final-state particles can be given by
(5)pT=(σc-(σc-σs)2-2R1C)-2lnR2,
where R1 and R2 are the random variable distributed in [0,1].
In our previous work [25], transverse momentum spectra of strange particles produced in high-energy collisions were investigated by the above thermalized cylinder model. Considering the simply cylinder shape, we obtained the emission source location dependence of the exciting degree specifically. It is shown that the single component distribution, that is, (5), is successful in describing the experimental results measured by the STAR and PHOBOS Collaborations. Moreover, we found that the single component distribution can describe a narrow transverse momentum range. To explain the wider transverse momentum spectra of identified particles produced in p+p collisions, we need to consider the relative importance of hard versus soft processes in the particle production mechanisms at different energies. Hard parton-parton scatterings with large momentum transfer occur on a short time scale and are governed by perturbative QCD. The bulk of particle production occurs via soft processes (with low momentum transfer and consequently longer time scales) which are described with phenomenological models.
In order to describe the invariant cross-section of hadrons as a function of pT over a wide range, Hagedorn proposed an empirical formula [28], which is described as “inspired by QCD” and is given by
(6)Ed3Ndp3=B(1+pT/λ)n≃B{(-npTλ)forpT→0,(λpT)nforpT→∞,
where B, λ, and n are fit parameters. The limiting cases show that the distribution is an exponential form at low transverse momenta, but it is a power law at large transverse momenta, which came from “QCD inspired” quark interchange model [29],
(7)Ed3Ndp3~(mT2)-4~1(pT)8.
The dominant contribution to the inclusive cross section is low pT particles. However, with center of mass energy increasing, the hardening of the pT distribution implies an increase of the transverse momentum. At both low pT and high pT, UA1 collaboration [30] gives a hybrid form,
(8)Ed3Ndp3=Aexp(-bmT),forpT<pTc,=B(1+pT/p0)n,forpT>pTc,
where pTc is a free parameter. In order to give a good description of pion spectra in wide pT range, the PHENIX Collaboration [18, 31, 32] found a single form referred as the modified Hagedorn formula. The modification is to better describe the π0 spectrum for wider pT range, in particular at high pT where the spectrum behaves close to a simple power law function. The single formula has been used to successfully describe the hadron spectra measured in p-p collisions at different energies [18, 31, 32]. The formula has been extensively used in terms of mT,
(9)Ed3Ndp3=B[exp(-amT-bmT2+mT/λ)]n,
where a and b are the fit parameters. The distribution function is close to an exponential form at low pT and a pure power law form at high pT. Similarly, our model should be improved to a two-component distribution,
(10)PT=k(σ1c-(σ1c-σ1s)2-2R1/C1)-2lnR2+(1-k)(σ2c-(σ2c-σ2s)2-2R3/C2)-2lnR4,
where R3 and R4 are the random variable distributed in [0,1], k(1-k) and C1(C2) are the contribution coefficient and normalization constant of the first (second) component respectively, and σ1s(σ2s) and σ1c(σ2c) are the momentum distribution widths of particles produced in side-surface region and central axis region of the first (second) component, respectively. Our calculation shows that C1≠C2. The first and second items in (10) correspond to the contributions of soft production and hard emission, respectively.
The basic model assumes the cylinders with the emitted momentum distribution width changing from the side-surface to the central axis. In fact, except for the most central collisions, the active (overlap) region is certainly not cylindrical. We could understand the “cylinder” as a nonperfect cylinder along the beam direction in the momentum space but not a perfect cylinder in the coordinate space. The emission sources with the same excitation degree stay at the same (sub)surface. This situation is similar to the equipotential surface in electromagnetism.
3. Comparison with PHENIX Results
The invariant differential cross-sections of neutral mesons produced in p-p collisions at nucleon-nucleon center-of-mass energy sNN=200 GeV in various decay modes are presented in Figure 1. The symbols represent the experimental data of the PHENIX Collaboration [18] and the curves are our calculated results by using (10). The χ2 testing provides statistical indication of the most probable value of corresponding parameters. The values of k, σ1c-σ1s, and σ2c-σ2s obtained by fitting the experimental data are given in Table 1 with χ2 per degree of freedom (dof). The maximum value of the observed pT reaches about 13.0 GeV. We see that the calculated results approximately agree with the experimental data for neutral mesons in the region. The values of k, σ1c, and σ2c are taken to be 99.98±0.02, 1.70±0.12 and 4.64±0.35, respectively. Both σ1s, and σ2s values do not change obviously.
Values of the parameters k, σ1s~σ1c, and σ2s~σ2c obtained by fitting the experimental pT spectra in Figures 1, 2, and 3. The units of σ1s, σ1c, σ2s, and σ2c are GeV/c.
Particle
k (%)
σ1s~σ1c (×0.1)
σ2s~σ2c (×0.1)
χ2/dof
ω
99.97±0.01
0.09±0.02~1.66±0.10
1.28±0.10~4.55±0.24
1.000
K
99.97±0.01
0.06±0.01~1.60±0.12
1.22±0.10~4.50±0.21
1.310
η′
99.98±0.01
0.07±0.01~1.68±0.10
1.28±0.10~4.65±0.27
0.472
ϕ
99.99±0.00
0.05±0.01~1.72±0.10
1.30±0.10~4.75±0.24
0.407
π
99.96±0.01
0.08±0.01~1.65±0.10
1.26±0.10~4.52±0.22
0.210
η
99.97±0.01
0.06±0.01~1.61±0.11
1.24±0.11~4.48±0.26
0.240
p+p-
100.00±0.01
0.10±0.02~1.78±0.10
—
0.278
J/ψ
100.00±0.01
0.16±0.03~2.25±0.15
—
0.447
ψ′
100.00±0.01
0.18±0.04~2.10±0.12
—
0.470
Invariant differential cross-section of neutral mesons measured in p+p collisions at sNN = 200 GeV in various decay modes. Experimental data measured by the PHENIX Collaboration [18] are shown by the scattered symbols. Our results calculated from the cylinder model are shown by the curves.
To compare these results to other particles, we show the invariant differential cross sections of different particles produced in p+p collisions at sNN = 200 GeV in various decay modes in Figure 2. The different symbols are the experimental data of the PHENIX Collaboration [18]. The solid curves are our results calculated by the model. By fitting the experimental data, the parameter values are given in Table 1 with the values of χ2/dof. Similar to Figure 1, the values of σ1s(σ2s) do not change significantly. We see again that the model describes approximately the invariant differential cross-sections of final-state mesons produced in p+p collisions at the highest RHIC energy. Compared with the case of Figure 1, the values of parameters do not change significantly.
Invariant differential cross-sections of different particles measured in p+p collisions at sNN = 200 GeV in various decay modes. Experimental data measured by the PHENIX Collaboration [18] are shown by the scattered symbols. Our results calculated from the cylinder model are shown by the curves.
In Figure 3, we present the invariant cross section at midrapidity for η in p+p collisions at sNN = 200 GeV. The symbols are the experimental data of the PHENIX Collaboration [33, 34]. The solid curves are our results calculated by the model. The dotted and dashed curves are our results corresponding to the contributions of the first and second items in (10). The parameter values and corresponding χ2/dof are given in Table 1. One can see the particle production of soft production and hard emission intuitively in the figure.
Invariant η cross section at midrapidity as a function of transverse momentum in p+p at sNN = 200 GeV. The error bars are the statistical uncertainties. Experimental data taken from the PHENIX Collaboration [33, 34] are shown by the scattered symbols. Our results calculated from the cylinder model are shown by the curves.
In order to testify the validity of the model, Figure 4 shows the invariant π0 cross section, Ed3σ/d3p, in p+p collisions at different RHIC energies. The symbols are the experimental data of the PHENIX Collaboration [35]. The solid curves are our results calculated by the model. The dotted and dashed curves are our results corresponding to the contributions of the first and second items in (10). The values of k, σ1s~σ1c, and σ2s~σ2c are given in Table 2 with the corresponding χ2/dof. The calculated results are in agreement with the experimental data in p-p collisions at the RHIC energies. The k values decrease with the increase of the collision energy, and σs has no obvious change, equals only about 0.60 GeV/c. The values of σ1c (0.99–1.81 GeV/c) and σ2c (3.94–4.83 GeV/c) increase with the collision energy increasing. As discussed in detail in [25], the excitation degree of emission source on the cylinder central axis becomes higher due to the higher collision energy.
Values of the parameters k, σ1s~σ1c, and σ2s~σ2c obtained by fitting the experimental pT spectra Figure 4. The units of σ1s, σ1c, σ2s, and σ2c are GeV/c.
Particle
Energy (GeV)
k (%)
σ1s~σ1c (×0.1)
σ2s~σ2c (×0.1)
χ2/dof
π0
22.4
100.00±0.01
0.05±0.01~1.05±0.06
—
0.600
62.4
99.95±0.01
0.05±0.01~1.26±0.09
1.10±0.07~4.14±0.20
1.104
200
99.69±0.02
0.07±0.01~1.69±0.12
1.29±0.10~4.58±0.25
0.422
Invariant π0 cross-section as a function of transverse momentum in p+p collisions at different RHIC energies. The error bars are the quadratic sum of statistical and systematic uncertainties. Experimental data taken from the PHENIX Collaboration [35] are shown by the scattered symbols. Our results calculated from the cylinder model are shown by the curves.
We would like to point out that the calculated results in Figures 1–4 are obtained by (10), which is a general representation in the Monte Carlo method. Equations (6)–(9) are only used to explain why the two components need to be considered in (10). The curves in the figures are numerical results from (10). Some of them are not smooth due to low statistics.
4. Discussion
The above comparison shows that the improved cylinder model can be used in the description of the meson production in the wider range of transverse momenta. The values of χ2/dof for all fits are shown in Tables 1 and 2. The maximum value is 1.420, and the minimum value is 0.246. We demonstrate that the calculated results are in good agreement with the available experimental data. In our calculation, k is the most important parameter which determines the distributive slope and range. In most cases, the second item cannot be neglected due to the contribution of the hard process. A large k (≃1) gives a precipitous and narrow distribution, whereas a smaller k (<1) gives a subdued and wide distribution. The second item in (10) is very small (<0.01) and does not contribute to the spectra at the lowest RHIC energy. In p+p collisions for different energies, we see that 1-k increases slightly with the increasing of collision energies. This means that the contribution of hard emission increases slightly. For the asymmetric system of d+Au collisions, it is predicted that the values of k do not change significantly due to the spatial asymmetry of the collision nucleons. Scattering processes at high energy hadron colliders can be classified as either hard or soft. Quantum Chromodynamics (QCD) is the underlying theory for all such processes, but, the approach and level of understanding arevery different for the two cases. For hard processes (e.g., Higgs boson or high pT jet production), the rates and event properties can be predicted by using perturbation theory. In addition, for soft processes (e.g., the total cross section, the underlying event, etc.), the rates and properties are dominated by nonperturbative QCD effects, which are less well understood. For many hard processes, soft interactions are occurring along with the hard interactions and their effects must be understood for comparisons to be made to perturbative predictions. An understanding of the rates and characteristics of predictions for hard processes, both signals and backgrounds, using perturbative QCD (pQCD) is crucial for both the Tevatron and LHC. We are in a position to evaluate the soft and hard contributions to the observed spectra by using k in the statistical model.
The particles produced in high-energy nuclear collisions have attracted much attention since people are trying to understand the properties of strongly interacting quark-gluon plasma by studying the possible production mechanisms [36, 37]. The final-state particles are emitted isotropically in the rest frame of emission sources with the different excitation degree in collisions. If we consider that the local emission source has a motion in the transverse direction, there are interactions among the emission sources, then the transverse flow (directed flow and elliptic flow) can be explained by the model. Thermal-statistical models have been successful in describing particle yields in various systems at different energies. The cylinder model is developed from the fireball model, which is suggested in heavy-ion collisions [20]. The excitation degree varies with location in the cylinder. In the previous work [25], we obtained the emission source location dependence of the exciting degree specifically. From central axis to side-surface of the cylinder, the excitation degree of the emission source decreases linearly with the direction of radius. The excitation degree can easily be characterized by the corresponding distribution width, which is written as (4) by using the Monte Carlo calculation. In the model, because of the influence of the emission source temperature, the values of σs are much smaller than those of the σc. It is consist with the conclusion of [15], where the temperature of the fireball decrease linearly with time.
In the present work, we improve the method by considering the difference between the soft and hard emission. The parameters σ1c(σ1s) and σ2c(σ2s) are used to describe the excitation degrees of emission sources closing to the side-surface and central axis of the cylinders, respectively. They are in fact reflections of the excitation degree of the soft emission process and the hard emission process, respectively. Our results agree well with the considered distributions of mesons with high transverse momenta produced in nuclear collisions at RHIC energies. At finite temperature, the stronger the collision strength, the larger the excitation degree. Therefore, on the central axis of the cylinder, the interaction between particles is strongest, and the excitation degree is highest. For different collision energies, our results in the improved cylinder model show that the parameters σ1c and σ2c slowly increase with the increasing of energy, and no obvious change can be observed in the values of σs1(σs2). The dependence of parameters on energies renders that intranuclear cascade collisions in the central axis regions play a more important role at higher energy. In the side-surface regions, a high multiplicity at high energy does not contribute to a high excitation degree due to weak intranuclear cascade collisions. In our opinion, high pT hadrons come primarily from the central region due to their large emission angles. And high momentum hadrons come primarily from the side-surface region due to their weak cascade collisions.
5. Conclusions
In summary, the transverse momentum distributions of mesons produced in p+p collisions at the RHIC energies have been studied in the framework of the improved cylinder model. The model is successful in description of the high transverse momentum meson production. Based on our phenomenological approach, we evaluated the soft and hard contributions to the observed spectra. At the same time, it can offer some information about soft and hard interactions in the collisions. In our previous work, a rudimentary investigation of the azimuthal anisotropy has been carried basing upon the cylinder model. Combined with the present work, the model can be used to describe uniformly the momentum distributions of final-state particles produced in high energy collisions.
Acknowledgments
This work is supported by the National Natural Science Foundation of China under Grant no. 11247250, no. 11005071, and no. 10975095, the National Fundamental Fund of Personnel Training under Grant no. J1103210, the Shanxi Provincial Natural Science Foundation under Grant no. 2013021006 and no. 2011011001, the Open Research Subject of the Chinese Academy of Sciences Large-Scale Scientific Facility under Grant no. 2060205, and the Shanxi Scholarship Council of China.
CainesH.Recent results from RHIC: first steps towards quantitative QGP measurements201214225AaijR.Abellan BetetaC.AdametzA.AdevaB.AdinolfiM.LHCb CollaborationMeasurement of prompt hadron production ratios in pp collisions at √s=0.9 and 7 TeV201272216810.1140/epjc/s10052-012-2168-xAamodtK.AbelN.AbeysekaraU.Abrahantes QuintanaA.AbramyanA.AdamovD.ALICE CollaborationTransverse momentum spectra of charged particles in proton-proton collisions at s=900 GeV with ALICE at the LHC20106932536810.1016/j.physletb.2010.08.026KhachatryanV.SirunyanA. M.TumasyanA.AdamW.BergauerT.DragicevicM.EroJ.CMS CollaborationTransverse-momentum and pseudorapidity distributions of charged hadrons in pp collisions at √s=7 TeV20101051402200210.1103/PhysRevLett.105.022002BackB. B.BakerM. D.BallintijnM.BartonD. S.BettsR. R.BickleyA. A.PHOBOS CollaborationCentrality dependence of charged hadron transverse momentum spectra in Au+Au collisions from √sNN=62.4 to 200 GeV2005944082304BackB. B.BakerM. D.BartonD. S.BettsR. R.BallintijnM.BickleyA. A.PHOBOS CollaborationCharged hadron transverse momentum distributions in Au+Au collisions at sNN=200 GeV200457829730310.1016/j.physletb.2003.10.101AdamsJ.AggarwalM. M.AhammedZ.AmonettJ.AndersonB. D.AndersonM.ArkhipkinD.STAR CollaborationScaling properties of hyperon production in Au+Au collisions at sNN=200 GeV200798606230110.1103/PhysRevLett.98.062301VasconcelosG. M. S.STAR CollaborationMulti-strange particle production in relativistic heavy ion collisions at sNN=62.4 GeV201037094034PreghenellaR.ALICE CollaborationMeasurement of the production cross section for Z/gamma* in association with jets in pp collisions at s=7 TeV with the ATLAS detector20128503200910.1103/PhysRevD.85.032009MillerM. L.ReygersK.SandersS. J.SteinbergP.Glauber modeling in high-energy nuclear collisions20075720524310.1146/annurev.nucl.57.090506.123020De VriesH.De JagerC. W.De VriesC.Nuclear charge-density-distribution parameters from elastic electron scattering19873634955362-s2.0-12444310071Braun-MunzingerP.MagestroD.RedlichK.StachelJ.Hadron production in Au-Au collisions at RHIC20015181–24146RafelskiJ.LetessierJ.Testing limits of statistical hadronization200371598c107c10.1016/S0375-9474(02)01418-5RoyV.ChaudhuriA. K.Transverse momentum spectra and elliptic flow in ideal hydrodynamics and geometric scaling2010814067901SongT.KoC. M.LeeS. H.XuJ.J/Ψ production and elliptic flow in relativistic heavy-ion collisions2011831301491410.1103/PhysRevC.83.014914ArmestoN.BorghiniN.JeonS.WiedemannU. A.AbreuS.AkkelinS. V.Heavy-ion collisions at the LHC—Last call for predictions200835054001BlumeC.NA49 CollaborationStrangeness production at the SPS200912385AdareA.AfanasievS.AidalaC.AjitanandN. N.AkibaY.Al-BatainehH.PHENIX CollaborationMeasurement of neutral mesons in p+p collisions at s=200 GeV and scaling properties of hadron production20118305200410.1103/PhysRevD.83.052004WernerK.Strings, pomerons and the VENUS model of hadronic interactions at ultrarelativistic energies19932322–58729910.1016/0370-1573(93)90078-RWestfallG. D.GossetJ.JohansenP. J.PoskanzerA. M.MeyerW. G.Nuclear fireball model for proton inclusive spectra from relativistic heavy-ion collisions1976371202120510.1103/PhysRevLett.37.1202LiuF. H.PanebratsevY. A.Baryon rapidity distribution in nucleus-nucleus collisions at ultrarelativistic energies1999591193119510.1103/PhysRevC.59.1193LiuF. H.PanebratsevY. A.Pseudorapidity distribution of shower particles in heavy ion induced interactions in nuclear emulsion at high energy1999593179818012-s2.0-0033244069LiuF. H.Particle production in Au—Au collisions at RHIC energies20045831–2687210.1016/j.physletb.2003.12.059LiuF. H.YuanY.DuanM. Y.Widths of charged particle pseudorapidity density distributions in high energy pp- and AA collisions20088013–415416810.1016/j.nuclphysa.2007.12.006LiB. C.FuY. Y.WangL. L.WangE. Q.LiuF. H.Transverse momentum distributions of strange hadrons produced in nucleus-nucleus collisions at sNN=62.4 and 200 GeV20123902500910.1088/0954-3899/39/2/025009LandauL.PomeranchukI.Limits of applicability of the theory of bremsstrahlung electrons and pair production at high-energies195392535536LandauL.PomeranchukI.Electron cascade process at very high-energies195392735738HagedornR.Multiplicities, P(t) distributions and the expected hadron19846N10150BalnkenbeclerR.BrodskyS. J.Unified description of inclusive and exclusive reactions at all momentum transfers19741092973299210.1103/PhysRevD.10.2973AlbajarC.AlbrowM. G.AllkoferO. C.AndrieuB.AnkoviakK.ApsimonR.UA1 CollaborationA study of the general characteristics of proton-antiproton collisions at s=0.2 to 0.9 TeV1990335261AfanasievA.AdareS.AidalaC.AjitanandN. N.AkibaY.Al-BatainehH.AlexanderJ.PHENIX CollaborationDetailed measurement of the e+e- pair continuum in p+p and Au+Au collisions at s=200 GeV and implications for direct photon production201081034911AdareA.PHENIX CollaborationHeavy-quark production in p + p and energy loss and flow of heavy quarks in Au + Au collisions at s=200 GeV20118404490510.1103/PhysRevC.81.034911
AdlerS. S.AfanasievS.AidalaC.AjitanandN. N.AkibaY.Al-JamelA.AlexanderJ.PHENIX CollaborationCentrality dependence of π0 and η production at large transverse momentum in s=200 GeV d+Au collisions200798172302AdareA.AfanasievS.AidalaC.AjitanandN. N.AkibaY.Al-BatainehH.AlexanderJ.PHENIX CollaborationTransverse momentum dependence of η meson suppression in Au+Au collisions at sNN=200 GeV201082701190210.1103/PhysRevC.82.011902AdareA.AfanasievS.AidalaC.AjitanandN. N.AkibaY.Al-BatainehH.AlexanderJ.PHENIX CollaborationOnset of π0 Suppression Studied in Cu+Cu Collisions at sNN=22.4, 62.4, and 200 GeV2008101162301AndronicA.Braun-MunzingerP.StachelJ.Hadron production in central nucleus-nucleus collisions at chemical freeze-out20067723-41671992-s2.0-3374453980310.1016/j.nuclphysa.2006.03.012AndronicA.Braun-MunzingerP.StachelJ.Thermal hadron production in relativistic nuclear collisions2009404100510122-s2.0-65149098609