function [newP,newT,sth] = checkPT(Q,P,T,th, peakLabels, troughLabels);
% check for same-height pts around each peak and trough found in newPT
% look at undecimated signal - less distortion

% Q the undecimated respiration signal
% P  the vector of peak indexes, on the decimated-by-5 signal
% T  the vector of trough indexes, on the decimated-by-5 signal
% th the threshold used by newPT: (max(Qd) - min(Qd)) * .1
% peakLabels & troughLabels assigned by classifyPeaks, then modified
% manually if necessary:  1 = valid 2 = invalid 3 = questionable
m = 0;  % counter for pts near each peak or trough
newP = []; % matrix of new peak elements
newT = []; % matrix of new trough elements
sth = th * .3;  % the small threshold
ssth = th * .01; % the super-small threshold

[n] = max(size(peakLabels));  % # of peaks found by newPT
for i = 1:n
    if peakLabels(i) == 1  % only find pause on valid peaks
        thisP = P(i) * 5;       % set to undecimated size
        
        m = 1;  % set counter for checking to the right of the peak
        while abs(Q(thisP + m) - Q(thisP)) < th & (m < 1000)
            m = m + 5;  % loop until end point is th higher than the peak
        end;
        littles = 0;
        stop = 0;  % set marker for moving back to end of pause
        while stop == 0 
            thisSlope = Q(thisP + m + 1) - Q(thisP + m);  % slope of last point
            if thisSlope > (-1 * ssth) & thisSlope < ssth
                littles = littles + 1;  % zero or not too steep slope
            else
                littles = 0;  % other - medium or steep slope
            end;
            if m <= 1 
                stop = 1;  % back to trough point - stop
                s = 'PR m <= 1'
            elseif Q(thisP + m) >= Q(thisP)
                stop = 1;  % height is same as peak point - stop
                s = 'PR back to start level'
            elseif littles > 9
                stop = 1;  % slope not changing much anymore
                m = m + 5; % move back to middle of the 3 littles
                s = 'PR littles > 9'
            else
                m = m - 1; % back one point
            end;
        end;
        newP = [newP, (thisP + m)];  % add in new point
        
        m = -1;  % set counter for checking to the left
        while abs(Q(thisP + m) - Q(thisP)) < th & (m > -1000)
            m = m - 5;  % loop until end point is th higher than the peak
        end;
        littles = 0;
        stop = 0;  % set marker for moving back to end of pause
        while stop == 0 
            thisSlope = Q(thisP + m - 1) - Q(thisP + m);  % slope of last point
            if thisSlope > (-1 * ssth) & thisSlope < ssth
                littles = littles + 1;  % zero or not too steep
            else
                littles = 0;  % other - medium or steep slope
            end;
            if m == -1  % back to peak point - stop
                stop = 1;
                s = 'PL m <= 1'
            elseif Q(thisP) <= Q(thisP + m)
                stop = 1;  % height is same as peak point - stop
                s = 'PL back to start level'
            elseif littles > 9
                stop = 1;  % slope not changing much anymore
                m = m - 5; % move back to middle of "little" stretch
                s = 'PL littles > 9'
            else
                m = m + 1;  % back one point
            end;
        end;
        newP = [newP, (thisP + m)];
    end;
end;

[n] = max(size(troughLabels));  % # of troughs found by newPT
for i = 1:n
    if troughLabels(i) == 1  % only find pause on valid peaks
        thisT = T(i) * 5;       % set to undecimated size
        
        m = 1;  % set counter for checking to the right of the trough
        while abs(Q(thisT + m) - Q(thisT)) < th & (m < 1000)
            m = m + 5;  % loop until end point is th higher than the trough
        end;
        littles = 0;
        stop = 0;  % set marker for moving back to end of pause
        while stop == 0 
            thisSlope = Q(thisT + m + 1) - Q(thisT + m);  % slope of last point
            if thisSlope > (-1 * ssth) & thisSlope < ssth
                littles = littles + 1;  % zero or not too steep slope
            else
                littles = 0;  % other - medium or steep slope
            end;
            if m <= 1 
                stop = 1;  % back to trough point - stop
                s = 'TR m <= 1'
            elseif Q(thisT + m) <= Q(thisT)
                stop = 1;  % height is same as trough point - stop
                s = 'TR back to start level'
            elseif littles > 9
                stop = 1;  % slope not changing much anymore
                m = m + 5; % move back to middle of the 3 littles
                s = 'TR littles > 9'
            else
                m = m - 1; % back one point
            end;
        end;
        newT = [newT, (thisT + m)];  % add in new point
        
        m = -1;  % set counter for checking to the left
        while abs(Q(thisT + m) - Q(thisT)) < th & (m > -1000)
            m = m - 5;  % loop until end point is th higher than the trough
        end;
        littles = 0;
        stop = 0;  % set marker for moving back to end of pause
        while stop == 0 
            thisSlope = Q(thisT + m - 1) - Q(thisT + m);  % slope of last point
            if thisSlope > (-1 * ssth) & thisSlope < ssth
                littles = littles + 1;  % zero or not too steep
            else
                littles = 0;  % other - medium or steep slope
            end;
            if m == -1  % back to trough point - stop
                stop = 1;
                s = 'TL m <= 1'
            elseif Q(thisT + m) <= Q(thisT)
                stop = 1;  % height is same as trough point - stop
                s = 'TL back to start level'
            elseif littles > 9
                stop = 1;  % slope not changing much anymore
                m = m - 5; % move back to middle of "little" stretch
                s = 'TL littles > 9'
            else
                m = m + 1;  % back one point
            end;
        end;
        newT = [newT, (thisT + m)];
    end;
end;

plot(Q); hold on;
n = max(size(P));  % n has number of peaks found
for i = 1:(n - 1)
    if peakLabels(i) == 1 | peakLabels(i) == 3 
        plot(P(i)*5, Q(P(i)*5), 'xg');
    end;
end;
plot(newP, Q(newP), 'oy', 'MarkerSize',10);

n = max(size(T)); % number of troughs
for i = 1:(n - 1)
    if troughLabels(i) == 1 | troughLabels(i) == 3 
        plot(T(i)*5, Q(T(i)*5), '+g');
    end;
end;
plot(newT, Q(newT), 'ow', 'MarkerSize',10);
hold off;

plot(Qd); hold on;
n = max(size(P));  % n has number of peaks found
for i = 1:(n - 1)
    if peakLabels(i) == 1
        plot(P(i), Qd(P(i)), 'xg');
    elseif peakLabels(i) == 2
        plot(P(i), Qd(P(i)), 'xr');
    elseif peakLabels(i) == 3
        plot(P(i), Qd(P(i)), 'xb');
    end;
    text(P(i), Qd(P(i)) + th*.5, ['p', int2str(i)], 'FontSize', 8);
    %plot((P(i)+200), Qd(P(i)+200), 'oy');  % end of window
    %plot((P(i)-200), Qd(P(i)-200), '.y');  % start of window
end;
n = max(size(T)); % number of troughs
for i = 1:(n - 1)
    if troughLabels(i) == 1
        plot(T(i), Qd(T(i)), '+g');
    elseif troughLabels(i) == 2
        plot(T(i), Qd(T(i)), '+r');
    elseif troughLabels(i) == 3
        plot(T(i), Qd(T(i)), '+b');
    elseif troughLabels(i) == 4
        plot(T(i), Qd(T(i)), '+y');
    end;
    text(T(i), Qd(T(i)) - th*.5, ['t', int2str(i)], 'FontSize', 8);
end;
hold off;