source: MML/trunk/machine/SOLEIL/StorageRing/diagnostics/fa_zoomer.m @ 17

% fa_zoomer([server])
%
% Investigate and zoom into archive data for the last day.
function fa_zoomer(server)
    % h is used to store all persistent state that is not part of the generated
    % data.
    h = {};

    if nargin == 0
        h.server = 'ac7b727-pc1';
    else
        h.server = server;
    end

    % Create figure with the standard toolbar but no menubar
    fig = figure('MenuBar', 'none', 'Toolbar', 'figure', ...
        'Position', [0 0 900 600]);
    global data;
    data = {};

    % Create the controls.
    global h_pos;
    h_pos = 10;
    h.bpm_list = control('edit', '4', 60, 'List of BPM FA ids');
    h.bpm_name_list = control('popup', ...
        dev2tangodev('BPMx',family2dev('BPMx')), 130, ...
        'Choose decimated data type to display', 'Value', 1, ...
        'Callback', protect(@bpm_name_callback));
    control('pushbutton', 'Back', 40, 'Return to previous zoom', ...
        'Callback', protect(@back_callback));
    control('pushbutton', 'Full', 40, 'View entire archive history', ...
        'Callback', protect(@full_archive_callback));
    control('pushbutton', '24h', 40, 'View last 24 hours', ...
        'Callback', protect(@last_day_callback));
    control('pushbutton', 'Zoom', 60, 'Update zoomed area from archive', ...
        'Callback', protect(@zoom_in_callback));
    control('pushbutton', 'Spectrogram', 100, 'Show as spectrogram', ...
        'Callback', protect(@spectrogram_callback));
    h.message = control('text', '', 150, ...
        'Error message or [bpm count] samples/decimation');
    h.maxpts = control('edit', num2str(1e6), 80, ...
        'Maximum number of sample points');
    h.ylim = control('checkbox', 'Zoomed', 70, ...
        'Limit vertical scale to +-100um', 'Value', 0, ...
        'Callback', protect(@reload_plot));
    h.data_type = control('popup', {'min/max', 'std', 'mean'}, 90, ...
        'Choose decimated data type to display', 'Value', 1, ...
        'Callback', protect(@reload_plot));
    control('pushbutton', 'Save', 40, 'Save Position data currently displayed', ...
        'Callback', protect(@save_data_callback));
    clear global h_pos;
    h.history = cell(0, 2);

    % Hang onto the controls we need to reference later
    guidata(fig, h);

    last_day_callback(fig, 0);
end


% Places control with specified style, value, width  and tooltip.
function result = control(style, value, width, tooltip, varargin)
    global h_pos;
    position = [h_pos 10 width 20];
    h_pos = h_pos + width + 5;
    result = uicontrol( ...
        'Style', style, 'String', value, 'Position', position, ...
        'TooltipString', tooltip, varargin{:});
end


function prot = protect(func)
    function protected(fig, event)
        try
            func(fig, event)
        catch
            message('Error: see console')
            rethrow(lasterror)
        end
    end
    prot = @protected;
end


function back_callback(fig, event)
    h = guidata(fig);
    if size(h.history, 1) > 1
        h.history = h.history(1:end-1, :);
        range = h.history{end, 1};
        type  = h.history{end, 2};
        guidata(fig, h)
        load_data(fig, range, type, false);
    end
end


function full_archive_callback(fig, event)
    load_data(fig, [now-2000 now], 'D', true);    % Go back as far as possible!
end


% Loads data for the last 24 hours
function last_day_callback(fig, event)
    load_data(fig, [now-1 now], 'D', true);
end


% Loads data enclosed by the current zoom selection, returned by xlim.
function zoom_in_callback(fig, event)
    h = guidata(fig);
    global data;

    maxdata = str2num(get(h.maxpts,   'String'));
    pvs     = str2num(get(h.bpm_list, 'String'));
    points = diff(xlim) * 24 * 3600 * 10079 * length(pvs);

    type = 'F';
    if points > maxdata
        type = 'd';
        points = points / 64;
        if points > maxdata
            type = 'D';
        end
    end

    load_data(fig, xlim + data.day, type, true);
end


function reload_plot(fig, event)
    h = guidata(fig);
    global data;
    plotfa(h, data);
end


% Loads the requested range of data.
function load_data(fig, range, type, save)
    h = guidata(fig);
    if save
        h.history(end+1, :) = {range, type};
        guidata(fig, h)
    end
    global data;

    pvs = str2num(get(h.bpm_list, 'String'));

    busy;
    data = fa_load(range, pvs, type, h.server);
    plotfa(h, data);
    describe;
end


function spectrogram_callback(fig, event)
    global data;
    %len=1024;
    len=20480;

    if length(size(data.data)) == 3
        busy;
        for n = 1:2
        subplot(2, 1, n)
            e = squeeze(data.data(n, 1, :));
            cols = floor(length(e)/len);
            sg = log10(abs(fft(reshape(e(1:(len*cols)), len, cols))));
            imagesc(data.t, [0 data.f_s/10], sg(1:(round(len/10)), :));

            set(gca, 'Ydir', 'normal');
            %caxis([2 6])
            caxis('auto');
            label_axis(n)
        end
        describe;
    else
        message('Decimated data');
    end
end

function save_data_callback(fig, event)
    h = guidata(fig);
    global data;
    if length(size(data.data)) == 4
            % decimated data can not be saved 
    else
        BPM_list=family2tangodev('BPMx',data.ids')
        bufferX=squeeze(data.data(1, :, :));
        bufferZ=squeeze(data.data(2, :, :));
        if size(BPM_list,1)==1
            bufferX=bufferX';
            bufferZ=bufferZ';
        end
        
        date=datestr(data.timestamp,30)
        date=[date(1:4),'_',date(5:6),'_',date(7:8),'_',date(10:11),'h_',date(12:13),'mn']
        filename=['/home/operateur/GrpDiagnostics/matlab/data/Enregistrement_',date,'.mat'];

        uisave({'BPM_list','bufferX','bufferZ'},filename)
    end
    
    
end

function plotfa(h, d)
    data_type = get(h.data_type, 'Value');
    if get(h.ylim, 'Value')
        if length(size(d.data)) == 4
            switch data_type
                case 1;     set_ylim = [-100 100];
                case 2;     set_ylim = [0 10];
                case 3;     set_ylim = [-10 10];
            end
        else
            set_ylim = [-100 100];
        end
    else
        set_ylim = [];
    end

    for n = 1:2
        subplot(2, 1, n)
        if length(size(d.data)) == 4
            switch data_type
                case 1
                    plot(d.t, 1e-3 * squeeze(d.data(n, 2, :, :))); hold on
                    plot(d.t, 1e-3 * squeeze(d.data(n, 3, :, :))); hold off
                case 2
                    plot(d.t, 1e-3 * squeeze(d.data(n, 4, :, :)));
                case 3
                    plot(d.t, 1e-3 * squeeze(d.data(n, 1, :, :)));
            end
        else
            plot(d.t, 1e-3 * squeeze(d.data(n, :, :)))
        end

        xlim([d.t(1) d.t(end)]);
        if length(set_ylim) > 0; ylim(set_ylim); end
        label_axis(n)
    end
end


function label_axis(n)
    axes = {'X'; 'Y'};
    global data;
    ylabel(gca, 'µm');
    if diff(data.t([1 end])) <= 2/(24*3600)
        title([datestr(data.timestamp) ' ' axes{n}])
        set(gca, 'XTickLabel', num2str( ...
            get(gca,'XTick').'*24*3600-60*floor(data.t(1)*24*60),'%.4f'))
    else
        title([datestr(data.day) ' ' axes{n}])
        datetick('keeplimits')
    end
end

function message(msg)
    h = guidata(gcf);
    set(h.message, 'String', msg);
end

function busy
    message('Busy');
    drawnow;
end

% Prints description of currently plotted data
function describe
    global data;
    message(sprintf('[%d] %d/%d', ...
        length(data.ids), length(data.data), data.decimation))
end

% Gives id from selected BPM Name
function bpm_name_callback(fig, event)
    h = guidata(gcf);
    index = get(h.bpm_name_list, 'Value');
    set(h.bpm_list, 'String', num2str(index));
end