* Fits g+p->pi+ n data on Delta to E1+/M1+
*
* Setup hardcopy output
fort/file 66 gpidelta.ps
metafile  66 -111
*
* Read in the data and create a list of "zeros"
* The data comes in two flavors: "Good" and "x='Other'"
vec/read theta,sig,error    gpidelta.dat
vec/create zeros(500)
sigma nfit =NCO(theta)
sigma xfit =COS(theta*3.14159/180)
*
* Set up the plot
opt nbox
opt nfil
opt date
opt soft
opt fit
set xsiz 20.0
set ysiz 26.0
opt nsta
set gsiz 0.50
set ygti 0.75
title '[g]p"5#n[p]^+! on the [D](1232)'
*
* Fit the data.
*
vec/create pars(3) R 7.5 0.0 0.0
hplot/null -1 1 0 25
vec/fit xfit sig error emrat.f SE 3 pars
hplot/errors xfit sig  zeros error  nfit  20 0.15
set hcol 1
function/plot emrat.f -1.0 1.0 CS
atitle 'cos([q])' 'd[s]/d[W] (mb/sr)'
*
* Cleanup
close 66