source: ETALON/CLIO/Analysis/Data24May/WideScan.m @ 719

clear all; close all; clc;
date=20170524;
load('data2405.mat');
%cut by electron signal
%cut on signal on 12's detector (zero signal); and 11's detector (DAQ problem) 
cut=(key_value(:,8)>1e+3)&(abs(data(:,11,2))<100)&(abs(data(:,12,2))<50);%&(sum(data(:,:,2)~=0,2));

position=key_value(:,4)*7.94e-6;
hour=rem(key_value(:,6),date);
diapW814=(hour>81800)&(hour<84200);
diapW906=(hour>121000)&(hour<123500);

for c=1:11
[Xpos,Ypos(c,:), ~]=makemean(position(cut&diapW814)-48.58-3.46,data(cut&diapW814,c,2))
end
figure(123)
% [hc hc]=contourf(data(cut&diapW814,:,2)./repmat(key_value(cut&diapW814,8)*1000,1,12))
[hc hc]=contourf(300./(6*(1-cos((48:7:118)/180*pi))),Xpos,Ypos',40)
colorbar;
set(gca,'fontsize',16)
set(hc,'LineStyle','none');
title('\phi_B=8.14')
xlabel('Frequency, [GHz]')
ylabel('BGS, [mm]')

clear Xpos Ypos
for c=1:11
[Xpos,Ypos(c,:), ~]=makemean(position(cut&diapW906)-48.58-3.46,data(cut&diapW906,c,2))
end
figure(321)
% [hc hc]=contourf(data(cut&diapW814,:,2)./repmat(key_value(cut&diapW814,8)*1000,1,12))
[hc hc]=contourf(300./(6*(1-cos((48:7:118)/180*pi))),Xpos,Ypos',40)
colorbar;
set(gca,'fontsize',16)
set(hc,'LineStyle','none');
title('\phi_B=9.06')
xlabel('Frequency, [GHz]')
ylabel('BGS, [mm]')


% figure(1)
% semilogy(position(cut&diapW814)-48.58-3.46,data(cut&diapW814,7,2)./key_value(cut&diapW814,8)*1000,'ok','MarkerFaceColor','k')
% hold on 
% semilogy(position(cut&diapW906)-48.58-3.46,data(cut&diapW906,7,2)./key_value(cut&diapW906,8)*1000,'^r','MarkerFaceColor','r')
% set(gca,'fontsize',16)
% legend('Location','northwest')
% legend('\phi_B=8.14','\phi_B=9.06')
% grid on
% xlabel('Beam-grating separation, [mm]')
% ylabel('Amplitude, [a.u.]')
% ylim([6 25])
% xlim([-55 0])
% title('Detector 90^o')
% 
% 
% 
% figure(11)
% [Xpos,Ypos, Yerr]=makemean(position(cut&diapW814)-48.58-3.46,(data(cut&diapW814,7,2)./key_value(cut&diapW814,8)*1000));
% set(gca,'YScale','log');
% errorbar(Xpos,Ypos, Yerr,'ok','MarkerFaceColor','k')
% hold on 
% [Xpos,Ypos, Yerr]=makemean(position(cut&diapW906)-48.58-3.46,(data(cut&diapW906,7,2)./key_value(cut&diapW906,8)*1000));
% errorbar(Xpos,Ypos, Yerr,'^r','MarkerFaceColor','r')
% set(gca,'YScale','log');
% set(gca,'fontsize',16)
% legend('Location','northwest')
% legend('\phi_B=8.14','\phi_B=9.06')
% grid on
% xlabel('Beam-grating separation, [mm]')
% ylabel('Amplitude, [a.u.]')
% ylim([5 25])
% xlim([-55 0])
% title('Detector 90^o')
% 
% 
% figure(2)
% semilogy(position(cut&diapW814)-48.58-3.46,data(cut&diapW814,11,2)./key_value(cut&diapW814,8)*1000,'ok','MarkerFaceColor','k')
% hold on 
% semilogy(position(cut&diapW906)-48.58-3.46,data(cut&diapW906,11,2)./key_value(cut&diapW906,8)*1000,'^r','MarkerFaceColor','r')
% set(gca,'fontsize',16)
% legend('Location','northwest')
% legend('\phi_B=8.14','\phi_B=9.06')
% grid on
% xlabel('Beam-grating separation, [mm]')
% ylabel('Amplitude, [a.u.]')
% ylim([5 25])
% xlim([-55 0])
% title('Detector 118^o')
% 
% 
% 
% figure(21)
% [Xpos,Ypos, Yerr]=makemean(position(cut&diapW814)-48.58-3.46,(data(cut&diapW814,11,2)./key_value(cut&diapW814,8)*1000));
% set(gca,'YScale','log');
% errorbar(Xpos,Ypos, Yerr,'ok','MarkerFaceColor','k')
% hold on 
% [Xpos,Ypos, Yerr]=makemean(position(cut&diapW906)-48.58-3.46,(data(cut&diapW906,11,2)./key_value(cut&diapW906,8)*1000));
% errorbar(Xpos,Ypos, Yerr,'^r','MarkerFaceColor','r')
% set(gca,'YScale','log');
% set(gca,'fontsize',16)
% legend('Location','northwest')
% legend('\phi_B=8.14','\phi_B=9.06')
% grid on
% xlabel('Beam-grating separation, [mm]')
% ylabel('Amplitude, [a.u.]')
% ylim([5 25])
% xlim([-55 0])
% title('Detector 118^o')
% 
% 
% figure(3)
% semilogy(position(cut&diapW814)-48.58-3.46,data(cut&diapW814,1,2)./key_value(cut&diapW814,8)*1000,'ok','MarkerFaceColor','k')
% hold on 
% semilogy(position(cut&diapW906)-48.58-3.46,data(cut&diapW906,1,2)./key_value(cut&diapW906,8)*1000,'^r','MarkerFaceColor','r')
% set(gca,'fontsize',16)
% legend('Location','northwest')
% legend('\phi_B=8.14','\phi_B=9.06')
% grid on
% xlabel('Beam-grating separation, [mm]')
% ylabel('Amplitude, [a.u.]')
% ylim([5 25])
% xlim([-55 0])
% title('Detector 48^o')
% 
% 
% 
% figure(31)
% [Xpos,Ypos, Yerr]=makemean(position(cut&diapW814)-48.58-3.46,(data(cut&diapW814,1,2)./key_value(cut&diapW814,8)*1000));
% set(gca,'YScale','log');
% errorbar(Xpos,Ypos, Yerr,'ok','MarkerFaceColor','k')
% hold on 
% [Xpos,Ypos, Yerr]=makemean(position(cut&diapW906)-48.58-3.46,(data(cut&diapW906,1,2)./key_value(cut&diapW906,8)*1000));
% errorbar(Xpos,Ypos, Yerr,'^r','MarkerFaceColor','r')
% set(gca,'YScale','log');
% set(gca,'fontsize',16)
% legend('Location','northwest')
% legend('\phi_B=8.14','\phi_B=9.06')
% grid on
% xlabel('Beam-grating separation, [mm]')
% ylabel('Amplitude, [a.u.]')
% ylim([5 25])
% xlim([-55 0])
% title('Detector 48^o')