function [EMG]=readEMGSig(pathToFile,subject,number,maxFile,windowEMG,windowShift,EMGChannels)
    EMGcell=cell(maxFile,1);
    
    for i=1:maxFile
        [sig, ~, params] = load_bcidat(strcat(pathToFile,...
            sprintf('%s/%s_B100%i/%s_B1S00%iR0%i.dat',subject,subject,number,subject,number,i)));
        EMGcell{i}=waveformLengthAll(sig(:,ismember(params.ChannelNames.Value,EMGChannels)),...
            params.SamplingRate.NumericValue,windowEMG,windowShift);
    end
    
    EMG=cell2mat(EMGcell);
end