Changeset 407 in ETALON for SPESO/ana2015

Timestamp:

Mar 24, 2016, 6:43:15 PM (8 years ago)

Author:

malovyts

Message:

Changed script to detect contrast and create plot on figure

File:

• SPESO/ana2015/analyse_dataContrast.m (modified) (1 diff)

SPESO/ana2015/analyse_dataContrast.m

@@ -1 @@
-function signal_is_Good=analyse_catchBadSignals(var_signal);
-        % determine if a signal is over the background
+function analyse_dataContrast();
 
-        global dvar
+        % analyse_contrast
+
+        global dvar;
+        global allData;
+        global alertVar;
+
         define_data_var();
+        define_alerts_var();
 
-        [valsused vals_signal] = data_get_vals(var_signal);
-        vals_signal = vals_signal(1,:);
-        [valsused vals_time] = data_get_vals(dvar.DATA_TIME_MINUTES);
-        vals_time=vals_time(1,:);
-        % ranges (in minutes) determined by hand form histograms        
-        signalFallMinTime = 20;
-        signalFallMaxTime = 27;
-        signalRiseMinTime =  0;
-        signalRiseMaxTime = 10;
+        data_types=[dvar.DATA1_VAL_LPM dvar.DATA2_VAL_LPM dvar.DATA1_VAL_SPM dvar.DATA2_VAL_SPM];
+        data_names={'LPM det.1','LPM det.2','SPM det.1','SPM det.2'};
 
-        % selecting signal in time intervals    
-        signalTime = find(vals_time             < signalFallMinTime);%in minutes
-        signalTime = find(vals_time(signalTime) > signalRiseMaxTime);%in minutes
-        bgTime     = find(vals_time             > signalFallMaxTime);%in minutes
+        signalMean =zeros(1,length(data_types));
+        bgMean     =zeros(1,length(data_types));
+        signalRMS  =zeros(1,length(data_types));
+        bgRMS      =zeros(1,length(data_types));
+        
+        for iloop=1:length(data_types)
+                
+                % get mean signal and bg for all day
 
-        signalValues=vals_signal(signalTime);
-        bgValues=vals_signal(bgTime);
-        
-        signalMean=mean(signalValues);
-        signalDeviationRMS=rms(signalValues-signalMean);
-        bgMean=mean(bgValues);
-        bgDeviationRMS=rms(bgValues-bgMean);
-        
-        % meanSignal-err < meanBG+err
-        if (abs(signalMean)-signalDeviationRMS<abs(bgMean)+bgDeviationRMS)
-                signal_is_Good=false;
-        else
-                signal_is_Good=true;
-        end;
+                %fih=figure('visible', 'off');
+                        bin=data_histc(dvar.DATA_TIME_MINUTES, data_types(iloop),...
+                                                                {{{dvar.DATA_TIME_30, '>', 6},{dvar.DATA_TIME_30, '<', 20}},...
+                                                                 {{dvar.DATA_TIME_30, '>', 25}}}, [0 24*60]);
+                %close(fih);
 
-        figure(1);
+                signalMean(iloop)=abs(bin.bins1(1,1));
+                bgMean(iloop)=abs(bin.bins2(1,1));
+                signalRMS(iloop)=abs(bin.bins1(3,1));
+                bgRMS(iloop)=abs(bin.bins2(3,1));
+        end; % iloop
+
+        % check if signal is within errors of bg 
+        if (length(find(signalMean-(bgMean+bgRMS)<0)))
+                %write alert into a file
+                alert_problem(alertVar.SIGNAL_BAD, '', '');
+        end % if bad signal
+
+        % figure signal and bg with errorbars for each detector
+                clf
                 hold on
-                errorbar(signalMean, signalDeviationRMS)
-                errorbar(bgMean, bgDeviationRMS, 'red');
+                        errorbar(signalMean, signalRMS, 'ro');
+                        errorbar(bgMean, bgRMS, 'bo' );
                 hold off
-        %print(figure(1), '../plots/signal_contrast.png', '-dpng');
-        %close(figure(1));
+                        legend('Signal', 'Background', 'Location', 'best')
+                set(gca, 'XTick', [1 2 3 4]);
+                set(gca, 'XTickLabel', data_names)
+                ylabel('Absolute value');
+                title({'Signal selection, background selection';...
+                                 'whole day';});
 
-        figure(2);
-                hold on
-                plot(vals_time, vals_signal)    
-                scatter(vals_time(signalTime),signalValues, 'r', '+');
-                scatter(vals_time(bgTime),bgValues, 'k', '+');
-                hold off
-        %print(figure(2), '../plots/signal_contrast2.png', '-dpng');
-        %close(figure(2));
-        
-end
+end% function