% Data format: angle(deg) count(in 60s)

data = [  0 15420;

          2 14806;

          4 10432;

          6 3420;

          8 501;

         -2 15424;

         -4 12841;

         -6 9054;

         -8 3484;

         -10 111/3];

errorbar(data(:,1),data(:,2),sqrt(data(:,2)),'o',...

            'markersize',2);

xlim([-12 12]);



dataold = [ 0 15202;

            4 13136;

           -3 15162;

           -5 12091;

           -6 10026;

           -8 4215;

          -10 131];

      

hold on;

errorbar(dataold(:,1),dataold(:,2),sqrt(dataold(:,2)),'rx',...

            'markersize',2);

xlabel('Incident angle');

ylabel('Count at detector for 60s');

title('Beam profile calibration: blue (10/1) red (9/26)');

grid on;