* Fits e+e- -> pi+pi- data to three interfering Breit-Wigners
*
* Setup hardcopy output
fort/file 66 fit3rhos.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 ecm,fpisq,error    fitrho.dat
vec/read ecmx,fpisqx,errorx fitrhox.dat
vec/create zeros(500)
sigma nfit =NCO(ecm)
sigma nfitx=NCO(ecmx)
*
* Set up the plot
opt nbox
opt nfil
opt date
opt soft
opt fit
set xsiz 20.0
set ysiz 26.0
opt nsta
opt logy
set gsiz 0.50
set ygti 0.75
title 'F?[p]^2! and Three [r] Resonances'
*
* Fit the data.  Starting values from Bisello, et al, Phys Lett B220(1989)321
*
vec/create pars(11) R 7E5 770 150 0 2E5 1770 224 3.2 2.4E5 1424 269
hplot/null 200 2200 1E-2 1E+2
vec/fit ecm fpisq error relbw3.f S  11 pars
hplot/errors ecm  fpisq  zeros error  nfit  20 0.15
hplot/errors ecmx fpisqx zeros errorx nfitx 27 0.15
set hcol 1
function/plot relbw3.f 290.0 2200.0 CS
atitle 'Center of Mass Energy (MeV)' 'Square of Timelike Pion Form Factor'
key 500 0.5 20 'Fitted data'
key 500 0.3 27 'Other data'
*
* Cleanup
close 66