function EEG = FilteringOpenBCI_EEGdata(EEG)
    % High-pass filter the data at 1-Hz. Note that EEGLAB uses pass-band edge, therefore 1/2 = 0.5 Hz.
    EEG = pop_eegfiltnew(EEG, 1, 0, 1650, 0, [], 0);
    % low-pass filter the data at 60-Hz.
    EEG = pop_eegfiltnew(EEG, 'hicutoff',60);

    % filtering step 2 erase 60 hz noise - necessary before conducting ERP & ICA?
    signal      = struct('data', EEG.data, 'srate', EEG.srate);
    lineNoiseIn = struct('lineNoiseMethod', 'clean', ...
                         'lineNoiseChannels', 1:EEG.nbchan,...
                         'Fs', EEG.srate, ...
                         'lineFrequencies', [60 120 180 240],...
                         'p', 0.01, ...
                         'fScanBandWidth', 2, ...
                         'taperBandWidth', 2, ...
                         'taperWindowSize', 4, ...
                         'taperWindowStep', 1, ...
                         'tau', 100, ...
                         'pad', 2, ...
                         'fPassBand', [0 EEG.srate/2], ...
                         'maximumIterations', 20);
    [clnOutput, ~] = cleanLineNoise(signal, lineNoiseIn); % ~: lineNoiseOut
    EEG.data = clnOutput.data;
end