The squared momentum transfer spectra of η and η 0 , produced in high-energy photon-proton (γp) → meson + proton process in electron-proton (ep) collisions performed at CEBAF, NINA, CEA, SLAC, DESY, and WLS are analyzed. The Monte Carlo calculations are used in the analysis of the squared momentum transfer spectra, where the transfer undergoes from the incident γ to emitted meson or equivalently from the target proton to emitted proton. In the calculations, the Erlang distribution and Tsallis-Levy function are used to describe the transverse momentum spectra of emitted particles. Our results show that the average transverse momentum ( pT ), the initial-state temperature (Ti), and the final-state temperature (T0) roughly decrease from lower center-of-mass energy (W ) to higher one, and at several special W , this tendency is broken.