source: ETALON/SPESO/ana2015/analyse_dataContrast.m @ 407

function analyse_dataContrast();

        % analyse_contrast

        global dvar;
        global allData;
        global alertVar;

        define_data_var();
        define_alerts_var();

        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'};

        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

                %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);

                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, signalRMS, 'ro');
                        errorbar(bgMean, bgRMS, 'bo' );
                hold off
                        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';});

end% function