Initial- and final-state temperatures of emission source from differential cross-section in squared momentum transfer in high energy collisions

The differential cross-section in squared momentum transfer of $\rho$, $\rho^0$, $\omega$, $\phi$, $f_{0}(980)$, $f_{1}(1285)$, $f_{0}(1370)$, $f_{1}(1420)$, $f_{0}(1500)$, and $J/\psi$ produced in high energy virtual photon-proton ($\gamma$$^{*} p$), photon-proton ($\gamma p$), and proton-proton ($pp$) collisions measured by the H1, ZEUS, and WA102 Collaborations are analyzed by the Monte Carlo calculations. In the calculations, the Erlang distribution, Tsallis distribution, and Hagedorn function are separately used to describe the transverse momentum spectra of the emitted particles. Our results show that the initial- and final-state temperatures increase from lower squared photon virtuality to higher one, and decrease with increasing of center-of-mass energy.


I. INTRODUCTION
In high-energy collisions, it is interesting for us to describe the excitation and equilibrium degrees of an interacting system because of the two degrees related to the reaction mechanism and evolution process of the collision system [1][2][3][4][5][6][7][8][9][10]. In the progress of describing the excitation degree and structure character of the system, temperature is an important quantity in physics in view of intuitiveness and representation. In high-energy collisions, different types of temperatures are used [11][12][13][14][15][16][17][18], which usually refer to the initial-state temperature T i , quarkhadron transition temperature T tr , chemical freeze-out temperature T ch , kinetic freeze-out or final-state temperature ("confinement" temperature) T kin or T 0 , and effective temperature T ef f or T , etc. In this work, we emphatically discuss the initial-and final-state temperatures, though other types of temperatures are also important.
The initial temperature T i is the temperature of emis- * qiwang-sxu@qq.com; 18303476022@163.com † Correspondence: fuhuliu@163.com; fuhuliu@sxu.edu.cn ‡ Correspondence: khkolimov@gmail.com; kh.olimov@uzsci.net sion source or interacting system when a projectile particle or nucleus and a target particle or nucleus undergo initial stage of a collision. It represents the excitation degree of emission source or that of an interacting system in the initial state of collisions, and it is usually meant as describing the interacting system after thermalization. The initial temperature T i can be extracted by fitting the transverse momentum p T spectra of particles by using some distributions such as the Erlang distribution [19][20][21], Tsallis distribution [22,23], Hagedorn function [24], Lévy-Tsallis function [25]. Here, both the names of distribution and function are used according to the various accepted terminology in the literature, though they represent the similar probability density function in fitting the particle spectra. Meanwhile, the average transverse momentum p T can be obtained from the same function.
The final-state temperature T 0 is usually known as the kinetic freeze-out temperature, which refers to the temperature of emission source when the inelastic collisions ceased and there are only elastic collisions among particles. In the last stage of collisions, the momentum distribution of particles is fixed and the transverse momentum spectra can be measured in experiments. The excitation degree of the system in the last stage can be described by the final-state temperature T 0 in which the influence of flow effect is excluded. The temperature or related main parameters used in the Erlang distribution [19][20][21], Tsallis distribution [22,23], Hagedorn function [24], and Lévy-Tsallis function [25] are not T 0 , but the effective temperature T in which the influence of flow effect is not excluded.
The Mandelstam variables [26] consist of the fourmomentum of particles in two-body reaction. Both the squared momentum transfer and transverse momentum can represent the kinetic character of particles. Let us use the squared momentum transfer to replace transverse momentum in fitting the particle spectra. Then, we can fit the squared momentum transfer spectra with the related distributions to obtain the initial temperature T i , average transverse momentum p T , and other quantities. Of course, in fitting the squared momentum transfer spectra, the above mentioned distributions cannot be used directly. In fact, we have to use the Monte Carlo method to obtain the concrete value of a transverse momentum for a given particle from the mentioned distributions. Then, the concrete value of squared momentum transfer can be obtained from the definition.
Except the temperature parameter, other parameters also describe partly the characters of the interacting system. For instance, the entropy index q which describes the degree of equilibrium can be extracted from the Tsallis distribution [22,23] considering the particle mass. Meanwhile, q can be extracted from the Hagedorn function [24] which is the same as the Lévy-Tsallis function [25] for a particle neglecting its mass. If there is relation between the Tsallis distribution and Hagedorn function, we may say that the former one covers the latter one in which the mass is neglected. Because the universality, similarity or common characteristics exist in high energy collisions [27][28][29][30][31][32][33][34][35][36], some distributions used in large collision system can be also used in small collision system.

II. FORMALISM AND METHOD
i) The Erlang distribution The Erlang distribution is the convolution of multiple exponential distributions. In the framework of multisource thermal model [19][20][21], we may think that more than one parton (or parton-like) contribute to the transverse momentum of considered particle. The j-th parton (or parton-like) is assumed to contribute to the transverse momentum to be p tj which obeys an exponential distribution with the average p t which is j-ordinal number independent. We have the probability density function obeyed by p tj to be The average p t reflects the excitation degree of contributor parton and can be regarded as the effective temperature T .
The contribution of all n s partons to p T is the sum of various p tj . The distribution of p T is then the convolution of n s exponential functions [19][20][21]. We have the p T distribution (the probability density function of p T ) of final-state particles to be the Erlang distribution where N denotes the number of all considered particles and p T has an average of p T = ∞ 0 p T f 1 (p T )dp T = n s p t . Eq. (2) is naturally normalized to be 1. In Eq.
(2), there are two free parameters, n s and p t .
ii) The Tsallis distribution The Tsallis distribution [22,23] has more than one form, which are widely used in the field of high energy collisions. Conveniently, we use the following form where C is the normalization constant, m T = p 2 T + m 2 0 is the transverse mass, m 0 is the rest mass, n = 1/(q −1), and q is the entropy index [22,23]. Eq. (3) is valid only at mid-rapidity (y ≈ 0) which results in cosh y ≈ 1 and the particle energy E = m T cosh y ≈ m T . In Eq. (3), a large n corresponds to a q that is close to 1, and the source or system approaches to equilibrium. The larger the parameter n is, the closer to 1 the entropy index q is, with the source or system being at a higher degree of equilibrium. There is not exactly the minimum n or maximum q [22][23][24][25] which is a limit for approximate equilibrium. Empirically, in the case of n ≥ 4 or q ≤ 1.25 which is 25% more than 1 (even n ≥ 5 or q ≤ 1.2 which is 20% more than 1), the source or system can be regarded as being in a state of approximate (local) equilibrium. Usually, in high energy collisions, the source or system is approximately in equilibrium due to n being large enough.
iii) The Hagedorn function The Hagedorn function [24] is an inverse power-law which has the probability density function of p T to be where A is the normalization constant, n 0 is a free parameter which is similar to n in the Tsallis distribution [22,23], and p 0 is a free parameter which is similar to the product of nT in the Tsallis distribution. Note here that it appears as p 0 = nT is a perfect liquid like relation; however, p 0 is transverse momentum, and n is a dimensionless number. This is not meant in a perfect liquid sense, but the letters are just randomly coinciding. It should be noted that the Hagedorn function is a special case of the Tsallis distribution in which m 0 can be neglected. Generally, at high p T , we may neglect m 0 , observing the two distributions being very similar to each other. At low p T , the two distributions have obvious differences due to non-ignorable m 0 . To build a connection with the entropy index q, we have n 0 ≈ 1/(q − 1). To build a connection with the effective temperature T , we have p 0 ≈ n 0 T ≈ T /(q − 1).
iv) The squared momentum transfer In the center-of-mass reference frame, in two-body reaction 2 + 1 → 4 + 3 (e.g. γ * p → ρp) or in two-body-like reaction, it is supposed that particle 1 is incident along the z direction and particle 2 is incident along the opposite direction. In addition, particle 3 is emitted with angle θ relative to the z direction and particle 4 is emitted along the opposite direction. According to ref. [26], three Mandelstam-variables are defined as where P 1 , P 2 , P 3 , and P 4 are four-momenta of particles 1, 2, 3, and 4, respectively.
In the Mandelstam-variables, slightly varying the form, √ s is the center-of-mass energy, −t is the squared momentum transfer between particles 1 and 3, and −u is the squared momentum transfer between particles 1 and 4. Conveniently, let |t| be the squared momentum transfer between particles 1 and 3. We have Here E 1 and E 3 , p 1 and p 3 , as well as m 1 and m 3 are energy, momentum, and rest mass of particles 1 and 3, respectively. In particular, p 3T is the transverse momentum of particle 3, which is referred to be perpendicular to the z direction. As the energy of incoming photon in the center-of-mass reference frame of the reaction, E 1 in Eq. (8) should be a fixed value. However, E 1 has a slight shift from the peak value due to different experiments and selections. To obtain a good fit, we treat E 1 as a parameter which is the same or has small difference in the same/similar reactions. p 3T obeys one of Eqs. (2)-(4) and θ obeys an isotropic assumption in the center-of-mass reference frame, which will be discussed later in this section. To obtain |t|, we may perform the Monte Carlo calculations. Note that we may calculate |t| from two particles, i.e., particles 1 and 3, but not from one particle. Instead, for one calculation, |t| means the squared momentum transfer in an event. For many calculations, |t| distribution can be obtained from the statistics. For convenience in the description, the transverse momentum and rest mass of particle 3 are also denoted by p T and m 0 respectively. Based on the experiments cited from literature [37][38][39][40][41][42][43][44], we have used two main selection factors for the data. 1) The squared photon virtuality Q 2 = −P 2 γ , where P γ denotes the four-momentum of photon. 2) The center-ofmass energy √ s or W , i.e., W = √ s = −(P 1 + P 2 ) 2 . Let x denote the Bjorken scaling variable, one has v) The initial-and final-state temperatures According to refs. [45][46][47], in a color string percolation approach, the initial temperature T i can be estimated as where p 2 T is the root-mean-square of p T and p 2 T = max 0 p 2 T f 1,2,3 (p T )dp T . In the expression of initial temperature, we have used a single string in the cluster for a given particle production [48], though more than two partons or partons-like take possibly part in the formation of the string. That is, we have used the color suppression factor F (ξ) to be 1 in the color string percolation model [48]. Other strings, even if they exist, do not affect noticeably the production of a given particle. If other strings are considered, i.e., if we take the minimum F (ξ) to be 0.6 [48], a higher T i can be obtained by multiplying a revised factor, 1/F (ξ) = 1.291, in Eq. (9).
The extraction of final-state temperature T 0 is more complex than that of the initial temperature T i . Generally, one may introduce the transverse flow velocity β T in the considered function and obtain T 0 and β T simultaneously [49][50][51][52][53][54][55][56][57], in which the effective temperature T no longer appears. Alternatively, the intercept in T versus m 0 is assumed to be T 0 [50,[58][59][60][61][62][63], and the slope in p T versus m is assumed to be β T [62][63][64][65][66], where m denotes the average energy. However, the alternative method using intercept and slope is not suitable for us due to the fact that the spectra of more than two types of particles (e.g., pions, kaons, and protons) are needed in the extraction which is not our case.
In γ * p, γp, and pp collisions discussed in the present work, the flow effect is not considered by us due to the collective effect being small in the two-body process. This means that T 0 ≈ T in the considered processes. Here, T appears as that in Eq. (3). Meanwhile, T can be also approximated by p t in Eq. (1) and p 0 /n 0 in Eq. (4). Generally, we may regard different distributions or functions as different "thermometers". Just like the Celsius thermometer and the Fahrenheit thermometer, different thermometers measure different temperatures, though they can be transformed from one to another according to conversion rules. Although we may approximately regard T in Eq. (3) as T 0 , a smaller T 0 can be obtained if the flow effect is considered.
As mentioned above, T = p t = p T /n s in Eq. (1) and Erlang distribution, and T i = pT , so this would mean that T is basically encoded in σ 2 pT , the squared variance of p T in the distribution. This also means that T i and T are related through p T . It is understandable, because they reflect the violent degrees of collisions at different stages. Generally, T i > T , this is natural.
Note that although we may use the final-state temperature, it is not a freeze-out temperature for the small system discussed in this paper. In particular, for γ * p and γp reactions, these are just a process describable in terms of perturbative quantum chromodynamics (pQCD) and factorization [67], but not a process in which de-confinement or freeze-out is involved. The meaning of final-state temperature for the large system such as heavy ion collisions or small system such as pp collisions with high multiplicity is somehow different from here. At least, for the large system, we may consider the de-confinement or freeze-out involved picture. Meanwhile, the flow effect in the large system cannot be neglected.
vi) The process of Monte Carlo calculations In analytical calculation, the functions Eqs. (2)-(4) on p T distribution are difficult to be used in Eq. (8) to obtain |t| distribution. Instead, we may perform the Monte Carlo calculations. Let R 1,2 and r 1,2,3,...,ns be random numbers distributed evenly in [0,1]. To use Eq. (8), we have to know changeable p 3T (i.e. p T ) and θ. Other quantities such as E 1 , m 1 , and m 3 in the equation are fixed, though E 1 is treated by us as a parameter with slight variety.
To obtain a concrete value of p T , we need one of Eqs. (2)-(4). Solving the equation where i = 1, 2, and 3, respectively, and δp T is a small shift relative to p T , we may obtain concrete p T . It seems that Eq. (10) directly means that the integral of f 1 (p T ), f 2 (p T ), and f 3 (p T ) is the same for the [0, p T ] interval, which essentially means that the three functions are equal (except for a null measure set). In fact, the three functions are different in forms because of Eqs. (2)-(4), and we need to distinguish them.
In particular, for f 1 (p T ), we have a simpler expression. Let us solve the equation We have due to Eq. (1) being used, where r j in Eq. (12) replaced 1 − r j because both of them are random numbers in [0,1]. The simpler expression is due to p T being the sum of n s random p tj .
To obtain a concrete value of θ, we need the function which is obeyed by θ under the assumption of isotropic emission in the center-of-mass reference frame. Solving the equation we have which is needed by us.
According to the concrete values of p T and θ, and using other quantities, the value of |t| can be obtained from Eq. (8). After repeating the calculations many times, the distribution of |t| is obtained statistically. Based on the method of least squares, the related parameters are obtained naturally. Meanwhile, T i can be obtained from Eq. (9). p T and p 2 T can be obtained from one of Eqs. (2)-(4) or from the statistics. The errors of parameters are obtained by the general method of statistical simulation.

III. RESULTS AND DISCUSSION
A. Comparison with data Figure 1 shows the differential cross-section in squared momentum transfer, dσ/d|t|, of (a) γ * p → ρp, (b) γ * p → ρY , and (c) γ * p → ρ 0 p produced in electron-proton (ep) collisions at photon-proton center-of-mass energy (a)(b) W = 75 GeV and (c) W = 90 GeV, where σ denotes the cross-section and Y in panel (b) denotes an "elastic" scattering proton or a diffractively excited "proton dissociation" [37]. The experimental data points from (a)(b) non-exclusive and (c) exclusive productions are measured by the H1 [37] and ZEUS Collaborations [39], respectively, with different average squared photon virtuality (a) Q 2 = 3.3, 6.6, 11.5, 17.4, and 33.0 GeV 2 , (b) Q 2 = 3.3, 6.6, and 15.8 GeV 2 , as well as (c) Q 2 = 2.7, 5.0, 7.8, 11.9, 19.7, and 41.0 GeV 2 . The data points are fitted by the Monte Carlo calculations with the Erlang distribution Eq. (2) (the solid curves), the Tsallis distribution Eq. (3) (the dashed curves), and the Hagedorn function Eq. (4) (the dotted curves) for p 3T in Eq. (8). Some data are scaled by different quantities marked in the panels for clear visibility. In the calculations, the method of least squares is used to obtain the parameter values. The values of E 1 , p t , n s , T i , T , n, p 0 , and n 0 are listed in Tables 1, 2, or 3 with χ 2 and number of degree of freedom (ndof). One can see that in most cases the calculations based on Eq. (8) with Eqs. (2)-(4) for p 3T can fit approximately the experimental data measured by the H1 and ZEUS Collaborations. Figure 2 presents the differential cross-section in squared momentum transfer, dσ/d|t|, of (a) γp → ωp, GeV, and at (d) proton-proton center-of-mass energy per nucleon pair √ s N N = 29.1 GeV. The experimental data points from (a)(c) exclusive, (b) non-exclusive, and (d) exclusive productions are measured by the ZEUS [40,41], H1 [37], and WA102 Collaborations [43,44], respectively, with different Q 2 for only panels (b) Q 2 = 3.3, 5, 6.6, and 15.8 GeV 2 and (c) Q 2 = 2.4, 3.6, 5.2, 6.9, 9.2, 12.6, and 19.7 GeV 2 . Similar to Figure 1, the data points are fitted by the Monte Carlo calculations based on Eq. (8). The values of parameters are listed in Tables 1, 2, or 3 with χ 2 /ndof. One can see that in most cases the calculations based on Eq. (8) with Eqs. (2)-(4) for p 3T can fit approximately the experimental data measured by the H1 and ZEUS Collaborations. From the above comparisons, we see that some fits have large χ 2 compared to ndof, corresponding to low confidence levels. The parameters obtained from these fits are not representing the data well. We would like to say here that these values are used only for the qualitative description of the data tendencies, but not the quantitative interpretation of the data size. In some cases, ndof < 1, Table 1. Values of E1, pt , ns, Ti, and χ 2 /ndof corresponding to the solid curves in Figures 1-3, where ns = 3 for Figure 1(a) and ns = 4 for other panels, which are not listed in the table to avoid trivialness. In some cases, ndof is less than 1, which is denoted by " − " in the last column, and the corresponding curve is only to guide the eyes. For Figure 3(a), the first and second Q 2 = 6.8 GeV 2 are averaged from the ranges of Q 2 = 2-100 and 5-10 GeV 2 , respectively.  Table 2. Values of E1, T , n, and χ 2 /ndof corresponding to the dashed curves in Figures 1-3, where " − " in the last column denotes the case of ndof < 1 and the corresponding curve is only to guide the eyes.  Table 3. Values of E1, p0, n0, and χ 2 /ndof corresponding to the dotted curves in Figures 1-3, where " − " in the last column denotes the case of ndof < 1 and the corresponding curve is only to guide the eyes. which means that there were at least as many parameters as data points. This means that a perfect fit should have been found. However, this was not the case here. The reason is that we have used given functions, but not any function such as a polynomial.

B. Tendency of parameters
The dependences of energy E 1 of particle 1 on rest mass m 0 of particle 3 for different two-body reactions are given in Figure 4, where panels (a)-(c) correspond to the results from the Erlang distribution, Tsallis distribution, and Hagedorn function, respectively. The types of reactions are marked in the panels. Different symbols represent the results from different reactions or collaborations. One can see that the production of particle 3 with larger m 0 does not need the participation of particle 1 with larger E 1 .
The tendency of E 1 versus m 0 presented in Figure 4 has less fluctuation due to the given collision energy. The results from the three distributions or functions are almost the same, if not equal to each other, due to the same experimental data considered. In fact, E 1 should be a fixed value for given reaction in the present work. However, because different selections such as different Q 2 and W are used in experiments, E 1 has a slight shift from the peak value. Thus, we may regard E 1 as a parameter  [40,41], H1 [37], and WA102 Collaborations [43,44], respectively, with different Q 2 for only panels (b) and (c). Similar to Figure 1, the data points are fitted by the Monte Carlo calculations based on Eq. (8).
and obtain it from the fits.
The dependences of (a) p T , (b) T i , (c) T 0 , (d) n, (e) p 0 , and (f) n 0 on average squared photon virtuality Q 2 for different two-body reactions are shown in Figure 5. The types of reactions are marked in the panels. Different symbols for different reactions represent the parameter values extracted from Figures 1-3 Tables  1-3, where the Erlang distribution, Tsallis distribution, and Hagedorn function in the ranges of available data are used. In particular, p T = n s p t from Table 1 and T 0 = T from Table 2. One can see that p T , T i , T 0 , and p 0 increase generally with an increases in Q 2 , and n and n 0 decrease significantly with increasing of Q 2 .

and listed in
Because of Q 2 being a reflection of hard scale of reaction, this is natural that a harder scale results in a higher excitation degree and then a larger p T , T i , and T 0 . In most cases, one can see a large enough n or n 0 . This means that q is close to 1 and the reaction systems stay in an approximate equilibrium state. At harder scale, the degree of equilibrium decreases due to more disturbance to the equilibrated residual partons in target particle. Then, one has larger q and smaller n or n 0 when compared with those at softer scale. ZEUS Collaboration γ * p→J/ψp W=90 GeV  Tables 1-3. Again, p T = n s p t from Table  1 and T 0 = T from Table 2. One can see that p T , T i , T 0 , and p 0 decrease with an increase in W , and n and n 0 increase with increasing of W .
In γp → J/ψp reactions, at higher center-of-mass energy, the incident photon has higher energy. Although the emitted J/ψ also has higher energy, it is more inclined to have smaller angle. As a comprehensive result, the transverse momentum of J/ψ is smaller, and then T i and T 0 , which are obtained from the transverse momentum are also smaller. In addition, larger n and n 0 at higher collision energy means more equilibrium due to shorter collision time, and then less disturbance to the equilibrated residual partons in target particle. This situation is different from nucleus-nucleus collisions in which cold or spectator nuclear effect have to be considered.
In fact, in nucleus-nucleus collisions, secondary cascade collisions may happen among produced particles and spectator nucleons. The secondary collisions may cause the emission angle to increase, and then the transverse momentum to increase. The effect of secondary collisions is more obvious or nearly saturated at higher energy. In nucleus-nucleus collisions at lower energy, the system approaches equilibrium more easily due to longer interaction time. Conversely, at higher energy, the system does not approach equilibrium more easily due to shorter interaction time for secondary collisions.

C. Further discussion
Before summary and conclusions, we would like to point out that the concept of temperature used in the present work is valid. Generally, the concept of temperature is used in a large system with multiple particles, which stays in an equilibrium state or approximate (local) equilibrium state. From the macroscopic point of view, the systems of γ * p, γp, and pp reactions are indeed small. However, we know that there are lots of events under the same condition in the experiments. These events obey the law of grand canonical ensemble in which the concept of temperature is applicable.
Because the same experimental condition is used in statistics, lots of events are in equilibrium if they consist of a large statistical system which can be described by the grand canonical ensemble. Particles in the large statistical system obey the same distribution law such as the same transverse momentum distribution. From the statistical point of view, particle productions in high energy collisions are a statistical behavior, and the temperature reflects the width of distribution. The higher the temperature is, the wider the distribution is.
The temperature is also a reflection of the average kinetic energy based on large statistical system or a single particle. For single particle, if the distribution law of kinetic energies or transverse momenta is known, the temperature of emission source or interacting system is known, where the source or system means the large ther-  Tables 1-3. Here, pT = ns pt from Table 1 and T0 = T from Table 2.  Tables 1-3. Here, pT = ns pt from Table 1 and T0 = T from Table 2. mal source from the ensemble. Generally, we say the temperature of source or system, not saying the temperature of a given particle, from the point of view of statistical significance of temperature. Based on the temperature, we may compare the experimental spectra of different particles in different experiments. However, different methods have used different distributions or functions, i.e. different "thermometers". To unify these "thermometers" or to find transformations among them, one has to perform quite extensive analysis. Although one may use as far as possible the standard distribution such as the Boltzmann, Fermi-Dirac, or Bose-Einstein distribution to fit the experimental spectra, it is regretful that a single standard distribution cannot fit the experimental spectra very well in general. Naturally, one may use a two-, three-, or even multi-component standard distribution to fit the experimental spectra, though more parameters are introduced.
In fact, the two-, three-, or multi-component standard distribution can be fitted satisfactorily by the Tsallis distribution with q > 1, because the standard distribution is narrower than the Tsallis distribution [68]. In particular, the standard distribution is equivalent to the Tsallis distribution with q = 1. It is natural to use the Tsallis distribution to replace the standard distribution. That is, one may use the Tsallis distribution with q > 1 to fit the experimental spectra and obtain the temperature, though the Tsallis temperature is less than the standard one.
As mentioned in the first section and discussed above, some distributions applied in large collision system can be also applied in small collision system due to the universality, similarity or common characteristics existing in high energy collisions [27][28][29][30][31][32][33][34][35][36]. Based on the same reason, some statistical or hydrodynamic models applied in large system should be also applied in small system. Of course, lots of events are needed in experiments and high statistics is needed in calculation if performing a Monte carlo code.

IV. SUMMARY AND CONCLUSIONS
In summary, the differential cross-section in squared momentum transfer of ρ, ρ 0 ω, φ, f 0 (980), f 1 (1285), f 0 (1370), f 1 (1420), f 0 (1500), and J/ψ produced in γ * p, γp, and pp collisions have been analyzed by the Monte Carlo calculations in which the Erlang distribution, Tsallis distribution, and Hagedorn function (inverse powerlaw) are separately used to describe the transverse momentum spectra of the emitted particles. In most cases, the model results are approximately in agreement with the experimental data measured by the H1, ZEUS, and WA102 Collaborations. In some cases, the fits show qualitatively the data tendencies. The values of the initial and final-state temperatures and other related parameters are extracted from the fitting process. The squared photon virtuality Q 2 and center-of-mass energy W dependent parameters are obtained.
With an increase in Q 2 , the quantities p T , T i , T 0 , and p 0 increase generally, and the quantities n and n 0 decrease significantly. Q 2 is a reflection of hard scale of reaction. A harder scale results in a higher excitation degree, and then a larger p T , T i , and T 0 . In most cases, the reaction system can be regarded as an equilibrium state. At harder scale (larger Q 2 ), the degree of equilibrium decreases due to more disturbance to the equilibrated residual partons in target particle, though the degree of excitation is high.
With increasing of W , the quantities p T , T i , T 0 , and p 0 decrease, and the quantities n and n 0 increase. In γp → J/ψp reactions at high energy, the emitted J/ψ is more inclined to have small angle and hence small p T , T i , and T 0 . In addition, the system stays in a state with higher degree of equilibrium at high energy due to less disturbance to the equilibrated residual partons in target particle. This situation is different from nucleus-nucleus collisions in which the influence of cold or spectator nuclear effect is existent.

Data Availability
This manuscript has no associated data or the data will not be deposited. (Authors' comment: The data used to support the findings of this study are included within the article and are cited at relevant places within the text as references.)

Ethical Approval
The authors declare that they are in compliance with ethical standards regarding the content of this paper.

Disclosure
The funding agencies have no role in the design of the study; in the collection, analysis, or interpretation of the data; in the writing of the manuscript; or in the decision to publish the results.

Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.