Multiparticle Production in High Energy CollisionsView this Special Issue
Research Article | Open Access
Provash Mali, Soumya Sarkar, Amitabha Mukhopadhyay, Gurmukh Singh, "Wavelet Analysis of Shower Track Distribution in High-Energy Nucleus-Nucleus Collisions", Advances in High Energy Physics, vol. 2013, Article ID 759176, 13 pages, 2013. https://doi.org/10.1155/2013/759176
Wavelet Analysis of Shower Track Distribution in High-Energy Nucleus-Nucleus Collisions
A continuous wavelet analysis is performed for pattern recognition of charged particle emission data in 28Si-Ag/Br interaction at 14.5A GeV and in 32S-Ag/Br interaction at 200A GeV. Making use of the event-wise local maxima present in the scalograms, we try to identify the collective behavior in multiparticle production, if there is any. For the first time, the wavelet results are compared with a model prediction based on the ultrarelativistic quantum molecular dynamics (UrQMD), where we adopt a charge reassignment algorithm to modify the UrQMD events to mimic the Bose-Einstein type of correlation among identical mesons—a feature known to be the most dominating factor responsible for local cluster formation. Statistically significant deviations between the experiment and the simulation are interpreted in terms of nontrivial dynamics of multiparticle production.
The primary objective of studying high-energy heavy-ion interactions is to compress and heat up the nuclear matter beyond the critical values of certain thermodynamic parameters in such a way that the boundaries of individual nucleons melt down to form a thermally and chemically equilibrated color deconfined state of quark-gluon plasma (QGP) [1–3]. As the collision process evolves in space-time, such an exotic state, if formed, subsequently expands and cools down to undergo a reverse transition to the usual hadronic state. In high-energy physics, the process is known as multiparticle production. Each nucleus-nucleus event has its own collision history that ultimately leads to large local fluctuations in the final state particle densities, apparently lacking any definite pattern. In different events, dense clusters of particles are formed at different locations and at different scales of phase-space variables. It is therefore, necessary to formulate a technique that can examine these clusters on an event by event basis. Often, the fluctuations are so large that they can not be explained simply in terms of statistical reasons. It is all very likely that nonstatistical (dynamical) components are present as well, but they are contaminated with trivial noise. With the help of suitable data analysis techniques, it is possible to filter out the genuine clusters of produced particles that in many high-energy interactions are found to scale self-similarly with the phase-space resolution size, approximately following a power law . Global analysis techniques such as the scaled factorial moment method [5, 6], the frequency moment method [7, 8], and the “-parameter” method  have extensively been used to characterize the particle correlations, and efforts are made to interpret the results in the framework of various mechanisms that are mostly speculative in nature. Formation of the QCD parton shower cascade , formation of the disoriented chiral condensate [11, 12], and collective phenomena like the emission of Cherenkov gluons and/or Mach shock wave in the nuclear/partonic medium [13, 14] are examples of some such speculative measures.
The wavelet analysis technique has found its application in many branches of physics [15, 16]. It is capable of revealing the local properties of particle distributions in individual events and at different scales. The technique is, therefore, very appropriate for pattern recognition in multiparticle distribution. In the present paper, we report some results on the wavelet analysis of the angular distribution of shower tracks that are caused by the singly charged produced particles moving with relativistic speed. Our data samples comprise 28Si-Ag/Br events at an incident energy of 14.5A GeV and 32S-Ag/Br events at an incident energy of 200A GeV. Nuclear emulsion technique has been used to collect the experimental data. Several works on the wavelet analysis of multiparticle production at –have so far been reported [17–21]. These works suffer from a common drawback in the sense that there has hardly been any comparison between the experiment and a proper simulation on interaction. It is, therefore, difficult to conclude whether the experimental observations are significant or they are consequences of mere statistical artifacts. We compare our results with the predictions of a microscopic transport model based on the ultrarelativistic quantum molecular dynamics (UrQMD) [22, 23]. It may also be noted that the UrQMD code does not incorporate the Bose-Einstein correlation (BEC) among the identical mesons, a phenomenon considered to be the most dominant factor behind particle cluster formation. Therefore, keeping the phase-space distribution of the produced particles (mostly -mesons) unaltered, we implement a charge reassignment algorithm [24–26] to each UrQMD generated event and thereby try to mimic the BEC into simulation. Any discrepancy between the experiment and the simulation should now be recognized as a genuine collective behavior of the final state particle emission, which has to be interpreted in terms of nontrivial dynamics. Thus, the present analysis on one hand allows us to compare experiments induced by very close projectile masses, while the corresponding values differ by an order of magnitude; on the other hand, it provides an opportunity to compare the experiment with such simulated data where the known dominant source(s) of cluster formation is taken into account. Our paper is organized according to the following sequence: in Section 2, we briefly describe the experiment and the gross characteristics of the data samples used in the paper; in Section 3, we summarily discuss the basic aspects of the UrQMD model and explain the charge reassignment algorithm; in Section 4, without claiming any originality we outline the method of wavelet analysis; in Section 5, we discuss our results—experimental as well as simulated; and in Section 6, we conclude with a critical assessment of these results.
Ilford G5 nuclear photoemulsion pellicles of size 16 cm × 10 cm × 600 m were horizontally irradiated with 28Si beam at an incident energy of 14.5A GeV from the alternating gradient synchrotron (AGS) of the Brookhaven National Laboratory (BNL). Similarly pellicles of size 18 cm × 7 cm × 600 m were irradiated with a 32S beam at an incident energy of 200A GeV from the super proton synchrotron (SPS) at CERN. The primary interactions (also called events/stars) within the emulsion plates are found by following individual projectile tracks, that is, tracks caused by the 28Si and 32S nuclei, along the forward as well as along the backward direction. The process known as line scanning was performed with Leitz microscopes under a total magnification of 300x. On the other hand, Koristka microscopes were utilized for the track counting and angle measurement purposes, for which a total magnification of 1500x was used. The secondary charged particles coming out of an event are categorized in the following way.(i)The shower tracks—caused by the singly charged produced particles most of which are mesons. In an event, their number is denoted by .(ii)The grey and black tracks—resulting from the fragments of the target (Ag/Br) nuclei. Their numbers are denoted, respectively, by and , and the total number denotes the number of target fragments in an event.(iii)The projectile fragments—caused by the spectator parts of the incident projectile (Si/S) nuclei. In an event, their number is denoted by . The details of emulsion experiments, track selection criteria, and data acquisition techniques are nicely elaborated in [27, 28]. To ensure that an interaction involves either an Ag or a Br nucleus as the target, in each event we impose a cut . Thus, altogether 331 28Si-Ag/Br events and 200 32S-Ag/Br events are selected for further analysis, which is confined only to the angular distribution of the shower tracks. The average shower track multiplicity for the 28Si-sample, and for the 32S-sample. The pseudorapidity variable is an approximation of the dimensionless boost parameter rapidity, and it is related to the emission angle of a track as . An accuracy of unit is achieved through the reference primary method of angle measurement. For each data set, the distribution can be crudely approximated by a Gaussian function, whereas the azimuthal angle distributions are in both cases more or less uniform between and . The Gaussian fit parameters for the -distribution in the 28Si-sample are the peak density , the centroid , and the width . For the 32S-sample they are , , and . Due to event averaging, the statistical noise as well as the nonstatistical components of the fluctuations present in individual events are simultaneously smoothed out in the overall distributions. Our basic task is, therefore, (i) to look for statistically significant unusual local structures in the particle distributions in individual events and (ii) to study systematic collective behaviour in large samples of events, if there is any.
To eliminate the background noise, we compare the experiment with the UrQMD (version 3.3p1) model [22, 23]. UrQMD itself does not incorporate any kind of particle correlation, and therefore, in this regard it can be utilized only to generate the statistical background. The rationale behind using a transport model like the UrQMD is that it treats the final freeze-out stage dynamically. It does not make any equilibrium assumption and describes the dynamics of a hadron gas system very well in and out of the chemical and/or thermal equilibrium. In the present case, neither the incident nuclei are too large nor the collision energies are extremely high. Hence, one can not be sure whether local thermal and/or chemical equilibrium are/is achieved. To describe such nonequilibrium many-body dynamics, a transport model is a natural choice. The UrQMD model is applicable over a wide range of energies starting from GeV and ending up at GeV. In this scheme, particle production at high-energy interactions is implemented by the color string fragmentation mechanism similar to that of the Lund model. The UrQMD code has been successfully used to reproduce the particle density distributions and the transverse momentum spectra of various particle species in proton-proton, proton-nucleus, and collisions. However, as mentioned above, the model does not incorporate the symmetry aspects of the fields associated with the produced particles.
It is well known that the Bose-Einstein correlation (BEC), an identical particle effect, dominates the origin of cluster formation. Due to the correlated emission of like sign and/or opposite sign mesons, the particle yield with small relative momenta may be enhanced, which is one of the reasons of large local densities in the final state particles in any high-energy interaction. The effect is quantum statistical in nature and it is not incorporated in the framework of a transport model like the UrQMD. Recently, a new algorithm has been developed [24, 25], where the BEC is introduced by reassigning the charges of produced mesons in such a way that the overall phase-space distribution in each simulated event remains unaltered. The event-wise particle multiplicities are not changed, and it looks as if the particles (mesons) are satisfying the BE statistics. The method of numerically modeling the BEC at the level of a so-called “after burner,” where the output of the UrQMD code is used, is very briefly described below. The UrQMD code provides the four coordinates and the four momenta of all particles. The particle information are contained in an ASCII file written in the OSCAR format. Each particle entry in an event contains a serial number, a particle ID, the particle four momentum , the particle mass , and the final freeze-out four coordinates .(i)In the first step, we arbitrarily choose a meson from an event, and irrespective of its original charge, assign a charge “sign” , − or to it with weight factor . Here, are the numbers, respectively, of the +ve, −ve and neutral mesons, and is the total number of mesons in the event. The chosen meson, say the th one, defines a phase-space cell.(ii)In the next step, we calculate the distances in the four momenta and the four coordinates between the already chosen meson (i.e., the th one) and all other mesons (indexed by ) that are not yet assigned any charge “sign.” Each th meson is associated with a weight factor  which characterizes the bunching probability of the particles in a given cell.(iii)Then, we start to generate uniformly distributed random numbers . If , we reassign the same charge “sign” to the th meson and put it in the same phase-space cell as the th one. We continue the process until either (a) exceeds or (b) all mesons in the event having same charge “sign” as the th one are exhausted.(iv)Now, we go back to our first step and again randomly choose a meson from the pool of the left over mesons for which the charge reassignment has not yet been done. Obviously, the weight factors will now be updated, as some of the particles present in the event are already used up.(v)The algorithm is then repeated until all mesons belonging to each charge variety in the event are used up, and then we move to the next event. Only the meson pairs with space-like separation are accepted, and appropriate checks are imposed so that does not exceed unity . Without changing the overall set of four momenta, four coordinates, or total mesonic charge of the system, we can in this way generate clusters of closely spaced identical charge states of mesons.
We use the UrQMD code in its default setting and generate the minimum bias event samples in the laboratory frame, separately for the Ag and the Br targets and, respectively, for the 28Si and 32S projectiles. For each projectile, the Ag and Br event samples are then mixed up. While doing so, the proportional abundances of these nuclei in the G5 emulsion  are maintained. Only the produced charged mesons are retained for subsequent analysis. From the minimum bias samples, we select subsamples in such a way as to match the respective experimental -distributions. For each projectile, the final sample of simulated events is five times as large as the corresponding experimental one, and the corresponding normalized and/or distributions can be approximately described by more or less the same set of parameters as quoted above for the respective experiment. The UrQMD events are then passed through the charge reassignment algorithm as mentioned above, and from now on these modified samples will be known as the UrQMD + BEC samples.
4. Wavelet Analysis
The wavelet method is used to analyze nonstationary as well as inhomogeneous signals that can be any ordered set of numerically recorded information on some processes, objects, functions, and so forth. A wavelet construction is based on a dilation and a translation parameter. By changing , the local characteristics of a signal are distinguished, while by doing the same with the whole range of a spectrum is analyzed. Unlike the Fourier transformation method which uses only two basis functions, the wavelet transformation method can in principle use an infinite set of discrete or continuous functions as the basis. However, a suitable choice of the basic wavelet is made only after looking at the basic features of the signal to be processed. In the present case, we use a continuous wavelet to find out the strongest fluctuations on an event by event basis that may exceed the expected statistical noise at a particular scale and at a particular point of the underlying phase-space variable, say . The wavelet transform of a function is its decomposition into an orthogonal functional family like where is the wavelet characterized by and as mentioned above, is a normalisation constant, and is the Fourier transform of . Derivatives of the Gaussian function are often used as mother wavelets. In particular, the second derivative: popularly known as the Mexican hat (MHAT) distribution, is customarily used to analyze multiparticle emission data. In Figure 1, we show the plots of and . In the present case, the phase-space variable is and the signal to be analyzed is the density function where is the number of shower tracks in the event sample considered and is the pseudorapidity of the th particle. may either be the value of a single event, or it may be the total number of shower tracks present in the entire event sample/subsample within the interval considered. The wavelet transform of , therefore, becomes is the contribution of to the spectrum in the sense that it represents the probability to find out a particle at some at the scale . A wavelet image at large scale shows only the coarse features, while the same at small reveals the more detailed and finer structures of the underlying distribution.
In Figure 2, we have presented the pseudorapidity spectra of the shower tracks coming out of all 331 28Si-Ag/Br events at an incident energy of 14.5A GeV at different scales (four different values). For comparison, the UrQMD and the UrQMD + BEC graphs are plotted along with those of the experiment. Though the overall multiplicity and the distributions of the simulated and experimental event samples are identical, we observe that the spectra of the experiment are quite different from those of the simulations. The fluctuations are more rapid in the experiment. In Figure 2(a), we can see peaks at and in the experimental distribution. These are the preferred values where particle clusters are formed, and one can relate them, respectively, to the target fragmentation, the central particle producing, and the projectile fragmentation regions. However, we also notice that the central particle producing peak around is well reproduced by the UrQMD + BEC plot. As expected, with increasing the fluctuations are smoothed out, and the distributions gradually converge to the mother wavelet . Needless to say that such plots do not reflect any unique structure of particle distribution in individual events. They would rather correspond to a systematic collective behaviour of the particle emission of the entire sample, if there is any. Similar plots for all the entire 32S-Ag/Br event sample at 200A GeV/c are presented in Figure 3. While the general features of Figures 2 and 3 are more or less similar, we notice that more peaks are present in the 32S-sample than in the 28Si-sample. There are at least 6 prominent peaks within in the experiment, out of which two very prominent peaks are lying within the central particle producing region , and the simulations can not replicate them. Even at a large scale , we find a hump to the left of the peak of the experimental distribution that refuses to be smoothed out, a feature that is absent in the 28Si case. The other peaks, one to the right and three to the left side of the central region, can be related, respectively, to the projectile and the target fragmentations.
The wavelet spectra can be generated for individual events at many different scales that can be used to simultaneously study the location and scale dependence of . We have obtained such distributions for two high multiplicity events, one for a 28Si-Ag/Br event and the other for a 32S-Ag/Br one . We have schematically presented the distributions, respectively, in Figures 4(a) and 4(b). The dark (white) regions in the graphs correspond to the low (high) values. As mentioned before, at the finest scales we only get information about individual particles while at large particles loose their individual identities to coalesce into a big cluster. It is, therefore, pointless to study any event at these two extreme but trivial scales. We see that in both diagrams several small and large clusters are present at . Looking at the 28Si-Ag/Br diagram, we recognize that two large groups of particles are present, one centered around and the other around . Similarly in the 32S-Ag/Br diagram again, there are two large groups, one at and the other at . Beside them several other smaller groups of particles are present, all belonging to the fragmentation regions.
Identification of the peculiarities in particle distribution in individual events from the 2-d energy spectrum is a difficult proposition. Instead we may concentrate on the scalogram defined as which represents the 1-d energy distribution with respect to the scale . A scalogram reflects some of the characteristic features of an event. As, for example, a minimum on it represents the average distance between the particle clusters, while a maximum represents the most compact groups of particles present. Two such scalograms, one for the 28Si-Ag/Br event and the other for the 32S-Ag/Br event considered above, are plotted in Figure 5. In each diagram, a peak or a small rise seen at the lowest scale represents individual particles, and they are of little significance. In the scalogram of the 28Si-Ag/Br event a peak at and a minimum at are seen. On the other hand, in the 32S-Ag/Br event there are a couple of maxima and minima. The maxima are located at and , while the minima are located at and . The simulation, either with or without BEC, fails to reproduce the experiment at the significant scales. It is now amply clear that the scales and the values at which the clusters are formed will vary from one event to the other. Most of the local maxima (minima) are found within , and in most of the events only a few such maxima (minima) are found. To check whether there exists any systematic behaviour of particle emission, or the clusters occur at random, we investigate the distributions of the extremum points over our entire event sample(s). In Figure 6, we plot the frequency distribution of the scales and at which, respectively, the maxima and the minima of the scalograms belonging to individual 28Si-Ag/Br events are graphically seen. The experiments as usual are plotted together with the simulations. Except in Figure 6(b), where the experiment slightly exceeds the simulation at the characteristic scale of , no significant difference between experiment and simulation is observed. In Figure 7, similar histograms for the 32S-Ag/Br events are plotted. In this case also no significant difference between the experiment and the corresponding simulation is seen.
The wavelet analysis is not complete unless we study the distributions of the locations , where the local maxima in are observed. We do this with different choices of scale intervals, cumulative as well as differential. In Figure 8, such distributions for the 28Si-Ag/Br sample (both experiment and simulation) are graphically presented at different cumulative scale windows. The common features of these diagrams are that at the lowest range the distributions are rapidly fluctuating, and as expected with increasing scale window size the fluctuations get reduced. In comparison with the experiment, the UrQMD distributions vary more smoothly. However, when the BEC is incorporated into the UrQMD, to some extent the fluctuating patterns are retrieved. The 32S-Ag/Br sample on the other hand behaves a little differently. The distributions are shown in Figure 9. In this case, the experiment is still more fluctuating than both the UrQMD and the UrQMD + BEC plots. However, incorporating BEC into UrQMD apparently has little effect on the respective distributions. In Figures 10 and 11, the distributions, respectively, for the 28Si-Ag/Br and 32S-Ag/Br samples are once again graphically shown, where we choose differential scale intervals to draw the histograms. For both sets of data, the basic features are more or less the same. As expected at the smallest scale most rapid fluctuations are seen, which are systematically smoothed out with increasing . The distributions for the 32S-Ag/Br interaction are slightly wider than those for the 28Si-Ag/Br interaction. It seems that the inclusion of BEC into the UrQMD in both interactions increases the heights of the local peaks to a small extent.
Pseudorapidity distributions of singly charged particles coming out with relativistic speeds from the 28Si-Ag/Br and 32S-Ag/Br interactions, respectively, at 14.5A GeV and 200A GeV, are analyzed by using the continuous wavelet transform technique. Compared to similar other such emulsion investigations [17–21], the target nuclei in the present case have less uncertainties.
For background noise elimination, the experiments are compared with a set of ordinary UrQMD simulated data, and also with the same set of UrQMD output that is modified by a mimicry of the Bose-Einstein type of correlation. The observed discrepancies between the experiment and the corresponding simulation should, therefore, result from nontrivial dynamics like collective flow of hadronic matter.
Irregularities in the wavelet pseudorapidity spectra, not reproducible by the simulation, are observed in individual events, and the cluster characteristics are not reproducible by the simulations. As far as systematic behavior in many events is concerned, we observed certain differences between experiment and simulation in the 28Si event sample under consideration. The differences with all probability are not a result of ordinary correlations among identical bosons. They should be interpreted in terms of certain nontrivial dynamical reason(s), which are not very much clear from the present analysis.
The present study can be extended to the azimuthal angle distribution of the irregularities and to the 2-d wavelet analysis with larger statistics, so that the impact parameter dependence of the observed irregularities can also be investigated.
Provash Mali acknowledges financial support from the University Grants Commission, Government of India; Project no. 40-510/2011 (SR). Soumya Sarkar is supported by the INSPIRE fellowship of the Department of Science & Technology, Government of India.
- E. V. Shuryak, “Quantum chromodynamics and the theory of superdense matter,” Physics Reports C, vol. 61, no. 2, pp. 71–158, 1980.
- K. Yagi, T. Hatsuda, and Y. Make, Quark-Gluon Plasma, Cambridge University Press, New York, NY, USA, 2005.
- Y.-G. Ma, E. Wang, X. Cai, H.-Z. Huang, X.-N. Wang, and Z.-Y. Zhu, “Quark matter 2006 conference,” Journal of Physics G, vol. 34, 2006.
- W. Kittel and E. A. de Wolf, Soft Multihadron Dynamics, World Scientific, Singapore, 2005.
- A. Bialas and R. Peschanski, “Moments of rapidity distributions as a measure of short-range fluctuations in high-energy collisions,” Nuclear Physics B, vol. 273, no. 3-4, pp. 703–718, 1986.
- A. Bialas and R. Peschanski, “Intermittency in multiparticle production at high energy,” Nuclear Physics B, vol. 308, no. 4, pp. 857–867, 1988.
- R. C. Hwa, “Fractal measures in multiparticle production,” Physical Review D, vol. 41, no. 5, pp. 1456–1462, 1990.
- F. Takagi, “Multifractal structure of multiplicity distributions in particle collisions at high energies,” Physical Review Letters, vol. 72, pp. 32–35, 1994.
- M. I. Adamovich, M. M. Aggarwal, Y. A. Alexandrov et al., “On the jet-like and ring-like substructure in distributions of produced particles in central heavy-ion collisions at ultra-relativistic energies,” Journal of Physics G, vol. 19, pp. 2035–2044, 1993.
- S. A. Bass, B. Müller, and D. K. Srivastava, “Parton rescattering and screening in Au + Au collisions at RHIC,” Physics Letters B, vol. 551, no. 3-4, pp. 277–283, 2003.
- L. Susskind, “Lattice fermions,” Physical Review D, vol. 16, pp. 3931–3039, 1977.
- H. Kluberg-Stern, A. Morel, O. Napoly, and B. Petersson, “Flavours of lagrangian Susskind fermions,” Nuclear Physics B, vol. 220, no. 4, pp. 447–470, 1983.
- I. M. Dremin, “Coherent hadron radiation at ultrahigh-energies,” Pisma v Zhurnal Eksperimentalnoi i Teoreticheskoi Fiziki, vol. 30, pp. 152–156, 1979, English translation Journal of Experimental and Theoretical Physics Letters, vol. 30, pp. 140–147, 1979.
- J. Hofmann, H. Stöcker, U. Heinz, W. Scheid, and W. Greiner, “Possibility of detecting density isomers in highdensity nuclear mach shock waves,” Physical Review Letters, vol. 36, pp. 88–91, 1976.
- C. K. Chui, An Introduction to Wavelets, Academic Press, New York, NY, USA, 1992.
- D.-W. Huang, “Wavelet analysis in multiplicity fluctuations,” Physical Review D, vol. 56, no. 7, pp. 3961–3969, 1997.
- I. M. Dremin, O. V. Ivanov, S. A. Kalinin, K. A. Kotelnikov, V. A. Nechitailo, and N. G. Polukhina, “Wavelet patterns in nucleus-nucleus collisions at 158A GeV,” Physics Letters B, vol. 499, no. 1-2, pp. 97–103, 2001.
- V. V. Uzhinsky, V. S. Navotny, G. A. Ososkov, A. Polanski, and M. M. Chernyavsky, “Wavelet analysis of angular distributions of secondary particles in high-energy nucleus-nucleus interactions: irregularity of particle pseudorapidity distributions,” Physics of Atomic Nuclei, vol. 67, pp. 156–162, 2004.
- J. Fedorišin and S. Vokál, “Search for the righ-like structures in the emission of secondary particles in central 197Au collisions with emulsion nuclei at 11.6A GeV/c,” JINR Preprint E1-2007-4, Dubna, Russia, 2007.
- J. Fedorisin and S. Vokal, “Wavelet analysis of angular spectra of relativistic particles in 208 Pb induced collisions with emulsion nuclei at 158A GeV/c,” JINR Preprint E1-2007-66, Joint Institute for Nuclear Research (JINR) publishing department, Dubna, Russia, 2007, http://www1.jinr.ru/Preprints/prepr_2007_eng.html.
- J. Fedorišin and S. Vokál, “Wavelet analysis of multiparticle correlations,” FIZIKA B, vol. 17, no. 2, pp. 273–278, 2008.
- S. A. Bass, M. Belkacem, M. Bleicher et al., “Microscopic models for ultrarelativistic heavy ion collisions,” Progress in Particle and Nuclear Physics, vol. 41, pp. 255–369, 1998.
- M. Bleicher, E. Zabrodin, C. Spieles et al., “Relativistic hadron-hadron collisions in the ultra-relativistic quantum molecular dynamics model,” Journal of Physics G, vol. 25, no. 9, pp. 1859–1896, 1999.
- O. V. Utyuzh, G. Wilk, and Z. Wlodarczyk, “Numerical modelling of Bose-Einstein correlations,” Physics Letters B, vol. 522, no. 3-4, pp. 273–279, 2001.
- O. Utyuzh, G. Wilk, and Z. Wlodarczyk, “Modeling Bose-Einstein correlations via elementary emitting cells,” Physical Review D, vol. 75, Article ID 074030, 2007.
- M. Bystesrky, “Modeling Bose-Einstein correlation in heavy ion collision at RHIC,” Nukleonika, vol. 49, supplement 2, pp. S37–S39, 2004.
- C. F. Powell, P. H. Fowler, and D. H. Perkins, The Study of Elementary Particles by Photographic Method, Pergamon Press, Oxford, UK, 1959.
- W. H. Barkas, Nuclear Research Emulsions Vol. I and II, Academic Press, New York, NY, USA, 1963.
Copyright © 2013 Provash Mali et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.