Changeset 408 in ETALON
- Timestamp:
- Mar 24, 2016, 6:44:09 PM (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
SPESO/ana2015/analyse_dayInfo.m
r401 r408 1 1 function analyse_dayInfo(date_str, figure_idx) 2 % analyse yyyymmdd.dayInfo file 2 3 global basedir; 3 4 global imageSizes; … … 21 22 22 23 % turn off figure windows 23 set(0, 'DefaultFigureVisible', 'off');24 %set(0, 'DefaultFigureVisible', 'off'); 24 25 25 26 data_types=[dvar.DATA1_VAL_LPM dvar.DATA2_VAL_LPM dvar.DATA1_VAL_SPM dvar.DATA2_VAL_SPM]; 26 27 data_names={'LPM det.1','LPM det.2','SPM det.1','SPM det.2'}; 27 28 28 load_analysis_data([date_str '.dayInfo']);29 if length(allData.selected)==030 display(['No events to analyse on ' date_str]);31 return;32 end;33 29 34 for iloop=1:length(data_types) 35 36 % get mean signal and bg for all day 30 % if there is no file to load 31 try 32 load_analysis_data([date_str '.dayInfo']); 33 catch 34 display(['No events to analyse on ' date_str]); 35 return; 36 end 37 37 38 fih=figure('visible', 'off');39 bin=data_histc(dvar.DATA_TIME_MINUTES, data_types(iloop),...40 {{{dvar.DATA_TIME_30, '>', 6},{dvar.DATA_TIME_30, '<', 20}},...41 {{dvar.DATA_TIME_30, '>', 25}}}, [0 24*60]);42 close(fih);43 44 signalMean(iloop)=abs(bin.bins1(1,1));45 bgMean(iloop)=abs(bin.bins2(1,1));46 signalRMS(iloop)=abs(bin.bins1(3,1));47 bgRMS(iloop)=abs(bin.bins2(3,1));48 end; % iloop49 50 % check if signal is within errors of bg51 if (length(find(signalMean-(bgMean+bgRMS)<0)))52 %write allert into a file53 alert_problem(alertVar.SIGNAL_BAD, date_str,...54 [get_var_name(data_type) '; ' ...55 'Between ' num2str(hours(hour_loop)*hours_multi)...56 'th and ' num2str((hours(hour_loop)+1)*hours_multi) 'th hour']);57 end % if bad signal58 59 % figure signal and bg with errorbars for each detector60 38 fih = figure(figure_idx); % scatter 61 hold on 62 errorbar(signalMean, signalRMS, 'ro'); 63 errorbar(bgMean, bgRMS, 'bo' ); 64 hold off 65 legend('Signal', 'Background', 'Location', 'best') 66 set(gca, 'XTick', [1 2 3 4]); 67 set(gca, 'XTickLabel', data_names) 68 ylabel('Absolute value'); 69 title({'Signal selection, background selection';... 70 'whole day';... 71 ['Data for ' date_str(1:4) '/' date_str(5:6) '/' date_str(7:8) char(10)]}); 39 clf 40 analyse_dataContrast(); 72 41 print(fih, [ basedirfigures 'analysis_figure_' num2str(figure_idx) '.png'], '-dpng'); 73 42 close(fih); … … 112 81 %-------------------------------------------------- 113 82 114 %Save selection115 selected_data=allData.selected;116 83 117 84 %====== Colored plot of data within 30 minutes ======= 118 85 fih = figure(dataloop); 119 86 clf 120 plot_handle = data_scatter2d(dvar.DATA_TIME_30, data_types(dataloop)); 121 set(plot_handle, 'MarkerEdgeColor','k'); 122 123 124 % selecting signal in time intervals 125 % > maximum time of signal rise, 126 % < minimum time of signal drop 127 128 data_select(dvar.DATA_TIME_30, '<', 20); 129 data_select(dvar.DATA_TIME_30, '>', 10); 130 hold on 131 plot_handle = data_scatter2d(dvar.DATA_TIME_30, data_types(dataloop)); 132 set(plot_handle, 'MarkerEdgeColor','r'); 133 set(plot_handle, 'MarkerFaceColor','r'); 134 hold off 135 136 % virtually reset selection 137 allData.selected=selected_data; 138 139 %selecting background in time intervals 140 % > maximum time of signal drop 141 % < end of range 142 143 data_select(dvar.DATA_TIME_30, '>', 27); 144 hold on 145 plot_handle = data_scatter2d(dvar.DATA_TIME_30, data_types(dataloop)); 146 set(plot_handle, 'MarkerEdgeColor','b'); 147 set(plot_handle, 'MarkerFaceColor','b'); 148 hold off 87 data_scatter30min(data_types(dataloop)); 149 88 print(fih,'-dpng',[ basedirfigures 'analysis_figure_' num2str(figure_idx+16+dataloop) '.png' ]); 150 89 close(fih);
Note: See TracChangeset
for help on using the changeset viewer.