  We have systematically estimated the possible temperatures obtained from an
analysis of recent data on $p_t$ distributions observed at RHIC experiments.
Using the fact that observed $p_t$ distributions cannot be described by the
original Hagedorn formula in the whole range of transverse momenta (in
particular above 6 GeV/c), we propose a modified Hagedorn formula including
temperature fluctuation. We show that by using it we can fit $p_t$
distributions in the whole range and can estimate consistently the relevant
temperatures, including their fluctuations.
