function[mean] = mean(session, varargin)
%% options :  '''Trials', [list of trials]}, ...
%              '''Type'', {raw|interpolated|filtered}'
%              '''BaselineCorrected''0/1'
%              '''Marker'', marker id'
%              '''Settings''' , settings object
params = struct('Trials', [], ...
                'Type', 'none', ...
                'BaselineCorrected', 1, ...
                'Marker', 'stimulus onset',...
                'Settings', []);
            
params = parseparams(varargin);


if isempty(params.Trials)
    params.Trials = 1:length(sessions.trials);
end

if isstruct(params.Settings)
    threshold = params.Settings.QualityThreshold;    
    quality_values = [ session.trials(Trials).quality ];    
	trials_to_kill = find(quality_values < threshold);
    
    %exclude trials that fail the quality measure
    params.Trials = setdiff(params.Trials, trials_to_kill); 
end


for t = 1:length(params.Trials)    
    trial_id = params.Trials(t);    
    trial = session.trials(trial_id);
    [times{t} traces{t}] = trial.getTrialDataShiftedForMarker(marker_id);    
end


    
if nargin == 1