Simulation of heavy ion collisions at {radical}s = 20--200 GeV

A new hadronic cascade code (LUCIFER) is introduced, for simulation of relativistic heavy ion collisions at CERN energies and up to RHIC. It is based on a simple, experimentally and theoretically motivated picture of hh interactions. Final state hadrons are produced by decay of intermediate state clusters, or lumps of excited hadronic matter. These are similar to resonances, but have a continuous mass distribution. Clusters are the objects that re-interact in the cascade. Single diffractive dissociation is used to fix the cluster properties. The model has just two parameters: {tau}{sub d}, the decay time of the clusters, and {tau}{sub f} the formation time of the clusters. Comparison is made with recent CERN data in the Pb + Pb system. The first consistent cascade simulation, of J/{psi} production/suppression is presented. It appears likely that a purely hadronic interpretation can be given to recent CERN data on apparently anomalous J/{psi} suppression in Pb + Pb