1 | close all; clear all; clc; |
---|
2 | load analyzeSP.mat |
---|
3 | |
---|
4 | |
---|
5 | %% Spectrum |
---|
6 | |
---|
7 | |
---|
8 | separation=-10; |
---|
9 | |
---|
10 | |
---|
11 | figure(9) |
---|
12 | hold on |
---|
13 | Spectrum=amp(6,:).*exp(lamb(6,:)*separation); |
---|
14 | errorbar(48:7:118,Spectrum,Spectrum.*sqrt((Ierrp(6,:,1)./amp(6,:)).^2+(separation*Ierrp(6,:,2)).^2),'k','linewidth',2) |
---|
15 | |
---|
16 | % plot(48:7:118, SP_spec_1(1:2:(end-1))*5e+16,'k','linewidth',2) |
---|
17 | load('GFW/Spec3.mat') |
---|
18 | plot(48:7:118, Spectr(1:11)./max(Spectr(1:11))*max(Spectrum),'--or','linewidth',2) |
---|
19 | load('GFW/Spec.mat') |
---|
20 | plot(48:7:118, Spectr(1:11)./max(Spectr(1:11))*max(Spectrum),'-^r','linewidth',2) |
---|
21 | load('GFW/Spec7.mat') |
---|
22 | plot(48:7:118, Spectr(1:11)./max(Spectr(1:11))*max(Spectrum),'-.*r','linewidth',2) |
---|
23 | |
---|
24 | set(gca,'fontsize',16) |
---|
25 | |
---|
26 | legend('Data','3ps','5ps','7ps','location','southeast') |
---|
27 | xlabel('Angle, [deg.]') |
---|
28 | ylabel('Amplitude, [a.u.]') |
---|
29 | title('Spectrum') |
---|
30 | grid on |
---|
31 | ylim([0 0.02]) |
---|
32 | load('GFW/SEY.mat') |
---|
33 | |
---|
34 | |
---|
35 | figure(10) |
---|
36 | [hc hc]=contourf (48:7:118, Bphase, amp.*exp(lamb.*separation)) |
---|
37 | set(hc,'LineStyle','none'); |
---|
38 | |
---|
39 | set(gca,'fontsize',16) |
---|
40 | xlabel('Angle, [deg.]') |
---|
41 | ylabel('\phi_S, [a.u.]') |
---|
42 | title('Spectrum map') |
---|
43 | |
---|
44 | |
---|
45 | load('GFW/SEYP5.mat') |
---|
46 | |
---|
47 | figure(10341) |
---|
48 | [hc hc]=contourf(300./(6*(1-cos((48:7:118)/180*pi))), Bphase, amp.*exp(lamb.*(-5))./repmat((Integral(1:11)),[11 1])*1e-17,40) |
---|
49 | colorbar |
---|
50 | set(hc,'LineStyle','none'); |
---|
51 | set(gca,'fontsize',16) |
---|
52 | xlabel('Frequency, [GHz]') |
---|
53 | ylabel('\phi_S, [a.u.]') |
---|
54 | title('Form Factor (BGS=5 mm)') |
---|
55 | |
---|
56 | |
---|
57 | load('GFW/SEY.mat') |
---|
58 | |
---|
59 | figure(10342) |
---|
60 | [hc hc]=contourf(300./(6*(1-cos((48:7:118)/180*pi))), Bphase, amp.*exp(lamb.*(-10))./repmat((Integral(1:11)),[11 1])*1e-17,40) |
---|
61 | colorbar |
---|
62 | set(hc,'LineStyle','none'); |
---|
63 | set(gca,'fontsize',16) |
---|
64 | xlabel('Frequency, [GHz]') |
---|
65 | ylabel('\phi_S, [a.u.]') |
---|
66 | title('Form Factor (BGS=10 mm)') |
---|
67 | |
---|
68 | |
---|
69 | load('GFW/SEYP15.mat') |
---|
70 | |
---|
71 | figure(10343) |
---|
72 | [hc hc]=contourf(300./(6*(1-cos((48:7:118)/180*pi))), Bphase, amp.*exp(lamb.*(-15))./repmat((Integral(1:11)),[11 1])*1e-17,40) |
---|
73 | colorbar |
---|
74 | set(hc,'LineStyle','none'); |
---|
75 | set(gca,'fontsize',16) |
---|
76 | xlabel('Frequency, [GHz]') |
---|
77 | ylabel('\phi_S, [a.u.]') |
---|
78 | title('Form Factor (BGS=15 mm)') |
---|
79 | % |
---|
80 | % figure(11) |
---|
81 | % [hc hc]=contourf(300./(6*(1-cos((48:7:118)/180*pi))), Bphase, amp.*exp(lamb.*separation),40) |
---|
82 | % |
---|
83 | % colorbar |
---|
84 | % set(hc,'LineStyle','none'); |
---|
85 | % set(gca,'fontsize',16) |
---|
86 | % xlabel('Frequency, [GHz]') |
---|
87 | % ylabel('\phi_B, [a.u.]') |
---|
88 | % title('Spectrum without background') |
---|
89 | % |
---|
90 | % |
---|
91 | |
---|
92 | figure(12) |
---|
93 | c=6; |
---|
94 | hold on |
---|
95 | load('GFW/SEYP5.mat') |
---|
96 | plot(300./(6*(1-cos((48:7:118)/180*pi))), ((amp(c,:).*exp(lamb(c,:)*separation))./Integral(1:11)*1e-25),'-^k','linewidth',2) |
---|
97 | load('GFW/SEY.mat') |
---|
98 | plot(300./(6*(1-cos((48:7:118)/180*pi))), ((amp(c,:).*exp(lamb(c,:)*separation))./Integral(1:11)*1e-25),'-^r','linewidth',2) |
---|
99 | load('GFW/SEYP15.mat') |
---|
100 | plot(300./(6*(1-cos((48:7:118)/180*pi))), ((amp(c,:).*exp(lamb(c,:)*separation))./Integral(1:11)*1e-25),'-^b','linewidth',2) |
---|
101 | set(gca,'fontsize',16) |
---|
102 | legend('BGS=5mm','BGS=10mm','BGS=15mm','Location','northwest') |
---|
103 | xlabel('Frequency, [GHz]') |
---|
104 | ylabel('Amplitude, [a.u.]') |
---|
105 | title('FormFactor') |
---|
106 | grid on |
---|
107 | |
---|
108 | load('GFW/SEY.mat') |
---|
109 | figure |
---|
110 | hold on |
---|
111 | plot(300./(6*(1-cos((48:7:118)/180*pi))), amp(2,:).*exp(lamb(2,:).*(-10))./(Integral(1:11))*1e-24,'-^b','linewidth',2) |
---|
112 | plot(300./(6*(1-cos((48:7:118)/180*pi))), amp(4,:).*exp(lamb(4,:).*(-10))./(Integral(1:11))*1e-24,'-^m','linewidth',2) |
---|
113 | plot(300./(6*(1-cos((48:7:118)/180*pi))), amp(6,:).*exp(lamb(6,:).*(-10))./(Integral(1:11))*1e-24,'-^k','linewidth',2) |
---|
114 | plot(300./(6*(1-cos((48:7:118)/180*pi))), amp(8,:).*exp(lamb(8,:).*(-10))./(Integral(1:11))*1e-24,'-^g','linewidth',2) |
---|
115 | plot(300./(6*(1-cos((48:7:118)/180*pi))), amp(10,:).*exp(lamb(10,:).*(-10))./(Integral(1:11))*1e-24,'-^r','linewidth',2) |
---|
116 | |
---|
117 | set(gca,'fontsize',16) |
---|
118 | xlabel('Frequency, [GHz]') |
---|
119 | legend('\phi_S=235','\phi_S=245','\phi_S=255','\phi_S=265','\phi_S=275') |
---|
120 | title('Form Factor (BGS=10 mm)') |
---|
121 | ylabel('Amplitude, [a.u.]') |
---|
122 | grid on |
---|
123 | Freq=300./(6*(1-cos((48:7:118)/180*pi))); |
---|
124 | for c=1:11 |
---|
125 | FF(c,:)=amp(c,:).*exp(lamb(c,:).*(-10))./(Integral(1:11))*1e-24; |
---|
126 | end |
---|
127 | save FFXSP.mat Freq FF |
---|