source: ETALON/SPESO/SPESO_Polariz/Analyse5bis.m @ 735

Last change on this file since 735 was 735, checked in by hodnevuc, 7 years ago
File size: 5.9 KB
Line 
1
2clear all; close all; clc
3% plotZ(-1.4e+5)
4% plotZ(-1.45e+5)
5% plotZ(-1.5e+5)
6% plotZ(-1.55e+5)
7% plotZ(-1.6e+5)
8load DATAdel.mat
9BTIME =BID(:,1);
10XBPM = BID(:,2);
11QICT1 =BID(:,3);
12QICT2 =BID(:,4);
13RESEAU_POS =BID(:,19);
14DETECT_XPOS=BID(:,15);
15DETECT_ZPOS=BID(:,17);
16DETECT_SPOS=BID(:,13);
17DETECT_RS1 =BID(:,5);
18DETECT_RS2 =BID(:,7);
19DETECT_RS3 =BID(:,9);
20DETECT_RS4 =BID(:,11);
21LPM =BID(:,21);
22%%% Rotation should be taken into account
23Dshift=+600000;
24H=345520;%distanse from grating to detector. Check if that is true!!!
25
26
27figure
28hold on
29cut=(LPM==0)&(BTIME>736614)&(BTIME<736628)&(AMP(:,3)~=1e+6)&(RESEAU_POS==1.45e+6)&(DETECT_SPOS<-5e+5)&(QICT1>0.2)&(QICT1<0.6)&(DETECT_ZPOS==-1.4e+5);%&(RESEAU_POS<=8000);%
30cut=cut&(~((DETECT_SPOS>-6.68e+5)&(DETECT_SPOS<-6.64e+5)&((AMP(:,1)./(QICT1.^2))<6e-3)));%delete bad data
31
32plot(DETECT_RS3(cut)+90,AMP(cut,2)./(QICT1(cut).^2),'ok','linew',2)
33plot(DETECT_RS3(cut)+90,AMP(cut,1)./(QICT1(cut).^2),'oy','linew',2)
34[Xpos,Ypos2, Yerr]=makemean(DETECT_RS3(cut)+90,AMP(cut,2)./(QICT1(cut).^2));
35plot(Xpos,Ypos2, '-r','linew',2)
36[Xpos,Ypos1, Yerr]=makemean(DETECT_RS3(cut)+90,AMP(cut,1)./(QICT1(cut).^2));
37plot(Xpos,Ypos1, '-b','linew',2)
38
39set(gca, 'fontsize', 16)
40xlabel('Angle, [deg]')
41ylabel('Signal/Q^2, [a.u.]')
42legend('Raw data, H', 'Raw data, V','Mean, H', 'Mean, V','location','northwest')
43title('Z=1.4')
44grid on
45  xlim([77 102])
46ylim([0 0.03])
47
48 TH1=Xpos;
49 DP1=(Ypos2-Ypos1)./(Ypos2+Ypos1);
50
51
52figure
53hold on
54cut=(LPM==0)&(BTIME>736616)&(BTIME<736617.6)&(AMP(:,3)~=1e+6)&(RESEAU_POS==1.45e+6)&(DETECT_SPOS<-5e+5)&(QICT1>0.2)&(QICT1<0.6)&(DETECT_ZPOS==-1.45e+5);%&(RESEAU_POS<=8000);%
55cut=cut&(~((DETECT_SPOS>-8e+5)&(DETECT_SPOS<-6.56e+5)&((AMP(:,1)./(QICT1.^2))<6e-3)));%delete bad data
56plot(DETECT_RS3(cut)+90,AMP(cut,2)./(QICT1(cut).^2),'ok','linew',2)
57plot(DETECT_RS3(cut)+90,AMP(cut,1)./(QICT1(cut).^2),'oy','linew',2)
58[Xpos,Ypos2, Yerr]=makemean(DETECT_RS3(cut)+90,AMP(cut,2)./(QICT1(cut).^2));
59plot(Xpos,Ypos2, '-r','linew',2)
60[Xpos,Ypos1, Yerr]=makemean(DETECT_RS3(cut)+90,AMP(cut,1)./(QICT1(cut).^2));
61plot(Xpos,Ypos1, '-b','linew',2)
62
63set(gca, 'fontsize', 16)
64xlabel('Angle, [deg]')
65ylabel('Signal/Q^2, [a.u.]')
66legend('Raw data, H', 'Raw data, V','Mean, H', 'Mean, V','location','northwest')
67title('Z=1.45')
68grid on
69  xlim([77 102])
70ylim([0 0.03])
71
72 TH2=Xpos;
73 DP2=(Ypos2-Ypos1)./(Ypos2+Ypos1);
74
75
76figure
77hold on
78 cut=(LPM==0)&(BTIME>736614)&(BTIME<736628)&(AMP(:,3)~=1e+6)&(RESEAU_POS==1.45e+6)&(DETECT_SPOS<-5e+5)&(QICT1>0.2)&(QICT1<0.6)&(DETECT_ZPOS==-1.5e+5);%&(RESEAU_POS<=8000);%
79cut=cut&(~((DETECT_SPOS>-8e+5)&(DETECT_SPOS<-6.56e+5)&((AMP(:,1)./(QICT1.^2))<7e-3)));%delete bad data
80plot(DETECT_RS3(cut)+90,AMP(cut,2)./(QICT1(cut).^2),'ok','linew',2)
81plot(DETECT_RS3(cut)+90,AMP(cut,1)./(QICT1(cut).^2),'oy','linew',2)
82[Xpos,Ypos2, Yerr]=makemean(DETECT_RS3(cut)+90,AMP(cut,2)./(QICT1(cut).^2));
83plot(Xpos,Ypos2, '-r','linew',2)
84[Xpos,Ypos1, Yerr]=makemean(DETECT_RS3(cut)+90,AMP(cut,1)./(QICT1(cut).^2));
85plot(Xpos,Ypos1, '-b','linew',2)
86
87set(gca, 'fontsize', 16)
88xlabel('Angle, [deg]')
89ylabel('Signal/Q^2, [a.u.]')
90legend('Raw data, H', 'Raw data, V','Mean, H', 'Mean, V','location','northwest')
91title('Z=1.5')
92grid on
93  xlim([77 102])
94ylim([0 0.03])
95
96 TH3=Xpos;
97 DP3=(Ypos2-Ypos1)./(Ypos2+Ypos1);
98
99
100figure
101hold on
102 cut=(LPM==0)&(BTIME>736614)&(BTIME<736628)&(AMP(:,3)~=1e+6)&(RESEAU_POS==1.45e+6)&(DETECT_SPOS<-5e+5)&(QICT1>0.2)&(QICT1<0.6)&(DETECT_ZPOS==-1.55e+5);%&(RESEAU_POS<=8000);%
103cut=cut&(~((DETECT_SPOS>-8e+5)&(DETECT_SPOS<-6.56e+5)&((AMP(:,1)./(QICT1.^2))<7e-3)));%delete bad data
104plot(DETECT_RS3(cut)+90,AMP(cut,2)./(QICT1(cut).^2),'ok','linew',2)
105plot(DETECT_RS3(cut)+90,AMP(cut,1)./(QICT1(cut).^2),'oy','linew',2)
106[Xpos,Ypos2, Yerr]=makemean(DETECT_RS3(cut)+90,AMP(cut,2)./(QICT1(cut).^2));
107plot(Xpos,Ypos2, '-r','linew',2)
108[Xpos,Ypos1, Yerr]=makemean(DETECT_RS3(cut)+90,AMP(cut,1)./(QICT1(cut).^2));
109plot(Xpos,Ypos1, '-b','linew',2)
110
111set(gca, 'fontsize', 16)
112xlabel('Angle, [deg]')
113ylabel('Signal/Q^2, [a.u.]')
114legend('Raw data, H', 'Raw data, V','Mean, H', 'Mean, V','location','northwest')
115title('Z=1.55')
116grid on
117  xlim([77 102])
118ylim([0 0.03])
119
120 TH4=Xpos;
121 DP4=(Ypos2-Ypos1)./(Ypos2+Ypos1);
122
123
124figure
125hold on
126 cut=(LPM==0)&(BTIME>736614)&(BTIME<736628)&(AMP(:,3)~=1e+6)&(RESEAU_POS==1.45e+6)&(DETECT_SPOS<-5e+5)&(QICT1>0.2)&(QICT1<0.6)&(DETECT_ZPOS==-1.6e+5);%&(RESEAU_POS<=8000);%
127cut=cut&(~((DETECT_SPOS>-8e+5)&(DETECT_SPOS<-6.56e+5)&((AMP(:,1)./(QICT1.^2))<7e-3)));%delete bad data
128plot(DETECT_RS3(cut)+90,AMP(cut,2)./(QICT1(cut).^2),'ok','linew',2)
129plot(DETECT_RS3(cut)+90,AMP(cut,1)./(QICT1(cut).^2),'oy','linew',2)
130[Xpos,Ypos2, Yerr]=makemean(DETECT_RS3(cut)+90,AMP(cut,2)./(QICT1(cut).^2));
131plot(Xpos,Ypos2, '-r','linew',2)
132[Xpos,Ypos1, Yerr]=makemean(DETECT_RS3(cut)+90,AMP(cut,1)./(QICT1(cut).^2));
133plot(Xpos,Ypos1, '-b','linew',2)
134
135set(gca, 'fontsize', 16)
136xlabel('Angle, [deg]')
137ylabel('Signal/Q^2, [a.u.]')
138legend('Raw data, H', 'Raw data, V','Mean, H', 'Mean, V','location','northwest')
139title('Z=1.6')
140grid on
141  xlim([77 102])
142ylim([0 0.03])
143
144 TH5=Xpos;
145 DP5=(Ypos2-Ypos1)./(Ypos2+Ypos1);
146 
147
148 
149 figure
150hold on
151plot(TH1,-DP1, '-.oy','linew',2)
152plot(TH2,-DP2, '-.ok','linew',2)
153plot(TH3,-DP3, '-.or','linew',2)
154plot(TH4,-DP4, '-.ob','linew',2)
155plot(TH5,-DP5, '-.og','linew',2)
156
157set(gca, 'fontsize', 16)
158xlabel('Angle, [deg]')
159ylabel('Degree of polarisation, [1]')
160 legend('Z=1.4','Z=1.45','Z=1.5','Z=1.55','Z=1.6','location','southwest')
161% title('Q_e=0.45\pm0.15 [nC]')
162grid on
163  xlim([73 102])
164  ylim([-1 1])
165
166% load polariz.mat TH1b DP1b
167%  figure
168% hold on
169% plot(TH1,-DP1, '-.oy','linew',2)
170% plot(TH2,-DP2, '-.ok','linew',2)
171% plot(TH3,-DP3, '-.or','linew',2)
172% plot(TH4,-DP4, '-.ob','linew',2)
173% plot(TH5,-DP5, '-.og','linew',2)
174% plot(TH1b,DP1b, '-.*k','linew',2)
175%
176% set(gca, 'fontsize', 16)
177% xlabel('Angle, [deg]')
178% ylabel('Degree of polarisation, [1]')
179%  legend('Z=1.4','Z=1.45','Z=1.5','Z=1.55','Z=1.6','Z=1.45*','location','southwest')
180% % title('Q_e=0.45\pm0.15 [nC]')
181% grid on
182%   xlim([73 102])
183%   ylim([-1 1])
Note: See TracBrowser for help on using the repository browser.