function [cohfreq] = eeg_coherence(filepath,patient,window,noverlap,nfft,fs)

% Die Funktion cohfreq berechnet die Kohärenzen der jeweiligen
% elektrodenpaarungen und schreibt die Ergebnisse der Berrechnung in ein
% dafür vorgesehenes Array!Die entsprechenden Frequenzen werden in einer
% weiteren structure mit dem Namen 'frequency' abgespeichert.
% filepath=('Y:\EEG\Innsbruck_alt\Ruhe-EEGs\KJ\mat-Files');

filepath=[filepath,patient,'\mat-Files\'];
[filepath_array,number_of_files]=coherence_counter(filepath);

for i=1:number_of_files
    
    act_filepath=char(filepath_array(i));
    Data =load_Data(act_filepath);

%Data = importdata([filepath,patient,'\mat-Files\','A','*.mat'],'\t'); %
%öffne alle mat-files, die mit 'A' beginnen

cohfreq = struct('coherence',[],'frequency',[]);
cohfreq.coherence = struct('T5_O1',[],'T6_O2',[],'T3_T5',[],'T4_T6',[],'F3_P3',[],'F4_C4',[],'C3_P3',[],'C4_P4',[]);
cohfreq.frequency = struct('F',[]);

F3=Data(4,:);
F4=Data(6,:);
C3=Data(9,:);
C4=Data(11,:);
P3=Data(14,:);
P4=Data(16,:);
T3=Data(8,:);
T4=Data(12,:);
T5=Data(13,:);
T6=Data(17,:);
O1=Data(18,:);
O2=Data(19,:);  

el=fieldnames(cohfreq.coherence);
for el_pairs=1:8;
tmp = char(el(el_pairs));
if tmp == 'T5_O1'
    [elname,F] = mscohere(T5,O1,hanning(nfft),noverlap,(nfft/2+1),fs);
elseif tmp == 'T6_O2'
    [elname,F] = mscohere(T6,O2,hanning(nfft),noverlap,(nfft/2+1),fs); 
elseif tmp == 'T3_T5'
    [elname,F] = mscohere(T3,T5,hanning(nfft),noverlap,(nfft/2+1),fs);
elseif tmp =='T4_T6'
    [elname,F] = mscohere(T4,T6,hanning(nfft),noverlap,(nfft/2+1),fs);
elseif tmp == 'F3_P3'
    [elname,F] = mscohere(F3,P3,hanning(nfft),noverlap,(nfft/2+1),fs);
elseif tmp =='F4_C4'
    [elname,F] = mscohere(F4,C4,hanning(nfft),noverlap,(nfft/2+1),fs);
elseif tmp == 'C3_P3'
    [elname,F] = mscohere(C3,P3,hanning(nfft),noverlap,(nfft/2+1),fs);
elseif tmp == 'C4_P4'
    [elname,F] = mscohere(C4,P4,hanning(nfft),noverlap,(nfft/2+1),fs);
    
end

cohfreq.coherence = setfield(cohfreq.coherence,tmp,elname);
cohfreq.frequency =  setfield(cohfreq.frequency,'F',F);

coh=getfield(cohfreq.coherence,tmp)
for j= 1:number_of_files
for k=1:noverlap
    all_Data(k,j)= coh
    coh=coh+noverlap
end
end

% [a b c] = fileparts(act_filepath);
%save([filepath,'coherence_',patient,'_',b,'.mat'], 'cohfreq');

end
end