Changeset 613 in ETALON for SPESO

Timestamp:

May 10, 2016, 2:59:31 PM (8 years ago)

Author:

malovyts

Message:

Modified offline anlysis

File:

• SPESO/ana2015/offline_analysis/offline_normalize_spectrum.m (modified) (4 diffs)

SPESO/ana2015/offline_analysis/offline_normalize_spectrum.m

@@ -169 +169 @@
 nu_comb=cell(length(Q)/2,1);
 
-det1_idx=7;
-det2_idx=8;
+det1_idx=15;
+det2_idx=16;
+
+bins_elem_num{det1_idx}
+bins_elem_num{det2_idx}
 
 Det2_val_raw=[values{det2_idx}]/3;
 Det1_val_raw=[values{det1_idx}];
-bins_elem_num{det1_idx}
-bins_elem_num{det2_idx}
 Det2_val=[corrected_values{det2_idx}]/3;
 Det1_val=[corrected_values{det1_idx}];
@@ -191 +192 @@
 Det1_ant_eff=[ant_eff{det1_idx}];
 
+% 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
+
+Det1_val_raw=Det1_val_raw.*det1_coeff;
+Det2_val_raw=Det2_val_raw.*det2_coeff;
+
 load('SC_RRR_10000fs.mat');
 
@@ -207 +217 @@
 end
 
-        coeff=1/20*10^12/310^2*data_corrFactor(10, sp_theta, 9, r0{7})...
-                /1000*1000*(1-exp(-50*10^6*20*10^-12)).*exp(d_eff2/10)*30;
+% smoothing experimental curve
+%{
+averaging_n=3;
+Det1_sel=fix(length(Det1_val_raw)/averaging_n)*averaging_n;
+Det2_sel=fix(length(Det2_val_raw)/averaging_n)*averaging_n;
+Det1_val_raw=mean(reshape(Det1_val_raw(end-Det1_sel+1:end),averaging_n,Det1_sel/averaging_n));
+Det2_val_raw=mean(reshape(Det2_val_raw(end-Det2_sel+1:end),averaging_n,Det2_sel/averaging_n));
+Det1_th=mean(reshape(Det1_th(end-Det1_sel+1:end),averaging_n,Det1_sel/averaging_n));
+Det2_th=mean(reshape(Det2_th(end-Det2_sel+1:end),averaging_n,Det2_sel/averaging_n));
+%}
+
+% removing edge points
+
+%Det1_val_raw=mean(reshape(Det1_val_raw,3,length));
+%Det2_val_raw=Det2_val_raw(1:end);
+%Det1_th=Det1_th(9:end);
+%Det2_th=Det2_th(1:end);
+%Det1_val_raw=Det1_val_raw(9:end);
+%Det2_val_raw=Det2_val_raw(1: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
                 
         fih=figure(1);
-        plot(Det1_th, Det1_val_raw-0.002, '--', 'LineWidth', 2, 'Color', [0 0.5 0]);
+        plot(Det1_th-5, Det1_val_raw, '--', 'LineWidth', 2, 'Color', [0 0.5 0]);
         xlim([40 140]);
         hold on
-        plot(Det2_th, Det2_val_raw*3-0.002, '--', 'LineWidth', 2, 'Color', [0.19 0.55 0.91]);
+        plot(Det2_th, 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]);
@@ -262 +296 @@
                          'Location', 'northwest');
         xlabel('theta, degree');
-        ylabel('magnitude, V')
+        ylabel('magnitude, W/mm^2')
         title('SPM Data. R=310mm')
    set(findall(gcf,'type','text'),'FontSize',20)%'fontWeight','bold')