The main objective of relativistic heavy ion collisions lies in the investigating of properties of nuclear matter consisting of strongly interacting particles at extreme condition of temperature and density. The most important theory from strongly interacting matter i.e. quantum chromodynamics (QCD)   suggests that at a critical temperature of approximately 200 - 250 MeV, a phase transition from hadronic matter to a novel state of matter that is quark gluon plasma (QGP) occurred which consists of the deconfined quarks and gluons with in a small volume. The energy density for this new phase transition should be a few GeV/fm3 (~2 GeV/fm3) compared with 0.17 GeV/fm3 of normal nuclear matter in its ground state   . Observations of such type of phase transition would confirm the predictions of QCD and it is also of prime importance in Astrophysics. According to Big Bang theory the initial state of matter consists of dense and hot plasma of unconfined quarks and gluons: after that cosmological expansion process takes place and the energy density decreases, the plasma cools down and transforms into colorless objects      . This process is called harmonization. Therefore one can expect that the observation of quarks-gluon plasma (QGP) would reproduce in laboratory  , the hypothetical state of matter that prevailed in early universe. According to phase diagram shown in Figure 1 of strongly interacting matter following points may be noted about the newly state of matter (QGP)  .
(i) Low temperature and high nuclear energy density is believed to have existed in the early Universe and today it may exist in the heart of neuron star. One can achieve it directly.
(ii) High temperature and low nuclear energy density, which is applicable only by relativistic heavy ion collisions experiments in laboratory which is called Little Bang i.e. Relativistic Heavy Ion Collisions (RHIC) and Large Hadron Collider (LHC) and so many future experiments.
Already now there are plenty of studies that describe different experiments. In references   showed that the high energy nucleus-nucleus (A-A) collisions provide a system in which the properties of hot, dense strongly interacting matter can be achieved. Also there have been numerous experimental results, which indicate that collisions of hadron-nucleus (h-A), nucleus-nucleus (A-A) interactions
Figure 1. A schematic phase diagram of strongly interacting matter  .
cannot be completely understood in term of superposition of nucleon-nucleon scattering    . There are many studies where such experimental results have been compared with the data generated with the computer code FRITIOF based on Lund Monte Carlo Model for high energy nucleus-nucleus collisions   . The measured phenomena, such as strangeness enhancement and J/y suppression show that extremely dense strongly interacting matter has been found  . But still of such research is not enough in order to completely understand the nature of heavy ion collisions at ultra-relativistic energies.
One significant area of research is the study of the dynamics of the chaotic behavior in nuclear collisions at ultra-high energies. Here it is necessary to understand such things: the accuracy of the estimates of the chaotic behavior in relativistic heavy ion collisions in different phase spaces and the connection of such assessments with each other. These questions and is dedicated to present study. To solve this problem we use wavelets ideology   . This is due to the fact that of the chaoticity and/or errticity pattern in ultra-relativistic heavy ion collisions can have quite difficult structure, contain local features of the various form and time extent.
Finally, In this article, an attempt has been made to study the wavelet methodology over a chaotic system to understand the underlying relation of the chaotic behavior of relativistic shower particles produced in 28 Si-emulsion interactions at 14.6 AGeV. This will help to extract the valuable information’s about the mechanism of multi-particle production in different phase spaces; h-space, f-space (in one dimension) and hf-space (in two dimensions). The merits of wavelet analysis over chaotic behavior were foundto be in good agreements.
2. Brief Detail of Experimental Data
The track detectors (so called nuclear emulsion detector) have been extensively used in high energy physics over many decades. This nuclear emulsion detector also called 4p detector and used for recording of secondary charged particles and is accompanied by the emergence of observable traces (track producing particles) corresponding to the elementary particle trajectories. Therefore in the present experiment, a FUJI film detector has been used to pick the data samples. In this experiment a beam(projectile) of 28Si nuclei at total energy » 409 GeV hit to the heterogeneous mixture/various nuclei (fixed target) of nuclear emulsion at Alternating Gradient Synchro-phasotron (AGS) of Brookhaven National Laboratory (BNL), NewYork, USA.The other relevant details about the present experiments and target identifications may be seen in our earlier publications  -  . The projectile spectator fragments, singly and multiple charged, emitted inside fragmentation cone, usually we determine the number of alpha particles (nα) and the number of fragments with Z > 2, separately  .
The polar (q) and azimuthal ( f) emission angles of all tracks have been measured, and the pseudo-rapidity has been calculated by for each shower particle.
This study has been carried out for the experimental data along with the theoretical prediction of Ultra-relativistic Quantum Molecular Dynamics model (UrQMD) and Monte-Carlo (RanMC) simulation and the total number of events/data points were 10000.
3. Mathematical Approach
3.1. Analytical Tool for Wavelet Methodology
We made an approach for the generalization of the constant wavelet transformation over time interval “t” which later on changed into an input time series such as; with a wavelet parent    :
where, the usual meanings to symbols such as: represents to normalization processes, “u” indicates to location parameter and “s” is a scale parameter.
The wavelet analysis transforms the original data to a classified structure by which we get a new set of wavelet coefficients. This is done using so-called multiresolution analysis method. The multiresolution analysis consists in splitting of an investigated number into two components―approximating and detailing, with their subsequent crushing for the purpose of change of level of expansion of a signal to the set level of expansion. The importance of wavelet analysis in the study of time series is to determine the fact that the method of wavelet analysis allows to discover the local features of the studied time series due to the decomposition of the input data   . And finally, this allows determining the presence of special characteristics of the analyzed data, as well as the point where these characteristics may arise  . On the other hand, the wavelet analysis allows conducting a detailed analysis for the original data.
In practice, it is quite spread the use of discrete wavelet transformation (DWT)  . It is connected with that application DWT becomes especially effective when the signal has high-frequency components of short duration and extensive low-frequency components. In the consent of time series study of the wavelet analysis that, it stimulates for the various methods of wavelet transformation and their application in different sectors such as: scaled analysis, cross wavelet transformation, wavelet coherence.
But the leading role of the wavelet transformation methods is to use for generalized cross-reference analysis between different time series is wavelet coherence. The wavelet coherency simultaneously assess how the co-movement and causalities between two variables vary across different frequencies involved and change over time in a time-frequency window. Moreover, co-movement of the time series is recognizable among different time scales, which the standard approaches, failed to perform. It allows to calculate local correlation of two time series (x and y) in a region of time-frequency. Hence, one can use the following formalized model; e.g. the wavelet coherence as the squared absolute value of the smoothed cross wavelet spectra , should normalized by the product of the smoothed individual wavelet power spectra of each series   :
where; is a smoothing operator.
In the present work, we used the Morlet wavelet that was a complex wavelet with good time-frequency localization, as a parent one   . The squared wavelet coherency coefficient was in the rangeof , and it was found that the values close to zero indicate weak correlation, while values close to one are evidences of strong correlation.Thus, wavelet coherency analysis enables interconnection between the studied time series and analyzes the frequency of such communications. The Monte Carlo methods are used here to properly find statistical level of significance of the wavelet coherence.
3.2. Analysis for the Chaotic Behavior
To eliminate the effect of non-uniform density distribution, we used the proposed method of Bialas and Gazdzicki  . In this method the original values of pseudo-rapidity (h) distribution is transformed into a new cumulative variable ( ). And this was demonstrated by following mathematical relation:
where the usual meanings of such above parameters were; is the single particle (h) distribution of the secondary charged pions and and are the two extreme points in the density distribution of . The effect of this analysis the corresponding region of investigation the interval of pseudo-rapidityDh, become changed in 0 to 1 order of a variable .
Further, for the study of event-to event multiplicity fluctuations in relativistic nuclear collisions, it is necessary to investigate the performance of the event factorial foment, in small bins of variable . The factorial moments (FMs) for the order “q” of an event was mathematically represented such as   :
where “M” is the partition number of variable , nm is the number of produced secondary charged particles into the mth bin and q = 2 - 6 is the order of the moment.
Since the fluctuations of from event-to-event, we get a distribution of denoted by after a large number of events. The normalized factorial moments of a particular event can now be defined such as:
It is worth mentioning here that the scaled factorial moments (SFMs) defined by Bialas and Paschanski   to study of intermittency in nuclear collisions is only an estimate of mean of the distribution . It should be realized that the averaging procedure, apart from its clear advantages, brings also a danger of losing some important information on spatial patterns from event-to-event. In particular, some interesting effects, if present only in a part of sample of events produced in high energy collisions, may be lost. A possible example of this kind is the quark gluon plasma (QGP), which is expected to be characterized by specific intermittency exponents   . It is therefore, essential to investigate the full shape of the distribution and the way it changes with the bin size.
It has been argued  -  that the event-to-event fluctuations can probe the dynamics of multiparticle production more deeply than the variables such as the multiplicity distribution and the average factorial moments. Therefore, we study the chaoticity or event-to-event fluctuations in the density of particle produced in 28Si-emulsion collisions.
Theerraticitymoments, that quantify these fluctuations have been determined     such as following relation in Equation(8). Expediency of use of this approach is discussed in the articles     .
where, the usual meaning of various symbols/and or parameters are such as: the “p” is any positive real number, ―a normalized factorial moment of a single event that is associated with the component in 1-D phase space (h-space and f- space) (     and references therein] of ―variable ( ―a variable that characterizes the number of relativistic heavy ion collisions).
4. Results and Discussions
4.1. The Chaoticity in Relativistic Heavy Ion Collisions
First we have calculated the Chaoticity, so called chaotic behavior from the present experimental data. For this task, the event factorial moments and the erraticity moments, were calculated for order of q = 2 - 4 and parameter “p” = p = 0.5, 0.9, 1.2, 1.4 and 1.6 by using the Equation (7) and Equation (8), the values of “M” was varied from 1 to 40. And the outcomes of all above these calculations were represented in Figure 2, for different phase spaces; h-space, f-space (in one dimension) and hf-space (in two dimensions) respectively. Hence the dependence of lnCp,q(M) as function of lnM in different phase spaces for the produced secondary charged particle in the collisions of 28Si (projectile) with fixed target of Emulsion nuclei at energy 14.6A GeV has been depicted in Figure 2. Some part of these results already published by the author M. Ayaz Ahmad et al., in reference  . From the Figure 2, it can be seen that the expe-
Figure 2. Variations of lnCp,q(M) as function of lnM in h-space (1D) (f-space (1D) and hf-space (2D)) in the collisions of 28Si+Em at energy 14.6A GeV.
rimental data of Cp,q lie well above those for the generated uncorrelated (Monte Carlo) data in the region of “M ³ 13. It means that the contribution of the statistical fluctuations to Cp,q values in present experimental data is small.
Shaoshun and Zhaomin  and Jinghua Fu et al.,  have studied the chaoticity in NA27 data on p-p collisions at 400 GeV/c. They pointed out that the observed chaoticity could be reproduced by the statistical fluctuations only. But the chaoticity observed in the present experimental work cannot be reproduced by the statistical fluctuations only. Thus the chaotic behavior observed in our data and has the dynamical origin.
4.2. Results from Wavelet Methodology
Thus, we have a few of time series. These time series should be compared using the ideology of wavelets. For this important task we will use the wavelet coherency. We consider the wavelet coherency for each triple of time series for the produced particles in the relativistic nuclear collisions for various phase spaces; h-space, f-space (in one dimension) and hf-space (in two dimensions) respectively, where the factorial moments order was q = 2 to 4, and the parameter “p” was of the order of, 0.5, 0.9, 1.2, 1.4 and 1.6.
The findings of this analysis were depicted in Figures 3-17. One can check the results of wavelet coherence between selected time series in these Figures 3-17. Each of the following figures indicated a separate group of time series that match each other. In this case the time scale (x-axis) is equal to the consistent change of values (these changes represented a sequence number ). The correspondence between a sequence number and value of is shown in Table 1. In Figures 3-17 (from left to right): the first picture (a) is wavelet coherence between h-space and f-space, the second picture (b)―wavelet coherence between h-space and hf-space, the third picture (c)―wavelet coherence between f-space and hf-space.
Further it has been find that on the vertical axis were the weighted features of the analyzed data series in the frequency space. Corresponding to each of the figures there was a significant scale and is obtainable as separate columns for reflections. The wavelet coherence illustrates the regions in the time-scale space where the central variables co-vary (but do not necessarily have high power). The maximum reflection is addressed about the full coherence between the data which were analyzed. By definition, the regions inside the black lines plotted in warmer colors indicate a strong interdependence between the investigated time series. The colder color indicates relatively weak co-movement between these variables. The defined lines were a sign of localization for individual irregularities within studied time series according to importance of irregularities. The
Table 1. The correspondence between a sequence number and value of .
(a) (b) (c)
Figure 3. Wavelet coherence for the experimental data for p = 0.5 and q = 2.
(a) (b) (c)
Figure 4. Wavelet coherence for the experimental data for p = 0.9 and q = 2.
(a) (b) (c)
Figure 5. Wavelet coherence for the experimental data for p = 1.2 and q = 2.
(a) (b) (c)
Figure 6. Wavelet coherence for the experimental data for p = 1.4 and q = 2.
(a) (b) (c)
Figure 7. Wavelet coherence for the experimental data for p = 1.6 and q = 2.
(a) (b) (c)
Figure 8. Wavelet coherence for the experimental data for p = 0.5 and q = 3.
(a) (b) (c)
Figure 9. Wavelet coherence for the experimental data for p = 0.9 and q = 3.
(a) (b) (c)
Figure 10. Wavelet coherence for the experimental data for p = 1.2 and q = 3.
(a) (b) (c)
Figure 11. Wavelet coherence for the experimental data for p = 1.4 and q = 3.
(a) (b) (c)
Figure 12. Wavelet coherence for the experimental data for p = 1.6 and q = 3.
(a) (b) (c)
Figure 13. Wavelet coherence for the experimental data for p = 0.5 and q = 4.
(a) (b) (c)
Figure 14. Wavelet coherence for the experimental data for p = 0.9 and q = 4.
(a) (b) (c)
Figure 15. Wavelet coherence for the experimental data for p = 1.2 and q = 4.
(a) (b) (c)
Figure 16. Wavelet coherence for the experimental data for p = 1.4 and q = 4.
(a) (b) (c)
Figure 17. Wavelet coherence for the experimental data for p = 1.6 and q = 4.
areas delimited by the black lines cover coherence values significant at the 5% level. In general, each point of wavelet reflects and clearly shown in Figure 3-17. One can measured their values in the time-frequency space, and it is calculated through wavelet transformation. The phase difference, indicated by arrows, that gives us details about delays of oscillation of the two examined time series. Arrows pointing to the right (left) when the time series are in-phase (anti-phase) or are positively (negatively) correlated. Arrow pointing up means that the first time series leads the second one, arrow pointing down indicates that the second time series leads the first one. In our case, these data suggest interrelations between different spaces (h-space, f-space (in 1D) and hf-space (in 2D) respectively) and the impact of this relationship on the chaotic behavior in relativistic heavy ion collisions. Thus this implies that wavelet plots enable to appropriately identify both frequency bands and time intervals where data series move together.
The Figures 3-7 show the results of wavelet coherence for q = 2. We see that between f-space and hf-space there is full consistency of the chaotic behavior in relativistic heavy ion collisions. Between h-space and f-space and between h-space and hf-space wavelet coherence varies depending on the value .
Thus, of the chaotic behavior in relativistic heavy ion collisions from the point of view of different spaces is changing. The smallest values the wavelet coherence is observed for the mean values . With increasing values p, violation of wavelet coherence to first decreases (Figures 3-5) and then increases (Figure 6 and Figure 7).
Figures 8-12 shows the results of wavelet coherence for q = 3.We see that Figures 3-7 coincide with Figures 8-12. But we can to observe and minor differences (Figure 6 and Figure 11). It is connected with an existing fluctuation of the original data (Figure 6). At the same time it confirms the validity of the results that were obtained. In general, can be seen that with increasing q (from q = 2 to q = 3) significant changes of the chaotic behavior in relativistic heavy ion collisions between different phase spaces does not occur.
Figures 13-17 shows the results of wavelet coherence for q = 4.
We see that between h-space and hf-space there is full consistency of the chaotic behavior in relativistic heavy ion collisions. Between h-space and f-space and between f-space and hf-space wavelet coherence varies depending on the value . Other characteristics of the wavelet coherency have not changed. Thus, with increasing q (from q = 3 to q = 4) we can observe changes in the interaction phase spaces in the study of the chaotic behavior in relativistic heavy ion collisions. Therefore with the change in value of q there is a change influence of the chaotic behavior in relativistic heavy ion collisions, which is reflected in the change the values of in the phase space.
5. Conclusions and Final Remark
The relativistic heavy ion collisions provide an experimental setting for the study of the exotic behavior of the matter. It has been suggested that the strongly interacting matter at high energy densities produced in these collisions may undergo a phase transition to quark gluon plasma (QGP). The produced particles in such nuclear collisions are believed to carry relevant information about the collision mechanism.
We have received different wavelet coherency for the chaotic behavior in relativistic heavy ion collisions between different phase spaces h-space, f-space (in one dimension) and hf-space (in two dimensions) respectively). We have shown changes in the wavelet coherence when changing the values of q and p. This must be considered when analyzing the data.
The ideology of wavelet analysis and/or methodology was found an important tool for the wavelet coherency between the studied data series and it can be implemented very well for the study of relativistic heavy ion collisions.
The authors would like to acknowledge the keen support for this work of the Department of Physics, Faculty of Science, University of Tabuk, Saudi Arabia and also the Department of Informatics, Kharkiv National University of Radio- Electronics, Kharkiv, Ukraine.
 Ahmad, S. and Ahmad, M.A. (2006) A Comparative Study of Multifractal Moments in Relativistic Heavy-Ion Collisions. Journal of Physics G Nuclear and Particle Physics, 32, 1279-1293.
 Mohery, M. and Abd-Allah, N.N. (2002) Systematic Comparison of the Experimental Data with the Fritiof Model in Nucleus-Nucleus Interactions with Light and Heavy Target Nuclei in Nuclear Emulsion at 4.5 A GeV/c. International Journal of Modern Physics E, 11, 161-175.
 Mohery, M. (2012) Characteristics of the Total Disintegration Events of Emulsion Heavy Target Nuclei Caused by O16 and 28Si Nuclei at High Energies. Canadian Journal of Physics, 90, 1267-1278.
 Andersson, B., Gustafson, G. and Nilsson-Almqvist, B. (1987) A Model for Low-pT Hadronic Reactions with Generalizations to Hadron-Nucleus and Nucleus-Nucleus Collisions. Nuclear Physics B, 281, 289-309.
 Lyashenko, V., Kobylin, O. and Ahmad, M.A. (2014) General Methodology for Implementation of Image Normalization Procedure Using Its Wavelet Transform. International Journal of Science and Research, 3, 2870-2877.
 Ahmad, S., Ayaz Ahmad, M., Irfan, M. and Zafar, M. (2006) Study of Non-Statistical Fluctuations in Relativistic Nuclear Collisions. Journal of the Physical Society of Japan, 75, Article ID: 064604.
 Tariq, M., Ayaz Ahmad, M., Ahmad, S. and Zafar, M. (2007) Analysis of High NS-Multiplicity Events Produced in Relativistic Heavy Ion Collisions at 4.5A GeV/c. Romanian Reports in Physics, 59, 773-790.
 Ahmad, S., Ayaz Ahmad, M., Tariq, M. and Zafar, M. (2009) Charged Multiplicity Distribution of Relativistic Charged Particles in Heavy Ion Collisions. International Journal of Modern Physics E, 18, 1929-1944.
 Ayaz Ahmad, M., Ahmad, S. and Zafar, M. (2010) Intermittent and Scaling Behaviour of Shower Particles Produced in the Collisions of 28Si-Em at 14.6A GeV. Indian Journal of Physics, 84, 1675-1681.
 Kingsbury, N. (1999) Image Processing with Complex Wavelets. Philosophical Transactions of the Royal Society of London A: Mathematical, Physical and Engineering Sciences, 357, 2543-2560.
 Delbeke, L. and Abry, P. (2000) Stochastic Integral Representation and Properties of the Wavelet Coefficients of Linear Fractional Stable Motion. Stochastic Processes and Their Applications, 86, 177-182.
 Veitch, D., Flandrin, P., Abry, P., Riedi, R. and Baraniuk, R. (2002) The Multiscale Nature of Network Traffic: Discovery, Analysis, and Modelling. IEEE Signal Processing Magazine, 19, 28-46.
 Grinsted, A., Moore, J.C. and Jevrejeva, S. (2004) Application of the Cross Wavelet Transform and Wavelet Coherence to Geophysical Time Series. Nonlinear Processes in Geophysics, 11, 561-566.
 Ayaz Ahmad, M., Rasool, M.H., Ahmad, S., Ahmad, N.A. and Madani, J.H. (2013) Chaotic Behaviour of Multiparticle Production in Relativistic Heavy Ion Collisions. International Journal of Enhanced Research in Science Technology & Engineering, 2, 77-89.