function [x,y,h] = pickdata(varargin)
%PICKDATA   Picks data from the active curve.
%   [X,Y] = PICKDATA returns the data (X,Y) from the active curve. If
%   several curves are present and none is selected, the first one is
%   picked (see FITPARAM to change this default setting).
%
%   If some data are selected using the "Brushing" tool of Matlab
%   (works only for versions >= 7.6), only those "brushed" data are
%   returned.
%
%   [X,Y] = PICKDATA(H) returns the data (X,Y) of the curve specified by
%   the handle H (eg, GCF, GCO etc). GCF, the current figure, is taken by
%   default.
%
%   [X,Y,H] = PICKDATA(...) also returns the handle to the current curve.
%
%   See also SHOWFIT, EZFIT.

%   F. Moisy, moisy_at_fast.u-psud.fr
%   Revision: 1.31,  Date: 2017/01/23
%   This function is part of the EzyFit Toolbox

% History:
% 2005/05/31: v1.00, first version.
% 2006/02/23: v1.01, better check if a current figure exists.
% 2006/03/07: v1.02, also works with histogram plots.
% 2006/07/06: v1.03, pick the first or last curve, according to the
%                    variable fp.whichpickdata in fitparam.m
% 2006/10/18: v1.10, new parameter fp ; now takes the mean for x for hist
% 2008/03/14: v1.20, uses the 'BrushData' field of Matlab 7.6 (for
%                    EzyFit >=2.30)
% 2016/04/28: v1.30, new behaviour of 'BrushData' since R2014b
% 2017/01/23: v1.31, now works with 'histogram' (thanks J. Valenzuela!)
% 2017/01/24: v1.32, modification for brushed data in 2016b (


% gr_dummy not defined for old versions
if verLessThan('matlab','8.4')
    eval('gr_dummy = 0;');
else
    gr_dummy = groot;
end

hf=get(gr_dummy,'CurrentFigure');
if isempty(hf)
    error('EzyFit:pickdata:InvalidInput', 'No active figure.');
end

h=gco;
if nargin==1
    if ishandle(varargin{1})
        h=varargin{1};
    elseif isstruct(varargin{1})
        fp=varargin{1};
    end
elseif nargin==2
    h=varargin{1};
    fp=varargin{2};
end

if ~exist('fp','var')
    % loads the default fit parameters:
    try
        fp=fitparam;
    catch
        error('No fitparam file found.');
    end
end

co=get(h); % search for active curve

% if no active curve, take the 1st curve of the active fig.
if ((isempty(co)) || (~isfield(co,'XData')))
    ha=get(gcf,'CurrentAxes');
    if ~isempty(ha)
        ca=get(ha);
        if ~isempty(ca.Children)
            if strcmp(fp.whichpickdata,'first') % new v1.03
                h=ca.Children(end);  % (end) is the first curve
            else
                h=ca.Children(1);  % (1) is the last curve
            end
            co=get(h);
        else
            error('EzyFit:pickdata:InvalidInput', 'No curve in the active figure.');
        end
    else
        error('EzyFit:pickdata:InvalidInput', 'No axis in the active figure.');
    end
end


try
    x=co.XData;
    y=co.YData;
catch    % for 'histogram' plots (v1.31)
    x=co.BinEdges(2:size(co.BinEdges,2));
    y=co.Values;
end

% % BrushData behavior before R2014b (v1.20)
% if isfield(co,'BrushData')
%     hb = get(co,'BrushData');
%     if any(hb)
%         x=x(hb==1);
%         y=y(hb==1);
%     end
% end

% BrushData behavior changed in R2014b (v1.30), modified for R2016b (v1.32)
try
    hb = get(h,'BrushData');
    if any(hb)
        x=x(hb==1);
        y=y(hb==1);
    end
end


%new v1.02: if the data originates from an histogram plot,
%keeps only the 2nd (or 3rd) row:
if size(x,1)~=1
    x = (x(2,:)+x(3,:))/2;  % changed v1.10
    y = y(2,:);
end

if nargout==0
    x
    y
    clear x y;
end