The clustering of events in 7-dimensional phase space is studied in the reactions z f p + pnf n* n-at 16 GeV/c. The Yang variables are used for locating eventa in the phase space. Clusters are defined and events are allocated by a novel iterative technique. This is baaed on a measure of the distance of each event from the cluster centre, weighted by the covariance matrix of the event density. The clusters found are then related to reaction mechanisms such as diffraotion dissociation of pion and proton, resonance production, etc. High mass enhancements in the (3z) system and in the ( p n f n -) system are identified and separated from other mechanisms. Cross sections are given for individual sub-channels and compared to results obtained by other techniques.