  The temporal fluctuations of produced hadron density in heavy ion collisions,
modelled by 2D Ising model at temperatures $T_c$ and below, are studied through
a recently developed wavelet based fluctuation analysis method. At $T_c$,
long-range correlated multifractal behavior, matching with the recently
observed Hurst exponent $H\simeq 1$, is found. Below $T_c$ uncorrelated
monofractal behavior is seen. The correlation behavior compares well with the
results obtained from continuous wavelet based average wavelet co-efficient
method, as well as with Fourier power spectral analysis.
