Line | |
---|
1 | clear all; close all; clc |
---|
2 | load DATAdel.mat |
---|
3 | BTIME =BID(:,1); |
---|
4 | XBPM = BID(:,2); |
---|
5 | QICT1 =BID(:,3); |
---|
6 | QICT2 =BID(:,4); |
---|
7 | RESEAU_POS =BID(:,19); |
---|
8 | DETECT_XPOS=BID(:,15); |
---|
9 | DETECT_ZPOS=BID(:,17); |
---|
10 | DETECT_SPOS=BID(:,13); |
---|
11 | DETECT_RS1 =BID(:,5); |
---|
12 | DETECT_RS2 =BID(:,7); |
---|
13 | DETECT_RS3 =BID(:,9); |
---|
14 | DETECT_RS4 =BID(:,11); |
---|
15 | LPM =BID(:,21); |
---|
16 | %7.36815e+5 to the end |
---|
17 | %7.36756e+5 to 7.36815e+5 |
---|
18 | %7.3655e+5 to 7.36756e+5 |
---|
19 | |
---|
20 | cut=(LPM==0)&(BTIME>7.3655e+5)&(AMP(:,3)~=1e+6)&(RESEAU_POS==1.45e+6);%&(RESEAU_POS<=8000);% |
---|
21 | |
---|
22 | %0.4 0.6 |
---|
23 | %3.1 3.52 |
---|
24 | % AMP=BASE; |
---|
25 | figure |
---|
26 | hold on |
---|
27 | plot(BTIME(cut),QICT1(cut),'.r') |
---|
28 | plot(BTIME(cut),QICT2(cut),'.k') |
---|
29 | set(gca,'fontsize',16) |
---|
30 | legend('QICT1', 'QICT2') |
---|
31 | grid on |
---|
32 | xlabel('Time') |
---|
33 | ylabel('Charge, [nC?]') |
---|
34 | |
---|
35 | %%How to deal with turns??? |
---|
36 | figure |
---|
37 | hold on |
---|
38 | plot(BTIME(cut),DETECT_RS1(cut),'.r')%2 |
---|
39 | plot(BTIME(cut),DETECT_RS2(cut),'.k')% 0 -3 |
---|
40 | plot(BTIME(cut),DETECT_RS3(cut),'.b') |
---|
41 | plot(BTIME(cut),DETECT_RS4(cut),'.g')%-45 |
---|
42 | |
---|
43 | |
---|
44 | figure |
---|
45 | hold on |
---|
46 | plot(BTIME(cut),DETECT_XPOS(cut),'.r') |
---|
47 | plot(BTIME(cut),DETECT_ZPOS(cut),'.k') |
---|
48 | plot(BTIME(cut),DETECT_SPOS(cut),'.b') |
---|
49 | |
---|
50 | figure |
---|
51 | hold on |
---|
52 | plot(BTIME(cut),DEL(cut),'.r') |
---|
53 | |
---|
Note: See
TracBrowser
for help on using the repository browser.