source: ETALON/reconstruction/clab/chiN6/plotsA.m @ 289

Last change on this file since 289 was 289, checked in by hodnevuc, 9 years ago
File size: 1.8 KB
Line 
1% % % figure(1)
2% % % plot(SP,SS,'or')
3% % % hold on
4% % % plot(aFT)
5% % % xlim([0 5000])
6clear all
7close all
8Ndet=[1 2 3 4 7 11 14 20 30 40 50 70];
9s=7
10load(['result/PrN' num2str(1)],'ChiVk','Corrk','ChiVh','Corrh','RPK','RPH','RSP','Profile');
11aFT=abs(fft(Profile(s,:)));
12DetAngl=35+110/(Ndet(s)+1)*(1:Ndet(s));
13
14samp1=0.3./(0.05.*(1-cosd(DetAngl)));
15samp2=0.3./(0.25.*(1-cosd(DetAngl)));
16samp3=0.3./(1.5.*(1-cosd(DetAngl)));
17samp=[samp1 samp2 samp3];
18SFr=sort(samp);
19SP=unique(round(2100.*SFr./max(25.6459)));%ACHTUNG!!! remove same points
20SP1=unique(round(2100.*sort(samp1)./max(25.6459)));
21SP2=unique(round(2100.*sort(samp2)./max(25.6459)));
22SP3=unique(round(2100.*sort(samp3)./max(25.6459)));
23SS=aFT(SP);
24figure(2)
25hold on
26plot(SP1,aFT(SP1),'og')
27plot(SP2,aFT(SP2),'om')
28plot(SP3,aFT(SP3),'oc')
29grid on
30hold on
31plot(aFT)
32plot(RSP(s,:),'k')
33xlim([0 5010])
34%
35% s=1
36% load(['result/PrN' num2str(1)],'ChiVk','Corrk','ChiVh','Corrh','RPK','RPH','RSP','Profile');
37% close all
38% figure(3)
39% plot(Profile(s,:),'k','linew',2)
40% hold on
41% %plot(RPK(s,:),'r','linew',2)
42%
43% plot(RPH(5,:),'r','linew',2)
44% plot(RPH(6,:),'g','linew',2)
45% plot(RPH(10,:),'b','linew',2)
46% grid on
47% xlim([32600 32950])
48% ch=mean(Chi,1);
49%
50% bar(ch(1:9))
51% set(gca,'YGrid','on')  % horizontal grid
52% set(gca,'XTickLabel',{'Lw1', 'Lw5', 'Lw10', 'Lw1X','Lw5X', 'Lw10X', 'Tw10', 'Tw1X', 'Tw5X'},'fontsize',10)
53% colormap([1,0,0]);
54%
55% ylabel('\chi^2')
56%
57% title('Comparison of different sample methods, Hilbert')
58% load SData.mat
59% [Profile,~]=PR_func(4, 2^(-16));
60% aFT=abs(fft(Profile));
61% SFr=sort(Lw1);
62% SP=unique(round(2100.*SFr./max(Tw10)));%ACHTUNG!!! remove same points
63% SS=aFT(SP);
64% figure(1)
65% plot(SP,SS,'or')
66% hold on
67% plot(aFT)
68% xlim([0 5000])
Note: See TracBrowser for help on using the repository browser.