Within the framewok of the VENUS model we propose a new approach to generate hadron-nucleus and nucleus-nucleus events at ultrarelativisitic energies according to a many body probability distribution using dynamical Monte-Carlo methods. For the primary stage of a collision of two nuclei we make the following ansatz: the probability of two colliding nuclei to form a configuration K=(IP_1,-,IP_n, h_1~p,-h_A~p, h_1~t,-h_b~t) of Pomerons IP_k, projectile remnants h_i~p and target remnants h_j~t is given as prob(A-B->K) propor.to OMEGA(K), with a OMEGA(K) being a probability distribution based on Gribov-glauber theory, which we will define here. We are going to employ Monte Carlo techniques to generate randomly configurations K distributed as OMEGA(K) Because of the high dimensional configuration space (central Pb-Pb: 500 collisions, 700 Pomerons, 200 nucleon remnants), with complicated boundary conditions ,due to energy-momentum conservation, direct Monte Carlo methods fail. Instead we construct a Markov process~5 by generating a sequence K_t: K_0, K_1-K_Ieq, with I_eq being sufficiently large to have reached equilibrium. Repeating this procedure many times we obtain configurations K_Ieq~(1), K_Ieq(2),-which are distributed according to OMEGA(K). These are then treated further in the simulation (fragmentation, secondary interactions).
展开▼