clear all
%
load fpi
hold off
semilogy(ecmOLYA,fpOLYA,'o','markersize',4,'markerFaceColor','black','markerEdgeColor','black')
hold on
semilogy(ecmDM2, fpDM2, 'v','markersize',4,'markerFaceColor','black','markerEdgeColor','black')
semilogy(ecmCMD, fpCMD, 's','markersize',4,'markerEdgeColor','black')
semilogy(ecmPREV,fpPREV,'o','markersize',4,'markerEdgeColor','black')
axis([200 2200 1E-2 1E2])
%
set(gca,'FontSize',15)
xlabel('Center of Mass Energy s^{1/2} (MeV)','fontsize',20)
ylabel('|F_\pi(s)|^2','fontsize',20)
%
plot(1400,40,'o','markersize',4,'markerFaceColor','black','markerEdgeColor','black')
text(1450,40,'OLYA','fontsize',20)
plot(1400,20,'v','markersize',4,'markerFaceColor','black','markerEdgeColor','black')
text(1450,20,'DM2' ,'fontsize',20)
plot(1400,10,'s','markersize',4,'markerEdgeColor','black')
text(1450,10,'CMD' ,'fontsize',20)
plot(1400, 5,'o','markersize',4,'markerEdgeColor','black')
text(1450, 5,'Other','fontsize',20)
%
global mpi
mpi=139.57;
bwpars=[0.6495E6,767.1,134.0,0.9074,0.1744E6,1745,230.6,2.898,0.5495E6,1404,820.6];
%
x=linspace(300.0,2175);
bwall=abs(relbwN(3,x.^2,bwpars)).^2;
plot(x,bwall,'-','color','red')
text(500,2,'Full fit','color','red','fontsize',30)
bw1  =abs(relbwN(1,x.^2,bwpars(1:3))).^2;
plot(x,bw1,'-','color','green')
text(750,1,'\rho(770)','color','green','fontsize',30)
bw2  =abs(relbwN(1,x.^2,bwpars(5:7))).^2;
plot(x,bw2,'-','color','blue')
text(800,0.1,'\rho(1450)','color','black','fontsize',30)
bw3  =abs(relbwN(1,x.^2,bwpars(9:11))).^2;
plot(x,bw3,'-','color','black')
text(1500,0.03,'\rho(1700)','color','blue','fontsize',30)
%
orient landscape
print -depsc fpiData.eps