The analysis of the experimental data in terms of multiplicity distributions for different emitted secondaries (slow and fast protons) is one of the main sources of information about the mechanism of particle production. Figures 1(a)-(c) shows the multiplicity distributions of black, grey and heavily ionizing particles from ³²S- Emulsion interactions at 200 AGeV along with the distribution obtained from ²⁸Si-Em [8] and ¹⁶O-Em [9] interactions at 14.6 AGeV and 200 AGeV respectively for comparison. It is observed from the figures that the peaks of the distributions appear in the lower values of N_b, N_g, and N_h. These distributions seem to be independent of incident energy as well as projectile mass within statistical errors up to lower values of N_b, N_g, N_h. This result is consistent with those obtained by other workers [8] - [13] . It may also be noticed from the figures that the per-

centage of events with large values of N_b, N_g, or N_h increases with projectile mass. Finally, it may be concluded from the multiplicity distributions of slow and fast protons produced in nucleus-nucleus interactions that no significant differences are observed regarding the mechanism of their production with energy.

Figures 2(a)-(c) and Figures 3(a)-(c) show the N_b, N_g, and N_h multiplicity distributions from ³²S-AgBr and ³²S-CNO interactions at 200 AGeV. The results obtained by other workers for ¹²C and ²⁸Si at 4.5 and 14.6 AGeV [8] respectively are incorporated in the same figure for comparison. From the figures, it may be observed that the characteristics feature of the distribution is similar in shape for all targets, but the multiplicity range or decay tail increases with target size. It has also been found that the distributions for ³²S-AgBr are broader than those for ³²S-CNO interactions. This may reflect the effect of the target mass number on the number of collisions of ³²S beam with target nuclei. Experimental results obtained by other workers [8] [10] reproduce qualitatively similar results.

The possibility of scaling, i.e., similarity in the multiplicity distributions of grey tracks produced in hadron-nucleus and nucleus-nucleus interactions has also been examined. In the present analysis, the events with N_g = 0 have been excluded because the coherent processes may also contribute to such events. Figure 4(a)

shows the -distribution from ³²S-Em interactions at 200 AGeV. A straight line of the form:

is found to represent the present data, N_ev denotes the number of events in the given bin and A and B are constants. The best fit to the data is given as:

For comparison the -distributions from ²⁸Si-Em at 14.6 AGeV and 4.5 AGeV and ¹²C-Em at 4.5 AGeV [14] are also plotted in Figures 4(b)-(d). The values of the slopes for ²⁸Si-Em at 14.6 AGeV and ²⁸Si-Em and ¹²C-Em interactions at 4.5 AGeV are found to be −0.88 ± 0.04, −0.89 ± 0.08 and −0.86 ± 0.06 respectively. The slope parameters are ~0.90, which are very much consistent with the value obtained for ³²S-Em at 200 AGeV. The constancy in the values of slopes for nucleus-nucleus collisions at different energies may be interpreted as existence of some kind of scaling for the production of grey tracks.

The mean multiplicity is the average number of charged particles produced in various types of high-energy heavy ion collisions. Multiplicity of different charged particles is helpful in understanding the mechanism of multiparticle production. The average values of number of black, grey and heavily ionizing particles produced in ³²S-Em interactions at 200 AGeV are displayed in Table 1. The values of, and measured in ³²S-Em interactions at 200 AGeV by S. Dhamija et al. [2] and A. Dabrowska et al. [15] along with other results at different energies [8] - [10] [12] [13] [16] - [18] are also listed in Table 1 for the sake of comparison. It is

seen that the values of obtained in ³²S-Em interactions at 200 AGeV are slightly larger than the values obtained by S. Dhamija et al. and A. Dabrowska et al. This discrepancy may be due to the different criteria of ionization chosen for black tracks by S. Dhamija and A. Dabrowska. The values of compare reasonably well. It can be observed from the results shown in the table that the value of depends weakly with increasing mass of the projectiles as well as energy of the projectiles, whereas the values of do not exhibit any such trends. The increasing trend is missing in the values of because slow particles are produced due to evaporation of excited residual nucleus. We also study the dependence of the average multiplicities on projectile mass, A_p, using the following power law:

where i = b, g, h, and represents the mean multiplicity. The values of the slope (b) for black, grey, and heavy tracks obtained from the least - squares fits are given as (−0.0054 ± 0.009), (0.1407 ± 0.008) and (0.1125

± 0.016) respectively. The dependence of, and on the projectile mass number, A_p, is shown in Figures 5(a)-(c). It has been reported that the values of b_g for is slightly higher than the cor-

responding values in case of hadron-nucleus interactions [18] [19] . The observation may indicate that the ³²S projectile at a given impact parameter are an extended object rather than a point object as in the case of a hadron beam. The variation of with A_p is shown in Figure 5(a), which is almost independent of the projectile mass number and its energy. This constancy of indicates that the average excitation of the residual target nucleus has reached its saturation at present projectile energy. The strong dependence of the total charged secondary particles on the masses of colliding nucleus are due to the increase in the overlapping region of the two interacting nuclei.

Table 2 presents the mean multiplicities of black, grey and heavily ionizing particles for ³²S-emulsion inelastic interactions at 200 AGeV for different N_h-intervals. It is found that mean multiplicities of slow particles in-

creases with increase of centrality of collisions. A regular pattern in the values of the ratio has been recorded from the table except for CNO (2 ≤ N_h ≤ 7) events which is slightly higher.

The average values of dispersion D(N_g), and for different N_h-intervals of ³²S-Em in-

teractions at 200 AGeV are given in Table 3. From the table it may be clearly noticed that as the impact parameter decreases (the degree of disintegration of the target nuclei increases), the ratio of the number of slow evaporated particles (N_b) to heavily ionizing particles (N_h) i.e. the is nearly constant. This result contradicts the results obtained by Antonchilk et al. [20] for 7 ≤ N_h ≤ 27 and N_h ≥ 28 respectively.

Multiplicity correlations among the heavily ionizing particles produced in nucleus-nucleus collisions have been widely studied which help to investigate the mechanism of particle production. In order to examine the behaviour of multiplicity correlations of secondary particles produced in nucleus- nucleus collisions, we have studied the correlations in the interactions of ³²S-Em at 200 AGeV. Generally, the experimental results have been analyzed by using linear fits of the type:

where N_i, N_j = N_g, N_b and N_h with i ≠ j. The values of inclination coefficients, a_ij and intercepts, b_ij are given in Table 4. The behaviour of multiplicity correlations of secondary particles produced in ³²S-Em interactions at 200 AGeV is shown in Figures 6(a)-(c) along with their linear fits using least square method. From these plots following conclusions may be drawn.

i) A clear saturation in the values of is observed in vs N_g plot for values of N_g beyond ~10. A similar result was obtained by Otterlund [18] for proton-nucleus interactions over a wide range of energy. This means correlation does not depend upon mass of the projectiles and the contribution of the recoiling nucleons towards the excitation energy of the residual nucleus is approximately the same for P- nucleus and nucleus-nucleus interactions.

ii) Variation of with N_b is similar to that of with N_g. However, the saturations are not statistically significant. Also displays a linear dependence on N_h in the whole range of N_h.

iii) The values of increases with the increase of N_b and N_g in the whole range of N_b and N_g in nucleus? nucleus collisions.

In order to see the dependence of the average multiplicities on the mass numbers of the target nuclei, the following relation has been used as:

where “j” stands for black, grey and heavy particles respectively.

The variation of average multiplicities on masses of the target nuclei is shown in Figure 7. The coefficients K and α, are determined from the least-square fit using the experimental data of the present study. The values of these coefficients are given in Table 5. The results depict that the multiplicities of black particles are nearly proportional to linear dimensions of target nuclei. The multiplicities of grey particles are characterized by extremely weak target size dependence.

Figure 8(a) and Figure 8(b) shows the normalized angular distributions for black and grey particles emitted in ³²S-Em collisions at 200 AGeV along with the results obtained by other workers [18] [21] for different projectiles and energies. It has been observed from the figure that the distributions are independent of projectile mass as well as the energy of the projectile. The production mechanism of the black particles, based on nuclear evaporation processes, are isotropically distributed in the rest frame of target nucleus and the distribution is influenced by a strong electromagnetic field of relativistic projectiles, which is almost symmetric about q ≈ 90˚ in the laboratory system. Whereas the production mechanism of grey particles is based on inter-nuclear cascade processes and these particles are mainly emitted in forward direction. The forward and backward hemispheres are defined as the regions where emission angles are less than 90˚ (q < 90˚) and greater than 90˚ (q > 90˚) respectively. The forward (F) to backward (B) ratio for these distributions are calculated and presented in Table 6. It is clearly seen from the table that the F/B ratio shows the same behavior and its value for grey particles is more than black particles. Also the F/B ratio shows a weak dependence on the mass of projectile. The differences in the angular distributions of black and grey particles point to the fact that these particles originate in two different processes, the subsequent evaporation and the initial interaction of target nucleus.

Recently [24] - [26] the investigation of nuclear fragments produced in high energy nucleus-nucleus collisions shows that the multiplicity distributions of target fragments; black and grey particles can be described by Koba- Nielsen-Olesen scaling [27] . Koba, Nielsen and Olesen have predicted that the multiplicity distributions of the produced particles in high-energy hadron-hadron collisions should obey a simple scaling law known as KNO scaling when expressed in terms of the scaling variable Z. If P_n(s) represents the probability for the production of n charged particles in an inelastic hadron-hadron collision at a centre of mass energy Ös, then the multiplicity distributions in high energy collision obey a scaling law:

where s_n(s) is the partial cross-section for the production of n charged particles, s_inel is the total inelastic cross-section and is the average number of charged particles produced. The KNO scaling thus implies that the multiplicity distribution is universal and ψ(z) is an energy independent function at sufficiently high energies when expressed in terms of scaling variable Z.

It has been found by various workers that the empirical expression for ψ(z) in hadron-hadron and hadron- nucleus interactions obeys the semi-inclusive KNO scaling starting from few GeV. It is desirable to make similar studies in nucleus-nucleus collisions as it is expected that nucleus-nucleus collisions (A-A) at these energies can be visualized as superposition of nucleon-nucleon collisions. Several workers [24] [28] - [30] have reported that the validity of KNO scaling holds for the projectile helium particle and black tracks in heavy-ion interactions at Dubna, Bevatron, CERN and AGS energies. It has been shown [26] [31] that the multiplicity distributions of produced black and grey fragments obtained from the events of different projectiles over a wide range of energies in nucleus-nucleus collisions can be described by a KNO scaling law. These distributions can be represented by a universal function of the following form:

where A and B are constants.

In the present work an attempt has been made to study the KNO scaling for the multiplicity distribution of slow and fast target associated particles produced in ³²S-emulsion collisions at 200 AGeV. A plot of ψ(z) as a function of the scaling variable Z for these medium energy target associated protons is shown in Figure 9(a) and Figure 9(b). The experimental points for ¹²C and ²⁸Si at 4.5 AGeV [8] , ²⁸Si at 14.6A GeV [8] and ¹⁶O at 3.7 AGeV and 60 AGeV [9] respectively are also shown in the same figure. The solid curve in the figure is well represented by Equation (6).

It is easily noticed from the figures that the multiplicity distributions of slow and fast target associated protons in nucleus-nucleus collisions at different energies are well described by Equation (6) for the different projectiles and seem to satisfy the scaling function. The best values of A and B used in Equation (6) are found to be 6.21 ± 0.85, 2.65 ± 0.10 and 5.26 ± 0.90, 2.89 ± 0.13 respectively for slow and fast target associated protons. The values of corresponding c²/DOF are found to be 0.63 and 0.79 ± 0.04 respectively for slow and fast target associated protons which indicates that the fitting is good for different projectiles at different energies in case of slow protons but a small deviation from exact scaling can be seen for fast target associated protons in Figure 9(b). It is difficult to give any physical explanation of the multiplicity scaling for slow and fast protons and hence can be regarded as an empirical observation.

The studies of multiplicity distribution in high energy heavy ion collisions have revealed some striking phenomenon which can be described by various models [27] [32] - [34] . Among them, the most spectacular and well known is the negative binomial distribution (NBD) model [34] . The NB law has been successful in describing the multiplicity results in recent high energy experiments. In Figures 10(a)-(c) the multiplicity distributions of black, grey and heavy particles in the interactions of ³²S-Em at 200 AGeV are shown with the corresponding best NB fits. The values of the free parameters of negative binomial distributions “k” and “n” and the c²/DOF obtained using the CERN MINUIT curve fitting program are given in Table 7. From the table it can be concluded that the multiplicity distributions of secondary particles produced in ³²S-Em at 200 AGeV interactions can be fitted quite satisfactorily by the NBD having the form

where n and represent respectively the multiplicity and mean multiplicity, the value of k is determined from: