%GETRECALL get frame-based chord symbol recall
%
%[recall, nondictionary] = getrecall(estimates, annotations,cardinality)
% 
% params.comparison = 'ov'; 
% params.cardinality = 3; 
% params.inversion = 0;
% params.dictionary = {'X:maj', 'X:min'};

function [recall, nondictionary] = getrecall3(estimates_labfile, annotations_labfile,params)


[at1, at2, as] = merge_repeats(annotations_labfile);

[et1, et2, es] = merge_repeats(estimates_labfile);


%indeces for E and A
i = 1;
j = 1;

ilength = length(et1);
jlength = length(at1);

% find length in Nt of annotated sequence 
%alength = length(annotations);

% find length in Nt of estimated sequence 
%elength = length(estimates);

% use lesser value for evaluation
%ilength = min(alength,elength);

Nc = 0; % number of successfully matched Nt

Nt = 0; % total number of evaluated Nt

Ne = 0;

seg_start= 0;
seg_end = 0;

Tc = 0;
Ti = 0;
Tt = 0;
    
while i <= ilength & j <= jlength

    segstart = max(at1(j), et1(i));
    segend = min(at2(j), et2(i));
    
    [match, success, errormessage] = comparechords2(char(es{i}),char(as{j}), params.comparison, params.cardinality, params.inversion);
           
    includeE = indictionary(char(es{i}), params.dictionary, params.comparison, params.cardinality);
    includeA = indictionary(char(as{j}), params.dictionary, params.comparison, params.cardinality);
 
    seglength = segend-segstart;
    
    Tc = Tc + (match * includeE * seglength);
    
    Ti = Ti + (includeA * seglength);
    
    Tt = Tt + seglength;
    
    if at2(j) > segend
        i = i+1;
    elseif et2(i) > segend
        j= j+1;
    else
        i = i+1;
        j = j+1;
    end
    
end

% recall is number of correct matches out of total evaluated Nt
recall = Tc/Ti;

% percentage of annotations that were non-dictionary chords so ignored
nondictionary = 1-(Ti/Tt);
