Multiplicity distributions in interaction models used for cascade simulations

In cascade simulations, conservation of energy is of paramount importance. It is pointed out that it is possible to construct algorithms in which energies of secondaries are chosen from physically correct distributions in such a way as to conserve energy, and give physically correct multiplicity distributions without further adjustment. The resulting non-Poissonian nature of the distribution and the approximate KNO scaling are related to the large fluctuations in the energy available for pion production and to the detailed shape of the inclusive momentum distribution. Two multiparticle production algorithms are described which have these features. One is the algorithm used in the simulation of C jets, the other is for a hadron-nucleus interaction model used in the calculation of muon and neutrino fluxes.