colour tube forms between the ¯ qq pair, and fragmentation of the strings in the tubes causes the formation of hadrons (see [154] for more details). The Lund theory is part of a more general class of parton fragmentation models used by event generators such as PYTHIA and JETSET. However these are complex multi-parameter Monte Carlo event generators that do not reflect the simplicity of QCD, and have little predictive power. A recent model which has had striking success is that of [3]. It assumes that each quark/antiquark from the ¯ qq pair forms a hadronic gas in thermo- dynamical equilibrium. The gas then cools and expands, until it is no longer in equilibrium, at which point the hadrons are frozen out, and these are then the particles that travel to the detector. So this model only requires three parameters: the temperature and volume of the gas at the point of freeze- out, and a strangeness suppression factor to account for the superior mass of the strange quark. The strangeness supression factor should, in principle, be derivable from the known mass of the strange quark. However, such an attempt is not made in [3], and γ$_{s}$ is left as a free-fit parameter. (Factors for the charmed and bottom quarks can also be introduced, but these quarks are so massive, their production rate is negligible). This model can repro- duce the results of LEP and PEP-PETRA to 34% and 38% respectively. All models of hadronization have two parts: predicting the initial yield of hadrons directly after annihilation, and then allowing for decays of those particles in transit towards the detector. Both of these processes are QCD processes, and so should be fully predictable from the theory. However just predicting one of these two processes is a major challenge: here we concentrate our efforts on the first stage. We will then use experimentally measured values for the branching ratios of all the hadrons [24] to predict 110