Changeset 614 in ETALON

Timestamp:

May 13, 2016, 5:44:43 PM (8 years ago)

Author:

malovyts

Message:

Changed spectrum norlmalizing script

Location:

SPESO/ana2015/offline_analysis

Files:

• offline_analyse_multi.m (modified) (1 diff)
• offline_normalize_spectrum.m (modified) (7 diffs)

SPESO/ana2015/offline_analysis/offline_analyse_multi.m

@@ -8 +8 @@
 % data_selection==2 after 2016/03/10 before 2016/04/09
 
+% 2015
+offline_analysis('2015', '2015????.dayFullInfo', 1001); % histo
+
 % 2016 03/09 -- 01/01 
-offline_analysis('2016JanFeb', '20160{30?,2[^0]?}.dayFullInfo', 1001); % histo
+%offline_analysis('2016JanFeb', '20160{30?,2[^0]?}.dayFullInfo', 1001); % histo
 %offline_analysis('2016JanFeb', '20160{30?,2??,1??}.dayFullInfo', 2001); % colorplots
 %offline_analysis('2016JanFebMarch', '20160{3[^1]?,31[^3],2??,1??}.dayFullInfo', 3002); % reseau analysis

SPESO/ana2015/offline_analysis/offline_normalize_spectrum.m

@@ -169 +169 @@
 nu_comb=cell(length(Q)/2,1);
 
-det1_idx=15;
-det2_idx=16;
+det1_idx=11;
+det2_idx=12;
 
 bins_elem_num{det1_idx}
@@ -193 +193 @@
 
 % detector+antenna corections
-% det sensitivity V/W, Boresight efficiency, surface of the waveguide
-
-det1_coeff=1/1800/25.2./exp(2.3);%./exp((Det1_ant_eff-max(Det1_ant_eff))/10); % W/mm^2
-det2_coeff=1/1200/16./exp(2.5);%./exp((Det2_ant_eff-max(Det2_ant_eff))/10); % W/mm^2
+% det sensitivity V/W, Boresight efficiency, surface of the waveguide power
+%  
+
+det1_coeff=1/1800/25.2./exp(Det1_ant_eff/10)*100; % W/mm^2
+det2_coeff=1/1200/16./exp(Det2_ant_eff/10)*100; % W/mm^2
 
 Det1_val_raw=Det1_val_raw.*det1_coeff;
@@ -234 +235 @@
 %Det1_th=Det1_th(9:end);
 %Det2_th=Det2_th(1:end);
-%Det1_val_raw=Det1_val_raw(9:end);
+Det1_val_raw=Det1_val_raw(9:end);
 %Det2_val_raw=Det2_val_raw(1:end);
-%Det1_th=Det1_th(9:end);
+Det1_th=Det1_th(9:end);
 %Det2_th=Det2_th(1:end);
 
         % theory coeff
         % power: N_e*J/Sr/signal_time/R^2*far_correction = W/mm^2
-        coeff=104/260*10^9/310^2*data_corrFactor(10, sp_theta, 9, r0{7})...
-                *(1-exp(-50*10^6*40*10^-12)); % capacitance correction
+        coeff=104/260*10^9/310^2*data_corrFactor(10, sp_theta, 9, r0{7});%...
+                %*(1-exp(-50*10^6*40*10^-12)); % capacitance correction
                 
         fih=figure(1);
@@ -248 +249 @@
         xlim([40 140]);
         hold on
-        plot(Det2_th, Det2_val_raw, '--', 'LineWidth', 2, 'Color', [0.19 0.55 0.91]);
+        plot(Det2_th-3, Det2_val_raw, '--', 'LineWidth', 2, 'Color', [0.19 0.55 0.91]);
    %plot(Det1_th, Det1_val, 'LineWidth', 2, 'Color', [0 0.75 0.24]);
         %plot(Det2_th, Det2_val*3, 'LineWidth', 2, 'Color', [0 0 1]);
         %plot(sp_theta, sp_val_sum_RRR_far.*4*10^2, 'LineWidth', 2, 'Color', [1 1 0]);
-%{
+
         load('SC_RRR_10000fs.mat')
         %sp_val_sum_SC=sp_val_sum_SC(1,1:15).*exp(d_eff(1,1:15)/10);
@@ -289 +290 @@
         %                                       data_corrFactor(10, sp_theta(1,1:15), 9, r0{7});
         plot(sp_theta, sp_val_sum_SC, '-^', 'LineWidth', 1, 'Color', [0 0 0]);
-%}
+
         hold off
         legend('Det1 raw', 'Det2 raw',...
@@ -297 +298 @@
         xlabel('theta, degree');
         ylabel('magnitude, W/mm^2')
-        title('SPM Data. R=310mm')
+        title('LPM Data. R=310mm')
    set(findall(gcf,'type','text'),'FontSize',20)%'fontWeight','bold')
     set(gca,'fontsize',20)
@@ -314 +315 @@
 
 
+        load('SC_RRR_10000fs.mat')
 
 save('NormalizedSpectrum_20160421.mat',...