- Timestamp:
- Mar 3, 2016, 10:03:13 AM (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
SPESO/ana2015/analyse_fast_dataContrast.m
r332 r351 10 10 define_data_var(); 11 11 12 % load data that will be analysed13 [valsused, vals_signal] = data_get_vals(var_signal);14 vals_signal_full = vals_signal(1,:);15 [valsused, vals_time] = data_get_vals(dvar.DATA_TIME_30);16 vals_time_full = vals_time(1,:);17 18 12 % ranges (in minutes) determined by hand form histograms 19 13 signalFallMinTime = 20; … … 40 34 data_select(dvar.DATA_TIME_30, '<', signalFallMinTime); 41 35 data_select(dvar.DATA_TIME_30, '>', signalRiseMaxTime); 42 allData.selection=[ 'Between ' num2str(signalRiseMaxTime) 'th and ' num2str(signalFallMinTime) 'th minute' ]; 36 allData.selection=[ 'Between ' num2str(signalRiseMaxTime)... 37 'th and ' num2str(signalFallMinTime) 'th minute' ]; 43 38 44 39 [valsused, vals_signal] = data_get_vals(var_signal); … … 56 51 data_select(dvar.DATA_TIME_30, '>', signalFallMaxTime); 57 52 data_select(dvar.DATA_TIME_30, '<', 31); 58 allData.selection=[ 'Between ' num2str(signalFallMaxTime) 'th and ' num2str(31) 'th minute' ]; 53 allData.selection=[ 'Between ' num2str(signalFallMaxTime)... 54 'th and ' num2str(31) 'th minute' ]; 59 55 60 56 [valsused, vals_signal] = data_get_vals(var_signal); … … 83 79 84 80 % Making plots 85 figure(1); 86 hold on 87 errorbar(signalMean, signalDeviationRMS, 'ro'); 88 errorbar(bgMean, bgDeviationRMS, 'bo' ); 89 legend('Signal', 'Background', 'Location', 'best') 90 hold off 91 ylabel(get_var_name(var_signal)); 92 title({'Signal selection, background selection'; get_var_name(var_signal); 'whole day'}); 81 % hold on 82 % errorbar(signalMean, signalDeviationRMS, 'ro'); 83 % errorbar(bgMean, bgDeviationRMS, 'bo' ); 84 % legend('Signal', 'Background', 'Location', 'best') 85 % hold off 86 % ylabel(get_var_name(var_signal)); 87 % title({'Signal selection, background selection'; get_var_name(var_signal); 'whole day'}); 93 88 %print(figure(1), '../plots/signal_contrast.png', '-dpng'); 94 89 %close(figure(1)); 95 90 96 figure(2);97 hold on98 scatter(vals_time_full, vals_signal_full, 'k')99 scatter(signalTime,signalValues, 'filled','r');100 scatter(bgTime,bgValues, 'filled','b');101 xlabel(get_var_name(dvar.DATA_TIME_30));102 ylabel(get_var_name(var_signal));103 hold off104 title({'Signal selection, background selection'; get_var_name(var_signal); 'whole day'});105 %print(figure(2), '../plots/signal_contrast2.png', '-dpng');106 %close(figure(2));107 91 108 92 allData.selected=selected_data;
Note: See TracChangeset
for help on using the changeset viewer.