Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: event shape variables are perhaps the most popular observables for testing qcd and for improving our understanding of its dynamics . event shape studies began in earnest towards the late seventies as a simple quantitative method to understand the nature of gluon bremsstrahlung @xcite . for instance , it was on the basis of comparisons to event shape data that one could first deduce that gluons were vector particles , since theoretical predictions employing scalar gluons did not agree with experiment  @xcite . quite generally , event shapes parametrise the geometrical properties of the energy - momentum flow of an event and their values are therefore directly related to the appearance of an event in the detector . in other words the value of a given event shape encodes in a continuous fashion , for example , the transition from pencil - like two - jet events with hadron flow prominently distributed along some axis , to planar three - jet events or events with a spherical distribution of hadron momenta . thus they provide more detailed information on the final state geometry than say a jet finding algorithm which would always classify an event as having a certain finite number of jets , even if the actual energy flow is uniformly distributed in the detector and there is no prominent jet structure present in the first place . event shapes are well suited to testing qcd mainly because , by construction , they are collinear and infrared safe observables . this means that one can safely compute them in perturbation theory and use the predictions as a means of extracting the strong coupling @xmath1 . they are also well suited for determinations of other parameters of the theory such as constraining the quark and gluon colour factors as well as the qcd beta function . additionally they have also been used in other studies for characterising the final state , such as investigations of jet and heavy quark multiplicities as functions of event shape variables  @xcite .    aside from testing the basic properties of qcd , event shape distributions are a powerful probe of our more detailed knowledge of qcd dynamics . all the commonly studied event shape variables have the property that in the region where the event shape value is small , one is sensitive primarily to gluon emission that is soft compared to the hard scale of the event and/or collinear to one of the hard partons . such small - transverse - momentum emissions have relatively large emission probabilities ( compared to their high transverse momentum counterparts ) due to logarithmic soft - collinear dynamical enhancements as well as the larger value of their coupling to the hard partons . predictions for event shape distributions therefore typically contain large logarithms in the region where the event - shape is small , which are a reflection of the importance of multiple soft and/or collinear emission . a successful prediction for an event shape in this region requires all - order resummed perturbative predictions or the use of a monte - carlo event generator , which contains correctly the appropriate dynamics governing multiple particle production . comparisons of these predictions to experimental data are therefore a stringent test of the understanding qcd dynamics that has been reached so far . one other feature of event shapes , that at first appeared as an obstacle to their use in extracting the fundamental parameters of qcd , is the presence of significant non - perturbative effects in the form of power corrections that vary as an inverse power of the hard scale , @xmath2 . for most event shapes , phenomenologically it was found that @xmath3 ( as discussed in  @xcite ) , and the resulting non - perturbative effects can be of comparable size to next - to - leading order perturbative predictions , as we shall discuss presently . the problem however , can be handled to a large extent by hadronisation models embedded in monte - carlo event generators @xcite , which model the conversion of partons into hadrons at the cost of introducing several parameters that need to be tuned to the data . once this is done , application of such hadronisation models leads to very successful comparisons of perturbative predictions with experimental data with , for example , values of @xmath1 consistent to those obtained from other methods and with relatively small errors . however since the mid - nineties attempts have also been made to obtain a better insight into the physics of hadronisation and power corrections in particular . theoretical models such as those based on renormalons have been developed to probe the non - perturbative domain that gives rise to power corrections ( for a review , see @xcite ) . since the power corrections are relatively large effects for event shape variables , scaling typically as @xmath4 ( no other class of observable shows such large effects ) , event shapes have become the most widely used means of investigating the validity of these ideas . in fact an entire phenomenology of event shapes has developed which is based on accounting for non - perturbative effects via such theoretically inspired models , and including them in fits to event shape data alongside the extraction of other standard qcd parameters . in this way event shapes also serve as a tool for understanding more quantitatively the role of confinement effects . the layout of this article is as follows . in the following section we list the definitions of several commonly studied event shapes in @xmath5 annihilation and deep inelastic scattering ( dis ) and discuss briefly some of their properties . then in section  [ sec : perturbative ] we review the state of the art for perturbative predictions of event - shape mean values and distributions . in particular , we examine the need for resummed predictions for distributions and clarify the nomenclature and notation used in that context , as well as the problem of matching these predictions to fixed order computations . we then turn , in sections  [ sec : exp ] and [ sec : distributions ] , to comparisons with experimental data and the issue of non - perturbative corrections , required in order to be able to apply parton level calculations to hadronic final state data . we discuss the various approaches used to estimate non - perturbative corrections , ranging from phenomenological hadronisation models to analytical approaches based on renormalons , and shape - functions . we also discuss methods where the standard perturbative results are modified by use of renormalisation group improvements or the use of ` dressed gluons ' . we study fits to event shape data , from different sources , in both @xmath5 annihilation and dis and compare the various approaches that are adopted to make theoretical predictions . in addition we display some of the results obtained by fitting data for various event shape variables for parameters such as the strong coupling , the qcd colour factors and the qcd beta function . lastly , in section  [ sec : outlook ] , we present an outlook on possible future developments in the field . we list here the most widely studied event shapes , concentrating on those that have received significant experimental and theoretical attention .      _ the _ canonical event shape is the thrust  @xcite , @xmath6 : @xmath7 where the numerator is maximised over directions of the unit vector @xmath8 and the sum is over all final - state hadron momenta @xmath9 ( whose three - vectors are @xmath10 and energies @xmath11 ) . the resulting @xmath8 is known as the thrust axis . in the limit of two narrow back - to - back jets @xmath12 , while its minimum value of @xmath13 corresponds to events with a uniform distribution of momentum flow in all directions . the infrared and collinear safety of the thrust ( and other event shapes ) is an essential consequence of its linearity in momenta . often it is @xmath14 ( also called @xmath15 ) that is referred to insofar as it is this that vanishes in the @xmath16-jet limit . a number of other commonly studied event shapes are constructed employing the thrust axis . amongst these are the invariant squared jet - mass  @xcite and the jet - broadening variables  @xcite , defined respectively as @xmath17 where the plane perpendicular to the thrust axis is used to separate the event into left and right hemisphere , @xmath18 and @xmath19 . given these definitions one can study the heavy - jet mass , @xmath20 and wide - jet broadening , @xmath21 ; analogously one defines the light - jet mass ( @xmath22 ) and narrow - jet broadening @xmath23 ; finally one defines also the sum of jet masses @xmath24 and the total jet broadening @xmath25 and their differences @xmath26 , and @xmath27 . like @xmath14 , all these variables ( and those that follow ) vanish in the two - jet limit ; @xmath22 and @xmath23 are special in that they also vanish in the limit of three narrow jets ( three final - state partons ; or in general for events with any number of jets in the heavy hemisphere , but only one in the light hemisphere ) . another set of observables @xcite making use of the thrust axis starts with the thrust major @xmath28 , @xmath29 where the maximisation is performed over all directions of the unit vector @xmath30 , such that @xmath31 . the thrust minor , @xmath32 , is given by @xmath33 and is sometimes also known @xcite as the ( normalised ) out - of - plane momentum @xmath34 . like @xmath22 and @xmath23 , it vanishes in the three - jet limit ( and , in general , for planar events ) . finally from @xmath28 and @xmath32 one constructs the oblateness , @xmath35 . an alternative axis is used for the spherocity @xmath36  @xcite , @xmath37 while there exist reliable ( though algorithmically slow ) methods of determining the thrust and thrust minor axes , the general properties of the spherocity axis are less well understood and consequently the spherocity has received less theoretical attention . an observable similar to the thrust minor ( in that it measures out - of - plane momentum ) , but defined in a manner analogous to the spherocity is the acoplanarity  @xcite , which minimises a ( squared ) projection perpendicular to a plane . it is also possible to define event shapes without reference to an explicit axis . the best known examples are the @xmath38 and @xmath39 parameters @xcite which are obtained from the momentum tensor @xcite . because of the quadratic dependence on particle momenta , these observable are collinear unsafe and so no longer widely studied . ] @xmath40 where @xmath41 is the @xmath42 component of the three vector @xmath43 . in terms of the eigenvalues @xmath44 , @xmath45 and @xmath46 of @xmath47 , the @xmath38 and @xmath39 parameters are given by @xmath48 the @xmath38-parameter is related also to one of a series of fox - wolfram observables  @xcite , @xmath49 , and is sometimes equivalently written as @xmath50 the @xmath39-parameter , like the thrust minor , vanishes for all final - states with up to 3 particles , and in general , for planar events .    the @xmath38-parameter can actually be considered ( in the limit of all particles being massless ) as the integral of a more differential observable , the energy - energy correlation ( eec ) @xcite , @xmath51 where the average in eq . ( [ eq : sigmachi ] ) is carried out over all events . another set of variables that characterise the shape of the final state are the @xmath52-jet resolution parameters that are generated by jet finding algorithms . examples of these are the jade  @xcite and durham  @xcite jet clustering algorithms . one introduces distance measures @xmath53 @xmath54 for each pair of particles @xmath55 and @xmath56 . the pair with the smallest @xmath53 is clustered ( by adding the four - momenta  the @xmath57 recombination scheme ) and replaced with a single pseudo - particle ; the @xmath53 are recalculated and the combination procedure repeated until all remaining @xmath53 are larger than some value @xmath58 . the event - shapes based on these jet algorithms are @xmath59 , defined as the maximum value of @xmath58 for which the event is clustered to 3 jets ; analogously one can define @xmath60 , and so forth . other clustering jet algorithms exist . most differ essentially in the definition of the distance measure @xmath53 and the recombination procedure . for example , there are e0 , p and p0 variants  @xcite of the jade algorithm , which differ in the details of the treatment of the difference between energy and the modulus of the 3-momentum . the geneva algorithm  @xcite is like the jade algorithm except that in the definition of the @xmath53 it is @xmath61 that appears in the denominator instead of the total energy . an algorithm that has been developed and adopted recently is the cambridge algorithm  @xcite , which uses the same distance measure as the durham algorithm , but with a different clustering sequence . we note that a number of variants of the above observables have recently been introduced  @xcite , which differ in their treatment of massive particles . these include the @xmath62-scheme where all occurrences of @xmath11 are replaced by @xmath63 and the @xmath57-scheme where each 3-momentum @xmath43 is rescaled by @xmath64 . the former leads to a difference between the total energy in the initial and final states , while the latter leads to a final - state with potentially non - zero overall @xmath65-momentum . while such schemes do have these small drawbacks , for certain observables , notably the jet masses , which are quite sensitive to the masses of the hadrons ( seldom identified in experimental event - shape studies ) , they tend to be both theoretically and experimentally cleaner . ref . @xcite also introduced a ` decay ' scheme ( an alternative is given in  @xcite ) where all hadrons are artificially decayed to massless particles . since a decay is by definition a stochastic process , this does not give a unique result on an event - by - event basis , but should rather be understood as providing a correction factor which is to be averaged over a large ensemble of events . it is to be kept in mind that these different schemes all lead to identical perturbative predictions ( with massless quarks ) and differ only at the non - perturbative level . a final question relating to event - shape definitions concerns the hadron level at which measurements are made . since shorter lived hadrons decay during the time of flight , one has to specify whether measurements were made at a stage before or after a given species of hadron decays . it is important therefore when experimental results are quoted , that they should specify which particles have been taken to be stable and which have not .      as well as the @xmath66 variables discussed above , it is also possible to define , by analogy , event shapes in dis . the frame in which dis event shapes can be made to most closely resemble those of @xmath66 is the breit frame @xcite . this is the frame in which @xmath67 , where @xmath68 is the incoming proton momentum and @xmath69 the virtual photon momentum . one defines two hemispheres , separated by the plane normal to the photon direction : the remnant hemisphere ( @xmath70 , containing the proton remnant ) , and the current hemisphere ( @xmath71 ) . at the level of the quark - parton model , @xmath71 is like one hemisphere of @xmath0 and it is therefore natural to define event shapes using only the momenta in this hemisphere . in contrast , any observable involving momenta in the remnant hemisphere must take care to limit its sensitivity to the proton remnant , whose fragmentation can not be reliably handled within perturbation theory . a possible alternative to studying just the current hemisphere , is to take all particles except those in a small cone around the proton direction  @xcite .    a feature that arises in dis is that there are two natural choices of axis . for example , for the thrust @xmath72 one can either choose the unit vector @xmath73 to be the photon ( @xmath74 ) axis , @xmath75 , or one can choose it to be the true thrust axis , that which maximises the sum , giving @xmath76 . similarly one defines two variants of the jet broadening , @xmath77 for the jet - mass and @xmath38-parameter the choice of axis does not enter into the definitions and we have @xmath78 all the above observables can also be defined with an alternative normalisation , @xmath79 replacing @xmath80 , in which case they are named @xmath81 and so on . we note that @xmath81 was originally proposed in  @xcite . for a reader used to @xmath0 event shapes the two normalisations might at first sight seem equivalent  however when considering a single hemisphere , as in dis , the equivalence is lost , and indeed there are even events in which the current hemisphere is empty . this is a problem for observables normalised to the sum of energies in @xmath71 , to the extent that to ensure infrared safety it is necessary to exclude all events in which the energy present in the current hemisphere is smaller than some not too small fraction of @xmath82 ( see for example @xcite ) . additionally there are studies of variables that vanish in the 2 + 1-parton limit for dis , in particular the @xmath34 defined in analogy with the thrust minor of @xmath0 , with the thrust axis replaced by the photon axis @xcite , or an azimuthal correlation observable  @xcite . rather than being examined just in @xmath71 , these observables use particles also in the remnant hemisphere ( for @xmath34 , all except those in a small cone around the proton ) .    as in @xmath0 , jet rates are studied also in dis . unlike the event - shapes described above , the jet - shapes make use of the momenta in both hemispheres of the breit frame . their definitions are quite similar to those of @xmath0 except that a clustering to the proton remnant ( beam jet ) is also included  @xcite . though only @xmath0 and dis are within the scope of this review , we take the opportunity here to note that related observables are being considered also for other processes . notably for drell - yan production an out - of - plane momentum measurement has been proposed in  @xcite and various thrust and thrust - minor type observables have been considered in hadron - hadron dijet production in  @xcite . the observables discussed above are all infrared and collinear safe  they do not change their value when an extra soft gluon is added or if a parton is split into two collinear partons . as emerges from the original discussion of sterman and weinberg  @xcite , this is a necessary condition for the cancellation of real and virtual divergences associated with such emissions , and therefore for making finite perturbative predictions .    for an event shape that vanishes in the @xmath52-jet limit ( which we shall generically refer as an @xmath52-jet observable ) , the leading perturbative contribution is of order @xmath83 . for example the thrust distribution in @xmath0 is given by @xmath84 + { { \cal o}\left({\alpha_s}^2\right ) } \nonumber\end{aligned}\ ] ] ( see e.g.  @xcite ) . when calculating perturbative predictions for mean values of event shapes , as well as higher moments of their distributions , one has integrals of the form @xmath85 we note that in general , fixed - order event - shape distributions diverge in the limit as @xmath86 goes to @xmath87 , cf . ( [ eq : thrustdist ] ) for @xmath88 . in eq . ( [ eq : meanv ] ) , the weighting with a power of @xmath86 is sufficient to render the singularity integrable , and the integral is dominated by large @xmath86 ( and so large transverse momenta of order @xmath82 , the hard scale ) . this dominance of a single scale ensures that the coefficients in the perturbative expansion for @xmath89 are well - behaved ( for example they are free of any enhancements associated with logarithms of ratios of disparate scales ) .    beyond leading order , perturbative calculations involve complex cancellations between soft and collinear real and virtual contributions and are nearly always left to general - purpose `` fixed - order monte carlo '' programs . the current state of the art is next - to - leading order ( nlo ) , and available programs include : for @xmath0 to 3 jets , event  @xcite , eerad  @xcite and event2 @xcite ; for @xmath0 to 4 jets , menlo parc  @xcite , mercutio @xcite and eerad2  @xcite ; for dis to @xmath90 jets , mepjet  @xcite , disent  @xcite and disaster++ @xcite ; and for all the above processes and additionally dis to @xmath91 jets and various hadron - hadron and photo - production processes , nlojet++  @xcite ( for photo - production , there is additionally jetvip  @xcite ) . in recent years much progress has been made towards nnlo calculations , though complete results remain to be obtained ( for a review of the current situation , see @xcite ) . for the event shape distributions themselves , however , fixed order perturbative estimates are only of use away from the @xmath92 region , due to the singular behaviour of the fixed order coefficients in the @xmath93 limit . in this limit , at higher orders each power of @xmath94 is accompanied by a coefficient which grows as @xmath95 . these problems arise because when @xmath86 is small one places a restriction on real emissions without any corresponding restriction on virtual contributions  the resulting incompleteness of cancellations between logarithmically divergent real and virtual contributions is the origin of the order - by - order logarithmic enhancement of the perturbative contributions . to obtain a meaningful answer it is therefore necessary to perform an all - orders _ resummation _ of logarithmically enhanced terms .      when discussing resummations it is convenient to refer to the integrated distribution , @xmath96 quite generally @xmath97 has a perturbative expansion of the form @xmath98 one convention is to refer to all terms @xmath99 as leading logarithms ( ll ) , terms @xmath100 as next - to - leading logarithms ( nll ) , etc . , and within this hierarchy a resummation may account for all ll terms , or all ll and nll terms and so forth . such a resummation gives a convergent answer up to values of @xmath101 , beyond which terms that are formally subleading can become as important as the leading terms ( for example if @xmath102 , then the nnll term @xmath103 is of the same order as the ll term @xmath104 ) . at its limit of validity , @xmath105 , a n@xmath106ll resummation neglects terms of relative accuracy @xmath107 . an important point though is that for nearly all observables that have been resummed there exists a property of exponentiation : @xmath108 in some cases @xmath109 is written in terms of a sum of such exponentiated contributions ; in certain other cases the exponentiation holds only for a suitable integral ( e.g.  fourier ) transform of the observable . the fundamental point of exponentiation is that the inner sum , over @xmath110 , now runs only up to @xmath111 , instead of @xmath112 as was the case in eq . ( [ eq : dumbsum ] ) . with the exponentiated form for the resummation , the nomenclature `` ( next - to)-leading - logarithmic '' acquires a different meaning  n@xmath106ll now refers to all terms in the exponent @xmath113 . to distinguish the two classification schemes we refer to them as n@xmath106ll@xmath114 and n@xmath106ll@xmath115 . the crucial difference between n@xmath106ll@xmath114 and n@xmath106ll@xmath115 resummations lies in the range of validity and their accuracy . a n@xmath106ll@xmath115 resummation remains convergent considerably further in @xmath116 , up to @xmath102 ( corresponding usually to the peak of the distribution @xmath117 ) ; at this limit , neglected terms are of relative order @xmath118 . consequently a n@xmath106ll@xmath115 resummation includes not only all n@xmath106ll@xmath114 terms , but also considerably more . from now when we use the term n@xmath106ll it should be taken to mean n@xmath106ll@xmath119 . two main ingredients are involved in the resummations as shown in eq . ( [ eq : expresum ] ) . firstly one needs to find a way of writing the observable as a factorised expression of terms for individual emissions . this is often achieved with the aid of one or more integral ( mellin , fourier ) transforms . secondly one approximates the multi - emission matrix element as a product of individual independent emission factors . it is to be kept in mind that not all observables exponentiate . the jade jet resolution is the best known such example  @xcite ( and it has not yet been resummed ) . the general criteria for exponentiation are discussed in  @xcite . resummed results to nll accuracy exist for a number of @xmath16-jet observables in @xmath120 jets : the thrust  @xcite , the heavy - jet mass  @xcite and the single - jet and light - jet masses @xcite , jet broadenings @xcite , @xmath38-parameter  @xcite , durham and cambridge jet resolutions @xcite , the thrust major and oblateness @xcite , eec  @xcite and event - shape / energy - flow correlations  @xcite ; for @xmath65-jet observables in @xmath121 jet events we have the thrust minor  @xcite and the @xmath39-parameter @xcite . for dis to @xmath122 jets the current - hemisphere observables are resummed in refs . @xcite and the jet rates in  @xcite ( though only to nll@xmath114 accuracy ) while the @xmath90 jet observables , @xmath34 and azimuthal correlations , have been resummed in  @xcite . methods for automated resummation of arbitrary observables are currently in development for a range of processes@xcite , and techniques have also being developed for dealing with arbitrary processes  @xcite . we also note that for the @xmath0 thrust and heavy - jet mass there have been investigations of certain classes of corrections beyond nll accuracy  @xcite . the above resummations all apply to events with only light quarks . investigations to nll@xmath123 accuracy have also been performed for @xmath0 jet rates with heavy quarks  @xcite . the above resummations are all for @xmath52-jet observables in the @xmath52-jet limit . for some @xmath65-jet observables it is useful to have resummations in the @xmath16-jet limit ( because it is here that most of the data lie ) , though currently little attention has been devoted to such resummations , with just nll@xmath123 results for jet - resolution parameters @xcite and nll@xmath115 results for the light - jet mass and narrow - jet broadening  @xcite . furthermore in addition to the resummation of the @xmath38-parameter for @xmath124 there exists a ll resummation  @xcite ( that could straightforward be extended to nll ) of a ` shoulder ' structure at @xmath125 ( related to a step function at lo ) , which corresponds to symmetric three jet events . some of the observables ( for example many single - jet observables ) have the property that they are sensitive to emissions only in some of the phase space . these are referred to as non - global observables and nll resummed predictions for them require that one account for coherent ensembles of energy - ordered large - angle gluons . this has so far been done only in the large-@xmath126 limit  @xcite ( for reviews see @xcite ) , though progress is being made in extending this to finite @xmath126  @xcite .     at @xmath127 . the resummed distribution has been matched ( @xmath128 scheme , see below ) so as to include the full lo and nlo contributions . the inset differs from the main plot only in the scales.,scaledwidth=60.0% ]    to illustrate the impact of resummations we show in figure  [ fig : thrlonlo ] the lo , nlo and resummed results for the thrust distribution . not only is the divergence in the lo distribution clearly visible , but one also observes a marked difference in the qualitative behaviour of the lo and nlo results , both of which show large differences relative to the nll resummed prediction at small @xmath14 .    for the practical use of resummed results an important step is that of matching with the fixed - order prediction . at the simplest level one may think of matching as simply adding the fixed - order and resummed predictions while subtracting out the doubly counted ( logarithmic ) terms . such a procedure is however too naive in that the fixed - order contribution generally contains terms that are subleading with respect to the terms included in the resummation , but which still diverge . for example when matching nll and nlo calculations , one finds that the nlo result has a term @xmath129 , which in the distribution diverges as @xmath130 . this is unphysical and the matching procedure must be sufficiently sophisticated so as to avoid this problem , and in fact it should nearly always ensure that the matched distribution goes to zero for @xmath131 . several procedures exist , such as @xmath128 and @xmath109 matching  @xcite and multiplicative matching @xcite and they differ from one - another only in terms that are nnlo and nnll , i.e.  formally beyond the state - of - the - art accuracy . we note the matching generally ensures that the resummation , in addition to being correct at nll accuracy ( our shorthand for nll@xmath115 ) , is also correct at nnll@xmath123 accuracy . matching also involves certain other subtleties , for example the ` modification of the logarithm ' , whereby @xmath132 is replaced with @xmath133 , where @xmath134 is the largest kinematically allowed value for the observable . this ensures that the logarithms go to zero at @xmath134 , rather than at some arbitrary point ( usually @xmath135 ) , which is necessary in order for the resulting matched distribution to vanish at @xmath134 . [ sec : means ] [ sec : meanpc ]    much of the interest in event shapes stems from the wealth of data that is available , covering a large range of centre of mass energies ( @xmath0 ) or photon virtualities ( dis ) . the @xmath0 data comes from the pre - lep experiments  @xcite , the four lep experiments @xcite @xcite @xcite and sld  @xcite . in addition , because many observables have been proposed only in the last ten or so years and owing to the particular interest ( see below ) in data at moderate centre of mass energies , the jade data have been re - analysed @xcite . have also been obtained from lep  1 @xcite by considering events with an isolated hard final - state photon and treating them as if they were pure qcd events whose centre of mass energy is that of the recoiling hadronic system . such a procedure is untrustworthy because it assumes that one can factorise the gluon production from the photon production . this is only the case when there is strong ordering in transverse momenta between the photon and gluon , and is not therefore applicable when both the photon and the gluon are hard . tests with monte carlo event generators which may indicate that any ` non - factorisation ' is small  @xcite are not reliable , because the event generator used almost certainly does not contain the full 1-photon , 1-gluon matrix element . it is also to be noted that the isolation cuts on the photon will bias the distribution of the event shapes , in some cases @xcite similarly to an event - shape energy - flow correlation  @xcite . therefore we would argue that before relying on these data , one should at the very least compare the factorisation approximation with exact lo calculations , which can be straightforwardly obtained using packages such as grace  @xcite , comphep  @xcite or amegic++ @xcite . ] results in dis come from both h1 @xcite ( mean values and distributions ) and zeus  @xcite ( means ) , @xcite ( means and distributions in rapidity gap events ) .     thrust , compared to lo , nlo and nlo+@xmath136 ; right : data for the durham jet resolution @xmath137 compared to the nlo prediction ( figure taken from @xcite).,title="fig:",scaledwidth=48.0% ]   thrust , compared to lo , nlo and nlo+@xmath136 ; right : data for the durham jet resolution @xmath137 compared to the nlo prediction ( figure taken from @xcite).,title="fig:",scaledwidth=48.0% ]    let us start our discussion of the data by examining mean values ( distributions are left until section  [ sec : distributions ] ) . as we have already mentioned , one of the most appealing properties of event shapes , in terms of testing qcd , is the fact that they are calculable in perturbation theory and so provide a direct method for the extraction of @xmath1 as well as testing other parameters of the theory such as the colour factors @xmath138 and @xmath139 by using a wealth of available experimental data . however one obstruction to the clean extraction of these fundamental parameters is the presence in many cases , of significant non - perturbative effects that typically fall as inverse powers @xmath140 of @xmath82 , the hard scale in the reaction ( the centre - of  mass energy in @xmath5 annihilation and the momentum transfer in dis ) . the importance of such _ power corrections _ varies from observable to observable . for example , it can be seen from fig . [ fig : meandata ] ( left ) that the data for the mean value of the thrust variable need a significant component @xmath141 in addition to the lo and nlo fixed order perturbative estimates , in order to be described . on the other hand the comparison ( right ) for the mean value of the durham jet resolution parameter @xmath137 with the nlo prediction alone is satisfactory , without the need for any substantial power correction term . the problem of non - perturbative corrections is so fundamental in event - shape studies that it is worthwhile giving a brief overview of the main approaches that are used . no statement can be made with standard perturbative methods about the size of these power corrections and hence the earliest methods adopted to quantify them were phenomenological hadronisation models embedded in monte carlo event generators  @xcite . one of the main issues involved in using such hadronisation models is the existence of several adjustable parameters . this essentially means that a satisfactory description of data can be obtained by tuning the parameters , which does not allow much insight into the physical origin of such power behaved terms which are themselves of intrinsic theoretical interest ( for a more detailed critique of hadronisation models see  @xcite ) . however since the mid - nineties there have been developments that have had a significant impact on the theoretical understanding of power corrections . perhaps the most popular method of estimating power corrections is based on the renormalon model . in this approach one examines high - order terms @xmath142 of the perturbative series that are enhanced by coefficients @xmath143 , where @xmath144 is the first coefficient of the @xmath145 function . the factorial divergence leads to an ambiguity in the sum of the series of order @xmath2 , where the value of @xmath62 depends on the speed with which the high - order terms diverge . this approach has seen far more applications than can possibly be described here and for a full discussion and further references on renormalons the reader is referred to @xcite . the most extensive phenomenological applications of renormalon based ideas have been in the estimation of power corrections for event - shape variables . the reason is that event shapes have much larger ( and so experimentally more visible ) power corrections ( typically @xmath146 ) than most other observables ( typically @xmath147 or smaller ) .    for event shape variables in @xmath5 annihilation in the two - jet limit , renormalon - inspired ( or related ) studies of power corrections have been carried out in refs . for the @xmath122 jet limit of dis , results on power corrections based on the renormalon approach can be found in @xcite . studies in the @xmath65-jet limit have been given in  @xcite . the approach that has been most commonly used in comparing theoretical predictions with experimental data is that initiated by dokshitzer and webber  @xcite . here the full result for the mean value , @xmath148 , is given by @xmath149 where @xmath150 is the perturbative prediction for the mean value . the hadronisation correction is included through the term @xmath151 , where @xmath62 is related to the speed of divergence of the renormalon series ( as discussed above ) , @xmath152 is an observable - dependent ( calculable ) coefficient and @xmath153 is a non - perturbative factor , scaling as @xmath154 which is hypothesised to be common across a whole class of observables with the same value of @xmath62 ( strictly speaking the full story is a little more complicated , see e.g.  @xcite ) . most event shapes have @xmath3 , implying a leading power correction scaling as @xmath146 , essentially a consequence of the observables linearity on soft momenta . for these observables with @xmath3 the form for @xmath155 that has become standard is  @xcite , @xmath156 where @xmath157 , @xmath158 and @xmath159 . in eq . ( [ eq : cpfin ] ) an arbitrary _ infrared matching scale _ @xmath160 has been introduced , intended to separate the perturbative and non - perturbative regions . it is usually taken to be @xmath16  gev ( and for systematic error estimates it is varied between @xmath161 and @xmath65  gev ) . the only truly non - perturbative ingredient in eq . ( [ eq : cpfin ] ) is @xmath162 , which can be interpreted as the average value of an infrared finite strong coupling for scales below @xmath160 ( such a concept was first applied phenomenologically in  @xcite ) . though one could imagine estimating it from lattice studies of the coupling ( such as  @xcite ) , one should keep in mind that the coupling in the infrared is not a uniquely defined quantity . in practice the phenomenological test of the renormalon approach will be that @xmath163 has a consistent value across all observables ( in the @xmath3 class ) . the terms with negative sign in the parenthesis in eq . ( [ eq : cpfin ] ) are a consequence of merging standard nlo perturbative results , which include the small spurious contributions from the infrared region ( scales up to @xmath160 ) , with the non - perturbative power correction that accounts correctly for scales from zero to @xmath160 . carrying out the subtraction of the perturbative terms that arise from the infrared region , below @xmath160 , amounts to inclusion of the negative sign terms in parenthesis above . note that this subtraction procedure must be carried out to the level of accuracy of the corresponding perturbative estimate . in other words if a perturbative estimate becomes available at @xmath164 then an additional subtraction term proportional to @xmath165 is required . with the subtraction procedure carried out to @xmath166 , as above , one expects a residual @xmath160 dependence @xmath167 . the scale @xmath168 is the renormalisation scale , which should be taken of the same order as the hard scale @xmath82 . the ` milan factor ' , @xmath169 in eq . ( [ eq : cpfin ] ) , is @xcite @xmath170 where the numerical result is given for @xmath171 , since only light flavours will be active at the relevant ( low ) scales . it accounts for the fact that the usual , ` naive ' , calculation for @xmath152 is carried out on the basis of a decaying virtual gluon ( i.e.  cutting a bubble in the chain of vacuum polarisation insertions that lead to the running of the coupling ) , but without fully accounting for the non - inclusiveness of that decay when dealing with the observable . it was pointed out in  @xcite that this is inconsistent and full ( two - loop ) calculations revealed  @xcite that if the ` naive ' @xmath152 coefficient is calculated in an appropriate scheme , then the factor @xmath169 comes out to be universal . as discussed in @xcite this is essentially a consequence of the fact that regardless of whether one accounts for the virtual gluon decay , the distribution of ` non - perturbative ' transverse momentum is independent of rapidity ( the ` tube - model ' of  @xcite , based essentially on boost invariance ) . while the factor @xmath169 corrects the @xmath152 at two - loop level , the question of even higher order corrections is still open , although one can argue that such effects will be suppressed by a factor @xmath172 relative to the leading power correction . this argument relies on the hope that the strong coupling remains moderate even in the infrared region . .coefficients of @xmath4 power corrections for @xmath5 event shapes ; @xmath173 . [ cols="^,^,^,^,^,^,^",options="header " , ]     tables 1 and 2 display the values of the @xmath152 coefficients obtained in the 2-jet @xmath5 and and @xmath122-jet dis cases respectively . note the different behaviour , proportional to @xmath174 , that arises in the case of the jet broadenings in both processes , a consequence of the fact that the broadening is more sensitive to recoil induced by perturbative radiation  @xcite . the improved theoretical understanding that led to this prediction was strongly stimulated by experimental analysis  @xcite of features of earlier predictions  @xcite that were inconsistent with the data . since the origin of the @xmath174 enhancement is a perturbative , non - perturbative interplay best explained with reference to the power correction for the distribution , we delay its discussion to section  [ sec : anldist ] . another consequence of this interplay of the broadening power correction with perturbative radiation is the fact that the dis broadening @xmath175 has a bjorken @xmath176 dependent term corresponding to dglap evolution of the parton densities , as indicated by the term proportional to @xmath177  @xcite , where @xmath178 denotes a parton density function . the coefficients for the other dis variables are @xmath176 independent since in those cases the power correction arises from soft emission alone . we also mention the result for the power correction to another interesting variable , the energy - energy correlation ( eec ) , in the back - to - back region . like the broadening this variable exhibits the impact of perturbative - non  perturbative interplay . the results are fractional power corrections that vary as @xmath179 ( for the quark - gluon correlation ) and @xmath180 for the hard quark - antiquark correlation  @xcite , instead of @xmath4 and @xmath181 contributions respectively , that would be obtained by considering np emission without the presence of harder perturbative emissions .    observables for which the power corrections have yet to be understood include the jet resolution parameters . renormalon based predictions were given in  @xcite suggesting a @xmath146 correction for @xmath182 and @xmath147 for @xmath183 . however it seems that in both cases there could be significant perturbative non - perturbative interplay which will complicate the picture . nevertheless , as we have seen in fig .  [ fig : meandata ] ( right ) any correction to @xmath184 is certainly small compared to that for other observables . a general point to be kept in mind is that the universality of @xmath163 can be broken by contributions of order @xmath185 , where @xmath186 is the mass scale for hadrons @xcite . whether this happens or not depends on the hadron mass scheme in which the observable is defined ( cf . section  [ sec : definitions ] ) . most observables are implicitly in the @xmath62-scheme , which generally involves a small negative breaking of universality ( which is almost observable - independent , so that an illusion of universality will persist ) . the jet masses are an exception and in their usual definitions they have a significant ( positive ) universality breaking component . the @xmath57-scheme is free of universality breaking terms , as is the decay scheme ( since hadron masses are zero ) . we also note that the power correction associated with hadron mass effects is enhanced by an anomalous dimension , @xmath187 with @xmath188 . the normal ` renormalon ' power correction would also be expected to have an anomalous dimension , but it has yet to be calculated .      as well as studies of mean values ( first moments ) of observables there has also been work on higher - moments , @xmath189 with @xmath190 . simple renormalon - inspired arguments suggest that the @xmath191 moments of event shapes should have as their leading power correction at most a @xmath192 contribution ( essentially since @xmath193 vanishes at @xmath194 for small @xmath195 ) . however as was pointed out by webber  @xcite , the fact ( to be discussed in section  [ sec : anldist ] ) that the @xmath146 power correction essentially corresponds to a shift of the distribution ( for @xmath6 , @xmath38 and @xmath196 , but not for the broadenings ) , means that to a first approximation we will have @xmath197 given that @xmath198 is of order @xmath94 , all higher moments of event shapes will receive power corrections of order @xmath199 , which is parametrically larger than @xmath192 . strictly speaking it is not possible to guarantee the coefficient of @xmath199 given in eq . ( [ eq : momentpc ] ) , since the shift approximation on which it is based holds only in the @xmath16-jet region , whereas the dominant contribution to the relevant integral comes from the @xmath65-jet region . for the @xmath38-parameter an alternative coefficient for the @xmath199 power correction has been proposed in  @xcite , however since it also is not based on a full calculation of the power correction in the three - jet region ( which does not yet exist ) , it is subject to precisely the same reservations as eq . ( [ eq : momentpc ] ) . we note also an interesting result by gardi  @xcite regarding the exact renormalon analysis of @xmath200 . while it is clear that the physical answer @xmath199 will not appear in a leading renormalon analysis ( the extra factor of @xmath94 means that it is associated with a subleading renormalon ) the calculation @xcite has the surprising result that the leading renormalon contribution is not @xmath147 as naively expected but rather @xmath201 . considerable progress has been made in recent years in the calculation of power corrections for three - jet observables ( those that vanish in the two and the three - jet limits ) . in the three - jet limit , there exist explicit results for the thrust minor  @xcite and the @xmath39-parameter  @xcite in @xmath202 , the out of plane momentum in drell - yan plus jet  @xcite and in @xmath90-jet dis  @xcite , as well as azimuthal correlations in dis @xcite . except for the @xmath39-parameter , all these calculations involve perturbative , non - perturbative interplay in a manner similar to that of the broadenings and back - to - back eec . the explicit forms for the results are rather complicated and so we refer the reader to the original publications . simpler ( though numerical ) results have been obtained  @xcite in the case of the @xmath39-parameter ( confirmed also by @xcite ) integrated over all 3-jet configurations , where one finds @xmath203 the power of @xmath94 comes from the matrix - element weighting of the @xmath65-jet configurations and the fact that the @xmath39-parameter vanishes in the @xmath16-jet limit . a final @xmath65-jet result of interest  @xcite is that for the @xmath38-parameter just above the sudakov shoulder , at @xmath125 . this is the only case where a proper calculation exists for the power correction to a two - jet observable in the three jet limit : @xmath204 this is slightly less than half the correction that appears in the two - jet limit ( @xmath205 ) . a point that emerges clearly from eqs . ( [ eq : dpower ] ) and ( [ eq : c34power ] ) , but that is relevant for all @xmath65-jet studies is that the power correction acquires a dependence on @xmath139 , i.e.  there is sensitivity to hadronisation from a gluon . in situations where one selects only three - jet events there is additionally non trivial dependence on the geometry , due to the coherence between the three jets . comparisons to data for such observables would therefore allow a powerful test of the renormalon - inspired picture . in particular , other models , such as the flux tube model  @xcite ( based essentially on boost invariance along the @xmath206 axis ) , which in the @xmath16-jet limit give the same predictions as renormalon - inspired approaches , can not as naturally be extended to the @xmath65-jet case . one of the most widespread ways of testing the dokshitzer - webber approach to hadronisation is to carry out simultaneous fits to the mean value data for @xmath1 and @xmath163 . figure  [ fig : meandata ] , where the term labelled @xmath4 is actually of the form @xmath207 , shows the good agreement that is obtained with the thrust data . the quality of agreement is similar for other observables ( see for example fig .  11 of  @xcite ) . the true test of the approach however lies in a verification of the universality hypothesis , namely that @xmath163 is the same for all observables and processes . , one may wonder why one also fits for @xmath1 . the two principal reasons are ( a ) that one is in any case interested in the value of @xmath94 and ( b ) that the data and perturbative prediction differ also by higher - order terms in @xmath94 and fixing @xmath94 would mean trying to fit these higher - order terms with a power correction , which would be misleading . ] while there are strong general reasons for expecting universality within @xmath0 ( boost invariance and the flux - tube model ) , universality across difference processes is less trivial  for example one could imagine dis hadronisation being modified by interactions between outgoing low - momentum gluons ( or hadrons ) and the ` cloud ' of partons that make up the proton remnant ( present even for a fast - moving proton , since the longitudinal size of the cloud is always of order @xmath208 ) . figure  [ fig : meanfits ] shows @xmath161-@xmath209 contours for fits of @xmath1 and @xmath163 in @xmath0  @xcite ( left , all experiments combined ) and dis ( right , h1  @xcite and zeus @xcite merged into one plot ) . aside from the different scales used , some care is required in reading the figures because of different treatments of the errors in the different plots . for example h1 include experimental systematic errors whereas zeus do not . in the @xmath0 fit , the systematics are included but treated as uncorrelated , which may underestimate the final error . in neither figure are theoretical uncertainties included : over the available range of @xmath0 energies , renormalisation scale dependence gives an uncertainty of about @xmath210 on @xmath94 , while a ` canonical ' variation of the milan factor @xmath211 by @xmath212 ( to allow for higher - order terms ) leads to an uncertainty of @xmath213 on @xmath163  @xcite . in dis the corresponding uncertainties are larger for @xmath94 ( @xmath214 , essentially because the fits are dominated by lower @xmath82 values , and they vary substantially for @xmath163 ( larger for photon - axis observables , smaller for the others ) . also to be kept in mind is that the dis figure includes the jet mass , @xmath215 in the default scheme , and in @xmath0 the origins of the arrows indicate the default - scheme results for @xmath216 and @xmath196 . since the default scheme for jet masses breaks universality ( cf . section  [ sec : dwdetails ] ) these results should not be compared directly to those for other observables . , @xmath163 fit to @xmath0 data on mean values of event shapes , with the jet masses converted to the @xmath62-scheme ( the arrow indicates the change in going from default to @xmath62-scheme ) @xcite . right : fits to dis mean event - shape data by h1 @xcite and zeus  @xcite ( curves taken from figures in  @xcite ) . contours indicate 1-@xmath209 uncertainties ( statistical and experimental systematic , except for the zeus results which are just statistical).,title="fig:",scaledwidth=48.0% ] , @xmath163 fit to @xmath0 data on mean values of event shapes , with the jet masses converted to the @xmath62-scheme ( the arrow indicates the change in going from default to @xmath62-scheme ) @xcite . right : fits to dis mean event - shape data by h1 @xcite and zeus  @xcite ( curves taken from figures in  @xcite ) . contours indicate 1-@xmath209 uncertainties ( statistical and experimental systematic , except for the zeus results which are just statistical).,title="fig:",scaledwidth=48.0% ]    taking into account the uncertainties , which are not necessarily correlated from one observable to another , fig . [ fig : meanfits ] indicates remarkable success for the renormalon - inspired picture . the results are in general consistent with a universal value for @xmath163 in the range @xmath217 to @xmath218 . furthermore the @xmath0 results for @xmath94 are in good agreement with the world average . on the other hand the dis results for @xmath94 seem to be somewhat larger than the world average . the discrepancy is ( just ) within the theoretical systematic errors , however it remains a little disturbing and one wonders whether it could be indicative of large higher - order corrections , or some other problem ( known issues include the fact that the low-@xmath82 data can bias the fits , for example because of heavy - quark effects and @xmath219 corrections ; and also at high-@xmath82 there is @xmath220 interference which is not usually accounted for in the fixed - order calculations ) . there is potentially also some worry about the @xmath221 result which has an anomalously low @xmath163 and large @xmath94 ( and additionally shows unexpectedly substantial @xmath176 dependence  @xcite ) . another point relates to the default scheme jet masses  though the universality - breaking effects should be purely non - perturbative , they have an effect also on the fit results for @xmath1 . this is a consequence of the anomalous dimension that accompanies @xmath146 hadron mass effects . similar variations of @xmath94 are seen when varying the set of particles that are considered stable  @xcite . though as discussed above some issues remain , they should not however be seen as detracting significantly from the overall success of the approach and the general consistency between @xmath0 and dis . we also note that the first moment of the coupling as extracted from studies of heavy - quark fragmentation @xcite is quite similar to the value found for event shapes . there have also been experimental measurements of higher moments of event shapes , notably in  @xcite . the parameters @xmath1 and @xmath163 are fixed from fits to mean values and then inserted into eq . ( [ eq : momentpc ] ) to get a prediction for the second moments . for @xmath222 this gives very good agreement with the data ( though the fact that @xmath196 is in the default scheme perhaps complicates the situation ) , while @xmath223 and @xmath224 show a need for a substantial extra correction . in  @xcite this extra contribution is shown to be compatible with a @xmath219 term , though it would be interesting to see if it is also compatible with a @xmath225 term with a modified coefficient . currently ( to our knowledge ) no fits have been performed for three - jet event shapes , though we understand that such fits are in progress @xcite for the mean value of the @xmath39-parameter . we look forward eagerly to the results . in the meantime it is possible to verify standard parameters for @xmath94 and @xmath163 against a single published point at @xmath127  @xcite and one finds reasonable agreement within the uncertainties . two other points of view that have also been used in analysing data on @xmath5 event shapes are the approach of _ dressed gluons _ employed by gardi and grunberg @xcite and the use of _ renormalisation group improved _ perturbation theory  @xcite as carried out by the delphi collaboration  @xcite .      in the gardi - grunberg approach one uses the basic concept of renormalons ( on which the dw model is also based ) , but one choses to treat the renormalon integral differently to the dw model . this renormalon integral is ill - defined due to the landau singularity in the running coupling and instead of assuming an infrared finite coupling below some matching scale @xmath160 as was the case in the dw model , gardi and grunberg define the renormalon integral by its principal value ( a discussion of the theoretical merits of different approaches has been give in  @xcite ) . in doing this they explicitly include higher order renormalon contributions in their perturbative result , rather than including them via power behaved corrections below scale @xmath160 that result from assuming an infrared finite coupling ( though they do also compare to something similar , which they call a ` cutoff ' approach ) . however since one is dealing with an ambiguous integral ( a prescription other than a principal value one would give a result differing by an amount proportional to a power correction ) one must still allow for a power correction term . gardi and grunberg studied the mean thrust in @xmath5 annihilation and used the following form for fitting to the data : @xmath226+\frac{\lambda}{q}\ ] ] where the subscript pv denotes the principal value of the renormalon integral @xmath227 for the thrust ( which includes the full lo contribution and parts of the higher - order contributions ) and @xmath228 is a piece that accounts for the difference between the true nlo coefficient for @xmath14 and the contribution included in @xmath227 . the best fit values obtained for @xmath229 and @xmath230 are respectively @xmath231 and @xmath232 . we note that @xmath233 in this approach corresponds to a smaller power correction than is required in the dw model . this is probably a consequence of the inclusion of pieces of higher perturbative orders via a principal value prescription , although it leads to a somewhat small value ( compared to the world average ) of the coupling at scale @xmath234 . next we turn to the renormalisation group improved ( rgi ) perturbative estimates  @xcite that have also been used to compare with event shape data for the mean value of different event shapes in @xmath5 annihilation . the basic idea behind this approach is to consider the dependence of the observable on the scale @xmath82 , which can be expressed using renormalisation group invariance as @xmath235 where in the above formula @xmath236 , @xmath237 , with @xmath238 being the mean value of a given event shape and @xmath239 being the coefficient of @xmath240 in its perturbative expansion .     from comparisons to @xmath5 event - shape data using the rgi approach . the band shows the mean value of @xmath94 and uncertainty , as obtained from the observables below the dashed line . figure taken from ref . @xcite.,scaledwidth=40.0% ]    thus the first term in the perturbative expansion of @xmath109 is simply @xmath241 . the @xmath242 are renormalisation scheme independent quantities . in particular the quantity @xmath243 is simply the ratio of the first two coefficients of the qcd @xmath145 function while @xmath244 additionally depends on the first three ( up to nnlo ) perturbative coefficients in the expansion of @xmath245 . current studies using this method are therefore restricted to an nlo analysis involving @xmath246 alone . with this simplification , the solution of eq . ( [ eq : rgi ] ) simply corresponds to the introduction of an observable - specific renormalisation scheme and associated scale parameter @xmath247 which is such that it sets the nlo perturbative term to be zero ( and neglects higher perturbative orders ) . its relation to the standard @xmath248 is easily obtained : @xmath249 where @xmath250 and @xmath251 , with @xmath252 being the nlo coefficient in the @xmath253 expansion for the observable . the terms in the bracket account for different definitions of the coupling in terms of @xmath254 as used in the @xmath109 and @xmath253 schemes .    in the specific case of fits to @xmath5 event shape mean values , the above approach has met with a considerable success @xcite . in particular the introduction of @xmath247 seems to remove any need for a significant power correction . the delphi collaboration  @xcite have performed a combined fit of the parameter @xmath247 and a parameter @xmath255 that quantifies the power correction  @xcite . the value of @xmath255 was found to be consistent with zero in most cases and rather small in all cases which indicates that there is no real need for a power behaved correction once the perturbative expansion is fixed through the rgi technique . the agreement between @xmath1 values extracted from different observables is impressive ( see figure  [ fig : rgi ] ) with a spread that is only about half as large as that obtained using the standard perturbative result supplemented by a power correction term , using the dw initiated model . while it is clear that the size of the power correction piece inevitably depends on how one choses to define the perturbative expansion ( in fact all the methods discussed thus far , including the dw model and the gardi - grunberg approach , account for this effect ) , it is nevertheless very interesting that for several different observables , with significantly different perturbative coefficients , one observes after the introduction of rgi , hardly any need for a power correction term . certainly it is not clear on the basis of any physical arguments , why the genuine non - perturbative power correction should be vanishingly small and that a perturbative result defined in a certain way should lead to a complete description of the data . the clarification of this issue is still awaited and the above findings are worth further attention and study .    a final point to be kept in mind about the rgi approach in the above form , is that it is valid only for inclusive observables that depend on just a single relevant scale parameter @xmath82 . this makes its applicability somewhat limited and for instance it is not currently clear how to extend the procedure so as to allow a meaningful study of event shape distributions in @xmath5 annihilation or dis event shape mean values and distributions ( which involve additional scales ) . so far our discussion of comparisons between theory and experiment has been limited to a study of mean values of the event shapes . as we have seen , there is some ambiguity in the interpretation of these comparisons , with a range of different approaches being able to fit the same data . this is in many respects an unavoidable limitation of studies of mean values , since the principal characteristic of the different models that is being tested is their @xmath82-dependence , which can also be influenced by a range of ( neglected ) higher - order contributions . in contrast , the full distributions of event shapes contain considerably more information and therefore have the potential to be more discriminatory . in section  [ sec : perturbative ] we discussed the perturbative calculation of distributions . as for mean values though , the comparison to data is complicated by the need for non - perturbative corrections . these corrections however involve many more degrees of freedom than for mean values , and there are a variety of ways of including them . accordingly we separate our discussion into two parts : in section  [ sec : evgen ] we shall consider studies in which hadronisation corrections are taken from monte carlo event generators , and where the main object of study is the perturbative distribution , with for example fits of the strong coupling . in section [ sec : anldist ] we shall then consider studies which involve analytical models for the non - perturbative corrections , and where it is as much the non - perturbative models , as the perturbative calculations that are under study .       compared to a variety of event generators  @xcite . the shaded band in the lower plot indicates the experimental statistical and systematic errors added in quadrature . figure taken from  @xcite.,title="fig:",scaledwidth=48.0% ]   compared to a variety of event generators  @xcite . the shaded band in the lower plot indicates the experimental statistical and systematic errors added in quadrature . figure taken from  @xcite.,title="fig:",scaledwidth=48.0% ]      before considering the purely perturbative studies mentioned above , let us recall that one major use of event shape distributions has been in the testing and tuning @xcite of monte carlo event generators such as herwig  @xcite , jetset / pythia  @xcite and ariadne  @xcite . figure  [ fig : tuning ] shows comparisons for two event shapes , the thrust and thrust major , and the agreement is remarkable testimony to the ability of the event generators to reproduce the data . the use of event generators to probe the details of qcd is unfortunately rather difficult , essentially because they contain a number of parameters affecting both the non - perturbative modelling and in some cases the treatment of the perturbative shower . furthermore though considerable progress is being made in matching to fixed order calculations ( see for example  @xcite ) , event generators are in general able to guarantee neither the nll accuracy nor the nlo accuracy of full matched nll - nlo resummed calculations . event - shapes ( global @xcite two - jet event shapes ) herwig @xcite is expected to be correct to nll accuracy . ] nevertheless , the good agreement of the event generators with the data suggests that the bulk of the dynamics is correct and in particular that a good model for the hadronisation corrections can be had by comparing parton and hadron ` levels ' of the generator . there has been a very widespread use of event generators in this way to complement the nll+nlo perturbative calculations . such a method has both advantages and drawbacks and it is worth devoting some space to them here . on one hand , the parton level of an event generator is not a theoretically well - defined concept  it is regularised with some effective parton mass or transverse momentum cutoff , which already embodies some amount of non - perturbative correction . in contrast the nlo+nll partonic prediction integrates down to zero momenta without any form ( or need ) of regularisation . this too implies some amount of non - perturbative contribution , but of a rather different nature from that included via a cutoff . the resulting difference between the event - generator and the purely perturbative nlo+nll parton levels means that the ` hadronisation ' that must be applied to correct them to hadron level is different in the two cases . there are nevertheless reasons why event - generators are still used for determining the hadronisation corrections . the simplest is perhaps that they give a very good description of the data ( cf . [ fig : tuning ] ) , which suggests that they make a reasonable job of approximating the underlying dynamics . furthermore the good description is obtained with a single set of parameters for all event shapes , whereas as we shall see , other approaches with a single ( common ) parameter are currently able to give equally good descriptions only for a limited number of event shapes at a time . additionally , the objection that the monte carlo parton level is ill - defined can , partially , be addressed by including a systematic error on the hadronisation : it is possible for example to change the internal cutoff on the monte carlo parton shower and at the same time retune the hadronisation parameters in such a way that the monte carlo description of the hadron - level data remains reasonable . in this way one allows for the fact that the connection between monte carlo and nll+nlo ` parton - levels ' is not understood . this is the procedure that has been used for example in @xcite , and the results are illustrated in figure  [ fig : sldhadr ] as a hadronisation correction factor @xmath256 with an uncertainty corresponding to the shaded band . , with uncertainties , as discussed in the text ) obtained from the ratio of jetset parton to hadron - level distributions . it is compared to nlo and nll+nlo predictions . also shown is the detector to hadron - level correction factor @xmath257 . figure taken from  @xcite . ] two further points are to be made regarding monte carlo hadronisation corrections . firstly , as is starting now to be well known , it is strongly recommended not to correct the data to parton - level , but rather to correct the theoretical perturbative prediction to hadron level . this is because the correction to parton level entails the addition of extra assumptions , which a few years later may be very difficult to deconvolve from the parton - level ` data ' . it is straightforward on the other hand to apply a new hadronisation model to a perturbative calculation . secondly it is best not to treat the hadronisation correction as a simple multiplicative factor . to understand why , one should look again at figure  [ fig : sldhadr ] ( which , dating from several years ago , corrects data to parton level ) . the multiplicative factor @xmath256 ( parton@xmath258hadron ) varies very rapidly close to the peak of the distribution and goes from above @xmath161 to below @xmath161 . this is because the distribution itself varies very rapidly in this region and the effect of hadronisation is ( to a first approximation , see below ) to shift the peak to larger values of the observable . but a multiplicative factor , rather than shifting the peak , suppresses it in one place and recreates it ( from non - peak - like structure ) in another . instead of applying a multiplicative correction , the best way to include a monte carlo hadronisation correction is to determine a transfer matrix @xmath259 which describes the fraction of events in bin @xmath56 at parton level that end up in in bin @xmath55 at hadron level . then for a binned perturbative distribution @xmath260 , the binned hadronised distribution @xmath261 is obtained by matrix multiplication , @xmath262 . this has been used for example in @xcite . having considered the basis and methods for including hadronisation corrections from event generators , let us now examine some of the perturbative studies made possible by this approach . the majority of them use nll+nlo perturbative calculations to fit for the strong coupling ( as discussed below some , e.g.  @xcite , just use nlo results , which exist for a wider range of observables ) and such fits can be said to be one of the principal results from event - shape studies . in general there is rather good agreement between the data and the resummed predictions ( cf . [ fig : sldhadr ] ) .     obtained by fitting nll+nlo event - shape distributions to l3 data at a range of centre of mass energies ( figure taken from  @xcite ) . only experimental uncertainties are shown.,scaledwidth=60.0% ]    an example of such a study is given in figure  [ fig : l3asevln ] , which shows results for fits of @xmath1 to l3 data at a range of energies ( including points below @xmath263 from radiative events which , as discussed earlier , are to be taken with some caution ) . one important result from such studies is an overall value for @xmath1 , which in this case is @xmath264 @xcite , where the first error is experimental and the second is theoretical . similar results have been given by other collaborations , e.g.  @xmath265 from aleph  @xcite . there is also good evidence for the evolution of the coupling , though leaving out the potentially doubtful radiative points leaves the situation somewhat less clearcut , due to the limited statistics at high energies ( another resummed study which gives good evidence for the running of the coupling uses jade and opal data  @xcite , however it is limited to jet rates ) .    these and similar results from the other lep collaborations and jade highlight two important points . firstly at higher energies there is a need to combine results from all the lep experiments in order to reduce the statistical errors and so improve the evidence for the running of the coupling in the region above @xmath127 . work in this direction is currently in progress @xcite . so far only preliminary results are available @xcite , and as well as improving the picture of the high - energy running of the coupling , they suggest a slightly smaller value of @xmath1 than that quoted above , more in accord with world averages  @xcite .    secondly , in the overall result for @xmath266 , by far the most important contribution to the error is that from theoretical uncertainties . however theoretical uncertainties are notoriously difficult to estimate , since they relate to unknown higher orders . systematic investigations of the various sources of uncertainty have been carried out in  @xcite and in particular @xcite proposes a standard for the set of sources of uncertainty that ought to be considered , together with an approach for combining the different sources into a single overall uncertainty on @xmath94 . one of the main principles behind the method is that while one may have numerous estimates of sources of theoretical uncertainty , it is not advisable to combine them in quadrature , as one might be tempted to do , because this is likely to lead to a double - counting of uncertainties . rather , one should examine the different sources of uncertainty across the whole range of the distribution and at each point , take the maximum of all sources to build up an uncertainty envelope or band , represented by the shaded area in figure  [ fig : lepqcdbands ] ( shown relative to a reference prediction for a ` standard ' theory ) . the overall uncertainty on the coupling ( rather than the distribution ) is given by the range of variation of @xmath94 such that the prediction remains within the band . it is to be noted that this kind of approach is of relevance not just to event - shape studies but also quite generally to any resummed matched calculation , which is inherently subject to many sources of arbitrariness . allowed such that the prediction stays within the uncertainty band.,scaledwidth=63.0% ]      one of the drawbacks of nll+nlo studies is that the nll resummed predictions exist for only a fraction of observables . recently there has been an extensive study by the delphi collaboration @xcite of 18 event - shape distributions compared to nlo calculations . since in most cases nlo calculations with a renormalisation scale @xmath267 tend to describe the data rather poorly , they have examined various options for estimating higher orders so as to improve the agreement . their main approach involves a simultaneous fit of the coupling and the renormalisation scale , i.e.   the renormalisation scale dependence is taken as a way of parameterising possible higher orders . in such a procedure the choice of fit - range for the observable is rather critical , since the true higher - orders may involve have quite a different structure from that parameterised by the @xmath168 dependence . the delphi analysis restricts the fit range to be that in which a good fit can be obtained . with these conditions they find relatively consistent values of @xmath94 across their whole range of observables ( finding a scatter of a few percent ) , with a relatively modest theoretical uncertainty , as estimated from the further variation of @xmath168 by a factor of @xmath13 to @xmath16 around the optimal scale . their overall result from all observables treated in this way is @xmath268 ( including both experimental and theoretical errors ) . they also compare this approach to a range of other ways of ` estimating ' higher orders . they find for example that some theoretically based methods for setting the scale ( principle of minimal sensitivity  @xcite , effective charge @xcite ) lead to scales that are quite correlated to the experimentally fitted ones , while another method ( blm @xcite ) is rather uncorrelated . they compare also to resummation approaches , though there is no equivalent renormalisation scale choice that can be made . instead the comparison involves rather the quality of fit to the distribution and the final spread of final @xmath1 values as a measure of neglected higher - order uncertainties . not surprisingly , they find that while the nll results work well in the two - jet region , in the three - jet region the combined nll+nlo predictions fare not that much better than the pure ( @xmath267 ) nlo results . they go on to argue that nll+nlo is somewhat disfavoured compared to nlo with an ` optimal ' scale . this statement is however to be treated with some caution , since the fit range was chosen specifically so as to obtain a good fit for nlo with the ` optimal ' scale  one could equally have chosen a fit range tuned to the nll+nlo prediction and one wonders whether the nlo with optimal scale choice would then have fared so well . regardless of this issue of fit range , it is interesting to note that over the set of observables which can be treated by both methods , the total spread in @xmath1 values is quite similar , being @xmath269 in the case of nlo with optimal scale choice and @xmath270 for nll+nlo . this is similar to the spread seen for fits to mean values with analytical hadronisation models ( fig . [ fig : meanfits ] ) , though the rgi fits of  @xcite have a somewhat smaller error . the jet rates ( especially the durham and cambridge algorithms and certain members of the jade family ) are among the few observables for which the pure nlo calculation gives a reasonable description of the distribution ( cf .  table  3 of  @xcite ) . one particularly interesting set of nlo studies makes use of the 3-jet rate as a function of @xmath271 ( this is just the integral of the distribution of @xmath137 ) in events with primary @xmath144-quarks as compared to light - quark events @xcite ( some of the analyses use other observables , such as the 4-jet rate or the thrust ) . using nlo calculations which account for massive quarks @xcite makes it possible , in such studies , to extract a value for the @xmath144 mass at a renormalisation scale of @xmath263 , giving first evidence of the ( expected ) running of the @xmath144-quark mass , since all the analyses find @xmath272 in the range @xmath273 to @xmath274  gev ( with rather variable estimates of the theoretical error ) . we have already seen , in section  [ sec : meanpc ] , that analytical models for hadronisation corrections , combined with normal perturbative predictions , can give a rather good description of mean values . there the inclusion of a non - perturbative ( hadronisation ) correction was a rather simple affair , eq . ( [ eq : meancorrections ] ) , since it was an additive procedure . in contrast the general relation between the full distribution @xmath275 for an observable @xmath276 and the perturbatively calculated distribution @xmath277 is more complicated , @xmath278 where @xmath279 encodes all the non - perturbative information . ( [ eq : shapegenconvoluion ] ) can be seen as a generalised convolution , where the shape of the convolution function depends on the value @xmath86 of the variable as well as the coupling and hard scale . such a general form contains far more information however than can currently be predicted , or even conveniently parameterised . accordingly analytical hadronisation approaches generally make a number of simplifying assumptions , intended to be valid for some restricted range of the observable .      the most radical simplification of eq . ( [ eq : shapegenconvoluion ] ) that can be made is to replace @xmath280 with a @xmath281-function , @xmath282 , leading to @xmath283 this was proposed and investigated in  @xcite , and the combination @xmath284 is the same that appears in the power correction for the mean value  @xcite , discussed in section  [ sec : meanpc ] . ( [ eq : npshift ] ) holds for observables with an @xmath1 independent power correction in tables  [ tab : cvs ] and [ tab : cvsdis ] , in the region @xmath285 , where the lower limit ensures that one can neglect the width of the convolution function @xmath286 , while the upper limit is the restriction that one be in the born limit ( in which the @xmath152 were originally calculated ) . for certain 3-jet observables a similar picture holds , but with a shift that depends on the kinematics of the 3-jet configuration  @xcite .    for the broadenings , the situation is more complicated because the power correction is enhanced by the rapidity over which the quark and thrust axes can be considered to coincide , @xmath287 , @xmath288 being the angle between thrust and quark axes  this angle is strongly correlated with the value of the broadening ( determined by perturbative radiation ) and accordingly the extent of the shift becomes @xmath252-dependent  @xcite . the simplest case is the @xmath0 wide - jet broadening , for which one has @xmath289 where the shift is written for the integrated distribution @xmath290 in order to simplify the expressions . the @xmath291 enhancement for the power correction to the mean broadening ( tables  [ tab : cvs ] and [ tab : cvsdis ] ) comes about simply because the average value of @xmath292 , after integration over @xmath86 with the resummed distribution , is of order @xmath293 . the full form for @xmath294 and analogous results for the total @xmath0 and dis broadenings have been given in  @xcite , with results existing also for the thrust minor  @xcite , and the dis and drell - yan out - of plane momenta  @xcite . yet subtler instances of perturbative , non - perturbative interplay arise for observables like the eec  @xcite and dis azimuthal correlation  @xcite , with the appearance of _ fractional _ powers of @xmath82 in the power correction , as mentioned before . -dependent non - perturbative shift in the jet broadening distribution for @xmath5 annihilation . the dashed curve is the perturbative prediction ( nll+nlo ) , while the solid curve includes a non - perturbative ( @xmath252-dependent ) shift . figure taken from  @xcite . right : nll+nlo distribution with a @xmath4 shift for @xmath295 in dis at a range of @xmath82 values ( round dots indicate points used in fits )  @xcite . the @xmath82 values range from 15  gev ( bottom ) to 81  gev ( top).,title="fig:",scaledwidth=44.0%,scaledwidth=48.0% ] -dependent non - perturbative shift in the jet broadening distribution for @xmath5 annihilation . the dashed curve is the perturbative prediction ( nll+nlo ) , while the solid curve includes a non - perturbative ( @xmath252-dependent ) shift . figure taken from  @xcite . right : nll+nlo distribution with a @xmath4 shift for @xmath295 in dis at a range of @xmath82 values ( round dots indicate points used in fits )  @xcite . the @xmath82 values range from 15  gev ( bottom ) to 81  gev ( top).,title="fig:",scaledwidth=48.0% ]    one might a priori have thought that the formal domain of validity of the shift approximation , @xmath296 , would be somewhat limited . figure  [ fig : broadpc ] ( left ) illustrates what happens in the case of @xmath297  quite remarkably the shift describes the data well over a large range of @xmath297 , with slight problems only in the extreme two - jet region ( peak of the distribution ) and in the four - jet region ( @xmath298 ) . similar features are seen in the right - hand plot of fig . [ fig : broadpc ] for @xmath299 in dis at a range of @xmath82 values : a large part of the distribution is described for all @xmath82 values , and problems appear only in the @xmath91 jet region ( @xmath300 ) , and in the @xmath122 jet region for low @xmath82 . this success is reproduced for quite a range of observables in @xmath0 and dis . -@xmath209 contours for simultaneous fits of @xmath94 and @xmath163 to @xmath0 event - shape distributions . the error contours account for statistic , systematical and theoretical uncertainties . the shaded band indicates the result for @xmath1 when using monte carlo hadronisation models . figure taken from @xcite . right : @xmath161-@xmath209 contours from fits to h1 dis event shape distributions  @xcite with statistical and systematic errors added in quadrature . figure taken from  @xcite.,title="fig:",scaledwidth=48.0% ] -@xmath209 contours for simultaneous fits of @xmath94 and @xmath163 to @xmath0 event - shape distributions . the error contours account for statistic , systematical and theoretical uncertainties . the shaded band indicates the result for @xmath1 when using monte carlo hadronisation models . figure taken from @xcite . right : @xmath161-@xmath209 contours from fits to h1 dis event shape distributions  @xcite with statistical and systematic errors added in quadrature . figure taken from  @xcite.,title="fig:",scaledwidth=48.0% ]    as was the case for mean values , a more systematic test involves carrying out a simultaneous fit of @xmath94 and @xmath163 for each observable and then checking for consistency between observables and with the corresponding results for mean values . results from @xmath0 and dis are shown in figure  [ fig : distcont ] , taken from @xcite and @xcite . other fits that have been carried out in recent years include  @xcite .- parameter in @xmath65-jet events . ]    just as for the mean values , one notes that ( with exception of @xmath301 and @xmath302 in @xmath0 , to be discussed shortly ) there is good consistency between observables , both for @xmath0 and dis . this statement holds holds also for the eec , fitted in @xcite , not shown in fig . [ fig : distcont ] . there is also good agreement with the results for the mean values , fig . [ fig : meanfits ] , except marginally for dis @xmath94 results , which for distributions are in better accord with the world average . the good agreement between distributions and mean values is important in the light of alternative approaches to fitting mean values , such as the rgi which , we recall , seems to suggest that the power correction for mean values can just as well be interpreted as higher order contributions . were this really the case then one would expect to see no relation between @xmath163 results for mean values and distributions .    despite the generally good agreement , some problems do persist . examining different groups results for the @xmath1-@xmath163 fits one finds that for some observables there is a large spread in the results . for example for the @xmath0 @xmath297 , aleph @xcite find @xmath303 while the jade result  @xcite shown in fig . [ fig : distcont ] corresponds to @xmath304 ( and this differs from a slightly earlier jade result @xmath305  @xcite ) . while the results all agree to within errors , those errors are largely theoretical and would be expected to be a common systematic on all results . in other words we would expect the results to be much closer together than @xmath161-@xmath209 . that this is not the case suggests that the fits might be unstable with respect to small changes in details of the fit , for example the fit range ( as has been observed when including data in the range @xmath306 in the fit corresponding to the right hand plot of figure  [ fig : distcont ]  @xcite ) . these differences can lead to contradictory conclusions regarding the success of the description of different observables and it would be of value for the different groups to work together to understand the origin of the differences , perhaps in a context such as the lep qcd working group . there do remain two cases however where there seems to be clear incompatibility between the data and theoretical prediction : the heavy - jet mass , @xmath196 ( also referred to as @xmath302 ) and the wide - jet broadening , @xmath301 . the heavy - jet mass is subject to @xmath146 hadron mass effects , as discussed already in the context of mean values , however even after accounting for this by considering @xmath62 or @xmath57-scheme data , as has been done by delphi  @xcite , the results appear to be rather inconsistent with other observables : @xmath1 is too small and @xmath163 too large , as for @xmath301 , with a fixed @xmath307 lead to a consistent value for @xmath308 , though the overall fit quality is poor . ] though less extreme . the common point between these two observables is that they both select the heavier / wider hemisphere and so are less inclusive than other observables . studies with event generators  @xcite suggest that whereas for more inclusive observables the 2-jet limit of the shift works well even into the 3-jet region , for heavy - wide observables the 3-jet limit behaves very differently from the 2-jet limit and these differences becomes relevant even for relatively low values of the observable ( some partially related discussion of the problem has been given also in  @xcite ) . we will come back to this and other potential explanations below when we discuss shape functions .    beyond the ` usual ' observables discussed above , there are some recent results also for the @xmath38-parameter shoulder at @xmath125 , where @xcite the shoulder position is seen to be consistent with a shift with the usual 2-jet coefficient @xmath205  this is in contradiction with the ( smaller ) result eq . ( [ eq : c34power ] ) for the shift to the shoulder ; however the analysis  @xcite does nt take into account resummation effects at the shoulder which may also contribute to an effective shift . a first analysis of corrections to the @xmath309 distribution has been given in @xcite which shows evidence for @xmath310 contributions in the @xmath309 distribution at the two lowest jade centre of mass energies ( 14 and 22  gev ) . one last interesting point to discuss regarding the shift approximation for hadronisation is an application to the measurement of the qcd colour factors . the fact that the hadronisation is encoded via just a single parameter , with all the dependence on colour factors available explicitly for the perturbative and non - perturbative contributions , means that it is feasible to carry out a simultaneous fit not only for @xmath311 and @xmath163 but also for one or more of the qcd colour factors @xmath139 , @xmath138 and @xmath312  @xcite . this would have been much more difficult using monte carlo hadronisation because all the hadronisation parameters would have had to have been retuned for each change of the colour factors . the results for this fit are summarised and compared to other approaches in figure  [ fig : colfact ]  @xcite , which shows @xmath16-@xmath209 contours for @xmath138 and @xmath139 . while other approaches constrain the ratio @xmath313 , the event - shape approach gives by far the best independent measurement of @xmath139 and @xmath138 . as a result the combination of the results gives rather tight limits on the colour factors , which are in good agreement with the qcd expectation . ( shift ) power corrections and a variety of other methods . figure taken from  @xcite . the ellipses show correlated measurements from studies of the angular distributions of 4-jet events @xcite and event shapes @xcite , while the lines are obtained from studies of the gluon fragmentation function  @xcite and from multiplicities in gluon jets  @xcite.,scaledwidth=51.0% ]      we have already mentioned that the approximation of a simple shift is valid only in the region @xmath296 . to better understand the origin of the lower limit , let us expand eq . ( [ eq : npshift ] ) , @xmath314 where @xmath315 is the @xmath316 derivative of @xmath317 with respect to @xmath86 . since @xmath315 goes as @xmath318 one sees that when the shift @xmath319 all terms in the series are of the same order . while the @xmath320 term can be related to the power correction to the mean value of @xmath276 , the shift merely provides an ansatz for the higher terms . the breakdown for small values of @xmath321 is clearly visible for the ( dis ) thrust in fig . [ fig : broadpc ] ( right ) .    a physically transparent way of dealing with the @xmath322 terms has been developed by korchemsky , sterman and collaborators @xcite . one approximates @xmath323 in eq . ( [ eq : shapegenconvoluion ] ) by a _ shape function _  @xcite @xmath324 ( rather than the @xmath281-function which corresponds to a pure shift ) , @xmath325 in eq . ( [ eq : expandedshift ] ) , @xmath326 is then replaced by @xmath327 where @xmath328 is the @xmath316 moment of @xmath329 . there is considerable freedom in one s choice of the form for the shape function . accordingly , success in fitting a given event - shape distribution for any single value of @xmath82 can not , alone , be considered strong evidence in favour of the shape - function picture . the true test of the approach comes by determining the shape function at one value of @xmath82 and then establishing whether it applies for all values of @xmath82 . this is illustrated for the @xmath0 thrust in fig . [ fig : thrustshape ] ( left ) , where the following shape - function has been used @xcite , @xmath330 with two free parameters , @xmath331 and @xmath254 ( it has been argued that the behaviour of this function for small @xmath176 and large @xmath176 is similar to that expected from a simple physical model for the dynamics at the origin of the shape function  @xcite ) . one notes the remarkable agreement , far into the two - jet region , for the whole range of @xmath82 values .    ) . shown for @xmath332 and @xmath333  gev . figure taken from  @xcite . right : the thrust distribution at the perturbative ( dge ) level , with a non - perturbative shift , and with a shape function . figure taken from  @xcite.,title="fig:",scaledwidth=48.0%,scaledwidth=48.0% ] ) . shown for @xmath332 and @xmath333  gev . figure taken from  @xcite . right : the thrust distribution at the perturbative ( dge ) level , with a non - perturbative shift , and with a shape function . figure taken from  @xcite.,title="fig:",scaledwidth=48.0%,scaledwidth=48.0% ]    this is a clear success for the shape - function approach as compared to the simple shift discussed above . however the drawback of shape - functions is that a priori one loses all predictivity ( save for the first moment ) when going from one observable to another . this predictivity was one of the main strengths of the shift approach . considerable work has therefore been devoted to understanding both the general properties of shape functions and the relations between shape functions of different observables . one line of investigation has been to relate the shape functions for the thrust and heavy - jet mass ( @xmath196 )  this can be done with the help of the observation  @xcite that in the two - jet limit the thrust is just the sum of the light and heavy - jet ( squared normalised ) masses . in  @xcite korchemsky and tafat introduced a shape function @xmath334 describing the distribution of non - perturbative radiation in each of the two hemispheres . the shape function for the thrust is then obtained by @xmath335 instead for the heavy - jet mass one calculates the full ( np ) double differential distribution for the two jet masses and from that derives the heavy - jet mass distribution ( this is not quite equivalent to a simple shape - function for @xmath196 ) . this approach gave a good fit to the thrust and heavy - jet mass data , on the condition that there was significant correlation between the non - perturbative contribution to the two hemispheres .- parameter ( with an overall rescaling to ensure the correct first moment ) , however the data were less well described than for the thrust . this is perhaps not surprising : owing to different sensitivities to large - angle radiation , the thrust and @xmath38-parameter shape functions are expected to have a different structure of higher moments  @xcite . ] it is difficult to interpret this result however , because as has already been mentioned , the default measurement scheme for the jet masses suffers from extra ( universality - breaking ) @xmath146 corrections , and it is possible that in  @xcite these are being mimicked by inter - hemisphere correlations . the thrust and heavy - jet mass have been investigated in great detail also by gardi and rathsman  @xcite . in analogy with the approach for the mean thrust  @xcite , they attempt not only to account for nll+nlo and power correction contributions , but also for all ( perturbative ) enhanced effects associated with the running of the coupling , within a formalism that they refer to as dressed gluon exponentiation ( dge ) . among the interesting theoretical results ( see also  @xcite ) is that the renormalon analysis suggests that the second moment of the shape function should be precisely that expected from a pure shift ( the same is true for the @xmath38-parameter @xcite ) . this observation seems to be consistent with fits to the data , and would explain why the shift approximation works so well , at least for the relevant observables . the right - hand plot of fig . [ fig : thrustshape ] shows the result of their analysis for the thrust , with both a shift and a shape - function . while the shape - function is clearly needed in the immediate vicinity of the peak and below , it is striking to see , above the peak , how well it is approximated by the shift . one of the other interesting results in  @xcite relates to the analysis of the heavy - jet mass . in contrast to  @xcite , a common hadron - mass scheme ( decay - scheme ) is used for the thrust and the heavy - jet mass , and there are no inter - hemisphere correlations in their shape function . using the same parameters for the thrust and @xmath196 leads to reasonable fits , suggesting that the true inter - hemisphere correlations may actually be small . ( a fit for @xmath196 alone does however prefer a lower value of @xmath94 than for the thrust , as is seen also with the usual shifted nll+nlo approach @xcite ) . a puzzle that emerges in the work of  @xcite is that all the fits lead to quite small values of the coupling , @xmath336 with a theoretical error of about @xmath210 . a similarly low result was obtained in the analogous analysis of the mean thrust  @xcite . in contrast a study of @xmath144-fragmentation using dressed gluon exponentiation finds more standard values , @xmath337  @xcite . let us close our discussion of shape functions with a recent result by berger & sterman  @xcite . given that there are no simple relations between the shape functions for most common event shapes , they instead consider a new class of event shapes  @xcite whose definition involves a parameter @xmath331 which allows the nature of the event shape to be continuously varied from a broadening ( @xmath338 ) to a thrust - type observable ( @xmath339 ) and beyond . their observation is that for all variants of the observable with @xmath340 , the logarithm of any moment of the shape functions scales as @xmath341 . in practice this scaling is found to hold in monte carlo studies of this class of observables  it would therefore be interesting to have data for these observables with which to make a comparison . conceptually , event shapes are rather simple observables , yet they are sensitive to a range of characteristics of qcd radiation . this combination of theoretical simplicity and experimental sensitivity is perhaps one of the main reasons why event shapes have found so many applications . as we have seen , they provide vital inputs in studies of the ingredients of the qcd lagrangian , such as the coupling and the colour factors ; they play an important role in the tuning and testing of event generators ; and they are at the heart of recent investigations into analytical approaches to understanding the dynamics of hadronisation .    there remain several directions in which progress may be expected ( or at the very least hoped for ) in coming years . all the experimental comparisons discussed here have been for event - shapes that vanish in the @xmath16-jet limit , yet there have been significant theoretical developments in recent years for observables that vanish also in the @xmath65-jet limit , examples being the @xmath39-parameter and thrust minor @xcite . studies with jet rates @xcite suggest that going to such ` three - jet ' observables could reduce theoretical errors , though as discussed in @xcite this needs to be investigated more systematically . additionally three - jet observables will allow much more stringent tests of analytical hadronisation approaches since for the first time they introduce sensitivity to predictions about hadronisation from a gluon  @xcite , not just from a quark . one of the frontiers of event shape studies is their extension to new processes . while @xmath0 is the traditional domain for event - shape studies , we have seen that studies in dis are in many respects competitive with @xmath0 results . additionally they provide important input on the question of universality of analytical hadronisation models , ruling out for example significant modifications of the hadronisation picture stemming from interactions with the proton remnant . though so far the question has received only limited study , one interesting application of dis studies would be in the use of event shapes for obtaining information on parton distribution functions , as suggested in  @xcite . currently an option that is attracting growing experimental  @xcite and theoretical @xcite interest is the development of event shapes in hadron - hadron collisions . though much work remains to be done on this subject , it seems that the sensitivity to parton distributions will be somewhat stronger than in most dis observables . furthermore one expects a rich hadronisation structure associated both with the four - jet structure of hadronic dijet events and with the properties of the underlying event . it is to be hoped that the availability of automated resummation methods  @xcite will make the theoretical aspects of the study of new observables and processes more straightforward than has been the case in the past . another frontier on which progress is expected in near future is with respect to the accuracy of theoretical predictions . nnlo perturbative calculations are making rapid progress ( for a recent review see @xcite ) with the main outstanding problem being that of a full subtraction procedure for combining zero , one and two - loop contributions and its implementation in a monte carlo integrator . progress on nnll resummations is also being made , albeit so far only in inclusive cases  @xcite . though technically more involved , an extension to event shape resummations is certainly conceivable in the near future . aside from the expected gains in accuracy that are the main motivation for these calculations , it will also be interesting to compare the exact higher - order calculations with predictions from approaches such as  @xcite which aim to identify the physically dominant higher - order contributions and in some cases  @xcite claim to significantly reduce the theoretical uncertainties on @xmath94 . less predictable is what development can be expected on the subject of hadronisation corrections . as we have seen , much has been learnt from event - shape studies . yet a variety of open questions remain . at a practical level it would for example be of interest to see an extension of the shape - function approach to a wider range of variables , including the quite subtle case of the broadenings and also to dis event shapes . the wealth of dis data at low @xmath82 , where shape functions are most important , makes the latter especially interesting . there remain also important issues that are poorly understood even at the conceptual level . for example predictions are currently limited to the domain in which the event - shape is close to vanishing ( e.g.  @xmath14 in the @xmath16-jet limit ) . it is tempting to suggest that problems that persist for a couple of ` recalcitrant ' observables ( specifically , @xmath196 , @xmath301 ) may actually be due to large differences between the power correction in the three - jet region ( a significant part of the fit - range ) and that in the well - understood two - jet limit ( see e.g.   fig .  3 of  @xcite ) . but techniques that would allow the calculation of the power correction to , say , the thrust or the broadening in the @xmath65-jet limit do not yet exist . another ( partially related ) issue is that of anomalous dimensions  these have so far only been calculated for @xmath146 effects that are associated with hadron masses  @xcite . yet , as has been pointed out also in @xcite , they are bound to be present for all classes of power suppressed contributions . one clear physical origin for them is that soft - gluons ( or hadrons ) with transverse momenta of order @xmath254 are radiated not just from the @xmath206 system but also from all soft perturbative gluons with transverse momenta between @xmath254 and @xmath82 . accounting for radiation from the latter will clearly lead to enhancements of power corrections by terms @xmath342 . and a final point , not to be forgotten , is that most of the analytical hadronisation approaches are essentially perturbative methods ` in disguise ' . an understanding of how they really relate to _ _ hadron__isation is very much lacking . let us close this review by remarking that one of the characteristics of event shape studies is that progress is made not merely through a better understanding of existing observables and data , but also through critical experimental studies of new observables with ( supposedly ! ) better theoretical properties . the dialogue between experimenters and theorists is crucial in this respect . we look forward to its being as fruitful in coming years as it has been up to now . it goes without saying that we are grateful to our collaborators and numerous other colleagues , including members of the lep qcd working group , for the many stimulating discussions that we have had about event shapes . for discussions on issues that arose specifically while writing this review , we thank in particular g.  heinrich , f.  krauss , c.  j.  maxwell , g.  rodrigo , m.  h.  seymour , h.  stenzel and g.  zanderighi . furthermore we are indebted to a.  banfi , yu .  l.  dokshitzer and g.  zanderighi for their careful reading of and helpful comments on the manuscript . one of us ( gps ) is grateful also to the cern th division for hospitality during the course of this work . 999    j.  r.  ellis , m.  k.  gaillard and g.  g.  ross , nucl . b * 111 * ( 1976 ) 253 [ erratum - ibid . b * 130 * ( 1977 ) 516 ] . e.  farhi , phys . * 39 * ( 1977 ) 1587 . h.  georgi and m.  machacek , phys . * 39 * ( 1977 ) 1237 . g.  c.  fox and s.  wolfram , phys . * 41 * ( 1978 ) 1581 ; ibid . b * 149 * ( 1979 ) 413 [ erratum - ibid . b * 157 * ( 1979 ) 543 ] . c.  l.  basham , l.  s.  brown , s.  d.  ellis and s.  t.  love , phys . * 41 * ( 1978 ) 1585 ; c.  l.  basham , l.  s.  brown , s.  d.  ellis and s.  t.  love , phys . d * 19 * ( 1979 ) 2018 . r.  k.  ellis , w.  j.  stirling and b.  r.  webber , _ qcd and collider physics _ , cambridge university press , 1996 . g.  marchesini , b.  r.  webber , g.  abbiendi , i.  g.  knowles , m.  h.  seymour and l.  stanco , comput . commun . * 67 * ( 1992 ) 465 . g.  corcella _ et al . _ , jhep * 0101 * ( 2001 ) 010 [ hep - ph/0011363 ] . t.  sjstrand , comput . commun . * 82 * ( 1994 ) 74 ; t.  sjstrand , p.  eden , c.  friberg , l.  lnnblad , g.  miu , s.  mrenna and e.  norrbin , hep - ph/0010017 . l.  lnnblad , comput . * 71 * ( 1992 ) 15 . m.  beneke , phys . * 317 * ( 1999 ) 1 [ hep - ph/9807443 ] . s.  brandt , c.  peyrou , r.  sosnowski and a.  wroblewski , phys . * 12 * ( 1964 ) 57 . l.  clavelli , phys . b * 85 * ( 1979 ) 111 ; t.  chandramohan and l.  clavelli , nucl . b * 184 * ( 1981 ) 365 ; l.  clavelli and d.  wyler , phys . b * 103 * ( 1981 ) 383 . s.  catani , g.  turnock and b.  r.  webber , phys . b * 295 * ( 1992 ) 269 . d.  p.  barber _ et al . _ [ mark - j collaboration collaboration ] , phys . * 63 * ( 1980 ) 337 . a.  banfi , g.  marchesini , y.  l.  dokshitzer and g.  zanderighi , jhep * 0007 * ( 2000 ) 002 [ hep - ph/0004027 ] . a.  banfi , y.  l.  dokshitzer , g.  marchesini and g.  zanderighi , phys . b * 508 * ( 2001 ) 269 [ hep - ph/0010267 ] . a.  de rujula , j.  r.  ellis , e.  g.  floratos and m.  k.  gaillard , nucl . b * 138 * ( 1978 ) 387 . r.  k.  ellis , d.  a.  ross and a.  e.  terrano , nucl . b * 178 * ( 1981 ) 421 . g.  parisi , phys . b * 74 * ( 1978 ) 65 ; j.  f.  donoghue , f.  e.  low and s.  y.  pi , phys .  rev . d * 20 * ( 1979 ) 2759 . j.  d.  bjorken and s.  j.  brodsky , phys . d * 1 * ( 1970 ) 1416 ; g.  hanson _ et al . _ , phys . * 35 * ( 1975 ) 1609 . c.  berger _ [ pluto collaboration ] , phys . b * 82 * ( 1979 ) 449 . w.  bartel _ et al . _ [ jade collaboration ] , z.  phys . c * 33 * ( 1986 ) 23 . s.  catani , yu .  l.  dokshitzer , m.  olsson , g.  turnock and b.  r.  webber , phys . b * 269 * ( 1991 ) 432 . s.  bethke , z.  kunszt , d.  e.  soper and w.  j.  stirling , nucl . b * 370 * ( 1992 ) 310 [ erratum - ibid . b * 523 * ( 1998 ) 681 ] . y.  l.  dokshitzer , g.  d.  leder , s.  moretti and b.  r.  webber , jhep * 9708 * ( 1997 ) 001 [ hep - ph/9707323 ] . g.  p.  salam and d.  wicke , jhep * 0105 * ( 2001 ) 061 [ hep - ph/0102343 ] . e.  gardi and j.  rathsman , nucl . b * 638 * ( 2002 ) 243 [ hep - ph/0201019 ] . r.  p.   feynman , _ photon hadron interactions _ , w.  a.  benjamin , new york , 1972 ; r.  d.  peccei and r.  ruckl , nucl . b * 162 * ( 1980 ) 125 ; l.  v.  gribov , y.  l.  dokshitzer , s.  i.  troian and v.  a.  khoze , sov . jetp * 67 * ( 1988 ) 1303 [ zh . fiz .   * 94 * ( 1988 ) 12 ] . k.  h.  streng , t.  f.  walsh and p.  m.  zerwas , z.  phys . c * 2 * ( 1979 ) 237 . a.  banfi , g.  marchesini , g.  smye and g.  zanderighi , jhep * 0111 * ( 2001 ) 066 [ hep - ph/0111157 ] . m.  dasgupta and g.  p.  salam , jhep * 0208 * ( 2002 ) 032 [ hep - ph/0208073 ] . a.  banfi , g.  marchesini and g.  smye , jhep * 0204 * ( 2002 ) 024 [ hep - ph/0203150 ] . s.  catani , y.  l.  dokshitzer and b.  r.  webber , phys . b * 285 * ( 1992 ) 291 . a.  banfi , g.  marchesini , g.  smye and g.  zanderighi , jhep * 0108 * ( 2001 ) 047 [ hep - ph/0106278 ] . i.  a.  bertram [ d0 collaboration ] , acta phys . b * 33 * ( 2002 ) 3141 . a.  banfi , g.  p.  salam and g.  zanderighi , hep - ph/0304148 ; _ ibid . _  in preparation . a.  banfi , hep - ph/0309161 . g.  sterman and s.  weinberg , phys . * 39 * ( 1977 ) 1436 . w.  t.  giele and e.  w.  glover , phys . d * 46 * ( 1992 ) 1980 . s.  catani and m.  h.  seymour , phys . b * 378 * ( 1996 ) 287 [ hep - ph/9602277 ] ; s.  catani and m.  h.  seymour , nucl . b * 485 * ( 1997 ) 291 [ erratum - imbid . b * 510 * ( 1997 ) 503 ] [ hep - ph/9605323 ] . l.  j.  dixon and a.  signer , phys . d * 56 * ( 1997 ) 4031 [ hep - ph/9706285 ] ; a.  signer , comput . commun . * 106 * ( 1997 ) 125 . s.  weinzierl and d.  a.  kosower , phys . d * 60 * ( 1999 ) 054028 [ hep - ph/9901277 ] . j.  m.  campbell , m.  a.  cullen and e.  w.  glover , eur . j.  c * 9 * ( 1999 ) 245 [ hep - ph/9809429 ] . z.  nagy and z.  trocsanyi , nucl . b * 486 * ( 1997 ) 189 [ hep - ph/9610498 ] ; z.  nagy and z.  trocsanyi , phys . * 79 * ( 1997 ) 3604 [ hep - ph/9707309 ] ; z.  nagy and z.  trocsanyi , phys . d * 59 * ( 1999 ) 014020 [ erratum - ibid . d * 62 * ( 2000 ) 099902 ] [ hep - ph/9806317 ] ; z.  nagy and z.  trocsanyi , phys . lett .   * 87 * ( 2001 ) 082001 [ hep - ph/0104315 ] ; z.  nagy , phys . * 88 * ( 2002 ) 122003 [ hep - ph/0110315 ] ; z.  nagy , hep - ph/0307268 . b.  potter , comput .  phys . commun . * 133 * ( 2000 ) 105 [ hep - ph/9911221 ] ; b.  potter , comput . commun . * 119 * ( 1999 ) 45 [ hep - ph/9806437 ] . t.  gehrmann , hep - ph/0310178 . m.  dasgupta and g.  p.  salam , phys . b * 512 * ( 2001 ) 323 [ hep - ph/0104277 ] . n.  brown and w.  j.  stirling , phys . b * 252 * ( 1990 ) 657 . g.  leder , nucl . b * 497 * ( 1997 ) 334 [ hep - ph/9610552 ] . s.  catani , g.  turnock , b.  r.  webber and l.  trentadue , phys . b * 263 * ( 1991 ) 491 . s.  catani , l.  trentadue , g.  turnock and b.  r.  webber , nucl . b * 407 * ( 1993 ) 3 . s.  catani , g.  turnock and b.  r.  webber , phys . b * 272 * ( 1991 ) 368 . s.  j.  burby and e.  w.  glover , jhep * 0104 * ( 2001 ) 029 [ hep - ph/0101226 ] . y.  l.  dokshitzer , a.  lucenti , g.  marchesini and g.  p.  salam , jhep * 9801 * ( 1998 ) 011 [ hep - ph/9801324 ] . s.  catani and b.  r.  webber , phys . b * 427 * ( 1998 ) 377 [ hep - ph/9801350 ] . g.  dissertori and m.  schmelling , phys . b * 361 * ( 1995 ) 167 . a.  banfi , g.  p.  salam and g.  zanderighi , jhep * 0201 * ( 2002 ) 018 [ hep - ph/0112156 ] . y.  l.  dokshitzer and d.  diakonov , phys . b * 84 * ( 1979 ) 234 . g.  parisi and r.  petronzio , nucl . b * 154 * ( 1979 ) 427 . j.  kodaira and l.  trentadue , phys . b * 123 * ( 1983 ) 335 . j.  c.  collins and d.  e.  soper , phys . * 48 * ( 1982 ) 655 . p.  e.  rakow and b.  r.  webber , nucl . b * 187 * ( 1981 ) 254 . y.  l.  dokshitzer , g.  marchesini and b.  r.  webber , jhep * 9907 * ( 1999 ) 012 [ hep - ph/9905339 ] . c.  f.  berger , t.  kucs and g.  sterman , phys . d * 68 * ( 2003 ) 014012 [ hep - ph/0303051 ] . a.  banfi , y.  l.  dokshitzer , g.  marchesini and g.  zanderighi , jhep * 0105 * ( 2001 ) 040 [ hep - ph/0104162 ] . v.  antonelli , m.  dasgupta and g.  p.  salam , jhep * 0002 * ( 2000 ) 001 [ hep - ph/9912488 ] . m.  dasgupta and g.  p.  salam , eur . j.  c * 24 * ( 2002 ) 213 [ hep - ph/0110213 ] . r.  bonciani , s.  catani , m.  l.  mangano and p.  nason , phys . b * 575 * ( 2003 ) 268 [ hep - ph/0307035 ] . e.  gardi and j.  rathsman , nucl . b * 609 * ( 2001 ) 123 [ hep - ph/0103217 ] . e.  gardi and l.  magnea , `` the c parameter distribution in e+ e- annihilation , '' jhep * 0308 * ( 2003 ) 030 [ hep - ph/0306094 ] . f.  krauss and g.  rodrigo , phys . b * 576 * ( 2003 ) 135 [ hep - ph/0303038 ] . s.  catani and b.  r.  webber , jhep * 9710 * ( 1997 ) 005 [ hep - ph/9710333 ] . m.  dasgupta and g.  p.  salam , jhep * 0203 * ( 2002 ) 017 [ hep - ph/0203009 ] ; c.  f.  berger , t.  kucs and g.  sterman , phys . d * 65 * ( 2002 ) 094031 [ hep - ph/0110004 ] ; a.  banfi , g.  marchesini and g.  smye , jhep * 0208 * ( 2002 ) 006 [ hep - ph/0206076 ] . r.  b.  appleby and m.  h.  seymour , jhep * 0212 * ( 2002 ) 063 [ hep - ph/0211426 ] ; r.  b.  appleby and m.  h.  seymour , jhep * 0309 * ( 2003 ) 056 [ hep - ph/0308086 ] . y.  l.  dokshitzer and g.  marchesini , jhep * 0303 * ( 2003 ) 040 [ hep - ph/0303101 ] . m.  dasgupta , hep - ph/0304086 ; r.  b.  appleby and g.  p.  salam , hep - ph/0305232 . h.  weigert , hep - ph/0312050 .    i.  h.  park _ et al . _ ( amy collaboration ) , phys . lett .  * 62 * ( 1989 ) 1713 ; + y.  k.  li _ et al . _ ( amy collaboration ) , phys . * d 41 * ( 1990 ) 2675 ; + h.  j.  behrend _ et al . _ ( cello collaboration ) , z.  phys . * c44 * ( 1989 ) 63 ; + d.  bender _ et al . _ ( hrs collaboration ) , phys . * d 31 * ( 1985 ) 1 ; + d.  p.  barber _ et al . _ ( mark j collaboration ) , phys . * 43 * ( 1979 ) 830 ; + a.  petersen _ et al . _ ( mark ii collaboration ) , phys . d 37 * ( 1988 ) 1 ; + s.  bethke _ et al . _ ( mark ii collaboration ) , z.  phys . * c43 * ( 1989 ) 325 ; + c.  berger _ ( pluto collaboration ) , z.  phys . * c12 * ( 1982 ) 297 ; + w.  braunschweig _ et al . _ ( tasso collaboration ) , phys . * b214 * ( 1988 ) 286 ; + w.  braunschweig _ et al . _ ( tasso collaboration ) , z.  phys . * c45 * ( 1989 ) 11 ; + w.  braunschweig _ et al . _ ( tasso collaboration ) , z.  phys .  * c47 * ( 1990 ) 187 ; + i.  adachi _ et al . _ ( topaz collaboration ) , phys . * b227 * ( 1989 ) 495 ; + k.  nagai _ et al . _ ( topaz collaboration ) , phys . * b278 * ( 1992 ) 506 ; + y.  ohnishi _ et al . _ ( topaz collaboration ) , phys . * b313 * ( 1993 ) 475 . d.  decamp _ et al . _ [ aleph collaboration ] , phys . b * 255 * , 623 ( 1991 ) . d.  decamp _ et al . _ [ aleph collaboration ] , phys . b * 284 * , 163 ( 1992 ) . d.  buskulic _ et al . _ [ aleph collaboration ] , z.  phys . c * 73 * ( 1997 ) 409 . r.  barate _ [ aleph collaboration ] , phys . * 294 * ( 1998 ) 1 . aleph collaboration , `` qcd measurements in @xmath0 annihilations at centre - of - mass energies between 189 and 202  gev '' , aleph 2000 - 044 conf 2000 - 027 . p.  abreu _ et al . _ [ delphi collaboration ] , z.  phys . c * 73 * , 229 ( 1997 ) . p.  abreu _ et al . _ [ delphi collaboration ] , phys . b * 456 * , 322 ( 1999 ) . p.  abreu _ et al . _ [ delphi collaboration ] , eur . j.  c * 14 * , 557 ( 2000 ) [ hep - ex/0002026 ] . j.  abdallah _ et al . _ [ delphi collaboration ] , eur . j.  c * 29 * ( 2003 ) 285 [ hep - ex/0307048 ] . o.  adriani _ et al . _ [ l3 collaboration ] , phys . b * 284 * , 471 ( 1992 ) . b.  adeva _ et al . _ [ l3 collaboration ] , z.  phys . c * 55 * ( 1992 ) 39 . o.  adriani _ et al . _ [ l3 collaboration ] , phys . rept .   * 236 * , 1 ( 1993 ) . m.  acciarri _ et al . _ [ l3 collaboration ] , phys . b * 411 * , 339 ( 1997 ) . m.  acciarri _ et al . _ [ l3 collaboration ] , phys . b * 371 * , 137 ( 1996 ) . m.  acciarri _ et al . _ [ l3 collaboration ] , phys . b * 404 * , 390 ( 1997 ) . m.  acciarri _ et al . _ [ l3 collaboration ] , phys . b * 444 * , 569 ( 1998 ) . m.  acciarri _ et al . _ [ l3 collaboration ] , phys . b * 489 * , 65 ( 2000 ) [ hep - ex/0005045 ] . p.  achard _ et al . _ [ l3 collaboration ] , phys . b * 536 * ( 2002 ) 217 [ hep - ex/0206052 ] .    p.  d.  acton _ et al . _ [ opal collaboration ] , z.  phys . c * 55 * ( 1992 ) 1 . d.  acton _ et al . _ [ opal collaboration ] , z.  phys . c * 59 * , 1 ( 1993 ) . g.  alexander _ et al . _ [ opal collaboration ] , z.  phys . c * 72 * , 191 ( 1996 ) . k.  ackerstaff _ et al . _ [ opal collaboration ] , z.  phys . c * 75 * , 193 ( 1997 ) . g.  abbiendi _ et al . _ [ opal collaboration ] , eur . j.  c * 16 * ( 2000 ) 185 [ hep - ex/0002012 ] . p.  a.  movilla fernandez , o.  biebel , s.  bethke , s.  kluth and p.  pfeifenschneider [ jade collaboration ] , eur . j.  c * 1 * , 461 ( 1998 ) [ hep - ex/9708034 ] . o.  biebel , p.  a.  movilla fernandez and s.  bethke [ jade collaboration ] , phys . b * 459 * , 326 ( 1999 ) [ hep - ex/9903009 ] . p.  pfeifenschneider _ et al . _ [ jade collaboration ] , eur . j.  c * 17 * , 19 ( 2000 ) [ hep - ex/0001055 ] . p.  a.  movilla fernandez , s.  bethke , o.  biebel and s.  kluth , eur . j.  c * 22 * ( 2001 ) 1 [ hep - ex/0105059 ] . s.  kluth , m.  blumenstengel , p.  a.  movilla fernandez , s.  bethke , o.  biebel , c.  pahl and j.  schieck [ jade collaboration ] , hep - ex/0305023 . j.  schieck [ for the opal collaboration ] , hep - ex/0312036 . f.  yuasa _ et al . _ , prog . * 138 * ( 2000 ) 18 [ hep - ph/0007053 ] . a.  pukhov _ et al . _ , hep - ph/9908288 . f.  krauss , r.  kuhn and g.  soff , jhep * 0202 * ( 2002 ) 044 [ hep - ph/0109036 ] . c.  adloff _ et al . _ [ h1 collaboration ] , phys . b * 406 * ( 1997 ) 256 [ hep - ex/9706002 ] . c.  adloff _ et al . _ [ h1 collaboration ] , eur . j.  c * 14 * ( 2000 ) 255 [ addendum - ibid . c * 18 * ( 2000 ) 417 ] [ hep - ex/9912052 ] . h.  u.  martyn [ h1 collaboration ] , hep - ex/0010046 . s.  chekanov _ et al . _ [ zeus collaboration ] , eur . j.  c * 27 * ( 2003 ) 531 [ hep - ex/0211040 ] . j.  breitweg _ et al . _ [ zeus collaboration ] , phys . b * 421 * ( 1998 ) 368 [ hep - ex/9710027 ] . b.  r.  webber , nucl . * 71 * ( 1999 ) 66 [ hep - ph/9712236 ] . a.  v.  manohar and m.  b.  wise , phys . b * 344 * ( 1995 ) 407 [ hep - ph/9406392 ] . g.  p.  korchemsky and g.  sterman , nucl . b * 437 * ( 1995 ) 415 [ hep - ph/9411211 ] . g.  p.  korchemsky and g.  sterman , hep - ph/9505391 . g.  p.  korchemsky and g.  sterman , nucl . b * 555 * ( 1999 ) 335 [ hep - ph/9902341 ] . g.  p.  korchemsky and s.  tafat , jhep * 0010 * ( 2000 ) 010 [ hep - ph/0007005 ] . b.  r.  webber , phys . b * 339 * ( 1994 ) 148 [ hep - ph/9408222 ] . y.  l.  dokshitzer and b.  r.  webber , phys . b * 352 * ( 1995 ) 451 [ hep - ph/9504219 ] . y.  l.  dokshitzer , g.  marchesini and b.  r.  webber , nucl . b * 469 * ( 1996 ) 93 [ hep - ph/9512336 ] . r.  akhoury and v.  i.  zakharov , phys . b * 357 * ( 1995 ) 646 [ hep - ph/9504248 ] . y.  l.  dokshitzer , a.  lucenti , g.  marchesini and g.  p.  salam , nucl . b * 511 * ( 1998 ) 396 [ erratum - ibid . b * 593 * ( 2001 ) 729 ] [ hep - ph/9707532 ] . y.  l.  dokshitzer , g.  marchesini and g.  p.  salam , eur . j.  directc * 1 * ( 1999 ) 3 [ hep - ph/9812487 ] . m.  dasgupta , l.  magnea and g.  smye , jhep * 9911 * ( 1999 ) 025 [ hep - ph/9911316 ] . g.  e.  smye , jhep * 0105 * ( 2001 ) 005 [ hep - ph/0101323 ] . e.  gardi and g.  grunberg , jhep * 9911 * ( 1999 ) 016 [ hep - ph/9908458 ] . e.  gardi , jhep * 0004 * ( 2000 ) 030 [ hep - ph/0003179 ] . c.  w.  bauer , a.  v.  manohar and m.  b.  wise , phys . * 91 * ( 2003 ) 122001 [ hep - ph/0212255 ] . z.  trocsanyi , jhep * 0001 * ( 2000 ) 014 [ hep - ph/9911353 ] . m.  dasgupta and b.  r.  webber , eur . j.  c * 1 * ( 1998 ) 539 [ hep - ph/9704297 ] . m.  dasgupta and b.  r.  webber , jhep * 9810 * ( 1998 ) 001 [ hep - ph/9809247 ] . g.  p.  salam and z.  trcsnyi , _ limited automated predictor of power suppressed effects ( lappse ) _ , unpublished work . b.  r.  webber , hep - ph/9510283 . y.  l.  dokshitzer , `` why it seems too early to report on status of qcd , '' _ given at international school of subnuclear physics : 31th course : from supersymmetry to the origin of space - time , erice , italy , 4 - 12 jul 1993 _ j.  m.  campbell , e.  w.  n.  glover and c.  j.  maxwell , phys . * 81 * ( 1998 ) 1568 [ hep - ph/9803254 ] . d.  t.  barclay , c.  j.  maxwell and m.  t.  reader , phys . d * 49 * ( 1994 ) 3480 . m.  beneke , phys . b * 307 * ( 1993 ) 154 . p.  abreu _ et al . _ ( delphi collaboration ) , z.  phys .  * c73 * ( 1996 ) 11 . d.  buskulic _ et al . _ ( aleph collaboration ) , z.  phys . * c55 * ( 1992 ) 209 . i.  g.  knowles _ , `` qcd event generators , '' hep - ph/9601212 . s.  frixione and b.  r.  webber , jhep * 0206 * ( 2002 ) 029 [ hep - ph/0204244 ] . s.  frixione , p.  nason and b.  r.  webber , jhep * 0308 * ( 2003 ) 007 [ hep - ph/0305252 ] . s.  gieseke , a.  ribon , m.  h.  seymour , p.  stephens and b.  webber , hep - ph/0311208 . lep qcd working group , work in progress . h.  stenzel , hep - ex/0309061 . s.  bethke , nucl . * 121 * ( 2003 ) 74 [ hep - ex/0211012 ] . k.  hagiwara _ et al . _ [ particle data group collaboration ] , phys .  rev . d * 66 * ( 2002 ) 010001 . r.  w.  l.  jones , m.  ford , g.  p.  salam , h.  stenzel and d.  wicke , jhep * 0312 * ( 2003 ) 007 [ hep - ph/0312016 ] . p.  m.  stevenson , phys . d * 23 * ( 1981 ) 2916 . g.  grunberg , phys . d * 29 * ( 1984 ) 2315 . s.  j.  brodsky , g.  p.  lepage and p.  b.  mackenzie , phys . d * 28 * ( 1983 ) 228 . p.  abreu _ et al . _ [ delphi collaboration ] , phys . b * 418 * ( 1998 ) 430 . g.  abbiendi _ et al . _ [ opal collaboration ] , eur . j.  c * 21 * ( 2001 ) 411 [ hep - ex/0105046 ] . a.  brandenburg , p.  n.  burrows , d.  muller , n.  oishi and p.  uwer , phys . b * 468 * ( 1999 ) 168 [ hep - ph/9905495 ] . r.  barate _ [ aleph collaboration ] , eur . j.  c * 18 * ( 2000 ) 1 [ hep - ex/0008013 ] . p.  bambade , m.  j.  coster , j.  fuster and p.  tortosa , delphi conference note 2003 - 022-conf-642 ( unpublished ) . w.  bernreuther , a.  brandenburg and p.  uwer , phys . * 79 * ( 1997 ) 189 [ hep - ph/9703305 ] . g.  rodrigo , a.  santamaria and m.  s.  bilenky , phys . * 79 * ( 1997 ) 193 [ hep - ph/9703358 ] . p.  nason and c.  oleari , nucl . b * 521 * ( 1998 ) 237 [ hep - ph/9709360 ] . p.  a.  movilla fernandez , hep - ex/0209022 . t.  kluge , hep - ex/0310063 . g.  p.  salam and g.  zanderighi , nucl . * 86 * ( 2000 ) 430 [ hep - ph/9909324 ] . a.  heister _ et al . _ [ aleph collaboration ] , eur . j.  c * 27 * ( 2003 ) 1 . g.  abbiendi _ et al . _ [ opal collaboration ] , eur . j.  c * 20 * ( 2001 ) 601 [ hep - ex/0101044 ] . p.  abreu _ et al . _ [ delphi collaboration ] , eur . j.  c * 13 * ( 2000 ) 573 . g.  abbiendi _ et al . _ [ opal collaboration ] , eur . j.  c * 23 * ( 2002 ) 597 [ hep - ex/0111013 ] . g.  p.  korchemsky , hep - ph/9806537 . g.  sterman , nucl . b * 281 * ( 1987 ) 310 ; + s.  catani and l.  trentadue , nucl . b * 327 * ( 1989 ) 323 ; + s.  catani and l.  trentadue , nucl . b * 353 * ( 1991 ) 183 ; + g.  p.  korchemsky , mod . lett .  a * 4 * ( 1989 ) 1257 ; + g.  p.  korchemsky and g.  marchesini , phys . b * 313 * ( 1993 ) 433 ; + i.  i.  y.  bigi , m.  a.  shifman , n.  g.  uraltsev and a.  i.  vainshtein , int . j.  mod . phys .  a * 9 * ( 1994 ) 2467 [ hep - ph/9312359 ] ; + r.  d.  dikeman , m.  a.  shifman and n.  g.  uraltsev , int .  j.  mod . a * 11 * ( 1996 ) 571 [ hep - ph/9505397 ] ; + m.  neubert , phys . d * 49 * ( 1994 ) 3392 [ hep - ph/9311325 ] ; + m.  neubert , phys . d * 49 * ( 1994 ) 4623 [ hep - ph/9312311 ] ; + g.  p.  korchemsky and g.  sterman , phys . b * 340 * ( 1994 ) 96 [ hep - ph/9407344 ] ; + a.  g.  grozin and g.  p.  korchemsky , phys . d * 53 * ( 1996 ) 1378 [ hep - ph/9411323 ] . m.  cacciari and e.  gardi , nucl . b * 664 * ( 2003 ) 299 [ hep - ph/0301047 ] . g.  bozzi , s.  catani , d.  de florian and m.  grazzini , phys . b * 564 * ( 2003 ) 65 [ hep - ph/0302104 ] . l.  magnea , hep - ph/0211013 .
Abstract:

Answer: this article reviews the status of event - shape studies in @xmath0 annihilation and dis . it includes discussions of perturbative calculations , of various approaches to modelling hadronisation and of comparisons to data . cern - th/2003 - 306 + lpthe-03 - 40 + hep - ph/0312283