% -------------------------------------------------------------------------------------------------
% SpettriEvB.m: 
%
% Copyright (C) 2004 Federico Cantini, Clinical Physiology Institute, CNR, Pisa, Italy
%
% This program is free software; you can redistribute it and/or modify it under the terms
% of the GNU General Public License as published by the Free Software Foundation; either
% version 2 of the License, or (at your option) any later version.
%
% This program is distributed "as is" and "as available" in the hope that it will be useful,
% but WITHOUT ANY WARRANTY of any kind; without even the implied warranty of MERCHANTABILITY
% or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License along with this program;
% if not, write to the Free Software Foundation, Inc.,
% 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
%
% For any comment or bug report, please send e-mail to: federico.cantini@ifc.cnr.it
% -------------------------------------------------------------------------------------------------

function score = SpettriEvB(mainPath_in,casiDir,casoT,casoS)

echo off

close all
debug=0;
graf=0;
extension = '.dat';
ns=8;

mainPath_in= [mainPath_in, casiDir];
fpname = fullfile(mainPath_in, [casoT extension]);
fprintf('Lettura del file: %s\n', fpname);
fid  = fopen(fpname,'r');
[datapT, count] = fread(fid,[ns,inf],'double');
fclose(fid);

fpname = fullfile(mainPath_in, [casoS extension]);
fprintf('Lettura del file: %s\n', fpname);
fid  = fopen(fpname,'r');
[datapS, count] = fread(fid,[ns,inf],'double');
fclose(fid);


fs = 1024;

Fmin = 3;
t1 = 56;
t2 = 58;
t4 = 60;


Ch = 7;

[ft2T,fp2] = specgram(datapT(Ch,fs*t1 + 1:fs*t2),length(datapT(Ch,fs*t1 + 1:fs*t2)),1024,ones(1,length(datapT(Ch,fs*t1 + 1:fs*t2))),0);
[ft3T,fp3] = specgram(datapT(Ch,fs*t2 + 1:fs*t4),length(datapT(Ch,fs*t2 + 1:fs*t4)),1024,ones(1,length(datapT(Ch,fs*t2 + 1:fs*t4))),0);

[ft2S,fp2] = specgram(datapS(Ch,fs*t1 + 1:fs*t2),length(datapS(Ch,fs*t1 + 1:fs*t2)),1024,ones(1,length(datapS(Ch,fs*t1 + 1:fs*t2))),0);
[ft3S,fp3] = specgram(datapS(Ch,fs*t2 + 1:fs*t4),length(datapS(Ch,fs*t2 + 1:fs*t4)),1024,ones(1,length(datapS(Ch,fs*t2 + 1:fs*t4))),0);



bi2 = min(find(fp2>Fmin));
ei2 = max(find(fp2<=10));

bi3 = min(find(fp3>Fmin));
ei3 = max(find(fp3<=10));

[maxT56_58, imaxT56_58] = max(abs(ft2T((bi2:ei2))).^2);
[maxT58_60, imaxT58_60] = max(abs(ft3T((bi3:ei3))).^2);
FmaxT56_58 = fp2(bi2+imaxT56_58);
FmaxT58_60 = fp3(bi3+imaxT58_60);

[maxS56_58, imaxS56_58] = max(abs(ft2S((bi2:ei2))).^2);
[maxS58_60, imaxS58_60] = max(abs(ft3S((bi3:ei3))).^2);
FmaxS56_58 = fp2(bi2+imaxS56_58);
FmaxS58_60 = fp3(bi3+imaxS58_60);


if (graf)

    hold off
    subplot(2,1,1)
    plot(fp2(bi2:ei2),28*abs(ft2S((bi2:ei2))).^2,'--b')
    hold on
    h=plot(fp3(bi3:ei3),28*abs(ft3S((bi3:ei3))).^2,'--r');
    title(casoS)
    hold off
    subplot(2,1,2)
    plot(fp2(bi2:ei2),28*abs(ft2T((bi2:ei2))).^2,'b')
    hold on
    h=plot(fp3(bi3:ei3),28*abs(ft3T((bi3:ei3))).^2,'r');
    title(casoT)
    xlabel('56-60 s')   
    figname = ['fig\' casiDir '\' casoT '_' casoS '_56-60'];
    %saveas(h,figname,'jpeg');
    hold off
end


score = 0;
if (FmaxT58_60 < FmaxT56_58) score = score+1; end
if (FmaxT58_60 > FmaxT56_58) score = score-1; end
if (FmaxS58_60 < FmaxS56_58) score = score-1; end
if (FmaxS58_60 > FmaxS56_58) score = score+1; end

if (score == 0)
    if (FmaxT58_60 < FmaxS58_60) score = score +1; end
    if (FmaxT58_60 > FmaxS58_60) score = score -1; end
end
score