* Fits e+e- -> pi+pi- data to three interfering Breit-Wigners
*
* Setup hardcopy output
fort/file 66 fitRho3.eps
metafile  66 -113
*
* 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 ndat
opt soft
opt fit
set xsiz 20.0
set ysiz 26.0
opt nsta
opt logy
set gsiz 0.50
set ygti 0.75
*
* 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
key 500 0.5 20 'Fitted data'
key 500 0.3 27 'Other data'
*
* Cleanup
close 66