function[obj] = deconvolve(obj, settings)
%% Perform pupil deconvolution for the entire session

if nargin ==1
    selected_trials = 1:length(obj.trials)
else
    selected_trials = settings.selection;
end


h = waitbar(0, 'Performing pupil deconvolution...');

for t = 1:length(selected_trials)
    selected_trial = selected_trials(t)
    obj.trials(selected_trial) = obj.trials(selected_trial).deconvolve(settings);

    waitbar(t/length(selected_trials), h);
end