1 | function []=checkparam(cut,plt) |
---|
2 | |
---|
3 | load DATAdel.mat |
---|
4 | BTIME =BID(:,1); |
---|
5 | XBPM = BID(:,2); |
---|
6 | QICT1 =BID(:,3); |
---|
7 | QICT2 =BID(:,4); |
---|
8 | RESEAU_POS =BID(:,19); |
---|
9 | DETECT_XPOS=BID(:,15); |
---|
10 | DETECT_ZPOS=BID(:,17); |
---|
11 | DETECT_SPOS=BID(:,13); |
---|
12 | DETECT_RS1 =BID(:,5); |
---|
13 | DETECT_RS2 =BID(:,7); |
---|
14 | DETECT_RS3 =BID(:,9); |
---|
15 | DETECT_RS4 =BID(:,11); |
---|
16 | LPM =BID(:,21); |
---|
17 | |
---|
18 | if (plt) |
---|
19 | |
---|
20 | |
---|
21 | figure(1001) |
---|
22 | hold on |
---|
23 | plot(BTIME(cut),RESEAU_POS(cut),'.r') |
---|
24 | title('RESEAU_POS') |
---|
25 | |
---|
26 | figure(1002) |
---|
27 | hold on |
---|
28 | plot(BTIME(cut),QICT1(cut),'.r') |
---|
29 | plot(BTIME(cut),QICT2(cut),'.k') |
---|
30 | set(gca,'fontsize',16) |
---|
31 | legend('QICT1', 'QICT2') |
---|
32 | grid on |
---|
33 | xlabel('Time') |
---|
34 | ylabel('Charge, [nC?]') |
---|
35 | |
---|
36 | %%How to deal with turns??? |
---|
37 | figure(1003) |
---|
38 | hold on |
---|
39 | plot(BTIME(cut),DETECT_RS1(cut),'.r')%2 |
---|
40 | plot(BTIME(cut),DETECT_RS2(cut),'.k')% 0 -3 |
---|
41 | plot(BTIME(cut),DETECT_RS3(cut),'.b') |
---|
42 | plot(BTIME(cut),DETECT_RS4(cut),'.g')%-45 |
---|
43 | legend('RS1','RS2','RS3','RS4') |
---|
44 | |
---|
45 | figure(1004) |
---|
46 | hold on |
---|
47 | plot(BTIME(cut),DETECT_XPOS(cut),'.r') |
---|
48 | plot(BTIME(cut),DETECT_ZPOS(cut),'.k') |
---|
49 | plot(BTIME(cut),DETECT_SPOS(cut),'.b') |
---|
50 | |
---|
51 | legend('xpos','zpos','spos') |
---|
52 | |
---|
53 | figure(1005) |
---|
54 | hold on |
---|
55 | plot(BTIME(cut),LPM(cut),'.r') |
---|
56 | title('LPM') |
---|
57 | end |
---|
58 | if mean(LPM(cut))==2 |
---|
59 | iword='Long pulse mode, '; |
---|
60 | else |
---|
61 | iword='Short pulse mode, '; |
---|
62 | end |
---|
63 | disp([iword, 'grating position = ' num2str(mean(RESEAU_POS(cut))) ... |
---|
64 | ', $QICT_1$ = $' num2str(mean(QICT1(cut))) '\pm' num2str(std(QICT1(cut))) '$' ... |
---|
65 | ', $QICT_2$ = $' num2str(mean(QICT2(cut))) '\pm' num2str(std(QICT2(cut))) '$' ... |
---|
66 | ', $RS_1$ = ' num2str(mean(DETECT_RS1(cut))) ', $RS_2$ = ' num2str(mean(DETECT_RS2(cut))) ... |
---|
67 | ', $RS_3$ = ' num2str(mean(DETECT_RS3(cut))) ', $RS_4$ = ' num2str(mean(DETECT_RS4(cut))) ... |
---|
68 | ', XPOS = ' num2str(mean(DETECT_XPOS(cut))) ', ZPOS = ' num2str(mean(DETECT_ZPOS(cut))) ', SPOS = ' num2str(mean(DETECT_SPOS(cut))) ... |
---|
69 | ]) |
---|
70 | |
---|
71 | |
---|