Particle physics is the science of the fundamental structure of matter. One of the directions of particle physics is the study of high-energy nuclear matter. The interest in the study of high-energy nuclear matter has increased many folds due to the possibility of studying unstable states of nuclear matter under extreme condition of high energy density and high temperature. The study of non-statistical fluctuations in relativistic nuclear collisions has recently also attracted a great deal of attention due to the possibility of extracting important information about the mechanism of multiparticle production in such collisions   . The multiplicity of charged particles in high energy nucleus-nucleus interactions is an important parameter which indicates how many particles are produced in that interaction. In the same time the multiplicity distributions of produced particles help in learning the interaction mechanism   .
But it is appropriate to study the field in terms of theoretical and the experimental points of view. This will help understand the relationship between theory and experiment, understand the complicated moments which may arise. Theory predicts phenomenon, which can be verified by experiments, and experiments very often provide new insight through unexpected results, which in turn lead to the improvement in theoretical description. In the present paper, an attempt has been made to find the relationship between the stability of a theoretical model and experiment to study the collisions of relativistic heavy ion. This exercise has been made to perform the study of (E-by-E) spatial fluctuations of relativistic shower particles produced in the collisions of 28Si + Em at energy 14.6A GeV in 1-D phase space of -variable (the case -space)  . The findings are compared with the predictions of Ultra-relativistic Quantum Molecular Dynamics (UrQMD) model   .
To solve this problem we use wavelets ideology   . This is due to the fact that various processes and the phenomenon can have quite difficult structure, and contain local features of the various form and time extent. Wavelet-analysis transforms the original data to hierarchical structure by means of the wavelet transformations which results to the set of wavelet coefficients. If the signal is discontinuous, only those wavelets will have high amplitudes, where the maximum value will appear near the discontinuity point, which will allow detecting image contour. At the same time, discontinuity point is a sharp intermittent transition during some process. Quantitatively, it can be estimated by the value of the first derivative of such process, taking into consideration that the first derivative of intermittent transitions is very high. If the transition is in the form of discontinuity point, then the first derivative tends to infinity. The sharper the transition, the higher the derivative value is. Smooth transitions will have small derivative values. This allows us to determine the presence of special characteristics of the analyzed one, as well as the point where these characteristics may arise   . Thus methodology of wavelet analysis allows taking into consideration the particular characteristics study by decomposing source data into a plurality of approximate and detailed coefficients. This allows conducting a detailed analysis of the original data.
2. Some Experimental Details
To obtain the data, FUJI nuclear emulsion pellicles were irradiated horizontally with a beam of 28Si nuclei at 14.6A GeV at Alternating Gradient Synchrophasotron (AGS) of Brookhaven National Laboratory (BNL), New York, USA have been used. The nuclear emulsion experiment is a versatile detector for the study of nuclear reactions in high energy heavy ion collisions. It has the ability to detect and identify the secondary charged particles in the outlet channel of nuclear reactions. The method of line scanning has been adopted to scan the stacks, which was carried out carefully using Japan made NIKON (LABOPHOT and Tc-BIOPHOT) high-resolution microscopes with 8 cm movable stage using 40× objectives and 10× eyepieces by two independent observers, so that the bias in the detection, counting and measurements can be minimized. The interactions due to beam tracks making an angle < 2˚ to the mean direction and lying in emulsion at depths > 35 mm from either surface of the pellicles were included in the final statistics   .
3. Mathematical Formalism for Obtaining Data
In order to perform a meaningful analysis of chaoticity, normalized cumulative variable () were used to reduce the effect of non-uniformity in single charged particle distributions. At the same time, we consider only one component in 1-D phase space of -variable (this is the case -space). The single charged particle density distribution is not flat in the analysis of the fluctuation in phase space variable. This non-uniformity of the particle spectra influences the scaling behaviour of scaled factorial moments. Bialas and Gazdzicki  proposed a method to construct a set of variables, which drastically reduces the distortion of intermittency due to the non-uniformity of single particle density distribution. According to them, the new scaled variable is related to the single particle density distribution by the following relation:
where, is the single particle pseudo rapidity density distribution of the shower particles and and are the two extreme points in the distribution ρ(η) (or, ,). Should also be noted that, in terms of new scaled variable () the single particle density distribution is always uniform in between X = 0 and 1.
Various experimental efforts have established the existence of the empirical phenomenon of “intermittency” in multiparticle production using normalized scaled factorial moments. On the basis of bin averaging the normalized scaled factorial moments of the order of q is defined in vertical form as    :
and its horizontal form is defined as    :
where, , and also bracket of Equation (3) indicates the average over all events in the whole data sample; is the number of relativistic charged particles in the mth bin, m can take values from 1 to M and N represents the total multiplicity of charged shower particles in a particular event in the pseudo-rapidity interval (or).
Recently, Cao and Hwa first introduced to measure the spatial pattern of particles in an event using normalized factorial moments associated with it. In contrast to the horizontally averaged vertical moments, and vertically averaged horizontal moments, of the qth order, they define event factorial moments as  :
where, M is the partition number in phase space, is the number of shower tracks producing particles falling into the mth bin and is the order of the moment.
The event factorial moments, , fluctuates from event-to-event, and the degree of fluctuation can be estimated from the probability distribution over all events. One can obtain a distribution for the whole sample of events. In the given situation, a normalized factorial moment of a single event is defined as   :
where, is the number of events in a sample and represents the event factorial moment describing the spatial pattern of an event. It is important to mention that the SFMs introduced to study the intermittency or fractality in multiparticle production is only an estimate of the mean of the probability distribution of event by event Scaled factorial moment (SFM),. It should be realized that the simple mean procedure, apart from its clear advantages, suppresses a lot of important information about the fluctuations of spatial patterns of final state of multiparticle production. In particular, some interesting effects present only in a part of sample of events produced in high-energy collisions, may be lost.
In order to quantify the degree of the fluctuations, a new normalized moment related to the chaotic nature of the system is defined as    :
where, p is any positive real number, it should not be negative, may vanish for some events, if p is negative.
The other relevant details may be seen in works   .
Thus, we examine the relationship between, and (this technique of erraticity moments  ) for order of moments, and for p = 0.5, 0.9, 1.2, 1.4 and 1.6 for the present experimental data of nucleus-nucleus collisions and with the predictions of Ultra-relativistic Quantum Molecular Dynamics (UrQMD) model. The findings in the forms of the pictorial graphs have been plotted between the natural log of normalized erraticity moments as a function of for -phase spaces respectively at energy 14.6A GeV. These data can be found in  (see Figure 1).
We can see (see Figure 1), which is experimental data and UrQMD data can be varied in a certain range. Thus, we have a few of time series. These time series should be compared using the ideology of wavelets. For this we will use the wavelet coherency.
4. Methods of Wavelet Analysis as Tool to Study Time Series
The basis of formal generalization of continuous wavelet transformation on the time interval is converting the input time series using wavelet parent    :
Figure 1. Variations of as function of in h-space (1D) in the collisions of 28Si + Em at energy 14.6A GeV.
where denotes a normalization, is a location parameter, is a scale parameter,
In this case, feasibility of wavelet analysis in time series study is determined by 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   . That is, the wavelet transformation defines the hierarchical structure of the input of the studied time series. Underlying this statement is the fact that information flows, generated by fractal time series have properties detected by wavelet transform, making it even more informative  . Then wavelet transform allows adding input characteristics of time series that are investigated.
In addition, the empowerment of study time series using wavelet analysis methodology promotes the use of various procedures of wavelet transformation: even-scaled analysis, cross wavelet transformation, wavelet coherence. Then wavelet analysis methodology has been widely used in the disclosure of dynamics of time series that define various data   .
One of the main wavelet transformation methods used for generalized cross-reference analysis between different time series is wavelet coherence. It allows calculating local correlation of two time series in a region of time-frequency. It uses the following formalized model: wavelet coherence as the squared absolute value of the smoothed cross wavelet spectra (time () and frequency () of the time series), normalized by the product of the smoothed individual wavelet power spectra of each series   :
where is a smoothing operator.
We use Morlet wavelet that is a complex wavelet with a good time-frequency localization, as a parent one  -  . The squared wavelet coherency coefficient is in the range, 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.
5. Result and Discussion
We consider the wavelet coherency for each pair of time series of the natural log (experimental data and UrQMD data) for order of moments, and for p = 0.5, 0.9, 1.2, 1.4 and 1.6.
In Figures 2-6 you can check the results of wavelet coherence between selected time series. Each of the following figures indicated a sepa-rate group of time series that match each other. This comparison is displayed in the time-frequency plane. In this case the time scale 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.
On the vertical axis are the weighted characteristics of the analyzed data series in frequency space. Along each of the figures importance scale is presented as separate columns for reflections. Maximum reflection is speaks about full coherence between the experimental data and UrQMD data. Defined lines are a manifestation of localization for individual irregularities within studied time series according to importance of irregularities. In general, each point of wavelet reflects shown in Figures 2-6 is their value in the time-frequency space, which is calculated through wavelet transformation.
The phase difference, indicated by arrows, 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 discrepancy of the experimental data according to the model.
We see that with increasing values q of unbalance occurs between the experimental data and UrQMD data. This is observed for all values p. With increasing the values of p and when is increased the values q imbalance between the experimental data and UrQMD data increasing. Especially it is typical for q = 4. Such
Table 1. The correspondence between a sequence number and value of.
Figure 2. Wavelet coherence between the experimental data and UrQMD data for p = 0.5. (a) q = 2; (b) q = 3; (c) q = 4.
Figure 3. Wavelet coherence between the experimental data and UrQMD data for p = 0.9. (a) q = 2; (b) q = 3; (c) q = 4.
Figure 4. Wavelet coherence between the experimental data and UrQMD data for p = 1.2. (a) q = 2; (b) q = 3; (c) q = 4.
Figure 5. Wavelet coherence between the experimental data and UrQMD data for p = 1.4. (a) q = 2; (b) q = 3; (c) q = 4.
Figure 6. Wavelet coherence between the experimental data and UrQMD data for p = 1.6. (a) q = 2; (b) q = 3; (c) q = 4.
an increase is observed for the mean values. At the same time imbalance between the experimental data and UrQMD data in the frequency domain is shifted from the lowest frequency to the average frequency  -  . Therefore need to have more experimental data to confirm their values. Results which we have obtained are needed to confirm the authenticity of the experiment data and confirm what these data meet certain models.
The results of the study can determine the appropriate application of wavelet analysis methodology as disclosure of wavelet coherency between the studied data series as a tool for the study collisions of relativistic heavy ion. Some comparative results between the experimental data and UrQMD data have been obtained from the analysis of event-by-event fluctuations of produced charged particles in heavy ion collisions at 14.6A GeV. We showed overall consistency between the experimental data and UrQMD data. But we also note that coherence decreases with increasing values of q. This result was obtained for the case -space in 1-D phase space of variable when collisions of relativistic heavy ion. This should be considered when choosing a model for the study and to evaluate the reliability of experimental data.
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  -  .