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: [ sec : intro ]    an impressive amount of data on backward hadron production in p+c interactions has been collected over the past four decades . a literature survey reveals no less than 19 experiments which have contributed a total amount of more than 3500 data points covering wide areas in projectile momentum , lab angle and lab momentum . looking at the physics motivation and at the distribution in time of these efforts , two distinct classes of experimental approaches become evident . 15 experiments cluster in a first period during the two decades between 1970 and 1990 . all these measurements have been motivated by the nuclear part of proton - nucleus collisions , in particular by the width of the momentum distributions in the nuclear rest system which reach far beyond the narrow limits expected from nuclear binding alone . these studies have ceased in the late 1980 s with the advent of relativistic heavy ion collisions and their promise of `` new '' phenomena beyond the realm of classic nuclear physics .    a second class of very recent measurements has appeared and is being pursued after the turn of the century , with publications starting about 2008 . here the motivation is totally different . it is driven by the necessity of obtaining hadronic reference data for the study of systematic effects in cosmic ray and neutrino physics , in particular concerning atmospheric and long base line experiments as well as eventual novel neutrino factories . the main aim of these studies is the comparison to and the improvement of hadronic production models  models which are to be considered as multi - parameter descriptions of the non - calculable sector of the strong interaction , with very limited predictive power . this new and exclusive aim has led to the strange situation that if all recent publications contain detailed comparisons to available production models , no comparison to existing data is attempted . it remains therefore unclear how these new results compare to the wealth of already available data and whether they in fact may over - ride and replace the existing results . in this environment the studies conducted since 15 years by the na49 experiment at the cern sps have a completely different aim . here it is attempted to trace a model - independent way from the basic hadron - nucleon interaction via hadron - nucleus to nucleus - nucleus collisions . this aim needs precision data from a large variety of projectile and target combinations as well as a maximum phase space coverage . as the acceptance of the na49 detector is limited to lab angles below 45 degrees , it is indicated to use existing backward data in the sps energy range in order to extend the acceptance coverage for the asymmetric proton - nucleus interactions . this requires a careful study of the dependence on cms energy and of the reliability of the results to be used .    in the course of this work it appeared useful and even mandatory to provide a survey of all available data over the full scale of interaction energies , the more so as no overview of the experimental situation is available to date . this means that the present study deals with projectile momenta from 1 to 400  gev / c , for a lab angle range from 10 to 180  degrees , and for lab momenta from 0.2 to 1.2  gev / c . [ sec : var_kin ]    most available data have been obtained as a function of the lab momentum @xmath8 ( or kinetic energy @xmath9 ) at constant lab angle @xmath10 . in this publication all given yields are transformed to the double differential invariant cross section    @xmath11    in this context the term `` backward '' needs a precise definition . one possibility would be to define as `` backward '' the region of lab angles @xmath12  90  degrees . the present paper uses instead a definition which refers to the cms frame with the basic variables feynman @xmath13 and transverse momentum @xmath14 , defining as `` backward '' the particle yields at @xmath15  0 . this allows a clear separation of the projectile fragmentation region at positive @xmath13 with a limited feed - over into negative @xmath13 and the target fragmentation region at negative @xmath13 with a limited feed - over into positive @xmath13 . at the same time the notion of `` kinematic limit '' in participant fragmentation is clearly brought out at @xmath13  =  @xmath16  1 and the contributions from intranuclear cascading may be clearly visualized and eventually separated . and @xmath10 in the cms frame spanned by feynman @xmath13 and @xmath14 for protons and pions at two different projectile momenta , a ) protons at 158  gev / c , b ) protons at 3  gev / c , c ) pions at 158  gev / c and d ) pions at 3  gev / c , width=510 ]    the correlation between the two pairs of variables is presented in fig .  [ fig : kinematics ] which shows lines of constant @xmath8 and @xmath10 in the coordinate frame of @xmath13 and @xmath14 for protons and pions for the two values of projectile momentum at 158 and 3  gev / c which are representative of the typical range of interaction energies discussed in this paper . several comments are due in this context . the definition of feynman @xmath13 has been modified from the standard one ,    @xmath17    to    @xmath18    with @xmath19 the proton mass . this takes care of baryon number conservation and regularises the kinematic borders at low interaction energies . the @xmath20 dependence in fig . [ fig : kinematics ] is small to negligible for lab angles above about 50  degrees both for pions and protons but becomes noticeable at small @xmath10 . if at sps energy the full range of lab momenta up to 1.4  gev / c and angles above 10  degrees is confined to the backward region both for protons and pions , the coverage for pions extends to positive @xmath13 at low lab angles and low beam momenta . another remark concerns the overlap between target fragmentation and nuclear cascading . for protons , at all lab angles above about 70  degrees the kinematic limit for fragmentation of a target nucleon at rest in the lab system is exceeded . for pions on the other hand this is not the case as their @xmath13 value for @xmath8  =  0 is at @xmath21    this means that over the full range of lab angles and up to large @xmath8 values the contribution from target participants mixes with the nuclear component . the separation of the two processes therefore becomes an important task , see sect . [ sec : separation ] of this paper . a last remark is due to the limits of experimental coverage . all existing experiments run out of statistics at cross section levels of about 10  @xmath22b , that is about 4 orders of magnitude below the maximum yields . as visible from the momentum ranges indicated in tables 1 and 2 , this corresponds to a typical upper momentum cut - off in the region of 1  gev / c . [ sec : exp_sit ]    the backward phase space coverage in p+c interactions is surprisingly complete if compared with the forward direction and even with the available data in the elementary p+p collisions . this is apparent from the list of experiments given in tables  [ tab : prot_data ] and [ tab : pion_data ] with their ranges in beam momentum , lab angle , and lab momentum . although some effort has been spent to pick up all published results , this list is not claimed to be exhaustive as some results given as `` private communication '' , in conference proceedings or unpublished internal reports might have escaped attention . .data sets for proton production in p+c and n+c collisions from seven experiments giving the ranges covered in projectile momentum , lab angle , and lab momentum , the number of measured data points and errors [ cols="^,^,^,^,^,^,^,^ " , ]     for secondary protons , table  [ tab : prot_data ] , the important amount of low energy n+c data by franz et al . @xcite has been added to the survey as the isospin factors for the transformation into p+c results have been studied and determined with some precision , see sect . [ sec : prot ] .    for secondary pions , table  [ tab : pion_data ] , the situation is somewhat complicated by the fact that two independent sets of results have been published by the harp - cdp @xcite and the harp @xcite groups , based on identical input data obtained with the same detector . an attempt to clarify this partially contradictory situation is presented in sect . [ sec : harp ] of this paper . unfortunately , no commonly agreed scale in the three basic variables @xmath10 , @xmath8 and @xmath23 of the double - differential cross sections has been defined by the different collaborations providing the data contained in tables  [ tab : prot_data ] and [ tab : pion_data ] . this leads to the fact that not a single couple out of the more than 3500 data points contained in these tables may be directly compared . the application of an interpolation scheme as described in sect . [ sec : comp ] is therefore an absolute necessity . ideally the thus obtained interpolated cross sections would form an internally consistent sample of results which would be coherent within the given experimental errors . as will become apparent in the following data comparison , this assumption is surprisingly well fulfilled for the majority of the experiments . only four of the 20 quoted groups of results fall significantly out of this comparison ; those will be discussed in sect . [ sec : dev_exp ] of this paper . in this sense the overall survey of the backward proton and pion production results in a powerful constraint for the comparison with any new data sample . [ sec : comp ]    as stated above the main problem in bringing the wealth of available data into a consistent picture is given by the generally disparate position in phase space and interaction energy of the different experiments . the triplet of lab variables given by the beam momentum @xmath23 , the lab momentum @xmath8 and the lab angle @xmath10 has been used for the establishment of the following interpolation scheme . in addition and of course , the statistical and systematic errors have to be taken into account in the data comparison . [ sec : errors ]    the last columns of tables  [ tab : prot_data ] and [ tab : pion_data ] contain some information about the statistical and systematic errors of the different experiments . the given numbers are to be regarded as mean values excluding some upward tails as they are inevitable at the limits of the covered phase space in particular for the statistical uncertainties . in some cases only rudimentary information about the systematic errors is available or the systematic and statistical errors are even combined into one quantity . in the latter cases these values are given in between the respective columns of table  [ tab : pion_data ] . inspection of these approximate error levels reveals a rather broad band of uncertainties ranging from about 4% to about 20% , the latter limit being generally defined by overall normalization errors . the presence of extensive data sets well below the 10% range of both statistical and systematic errors gives however some hope that a resulting overall consistency on this level might become attainable by the extensive use of data interpolation . the term `` interpolation '' is to be regarded in this context as a smooth interconnection of the data points in any of the three phase space variables defined above . this interconnection is generally done by eyeball fits which offer , within the error limits shown above , sufficient accuracy . if the distributions in @xmath10 and interaction energy are anyway not describable by straight - forward arithmetic parametrization , the @xmath8 dependences are , as discussed in sect . [ sec : pldep ] below , in a majority of cases approximately exponential . in these cases exponential fits have been used if applicable .    as additional constraint physics asks of course for smoothness and continuity in all three variables simultaneously . therefore the resulting overall data interpolation has to attempt a three - dimensional consistency . if the data interpolation helps , by the inter - correlation of data points , to reduce the local statistical fluctuations , it does of course not reduce the systematic uncertainties . it is rather on the level of systematic deviations that the consistency of different experimental results is to be judged . it will become apparent from the detailed discussion described below that the majority of the quoted experiments allows for the establishment of a surprisingly consistent overall description in all three variables . [ sec : sdep ]    as the data discussed here span an extremely wide range of cms energy from close to production threshold to the upper range of fermilab energies , a suitable compression of the energy scale has been introduced in order to be able to present the results in a close - to - equidistant fashion against energy . the form chosen here is the variable @xmath24 . this choice is suggested by the considerable amount of work invested in studying the approach of hadronic cross sections to the scaling limit at high energy in the 1970 s @xcite . in fact the regge parametrization suggested a smooth dependence of the cross sections as @xmath25 , with @xmath26  =  0.250.5 depending on the choice of trajectories involved . such behaviour was indeed found experimentally . in the present study the cross sections turn out to have only a mild @xmath24 dependence for @xmath27  5  gev , a dependence which is however different for pions and protons . this dependence is strongly modified below @xmath28  2.5  gev due to threshold effects . [ sec : thdep ]    a convenient and often used scale for the lab angle dependence is given by @xmath29 . this scale has the advantage of producing shapes that are again to zero order exponential . of course , continuity through @xmath10  =  180  degrees imposes an approach to 180  degrees with tangent zero . as the data samples are generally not measured at common values of @xmath10 , a fixed grid of angles has been defined based on the @xmath10 values of the harp - cdp experiment @xcite dominating the range from 25 to 112  degrees . measured values down to 10  degrees and in the higher angular range at 137 , 160 , and 180  degrees have been added . measurements not corresponding to these grid values are interpolated using the @xmath29 distributions specified below . [ sec : pldep ]    all data discussed here have been transformed into invariant cross sections ( [ eq : ddcs ] ) . this facilitates the presentation in different coordinate systems and eliminates the trivial approach of the phase space element to zero with decreasing momentum . in addition , most of the invariant @xmath8 distributions are close to exponential within the measured @xmath8 range . there are notable deviations mostly at low momentum and in the lower ( higher ) range of lab angles for pions and protons , respectively , as well as in the approach to threshold . in these cases an eyeball fit has been used which can be reliably performed within the error margins indicated above .    at low lab momenta physics requires a deviation from the exponential shape as the invariant cross sections must approach @xmath8  =  0 with tangent zero . this limit appears in general at @xmath30  0.2  gev / c for pions and @xmath30  0.5  gev / c for protons . the data presented here fall practically all above these momentum limits . only the harp experiment @xcite gives results at @xmath8  =  0.125  gev / c for pions where indeed a substantial deviation from the exponential shape is visible . this is shown in fig . [ fig : first_point ] where the deviation from exponential fits at this @xmath8 is given in percent for all angles and beam momenta together with similar deviations observed in p+p interactions @xcite . =  0.125  gev / c from the exponential fits for @xmath31 and @xmath32 at all angles and beam momenta , b ) deviation of @xmath31 cross sections at @xmath14  =  0.125  gev / c from exponential fits to the higher @xmath14 region in p+p interactions as a function of @xmath13,width=453 ]    a number of examples of momentum distributions for protons and pions is given in the following figs . [ fig : pldist_harpcdp ] and [ fig : plabex ] which show the invariant cross sections as a function of @xmath8 and the corresponding exponential fits    @xmath33    which are , whenever necessary , supplemented by hand interpolations into the non - exponential regions . a first group of distributions in the medium angular range at 45 and 97  degrees is presented in fig .  [ fig : pldist_harpcdp ] for the harp - cdp data concerning protons and pions , including exponential fits .     and @xmath32 as a function of @xmath8 at @xmath10  =  45 and 97  degrees . full lines : exponential fits . broken line : hand - interpolation into the non - exponential region , width=529 ]    evidently the exponential shape is within errors in general a good approximation to the momentum dependence . more quantitative information is contained in the normalized residual distributions of the data points ,    @xmath34    where @xmath35 is the difference between data and fit and sigma the statistical error of the given data point . should the fit describe the physics and should systematic effects be negligible , the distribution of @xmath36 is expected to be gaussian with rms equal to unity . the @xmath36 distributions are given for the totality of the harp - cdp data in fig . [ fig : diff_dist ] . these distributions are well described by centred gaussians . the resulting rms values are however somewhat bigger than one signalling systematic experimental effects or a deviation of physics from the simple exponential parametrization . in view of the statistical errors of 4% to 6% given by harp - cdp ( tables  [ tab : prot_data ] and [ tab : pion_data ] ) these deviations are on the level of a few percent which is below the error margin to be anyway expected from the present general data survey . further examples of @xmath8 distributions from other experiments are given in fig . [ fig : plabex ] for a selection of particle type , beam momenta and angles .     for a variety of particle type , lab angle and beam momenta including exponential fits ( full lines ) and , when necessary , eyeball fits into the non - exponential regions of @xmath8 ( broken lines).,width=529 ]    again the basically exponential shape of these distributions is evident . characterizing the exponential fits by their inverse slopes @xmath37 a smooth and distinct dependence on lab angle and beam momentum becomes visible as shown in fig . [ fig : inverse ] .     as a function of @xmath10 for four beam momenta from 3 to 158  gev / c , a ) for protons , b ) for @xmath31 and c ) for @xmath32 . the full lines are drawn to guide the eye ]    compared to the strong dependence of b on @xmath10 which ranges from 0.3 to 0.05  gev / c , the only modest dependence on @xmath23 of @xmath380.03  gev / c for beam momenta from 3 to 158  gev / c is noticeable .    following the above data parametrization a generalized grid of @xmath8 values between 0.2 and 1.2  gev / c , in steps of 0.1  gev / c , may now be established . concerning the lower and upper limits of this grid , an extrapolation beyond the limits given by the experimental values has been performed in some cases . this extrapolation does not exceed the bin width of the respective data lists and is therefore defendable in view of the generally smooth , gentle and well - defined @xmath8 dependences . [ sec : ph ]    in the absence of theoretical predictability in the soft sector of the strong interaction , any attempt at bringing a multitude of experimental results into a common and consistent picture has to rely on a minimal set of model - independent physics constraints . in fact a `` democratic '' averaging of eventually contradictory data sets would only add confusion instead of clarity . [ sec : ph1 ]    two examples of the continuity constraint have already been mentioned above : invariant @xmath8 distributions have to approach zero momentum horizontally that is with tangent zero . the same is true for angular distributions in their approach to 180  degrees . [ sec : ph2 ]    it is a matter of experimental experience in the realm of soft hadronic interactions that in general distributions in any kind of kinematic variable tend to be `` smooth '' in the sense of absence of abrupt local upwards or downwards variations . the widespread use of simple algebraic parametrizations has its origin in this fact , specifically in the absence of local maxima and minima , with the eventual exception of threshold behaviour of which some examples will become visible below .      [ sec : ph3 ]    charge conservation has of course to be fulfilled by any type of experimental result . this means for instance that for the interaction of a positively charged projectile ( proton ) with an isoscalar nucleus ( carbon ) the @xmath39 ratio has to be greater or equal to unity over the full phase space invoking isospin symmetry ( and of course the experience from a wide range of experimental results ) . the presence of data with @xmath40  1 therefore immediately indicates experimental problems . the inspection of @xmath39 ratios has the further advantage that a large part of the systematic uncertainties , notably the overall normalization errors , cancel in this ratio . [ sec : ph4 ]    it has been shown that in proton induced nuclear collisions the yields of the secondary protons and neutrons are related by a constant factor of about 2 which is in turn related to the ratio of the basic nucleon - nucleon interaction @xcite . similarly , when rotating the projectile isospin from proton to neutron , it has been predicted that the yield ratio of secondary protons from proton and neutron projectiles    @xmath41    should be constant and equal to 2.5 for light nuclei @xcite . the extensive and precise low - energy data set of franz et al . @xcite from n+c interactions has therefore been included in the present survey . these data present a welcome extension of the @xmath24 scale into the region 0.47 to 0.49 which is not covered for most of the angular range with proton projectiles . as shown below , these data fit indeed very well , after re - normalization , into the general @xmath24 dependence of secondary protons where the low energy data by frankel et al . @xcite and komarov et al . @xcite at angles between 112 and 180  degrees provide an independent control of the normalization .      [ sec : ph5 ]    with these constraints in mind , and having established the parametrization and interpolation of the @xmath8 distributions as discussed above , one may now proceed to the attempt at sorting the 19 available experiments into a consistent global data set . it would of course be rather surprising if all experiments would fit into this global picture within their respective error limits . in fact it turns out that this procedure establishes a very strong constraint for possible deviations , as a large majority of results is creating a perfectly consistent picture both for protons and for pions . only four of the 19 data sets can not be brought into consistency with all other experiments without gravely affecting and contradicting the above constraints . these data are not included in the following global interpolation scheme . they will be discussed separately in sect . [ sec : dev_exp ] below . [ sec : prot ]      [ sec : prot_sdep ]    the invariant proton cross sections are shown in fig .  [ fig : oneoversq ] as a function of @xmath24 for a grid of ten lab angles between 25 and 180  degrees and constant lab momenta between 0.3 and 1.2  gev / c . the interpolated data points in each panel are identified by symbols corresponding to the different experiments .                                     at fixed @xmath8 and @xmath10 . the interpolated data points are indicated by symbols corresponding to the respective experiments in each panel . the solid lines represent the interpolation of the data , width=268 ]     at fixed @xmath8 and @xmath10 . the interpolated data points are indicated by symbols corresponding to the respective experiments in each panel . the solid lines represent the interpolation of the data , width=268 ]    the solid lines are eyeball interpolations through the data points . a first remark concerning this figure concerns the smoothness and continuity of the @xmath24 dependences . the achieved overall consistency of all data is rather impressive even if single points are deviating in some areas of phase space . the salient features of the physics contained in these plots may be summarized as follows :    * a strong yield suppression between 1/@xmath28  0.45 and the elastic limit at @xmath24  =  0.53 is evident . * the n+c data @xcite are well consistent with the p+c results in the overlap regions ; they define a broad maximum of the cross sections at 1/@xmath28  0.46 at medium angles and low @xmath8 . * there is a well - defined asymptotic behaviour of the cross sections for @xmath24 below about 0.2 or beam momenta above about 12  gev / c . * for the lower @xmath10 region and/or low @xmath8 the asymptotic region is approached from above . the latter point is reminiscent of the behaviour of the proton yields in p+p interactions , as shown in fig . [ fig : oneoversq_pp ] .     for p+p interactions at the four lab angles a ) 25 , b ) 35 , c ) 45 and d ) 55  degrees for @xmath8 values from 0.4 to 1  gev / c . the data are interpolated from blobel @xcite and na49 @xcite . the lines are drawn to guide the eye ]    another feature of fig . [ fig : oneoversq ] is the systematic droop of the cross sections from harp - cdp at their highest beam momentum of 15  gev / c or @xmath24  =  0.18 , demonstrating the discriminative power of the approach . this decrease is quantified in fig . [ fig : harp15_prot ] where the ratio @xmath43 between the measured invariant cross sections and the data interpolation is shown as a function of @xmath8 for the complete angular range from 25 to 97  degrees . here deviations of up to 50% are visible .     between the interpolated invariant proton cross sections from harp - cdp @xcite and the global interpolation as a function of @xmath8 for the angular range 25  @xmath44  97  degrees ] the abruptness of this decrease would necessitate a rather violent variation of the cross sections with increasing energy including a minimum between ps and sps energies . a final clarification of this situation is given by the proton data from serpukhov @xcite which , although suffering from a different and independent problem , at least exclude such variations in the region between 17 and 67  gev / c beam momentum , see sect . [ sec : belyaev ] below . [ sec : prot_thdep ]    in addition to the description of the energy dependence , the global interpolation has of course also to result in a smooth and continuous verification of the angular dependence presenting the third dimension of the present study . this constraint has to be fulfilled at any value of @xmath24 .    in a first example the situation at @xmath24  =  0.05 is shown in fig . [ fig : pc_prot_theta ] . this value lies in between the fermilab @xcite and na49 @xcite data in the region of negligible @xmath20-dependence . it therefore allows for the direct comparison of the two experiments in their respective angular regions which have no overlap . =  0.05 as a function of @xmath46 combining the fermilab and na49 data for @xmath8 between 0.2 and 1.4  gev / c . the global interpolation is shown as full lines . the measured cross sections in the angular ranges from 70 to 160  degrees @xcite and from 10 to 40  degrees @xcite are given on the vertical broken lines , width=472 ]    several observations are in place here :    * the two experimental results connect perfectly through the gap between the na49 ( @xmath47  40  degrees ) and the fermilab ( @xmath48  70  degrees ) data . * there is at most a few percent variation of the cross sections between the angles of 160 and 180  degrees taking into account the constraint of continuity through 180  degrees discussed in sect . [ sec : comp ] above . this allows the combination of results in this angular region as it is applied in the determination of the @xmath24 dependence , fig . [ fig : oneoversq ] . * the angular distributions are smooth and close to exponential in shape . in particular , no instability in the region around 90  degrees is visible where an eventual diffractive peak from target fragmentation would appear , see also @xcite . further angular distributions at four @xmath24 values between 0.1 and 0.4  gev@xmath49 are given in fig . [ fig : costheta ] . in fact such distributions at arbitrary values of @xmath24 may be obtained from the global interpolation as it is presented in numerical form at the na49 web page @xcite .     for four values of @xmath24 : a ) 0.1 , b ) 0.2 , c ) 0.3 , d ) 0.4  gev@xmath49 and for @xmath8 values between 0.4 and 1.2  gev / c . the standard grid of 10 angles , fig . [ fig : oneoversq ] , is indicated by the vertical broken lines ]    evidently the angular distributions maintain their smooth and continuous shape , specifically through 90  degrees , at all interaction energies . with the approach to low beam momenta however , a progressive rounding of the shape towards higher lab angles manifests itself . [ sec : pip ]    the global interpolation of the @xmath31 data is presented in this section in close analogy to the preceding section for protons .      [ sec : pip_sdep ]    the invariant @xmath31 cross sections are shown in fig .  [ fig : sqs_pip ] as a function of @xmath24 for the standard grid of ten lab angles between 25 and 180  degrees and for constant lab momenta between 0.2 and 1.2  gev / c . the interpolated data points in each panel are identified by symbols corresponding to the different experiments .                     in p+c collisions as a function of @xmath24 at fixed @xmath8 and @xmath50 . the interpolated data points are indicated by symbols corresponding to the respective experiments in each panel . the solid lines represent the global data interpolation , width=257 ]     in p+c collisions as a function of @xmath24 at fixed @xmath8 and @xmath50 . the interpolated data points are indicated by symbols corresponding to the respective experiments in each panel . the solid lines represent the global data interpolation , width=257 ]     in p+c collisions as a function of @xmath24 at fixed @xmath8 and @xmath50 . the interpolated data points are indicated by symbols corresponding to the respective experiments in each panel . the solid lines represent the global data interpolation , width=257 ]     in p+c collisions as a function of @xmath24 at fixed @xmath8 and @xmath50 . the interpolated data points are indicated by symbols corresponding to the respective experiments in each panel . the solid lines represent the global data interpolation , width=257 ]     in p+c collisions as a function of @xmath24 at fixed @xmath8 and @xmath50 . the interpolated data points are indicated by symbols corresponding to the respective experiments in each panel . the solid lines represent the global data interpolation , width=257 ]     in p+c collisions as a function of @xmath24 at fixed @xmath8 and @xmath50 . the interpolated data points are indicated by symbols corresponding to the respective experiments in each panel . the solid lines represent the global data interpolation , width=257 ]    the solid lines represent the global interpolation by eyeball fits of both the energy and the angular dependences . again the @xmath24 dependence is in general smooth and continuous , with an impressive overall consistency of all data with only few exceptions discussed below . there are some general trends to be pointed out :    * at the lowest lab momentum , the pion cross sections are practically @xmath20-independent , with variations of only 1020% in the range from 1 to 400  gev / c beam momentum . * this fact suggests @xmath31 production at low momentum transfer in the nuclear cascade . * for all lab momenta , the approach to high energies is very flat for @xmath51  0.2 or beam momenta above 12  gev / c . * the high energy cross sections are approached for all angles and beam momenta from below . there are two areas of deviation from the global interpolation which are both connected to the harp - cdp data @xcite . at their lowest angle of 25  degrees , the cross sections are systematically low by up to a factor of two below @xmath52  0.5  gev / c and @xmath24 above 0.2 . this is in contradiction to the available low energy data from other experiments also shown in fig . [ fig : sqs_pip ] . the second area concerns , as for the protons , the data at 15  gev / c beam momentum where a characteristic pattern of deviations is visible : at low angles and low @xmath8 , the data tend to overshoot the interpolation , whereas at angles above 45  degrees a progressive droop with increasing lab momentum is evident . this is quantified by the ratio @xmath43 between the harp - cdp data and the global interpolation shown in fig . [ fig : harp15_pip ] .     between the interpolated invariant @xmath31 cross sections from harp - cdp @xcite and the global interpolation as a function of @xmath8 for the angular range 25  @xmath53  112  degrees ]    these deviations are rather consistent with the ones found for protons . also in this case a rapid variation of the cross sections with increasing beam momentum can be excluded by the comparison with the pion data from the serpukhov experiment @xcite between 17 and 67  gev / c beam momentum , see sect . [ sec : belyaev ] below . [ sec : pip_thdep ]    as already shown in sect . [ sec : prot_thdep ] for protons , the angular distributions at @xmath24  =  0.05 , in between the fermilab @xcite and na49 @xcite energies , are presented in fig . [ fig : pc_pip_theta ] . this allows the comparison of the two data sets and their connection across the gap in lab angles between 40 and 70  degrees which represent the upper and lower limit of the respective experiment . cross sections at @xmath24  =  0.05 as a function of @xmath46 combining the fermilab @xcite and na49 @xcite data for @xmath8 between 0.2 and 1.2  gev / c . the global interpolation is shown as full lines . the measured cross sections in the angular ranges from 70 to 160  degrees ( @xcite ) and from 10 to 40  degrees ( @xcite ) are given on the vertical broken lines , width=472 ]    further angular distributions at four @xmath24 values between 0.1 and 0.4  gev@xmath49 are given in fig . [ fig : costheta_pip ] . cross sections as a function of @xmath46 for four values of @xmath24 : a ) 0.1 , b ) 0.2 , c ) 0.3 , d ) 0.4  gev@xmath49 and for @xmath8 values between 0.2 and 1.2  gev / c . the standard grid of 10 angles , fig . [ fig : oneoversq ] , is indicated by the vertical broken lines ]    the angular distributions are characterized by a smooth , close to exponential shape . at backward angles , the @xmath8 dependence is very steep with four orders of magnitude already between @xmath8  =  0.2 and 0.8  gev / c . in forward direction this dependence is much reduced with less than one order of magnitude between @xmath8  =  0.2 and 1.2  gev / c . this is due to the prevailance of target fragmentation in this region , see sect .  [ sec : separation ] for a quantitative study of this phenomenology . [ sec : pin ]    this section follows closely the discussion of the @xmath31 cross sections in the preceding section . [ sec : pin_sdep ]    the invariant @xmath32 cross sections are shown in fig .  [ fig : sqs_pin ] as a function of @xmath24 for the standard grid of ten lab angles between 25 and 180  degrees and for constant lab momenta between 0.2 and 1.2  gev / c . the interpolated data points in each panel are identified by symbols corresponding to the different experiments .                             in p+c collisions as a function of @xmath24 at fixed @xmath8 and @xmath50 . the interpolated data points are indicated by symbols corresponding to the respective experiments in each panel . the solid lines represent the global data interpolation , width=268 ]     in p+c collisions as a function of @xmath24 at fixed @xmath8 and @xmath50 . the interpolated data points are indicated by symbols corresponding to the respective experiments in each panel . the solid lines represent the global data interpolation , width=268 ]     in p+c collisions as a function of @xmath24 at fixed @xmath8 and @xmath50 . the interpolated data points are indicated by symbols corresponding to the respective experiments in each panel . the solid lines represent the global data interpolation , width=268 ]     in p+c collisions as a function of @xmath24 at fixed @xmath8 and @xmath50 . the interpolated data points are indicated by symbols corresponding to the respective experiments in each panel . the solid lines represent the global data interpolation , width=268 ]    the solid lines represent the global interpolation by eyeball fits to the data , with several features which are worth noticing :    * all the different data sets form a consistent ensemble without the systematic deviations visible in some regions of the proton and @xmath31 results . * the approach to large beam momenta happens from below for all @xmath8 . * the @xmath20-dependence is in general stronger than for @xmath31 , fig . [ fig : sqs_pip ] . if it is again flat up to @xmath54  0.2 at low @xmath8 , it becomes more pronounced both towards higher @xmath8 and in the approach to the production threshold at large @xmath24 indicating a marked increase of the @xmath39 ratio . * this effect has as physics origin the progressive change of the production mechanism from pion exchange at low energy to gluon or pomeron exchange at sps energy . this will be discussed in relation to the charge ratios in sect . [ sec : p2n ] . it is again interesting to compare the energy dependence to the one observed in p+p interactions as presented in fig . [ fig : oneoversq_pp_piminus ] . cross sections as a function of @xmath24 for p+p interactions at the two lab angles a ) 45 , and b ) 55  degrees for @xmath8 values from 0.2 to 1.2  gev / c . the data are interpolated from blobel @xcite and na49 @xcite . the lines are drawn to guide the eye ]    although for both reactions the asymptotic high energy region is approached from below , this comparison shows a stronger @xmath20-dependence , at the same lab angle , in p+c than in p+p collisions . this is due to the component of nuclear cascading which contributes , in the given angular range , with equal strength than the target fragmentation to the total yield ( see sect .  [ sec : separation ] below ) . [ sec : pin_thdep ]    as for protons and @xmath31 in figs . [ fig : pc_prot_theta ] and [ fig : pc_pip_theta ] , the @xmath32 cross sections from the fermilab @xcite and na49 @xcite experiments are compared and combined as a function of @xmath46 in fig .  [ fig : pc_pin_theta ] . cross sections at @xmath24=  0.05 as a function of @xmath46 combining the fermilab @xcite and na49 @xcite data for @xmath8 between 0.2 and 1.2  gev / c . the global interpolation is shown as full lines . the measured cross sections in the angular ranges from 70 to 160  degrees ( @xcite ) and from 10 to 40  degrees ( @xcite ) are given on the vertical broken lines , width=415 ]    further angular distributions at four @xmath24 values between 0.1 and 0.4  gev@xmath49 are given in fig . [ fig : costheta_pin ] . cross sections as a function of @xmath46 for four values of @xmath24 : a ) 0.1 , b ) 0.2 , c ) 0.3 , d ) 0.4  gev@xmath49 and for @xmath8 values between 0.2 and 1.2  gev / c . the standard grid of 10 angles , fig . [ fig : oneoversq ] , is indicated by the vertical broken lines , width=529 ]    concerning smoothness and continuity these distributions are similar to the @xmath31 data , including the large asymmetry between the forward and backward directions . the reduction of the cross sections for @xmath32 with respect to @xmath31 with increasing @xmath24 is however very apparent . this will be quantified in the following section on @xmath39 ratios . [ sec : p2n ]    as already evoked in sect . [ sec : ph3 ] above , the study of @xmath39 ratios has two main advantages . firstly , in this ratio a major fraction of the experimental systematic uncertainties cancels . secondly , the ratio is constrained by very fundamental and model - independent physics arguments like charge conservation and isospin symmetry . in addition , its @xmath20-dependence is governed by the hadronic meson exchange process which leads to a power - law behaviour that will be shown to be common to a wide range of interactions . in the following argumentation the ratio between the global data interpolation for @xmath31 and @xmath32 as described in the preceding sects . [ sec : pip ] and [ sec : pin ] will be used :    @xmath56    as a by - product , the fluctuation of this ratio as a function of angle and interaction energy will allow for the estimation of the local precision of the interpolation procedure . [ sec : p2n_high ]    it has been established by numerous experimental results that at collision energies in the sps / fermilab range and above the hadronic interactions are characterized by the absence of charge and flavour exchange . it has also been shown that the feed - over of pions from the projectile hemisphere into the backward region of @xmath13 is sharply limited to the range of @xmath57  -0.05 , see @xcite for a detailed discussion . this range is outside the coverage in @xmath10 and @xmath8 considered in this publication . it is therefore to be expected that the backward production of pions off an isoscalar nucleus should be charge - symmetric at high energy . this is indeed verified by the results on pion production shown in the preceding sections . it is quantified in fig . [ fig : costheta_pin ] which shows the @xmath39 ratio at @xmath24  =  0.04 or 330  gev / c beam momentum for all lab angles and lab momenta treated in this publication . this number distribution has a mean value of 1.0125 with an rms deviation of 3.2% . this rms value may be seen as a first estimate of the local precision of the three - dimensional interpolation scheme at this energy which has been established independently for both pion charges .    /@xmath32 ratio @xmath58 at @xmath24  =  0.04  gev@xmath49 for 25  @xmath59  162  degrees and 0.2  @xmath60  1.2  gev / c ]      [ sec : p2n_endep ]    with decreasing interaction energy or increasing @xmath24 the @xmath39 ratio develops a strong increase at all lab momenta and lab angles . this is shown in fig . [ fig : rpm_plab ] which gives @xmath61 as a function of @xmath24 for four lab momenta . the ratio of the global data interpolation is given in steps of 0.02 in @xmath24 . at each value of @xmath24 the number of points corresponds to the standard grid of angles available at this energy .     as a function of @xmath24 for four values of @xmath8 . the dots represent the ratio of the global data interpolation for @xmath31 and @xmath32 in steps of 0.02 in @xmath24 and for the chosen grid of 10 angles . a ) @xmath8  =  0.2 , b ) @xmath8  =  0.4 , c ) @xmath8  =  0.6 , d ) @xmath8  =  0.8  gev / c ]    several features of fig . [ fig : rpm_plab ] are noteworthy :    * considering the wide range of lab angles , @xmath61 is at each value of @xmath24 confined to a narrow band indicating an approximative angle independence . * large @xmath61 values in excess of 5 are reached at the upper limit of the available scale in @xmath24 . * there is a systematic increase of @xmath61 with @xmath8 . [ sec : p2n_high1 ]    the features pointed out above may be quantified and at the same time the local systematic fluctuations of the interpolation may be estimated by establishing the mean values @xmath62 averaged over the angular range at each @xmath24 . these mean values are well defined as shown in fig . [ fig : p2n_angdiff ] which presents the normalized distribution of the point - by - point deviations from the mean in percent ,    @xmath63    for four values of @xmath24 , summing the four @xmath8 values used in fig . [ fig : rpm_plab ] .     from the mean over the angular range , a ) @xmath24  =  0.1 , b ) @xmath24  =  0.2 , c ) @xmath24  =  0.3 , d)@xmath24  =  0.4 . the four @xmath8 values shown in fig .  [ fig : rpm_plab ] are summed up , width=283 ]    these distributions are of gaussian shape with an rms which increases with @xmath24 as indicated in fig . [ fig : sigsqs ] .     as a function of @xmath24 . the full line represents a hand interpolation . broken line : corresponding error margin of the mean values @xmath62 . dotted lines : corresponding errors for @xmath31 and @xmath32 separately ]    the observed energy dependence of the rms deviations is due to the fact that the invariant pion cross sections decrease , after a relatively flat behaviour up to @xmath54  0.15 , progressively steeper towards the production threshold , see figs . [ fig : sqs_pip ] and [ fig : sqs_pin ] . this leads inevitably to larger variations in the corresponding energy interpolation .    from the rms values given in fig . [ fig : sigsqs ] the error of @xmath62 may be derived which varies between 1% and 5% for the highest and lowest interaction energy , respectively ( broken line ) . also the corresponding error margins for the mean pion yields may be extracted as indicated by the dotted lines in fig . [ fig : sigsqs ] . from these plots it appears that the global interpolation induces fluctuations which increase from a few percent in the high @xmath20 region to about 10% in the approach to the pion threshold . [ sec : p2n_high2 ]    the dependence of @xmath61 on @xmath10 is shown in fig . [ fig : angdep ] for four values of @xmath24 and four values of @xmath8 .     as a function of @xmath10 for four values of @xmath24 , a ) @xmath8  =  0.2 , b ) @xmath8  =  0.4 , c ) @xmath8  =  0.6 , and d ) @xmath8  =  0.8  gev / c . the mean values @xmath62 are indicated as the horizontal broken lines in each panel , width=510 ]    evidently no systematic @xmath10 dependence is visible over the complete angular range within the quoted errors .      [ sec : p2n_high3 ]    in the absence of angular dependence of @xmath61 as shown above , the mean values @xmath62 may now be used in order to establish a precise view of the @xmath24 dependence for different @xmath8 values . this dependence is presented in fig . [ fig : p2n_meanang ] .     as a function of @xmath24 for five values of @xmath8 between 0.2 and 1.0  gev / c . the full lines are hand interpolations through the data points ]    within the errors of @xmath62 extracted above , a clear @xmath8 dependence is evident superposing itself to the strong common increase of @xmath62 with @xmath24 . this increase may be parametrized up to @xmath54  0.3 by the functional from @xmath64 which is , as discussed below , typical of meson exchange processes . indeed the exponent beta varies from 2 to 1.2 for @xmath8 increasing from 0.2 to 0.8  gev / c . [ sec : p2n_epdep ]    the strong increase of @xmath61 with @xmath24 merits a detailed study as it is directly connected to the basic hadronic production mechanisms in p+a interactions . the fact that the pion yields in the complete backward fragmentation region of an isoscalar nucleus remember the isospin of the projectile is clearly incompatible with charge and flavour independent exchange processes . instead a meson exchange mechanism may be invoked which has indeed been used successfully in a wide range of work at low projectile momenta , see for instance @xcite and references therein . close to the pion production threshold in the nuclear hemisphere , single excitation processes via pion exchange of the type    @xmath65    only allow @xmath31 and @xmath66 production , whereas @xmath32 production needs double excitation like    @xmath67    with in general an additional penalty for @xmath32 due to the isospin clebsch  gordan coefficients . all meson exchange mechanisms are characterized by a strong decrease with projectile energy . this energy dependence and its interplay with processes governing the high energy sector is studied here for the first time in p+a collisions using the @xmath39 ratio .    in this context it seems mandatory to first refer to the study of exclusive charge exchange reactions in elementary nucleon - nucleon collisions as the complete energy range discussed here has been covered there by a number of experiments @xcite . [ sec : p2n_epdep1 ]    charge exchange processes may be cleanly isolated experimentally in nucleon - nucleon interactions by studying the following exclusive channels :    * charge exchange scattering of the elastic type + & + + & * single dissociation with pion production + [ eq : single ] & + + ^++ + ( + ^+ ) & * double dissociation with pion production + [ eq : double ] & + ( + ^- ) + ( + ^+ ) &    these channels are characterized by a very steep energy dependence . this is to be confronted with non - charge - exchange exclusive channels like :    * elastic scattering + & + + & * single dissociation + [ eq : single0 ] & + + ( + ^+ + ^- ) & * double dissociation + [ eq : double0 ] & + ( + ^+ + ^- ) + ( + ^+ + ^- ) &    which show a constant or logarithmically increasing @xmath20-dependence . charge exchange scattering has been measured by five experiments in the range of neutron beam momenta from 3 to 300  gev / c . this is exactly covering the energy range discussed in this paper . the single and double dissociation has been studied at the cern isr by two experiments @xcite extending the energy scale to @xmath20  =  3700  gev@xmath68 . the two isr experiments may be directly compared to the charge exchange measurements after appropriate re - normalization of the cross sections in the overlap region at the lowest isr energy . the resulting @xmath20-dependence at a momentum transfer @xmath69  =  0.032  gev@xmath68 is presented in fig . [ fig : chargeexch ] .     at a momentum transfer @xmath69  =  0.032  gev@xmath68 . the full line represents an interpolation of the data points . the insert gives the local slope @xmath70 in the parametrization @xmath71 as a function of @xmath20,width=491 ]    several features of fig . [ fig : chargeexch ] are of interest :    * there is a decrease of about 4 orders of magnitude in cross section between the lowest and highest @xmath20 value . this decrease is to be compared to the constant or logarithmically increasing elastic and nucleon diffraction cross sections . the charge exchange contribution is therefore negligible compared to the inclusive baryon yields already at sps energy . * there is a steady decrease of the local slope @xmath72 with energy , from about 3.6 at 3  gev / c to about 1.1 above 80  gev / c beam momentum . * a characteristic change of slope manifests itself at around 30  gev / c beam momentum . these features have been interpreted in the 1970 s when the relevant experiments were performed , in the framework of regge theory which predicts an @xmath20-dependence of the form    @xmath73    where @xmath26 is the intercept of the leading trajectory . this should in the case of one - meson exchange at low energy be given by the pion trajectory with zero intercept . the actual beta values above 3 at low @xmath20 seem to contradict however this expectation . here threshold effects may play a role which are not included in the parametrization ( [ eq : regge ] ) . with increasing energy the slopes move through the region of pion exchange with @xmath74  2 down to values of about 1.1 at high energy which could be connected to @xmath75 and a@xmath76 exchange with correspondingly higher intercepts @xmath26 in the region of 0.5 . at isr energy the ratio of @xmath77 contributions has indeed been estimated to be about 2 @xcite . anyway the simple parametrization given by ( [ eq : regge ] ) should not be expected to hold over the full energy scale . what is interesting here is rather the strong decline of the charge exchange cross sections with energy and the experimentally rather precisely determined slope variation . [ sec : p2n_epdep2 ]    the single ( [ eq : single ] ) and double ( [ eq : double ] ) dissociation processes defined above are determined by the formation of @xmath35 resonances in the final states . they therefore constitute a source of direct @xmath35 production in nucleon - nucleon interactions . these channel cross sections decrease rapidly to the @xmath22barn level at sps energies . in contrast , the non - charge exchange channels like ( [ eq : single0 ] ) and ( [ eq : double0 ] ) have no @xmath20-dependence and stay on the mb level of cross sections . their final states have been shown to be governed by n@xmath78 resonances @xcite which may be excited by pomeron exchange . moreover , the p+@xmath31 combination of the p+@xmath31+@xmath32 final states has been shown to be dominated by @xmath79 @xcite . this is an indirect source of @xmath35 resonances as a decay product of n@xmath78 states which have large decay branching fractions into @xmath35+@xmath80 and @xmath35+@xmath75 . it is therefore questionable if , at sps energies and above , any direct @xmath35 production is persisting . this is an interesting question for the majority of microscopic models which produce final states by string fragmentation . in the baryonic sector , diquark fragmentation is generally invoked with a prevailing direct production of @xmath35 resonances which by isospin counting will dominate over n@xmath78 . indeed in practically all such models there is no or only negligible n@xmath78 production . as shown below , the decrease of charge exchange processes can be traced well into the non - diffractive , inelastic region of particle production . the multi - step , cascading decay of primordial n@xmath78 resonances into @xmath35 resonances and final state baryons should therefore be seriously considered , in particular also concerning the consequences for the evolution of final state energy densities with time . [ sec : p2n_epdep3 ]    the very characteristic decrease of @xmath62 with increasing @xmath20 derived from the global data interpolation , fig . [ fig : p2n_meanang ] , offers a tempting possibility of comparison to the phenomenology discussed above for the elementary nucleon - nucleon sector . indeed , two components should contribute to the observed @xmath39 ratios : at high energy this ratio should approach unity due to the absence of charge and flavour exchange in this region . at low energy on the contrary it should be governed by meson exchange with its strong @xmath20-dependence . these two components may be tentatively separated by using instead of @xmath62 the quantity    @xmath81    in order to extract the meson exchange contribution . this quantity is plotted in fig . [ fig : p2n_meanang1 ] as a function of @xmath20 for four @xmath8 values from 0.2 to 0.8  gev / c .     as a function of @xmath20 for different values of @xmath8 , a ) @xmath8  =  0.2 , b ) @xmath8  =  0.4 , c ) @xmath8  =  0.6 , and d ) @xmath8  =  0.8  gev / c . the elastic limit is indicated by the arrows ]    a very characteristic pattern emerges which resembles the @xmath20-dependence for the charge exchange in elementary interactions described above , see fig . [ fig : chargeexch ] . in general @xmath82 follows a power law dependence on @xmath20    @xmath83    with local slopes @xmath84 which are in turn a function of @xmath20 . three different regions with distinct local slopes can be identified in fig . [ fig : p2n_meanang1 ] :    * a first region with large slopes is located at @xmath20 below about 6  gev@xmath68 . this region is strongly influenced by threshold effects as the threshold for inelastic production is placed at the elastic limit @xmath85 indicated in fig . [ fig : p2n_meanang1 ] . in the approach to pion threshold the @xmath39 ratio has to diverge as @xmath32 is progressively suppressed , see above . with increasing @xmath8 this suppression will of course be more pronounced . * an intermediate region between about 8 and 40  gev@xmath68 with an @xmath20 dependence decreasing with increasing @xmath8 . * a third region with flattening @xmath20-dependence above about 40  gev@xmath68 .    at the lowest @xmath8 value of 0.2 gev / c corresponding to the lowest momentum transfer , the similarity to the charge exchange process in nucleon - nucleon interactions , fig . [ fig : chargeexch ] , is absolutely striking . this concerns both the detailed shape and the overall suppression factors . with increasing @xmath8 , the @xmath20 dependence is modified in a systematic way by a general reduction of slopes , with the exception of the threshold enhancement this is quantified in fig . [ fig : beta ] which shows the local slopes as a function of @xmath20 for @xmath8 values between 0.2 and 1  gev / c .     of the @xmath20-dependence of @xmath82 as a function of s : a ) @xmath8  =  0.2 , b ) @xmath8  =  0.4 , c ) @xmath8  =  0.6 , d ) @xmath8  =  0.8 , and e ) @xmath8  =  1.0  gev / c . the shaded regions mark the error margins ]    with the exception of the threshold region , the slopes are confined to the region between 2 and 1 typical of meson exchange processes . the dependence on @xmath8 is given in fig . [ fig : beta_pldep ] where the slopes in the three regions of @xmath20 specified above are presented . 6  gev@xmath68 ( broken line ) , 6  @xmath86  40  gev@xmath68 ( full line ) and @xmath87  40  gev@xmath68 ( dotted line ) as a function of @xmath8 between 0.2 and 1  gev / c ]    this figure shows clearly the different nature of the low @xmath20 enhancement where the slopes increase strongly with @xmath8 . the two other regions , full and dotted lines , are compatible with a regge parametrization with trajectory intercepts which increase with @xmath8 . this is insofar interesting as the region of measurements regarded here covers the complete backward angular range and the corresponding interactions are by no means confined to diffractive or low momentum transfer collisions . it is shown in sect . [ sec : separation ] of this paper that in the backward hemisphere the pion yields from nuclear cascading and target fragmentation are comparable . if the nuclear component is characterized by low momentum transfer reactions @xcite the target fragmentation is manifestly inelastic and non - diffractive . it governs the total yield at all angles below about 70  degrees .    in conclusion of this study of @xmath39 ratios in p+c interactions the following points should be stressed :    * the global data interpolation leads to a precise and consistent description of the behaviour of the @xmath39 ratios in the full backward hemisphere , thus offering an additional tool for the discrimination of experimental deviations . * the inspection of the detailed @xmath20-dependence of the ratios opens a new window on the underlying exchange processes . * in particular the comparison to the elementary nucleon - nucleon collisions establishes a close relation between apparently disjoint sectors of the different hadronic interactions . [ sec : dev_exp ]    as mentioned in sect . [ sec : ph ] four of the 19 investigated data sets are incompatible with the attempt at generating an overall consistent description of the experimental situation . these data will be shortly discussed below .      [ sec : geaga ]    these data have been obtained at the bevalac using beam momenta of 1.75 , 2.89 and 5.89  gev / c , spanning a lab momentum range from 0.3 to 0.9  gev / c at a lab angle of 180  degrees . the resulting cross sections trace the shape of the @xmath24 dependence rather precisely but are consistently about a factor of two below the global interpolation as shown in fig . [ fig : bevalac ] . here the full lines correspond to the global interpolation and the broken lines give the interpolation divided by a factor of two .         for three values of @xmath24 ; a ) 0.427 , b ) 0.366 and c ) 0.281 gev@xmath49 . the full lines give the global data interpolation , the open circles the data from @xcite , width=585 ]    as the angular bin from 160 to 180  degrees is mostly covered by data around 160162  degrees , a steep angular dependence in this region can not a priori be excluded . the smooth and gentle angular dependence of the interpolated data shown in fig . [ fig : geaga ] for the angular range from 82 to 180  degrees and for the three @xmath24 values of ref . @xcite , together with the constraint of the approach to 180  degrees with tangent zero , excludes however a drop of the cross sections by a factor of two between 160 and 180  degrees . [ sec : belyaev ]    a sizeable set of data on proton @xcite and pion @xcite production has been obtained at the serpukhov accelerator spanning the range of beam momenta between 17 and 57  gev / c . this fills the gap between the ps and sps energies where no other data are available . the data cover the @xmath8 range from 0.25 to 1.2  gev / c at @xmath10  =  159  degrees . they are presented in fig . [ fig : serp_sdep ] in comparison to the global data interpolation .     at @xmath10 =  159  degrees ( open circles ) in comparison to the global data interpolation at 160 degrees ( full lines),width=604 ]    several features of this comparison are noteworthy :    * the shape of the @xmath24 dependences complies precisely with the global interpolation . this is compatible with the absence of rapid variations of the cross sections with energy in the region between ps and sps . * there is a pronounced suppression of these data with respect to the interpolation with increasing @xmath8 reaching factors of three at the upper ranges for protons and pions . * the @xmath31 and @xmath32 data show an identical behaviour . * the proton data are tracing the interpolation up to @xmath8  =  0.4  gev / c whereas the pion data are already suppressed in this @xmath8 range . * the suppression factors are generally bigger for the pions at equal @xmath8 . * these features might be compatible with a momentum scale error . in addition to the reproduction of the shape of the @xmath24 dependence , also the @xmath39 ratio complies exactly with the one extracted from the global interpolation , fig . [ fig : serp_p2n ] , up to @xmath8  =  0.7  gev / c . above this value there is a sharp drop of @xmath61 reaching unphysical values at the upper limit of @xmath8 . this drop of about 20% has however to be compared to a drop of 300% of the invariant cross sections at this limit .     from @xcite as a function of @xmath8 averaged over beam momentum between 17 and 57  gev / c . the full line gives the result of the global interpolation averaged over the same beam momentum scale , sect . [ sec : p2n],width=340 ]      [ sec : harp ]    these results cover a range from 3 to 12  gev / c beam momentum at @xmath10 between 25 and 117  degrees and 0.125  @xmath60  0.75  gev / c . they are thus directly comparable to the ones from @xcite which are part of the global data interpolation . their differences to this interpolation are presented in fig . [ fig : harpdiffdist ] for all beam momenta and the standard grid of @xmath8 and @xmath10 values . , panel b ) @xmath31 ]    if the mean values of the differences are close to zero , their number distributions show wide spreads especially for @xmath31 . this is exemplified in fig . [ fig : pion_sdep ] where a typical comparison to the global interpolation ( full lines ) is given as a function of @xmath24 at @xmath10  =  67  degrees for four @xmath8 values . =  67  degrees and four @xmath8 values as a function of @xmath24 , ( open circles ) in comparison to the global data interpolation ( full lines ) . panel a ) @xmath31 , panel b ) @xmath32,width=585 ]    a comparison of @xmath39 ratios is given in fig . [ fig : harp_p2n_5 ] as a function of @xmath8 for four values of @xmath10 at @xmath23  =  5  gev / c .     from @xcite as a function of @xmath8 for four lab angles and a beam momentum of 5  gev / c . the full lines give the results of the data interpolation from @xcite , the full circles correspond to the ratios of the measured cross sections . the broken lines give the result of the global interpolation ]    this figure demonstrates the importance of using , in addition to the invariant cross section proper , the particle ratios which are strongly constrained by physical arguments , see sect . [ sec : p2n ] . [ sec : shine ]    these data have been obtained at a beam momentum of 31  gev / c in a @xmath10 range from 0.6 to 22.3  degrees and @xmath8 from 0.2 to 18  gev / c . if a large part of the given angular and momentum coverage falls outside the backward region regarded here , the low momentum range up to @xmath52  0.5  gev / c for all angles and the range 0.6  @xmath60  1  gev / c for angles above about 9  degrees corresponds to negative @xmath13 and can therefore be considered here . the complete @xmath24 dependence established in the preceding sections has a lower angular limit at 25  degrees corresponding to the lowest value of the standard grid of angles . this angle is close to the highest angle of @xcite at 22.3  degrees allowing for a safe interpolation . this is shown in fig . [ fig : na61 ] where the global interpolation is compared to the data of @xcite at their highest angles between 12 and 22  degrees for two @xmath8 values both for @xmath31 and for @xmath32 .     at @xmath24 =  0.13  gev@xmath49 in the @xmath10 region from 12 to 35  degrees for @xmath8  =  0.4 and 0.6 full lines : global data interpolation . full dots and broken lines : data from @xcite ]    as the global interpolation is limited to @xmath12  25  degrees , another way of comparison is offered by the combined na49 and fermilab results at 158/400  gev / c where the former data cover the complete angular range of @xcite . the ratio of the available high energy data to the results at 31  gev / c beam momentum is shown , as a function of @xmath10 , in fig .  [ fig : na61_ratio ] .     for @xmath8 =  0.4 and 0.6  gev / c . the full circles correspond to the global data interpolation , the open circles to the direct ratio @xcite/@xcite . the broken lines give the mean ratio over the complete range of the global data interpolation ]    apparently this cross section ratio is within errors angle independent over the full range of the global data survey , with well defined averages below 1.05 for @xmath31 and 1.1 for @xmath32 . in contrast , the cross section ratio between na49 and ref . @xcite shows values in the region of 1.4 increasing with decreasing @xmath10 . in conclusion to this section it may be stated that the global data interpolation between 15 different experiments attempted in this paper proves to be a useful tool for the detection of deviating data sets . further details concerning the above comparisons can be found in an internal report on ref . [ sec : separation ]    hadronic production in the backward direction of p+a collisions has two components : the fragmentation of the target nucleons which have been hit by the projectile proton , and the propagation of momentum transfer into the nucleus by secondary nucleon - nucleon interaction which follow , on a longer time scale , the initial excitation process . both processes are governed by the mean number of collisions @xmath89 suffered by the projectile on his trajectory through the nucleus .    as only the sum of these two separate mechanisms is experimentally accessible , a minimum assumption about the fragmentation of the target nucleons is needed in order to allow the separation of the components in an otherwise model - independent fashion . this minimal assumption consists in assuming that the fragmentation process of the hit nucleons is equal to the basic nucleon - nucleon interaction , taking full account of course of isospin symmetry . in addition and only valid for the relatively small value of @xmath89 in the carbon nucleus , it will be assumed that successive collisions result in hadronization at full interaction energy of the corresponding elementary interactions . as far as the value of @xmath89 is concerned , this has been determined for pion production in some detail in @xcite using the forward and the backward region at @xmath90  -0.1 where no intra - nuclear cascading is present , see below . this determination used three independent approaches :    * a monte - carlo calculation using the measured nuclear density distributions . * the relation between the inelastic cross sections of p+p and p+c interactions . * the approach to @xmath13  =  -0.1 of the ratio of pion densities in p+c and p+p collisions . the two former methods have to make the assumption that the inelastic interaction cross sections are independent of the number of subsequent collisions @xmath91 .    in @xcite a similar approach is used concerning the production of protons and anti - protons , again in the regions where there is no contribution from nuclear cascading as well as in the full backward hemisphere . all methods mentioned above result in a consistent estimate of @xmath89  =  1.6 in p+c collisions , with a relative systematic uncertainty of the order of a few percent .    in the following argumentation a prediction of the mean pion density of target fragmentation in the backward hemisphere at @xmath92  =  17.2  gev will be used which is relying on the published pion data from na49 @xcite and the estimated mean number of collisions , @xmath89 . the invariant pion cross sections are divided by the inelastic cross section to yield the quantity    @xmath93    per inelastic event which establishes isospin symmetry , and    @xmath94    this prediction is transformed into the appropriate coordinates @xmath8 and @xmath10 and divided by the measured invariant p+c cross sections @xmath95 per inelastic event yielding the ratio    @xmath96    this ratio is shown in fig . [ fig : rpred_plab ] as a function of @xmath8 for the lab angles 10 , 20 , 30 , 40 and 45  degrees .     as a function of @xmath8 for the five angles 10 , 20 , 30 , 40 and 45  degrees ]    it is evident that the ratio is close to one for the three lowest angles at all @xmath8 and for the region below 0.8 gev / c for 40 and 45  degrees . this is quantified in fig . [ fig : rpred_hist ] which gives the distribution of the ratio for the mentioned @xmath8 ranges .     for the lab angles 10 , 20 , 30  degrees and 40 and 45  degrees at @xmath30  0.8  gev / c . full line : gauss fit to the distribution with a mean value at 1.01 and a relative rms of 4.3% ]    the results show that indeed the measured pion cross sections correspond for lab angles up to 45  degrees precisely to the prediction from elementary collisions . this indicates that there is no contribution from intra - nuclear cascading in this region , in accordance with the results of @xcite . a drop of the ratio becomes however visible in the higher @xmath8 range at 40 and 45  degrees . this marks the onset of a nuclear component which becomes clearly visible in the ratios at larger angles shown in fig . [ fig : rpred_plab_high ] .     as a function of @xmath8 for the angles of 45 , 55 , 70 , 90 , 112 , 137 and 160 the full lines are local interpolations , width=396 ]    it is interesting to note that the target fragmentation governs the pion density up to the highest lab angles at low @xmath8 , with @xmath97 values of more than 50% . the ratio decreases however steadily with increasing @xmath8 and reaches zero at distinct momentum values indicating the approach to @xmath13  =  -1 in the plots of fig . [ fig : kinematics ] . this kinematic effect is more clearly brought out in fig . [ fig : rpred_kin ] showing that the fraction of target fragmentation is essentially a function of @xmath13 and rather independent on lab angle .     and @xmath10 for constant values of @xmath98 between 0.1 and 0.9 , b ) correlation between @xmath8 and @xmath10 for fixed values of @xmath13 between -0.1 and -1.0 and c ) @xmath98 as a function of @xmath13 ]    the correlation between @xmath8 and @xmath10 for fixed values of @xmath97 shown in panel a traces rather exactly the kinematic correlation between the same variables for fixed values of @xmath13 , panel b. this allows to establish a direct dependence of @xmath97 on @xmath13 which is to first order angle - independent , panel c.    the invariant densities @xmath100 per inelastic event as predicted from the fragmentation of the participant target nucleons is presented in fig . [ fig : fpred_theta ] .     per inelastic event as a function of @xmath29 for fixed values of @xmath8 between 0.1 and 1.2  gev / c . the full lines represent data interpolations , width=415 ]    this density may be subtracted from the pion density @xmath101 measured in p+c interactions which is within errors equal for @xmath31 and @xmath32 , see figs . [ fig : pc_pip_theta ] and [ fig : pc_pin_theta ] . the resulting invariant density    @xmath102    is shown in fig . [ fig : nucl_theta ] .     from intra - nuclear cascading as a function of @xmath29 for fixed values of @xmath8 between 0.1 and 1.2  gev / c . the full lines represent data interpolations , width=415 ]    this subtraction procedure becomes of course uncertain in the small angle region where the nuclear component is on the few percent level and below with respect to the target fragmentation , see figs . [ fig : rpred_plab ] and [ fig : rpred_plab_high ] . the invariant angular distributions shown in figs .  [ fig : fpred_theta ] and [ fig : nucl_theta ] may be converted into number distributions following :    @xmath103    and    @xmath104    integrating these distributions over @xmath8 , the number distributions    @xmath105    and    @xmath106    are obtained which are shown in fig . [ fig : dndcostheta ] together with the ratio    @xmath107     as a function of @xmath108 ( full line ) , @xmath109 as a function of @xmath108 ( broken line ) , b ) the ratio @xmath110 as a function of @xmath108,width=453 ]    evidently the nuclear component of pion production stays comparable to the target fragmentation in the full backward hemisphere of @xmath10 . it decreases rapidly for @xmath10 below about 60  degrees and vanishes below @xmath10  25  degrees . integration of @xmath111 over @xmath108 results in the total single pion yield from nuclear cascading    @xmath112    per inelastic event . the predicted integrated yield from target fragmentation is    @xmath113    with    @xmath114    and    @xmath115    from p+p interactions as measured by na49 , @xcite . this means that for p+c interactions the nuclear component of pion production amounts to 4.9% of the pions originating from the fragmentation of the hit target nucleons . applying isospin symmetry on the isoscalar c nucleus with    @xmath116    the total pion yields are 6.45 from target fragmentation and 0.315 from nuclear cascading . making use of the kinematic relation between the coordinate pairs @xmath8 , @xmath10 and @xmath13 , @xmath14 , see fig . [ fig : kinematics]c , the double differential yields for the nuclear component as functions of @xmath13 and @xmath14    @xmath117    may be obtained where @xmath118 , ( [ eq : xf ] ) , and @xmath119 are cms quantities . the resulting pion density distributions are shown in fig . [ fig : dndxfdpt ] as a function of @xmath13 for @xmath14 values from 0.05 to 0.7  gev / c .     as a function of @xmath13 for a ) 0.05  @xmath120  0.3  gev / c and b ) 0.4  @xmath120  0.7  gev / c , width=453 ]    a peak at low @xmath14 and @xmath13  =  -0.15 is apparent which corresponds to the location of pions with small lab momentum , see fig . [ fig : kinematics ] . with increasing @xmath14 the maximum density decreases and shifts in @xmath13 to lower values which is again in accordance with the kinematic correlation visible in fig . [ fig : kinematics ] . integration over @xmath14 results in the single differential density @xmath121 shown in fig . [ fig : dndxf ] together with the predicted density distribution @xmath122 from target fragmentation and with the ratio of the two densities .     as a function of @xmath13 ( full line ) . the predicted density distribution from target fragmentation @xmath123 is shown as the broken line ; b ) ratio @xmath124 as a function of @xmath13,width=453 ]    the @xmath14 integrated pion density @xmath121 shows a peak at @xmath125  -0.2 and vanishes at @xmath125  -0.08 . as shown by the density ratio with the predicted target fragmentation @xmath122 in fig . [ fig : dndxf]b , the nuclear component reaches  10% of the target fragmentation at @xmath13  =  -0.15 and exceeds this contribution for @xmath15  -0.55 . the nuclear pion component extracted above is used in @xcite in conjunction with the complementary nuclear proton component to obtain the percentage of cascading protons which are accompanied by pion emission . [ sec : conclusion ]    this paper presents a survey of available data concerning backward proton and pion production in minimum bias p+c interactions , including new and extensive data sets obtained at the cern ps and sps . the backward direction being defined as the complete phase space at negative feynman @xmath13 , the data cover , for projectile momenta from 1 to 400  gev / c , the ranges from 0.2 to 1.2  gev / c in lab momentum @xmath8 and from 10 to 180  degrees in lab angle @xmath10 . the paper attempts an interconnection of the different data sets by a detailed three - dimensional interpolation scheme in the variables @xmath24 , @xmath8 , and @xmath108 . this attempt allows a precise control of the internal data consistency as well as the study of the evolution of the invariant inclusive cross sections in all three variables . a literature search has provided a set of 19 different experiments with a total of more than 3500 data points . these measurements were obtained over 40 years of experimentation by collaborations employing widely different experimental techniques . in this respect it may be stated as a first positive result that the majority of the data may be combined into a surprisingly self - consistent ensemble . this global interpolation scheme results in a considerable discriminative power against the systematic deviation of particular data sets . only 4 of the 19 quoted experiments show in fact deviations which clearly mark them as systematically diverging . these experiments are inspected in detail one by one in an attempt to clearly bring out the discrepancies . in some of the cases , possible experimental error sources are pointed out .    the underlying physics provides for additional constraints concerning basic quantities like charge conservation and isospin symmetry as well as the necessity of smoothness and continuity of the observed cross sections . whenever possible , contact to the complementary elementary nucleon - nucleon interactions is established . this concerns in particular the evocation of mesonic exchange processes for the description of @xmath39 ratios and the prediction of the target fragmentation from elementary interactions and its separation from the component of nuclear cascading . as far as the dependences of the invariant cross sections on the three basic variables @xmath8 , @xmath10 and @xmath24 is concerned , a well constrained phenomenology emerges . the @xmath8 dependences are exponential or close to exponential over a major part of the phase space with some exceptions mostly towards low interaction energies . this fact results in an important constraint for the data interpolation . the @xmath108 dependences are not far from exponential and smooth and continuous through all lab angles . in particular there is no indication of an instability around 90  degrees for the proton yields . the @xmath24 dependences converge , after strong variations close to production threshold , smoothly to asymptotic behaviour in the sps energy range . this region is approached from above by the protons and from below for the pions . this convergence is confirmed by the @xmath39 ratios which show , being governed by meson exchange at low @xmath92 with large values marked by the projectile isospin , a smooth decline with energy towards unity as expected from the underlying elementary exchange processes . this work was supported by the polish state committee for scientific research ( p03b00630 ) , the polish national science centre ( on the basis of decision no . dec-2011/03/b / st2/02634 ) the bulgarian national science fund ( ph-09/05 ) , the eu fp6 hrm marie curie intra - european fellowship program , the hungarian scientific research fund otka ( t68506 ) and the hungarian otka / nkth a08 - 77719 and a08 - 77815 grants . [ 20 ] y.  d.  bayukov et al . , phys . rev . * c20 * , ( 1979 ) 764 b.  baatar et al . , arxiv:1207.6520v1 [ hep - ex ] i.  m.  belyaev et al . , * jinr - p1 - 90 - 551 * ( 1990 ) ( in russian ) a.  bolshakova et al . , j. * c70 * ( 2010 ) 573 n.  a.  burgov et al . phys . * 30 * ( 1979 ) 371 , yad . fiz . * 30 * ( 1979 ) 720 y.  d.  bayukov et al . . phys . * 18 * ( 1974 ) 639 + y.  d.  bayukov et al . * 19 * ( 1974 ) 648 j.  v.  geaga et al . * 45 * ( 1980 ) 1993 s.  frankel et al . lett . * 36 * ( 1976 ) 642 v.  i.  komarov et al . , phys . lett . * 69b * ( 1977 ) 37 j.  franz et al . , nucl . phys . * a472 * ( 1987 ) 733 n.  a.  nikiforov et al . , phys . rev . * c22 * ( 1980 ) 700 c.  alt et al . j. * c49 * ( 2007 ) 897 i.  m.  belyaev et al . , sov . * 49 * ( 1989 ) 295 n.  abgrall et al . , phys . rev . * c84 * ( 2011 ) 034604 m.  g.  catanesi er al . , phys . rev . * c77 * ( 2008 ) 055207 n.  a.  burgov et al . * 32 * ( 1980 ) 219 , yad * 32 * ( 1980 ) 423 a.  m.  baldin et al . , * 20 * ( 1975 ) 629 d.  r.  f.  cochran , phys * d6 * ( 1972 ) 3085 j.  f.  crawford , phys . rev . * c22 * ( 1980 ) 1184 j.  whitmore et al . , phys . lett * 60b * ( 1976 ) 211 + t.  ferbel , phys . * 29 * ( 1972 ) 448 + a.  h.  mueller , phys d2 * ( 1970 ) 2963 c.  alt et al . j. * c45 * ( 2006 ) 343 b.  d.  anderson et al . lett . * 46 * ( 1981 ) 226 d.  h.  boal , phys . rev . * c25 * ( 1982 ) 3068 v.  blobel et al . , nucl .phys . * b135 * ( 1978 ) 379 t.  anticic et al . j. * c65 * ( 2010 ) 9 http://cern.ch/spshadrons v.  blobel et al . , nucl . phys . * b69 * ( 1974 ) 454 g.  barr et al . j. * c49 * ( 2007 ) 919 e.  l.  miller et al . * 26 * ( 1971 ) 984 m.  n.  kreisler et al . , nucl . phys . * b84 * ( 1975 ) 3 v.  bhmer et al . * b110 * ( 1976 ) 205 a.  babaev et al . , nucl . phys . * b110 * ( 1976 ) 189 h.  r.  barton et al . , phys . * 37 * ( 1976 ) 1656 h.  de  kerret et al . * 69b * ( 1977 ) 372 g.  goggi et al . , phys . lett . * 72b * ( 1977 ) 265 g.  goggi et al . , phys . lett . * 79b * ( 1978 ) 165 c.  conta et al . , nucl . phys . * b175 * ( 1980 ) 97
Abstract:

Answer: new data on proton and pion production in p+c interactions from the cern ps and sps accelerators are used in conjunction with other available data sets to perform a comprehensive survey of backward hadronic cross sections . this survey covers the complete backward hemisphere in the range of lab angles from 10 to 180 degrees , from 0.2 to 1.4  gev / c in lab momentum and from 1 to 400  gev / c in projectile momentum . using the constraints of continuity and smoothness of the angular , momentum and energy dependences a consistent description of the inclusive cross sections is established which allows the control of the internal consistency of the nineteen available data sets . o.  chvala@xmath0 , h.  g.  fischer@xmath1 , m.  makariev@xmath2 , a.  rybicki@xmath3 , d.  varga@xmath4 , s.  wenig@xmath5    @xmath4etvs lornd university , budapest , hungary + @xmath3h . niewodniczaski institute of nuclear physics , polish academy of sciences , cracow , poland + @xmath5cern , geneva , switzerland + @xmath6charles university , faculty of mathematics and physics , institute of particle and nuclear physics , prague , czech republic + @xmath2institute for nuclear research and nuclear energy , bas , sofia , bulgaria + @xmath7now at university of tennessee , knoxville , tn , usa +