*******************************************************************
*
* Kumac to plot the E791 D^0 production cross section versus pt^2.
* 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 +11-15% normalization uncertainty.
* (5) The data are integrated over xf>0.
* (6) A cross section with a linear dependence on the atomic number
*     is assumed.
*
******************************************************************

* The limits of the pt^2 bins
VEC/CRE PTBINS(21) r_
  0.00  0.10  0.20  0.30  0.40  0.50  0.75 _
  1.00  1.50  2.00  2.50  3.00  4.00  5.00 _
  6.00  8.00  10.0  12.0  14.0  16.0  18.0

* The midpoints of the pt^2 bins.
vec/cre ptmids(20) R _
  0.05  0.15  0.25  0.35  0.45  0.625 0.875 _
  1.25  1.75  2.25  2.75  3.50  4.50  5.50 _
  7.00  9.00  11.0  13.0  15.0  17.0

* The sizes of each pt^2 bin
vec/cre ptdif(20) R 5*0.05 2*0.125 4*0.25 3*0.5 6*1.0

*title 'p^2!&?T! distribution of D^0! (K[p] + K[ppp])'

* The central values
vec/cre d0pt(20) R  13.97   12.20   11.25    9.44    8.887   7.386  _
     5.802   3.888   2.498   1.7780  1.2017  0.7278  0.3672  0.2158 _
     0.09826 0.03563 0.01668 0.00976 0.00539 0.00222
* The positive error
vec/cre d0ptep(20) R 1.67    1.40    1.32    1.09    0.419   0.390  _
     0.244   0.147   0.104   0.0799  0.0702  0.0648  0.0460  0.0317 _
     0.01185 0.00535 0.00409 0.00398 0.00126 0.00188  
* The negative error
vec/cre d0ptem(20) R 1.74    1.56    1.38    1.15    0.466   0.327  _
     0.206   0.227   0.089   0.0640  0.0451  0.0269  0.0162  0.0109 _
     0.00739 0.00425 0.00264 0.00156 0.00303 0.00062

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

opt utit
title '   p^2!&?T!' '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
set ywin 1.7
set yhti 0.60

* Two plots (one for pt^2 < 4 GeV^2 and another for pt^2 < 18 GeV^2)
zone 1 2

* Plot the pt^2 < 4 GeV^2 data
opt logy
null 0.0 4.0 0.29999 20.001
atitle '(GeV/c)^2!      [ ]'
hplot/aerror ptmids d0pt ptdif ptdif d0ptem d0ptep 13 20 0.20 0

* Print text and label symbol
opt liny
null 0. 1. 0. 1. 'sa'
set txfp 130
set txal 20
set chhe 0.30
set tang 90
itx -0.07 0.5 'd[s] / dp^2!&?T!      ([m]barn /  (nucleon  (GeV/c)^2!))' 
set tang 0
set txal 10
set chhe 0.40
vec/cre xsym(1) R 0.42
vec/cre ysym(1) R 0.92
hplot/symbol xsym ysym 1 20 0.20
itx 0.455 0.902 'E791 D^0! data (K[p] + K[ppp])'

* Plot the pt^2 < 18 GeV^2 data
opt logy
null 0.0 18.0 0.0009999 20.001
atitle '(GeV/c)^2!      [ ]'
hplot/aerror ptmids d0pt ptdif ptdif d0ptem d0ptep 20 20 0.20 0

* Print text and label symbol
opt liny
null 0. 1. 0. 1. 'sa'
set txfp 130
set txal 20
set chhe 0.30
set tang 90
itx -0.07 0.5 'd[s] / dp^2!&?T!      ([m]barn /  (nucleon  (GeV/c)^2!))' 
set tang 0
set chhe 0.40
set txal 10
vec/cre xsym(1) R 0.42
vec/cre ysym(1) R 0.92
hplot/symbol xsym ysym 1 20 0.20
itx 0.455 0.902 'E791 D^0! data (K[p] + K[ppp])'