Proposition of numerical modelling of BEC

We propose an extension of the numerical method to model effect of Bose-Einstein correlations (BEC) observed in hadronization processes which allows for calculations not only the correlation functions C2(Qinv) (onedimensional), but also the corresponding C2(Qx,y,z) (i.e., three-dimensional). The method is based on the bunching of identical bosonic particles in elementary emitting cells (EEC) in phase space in a manner leading to proper Bose- Einstein form of distribution of energy (this was enough to calculate C2(Qinv)). To obtain C2(Qx,y,z) as well as one has to add to it the symmetrization of the multiparticle wave function to properly correlate space-time locations of produced particles with their energy-momentum characteristics.