source: ETALON/CLIO/2017_CLIO_report/scripts/Grating_correction.m @ 728

clear all; close all;clc
alpha=45/180*pi;
L=180;
D=25;
XposA=-25:1:25;
alphaA=(-42:7:35)/180*pi;
for c=1:51
    for k=1:12
    Xpos=XposA(c);
    alpha=alphaA(k);
x=(-L/2+Xpos*tan(alpha)):(L/1000):(L/2+Xpos*tan(alpha));
Norm(c,k)=sum(real(8*sqrt((D/2)^2-x.^2*cos(alpha)^2))/(pi*D^2)*(x(2)-x(1)));
    end
end
surf(XposA',alphaA'./pi*180+90,Norm')
ylim([40 140])
xlabel('Xpos, [mm]')
ylabel('Angle, [deg]')