*******************************************************************
*
* Kumac to plot the E791 D^0 production cross section versus xf.
* The paper using this data has been submitted to Physics Letters B
* and can be found at Los Alamos as hep-ex/9906034 and at Fermilab
* as FERMILAB-Pub-99/185-E.
*
* Kumac author: Kevin Stenson, stenson@fnal.gov.
* Date:         June 26, 1999.
*
* Information about this data:   
* (1) The data are fully acceptance corrected.
* (2) The errors include statistical and systematic errors.
* (3) The systematic errors dominate.
* (4) The error bars do not include a +10-11% normalization uncertainty.
* (5) The data are integrated over all pt^2.
* (6) A cross section with a linear dependence on the atomic number
*     is assumed.
*
******************************************************************

* The limits of the xf bins
VEC/CRE XFBINS(21) r_
 -0.125 -0.10  -0.075 -0.05  -0.025  0.00   0.025 _
  0.05   0.075  0.10   0.125  0.15   0.175   0.20 _
  0.25   0.30   0.35   0.40   0.50   0.60   0.80  

* The midpoints of the xf bins.
vec/cre xfmids(20) R _
  -0.1125 -0.0875 -0.0625 -0.0375 -0.0125  0.0125 _
   0.0375  0.0625  0.0875  0.1125  0.1375  0.1625  0.1875 _
   0.2250  0.2750  0.3250  0.3750  0.4500  0.5500  0.7000

* The sizes of each xf bin
vec/cre xfdif(20) R 13*0.0125 4*0.025 2*0.05 0.10

*title 'x?F! distribution of D^0! (K[p] + K[ppp])'

* The central values
vec/cre d0xf(20) R   29.05 45.43 50.92 59.22 60.52 64.19 62.40 59.36 _
   51.78 44.87 39.81 33.73 29.88 24.75 17.44 13.32 11.14  6.79  4.03  4.11
* The positive error
vec/cre d0xfep(20) R 11.55  6.39  3.65  3.24  2.55  1.89  1.46  1.40 _
    1.46  1.33  2.32  1.40  1.49  0.95  2.73  2.08  1.80  1.39  1.21  0.94
* The negative error
vec/cre d0xfem(20) R 13.74  8.19  2.93  2.15  1.95  1.43  1.31  1.42 _
    1.32  1.28  1.05  1.06  1.02  1.08  2.62  2.33  2.12  2.69  1.65  2.03

* N.B. The errors do not include a +10-11% normalization uncertainty.

opt utit
title 'x?F!' 'U'
title ' '
opt ndat
opt nsta
set hwid 4
set fwid 4
set lwid 4
set pmci 1
set hcol 1
set *fon 130
set tsiz 0.50

* Plotting commands
opt logy
null -0.125 0.8 1.999 80.001
hplot/aerror xfmids d0xf xfdif xfdif d0xfem d0xfep 20 20 0.20 0

* Text printing and symbol labeling
opt liny
null 0. 1. 0. 1. 'sa'
set txfp 130
set txal 20
set chhe 0.45
set tang 90
itx -0.07 0.5 'd[s] / dx?F!      ([m]barn / nucleon)' 
set tang 0
set txal 10
set chhe 0.45
vec/cre xsym(1) R 0.39
vec/cre ysym(1) R 0.96
hplot/symbol xsym ysym 1 20 0.22
itx 0.435 0.951 'E791 D^0! data (K[p] + K[ppp])'