clear all
%
load gpi1232.dat
ifit=1:21;
idat=22:length(gpi1232);
xfit=gpi1232(ifit,1);
yfit=gpi1232(ifit,2);
xdat=cos(pi*gpi1232(idat,1)/180);
ydat=gpi1232(idat,2);
edat=gpi1232(idat,3);
%
hold off
subplot(1,2,1)
bardat=sum(ydat)/length(ydat);
yA=bardat*(5-3*xfit.^2)./4;
yB=bardat*3*(1+xfit.^2)./4;
plot(xfit,yfit,'-k',xfit,yA,'--k',xfit,yB,':k',xdat,ydat,'ok','MarkerFaceColor','black')
adderror(xdat,ydat,edat)
axis([-1 1 0 25])
legend('Full Fit','M_{1+}','E_{1+}',4)
set(gca,'FontSize',20)
%
xlabel('cos\theta')
ylabel('d\sigma/d\Omega (\mu b/sr)')
text(-0.50,23,'E_\gamma\approx340 MeV','FontSize',20)
%
load gpi1530.dat
ifit=1:21;
idat=22:length(gpi1530);
xfit=gpi1530(ifit,1);
yfit=gpi1530(ifit,2);
xdat=cos(pi*gpi1530(idat,1)/180);
ydat=gpi1530(idat,2);
edat=gpi1530(idat,3);
%
subplot(1,2,2)
bardat=sum(ydat)/length(ydat);
yA=bardat*(5-3*xfit.^2)./4;
yB=bardat*3*(1+xfit.^2)./4;
yC=bardat*ones(length(xfit));
plot(xfit,yfit,'-k',xfit,yA,'--k',xfit,yB,':k',xfit,yC,'-.k',xdat,ydat,'ok','MarkerFaceColor','black')
legend('Full Fit','E_{2-}','M_{2-}','E_{0+}',4)
adderror(xdat,ydat,edat)
axis([-1 1 0 15])
set(gca,'FontSize',20)
%
xlabel('cos\theta')
text(-0.50,14,'E_\gamma\approx780 MeV','FontSize',20)
%
orient portrait
print -deps gamppip.eps