%function [faces] = getFaces(le,re,mo,thres);
function [faces] = getFaces(le,re,mo,thres);

nLe = length(le(1,:));
nRe = length(re(1,:));
nMo = length(mo(1,:));
nFace=0;
for j=1:nLe
  for k=1:nRe
     for t=1:nMo
        nFace = nFace + 1;        
        [p, geo] = eyesMouthLikelihood(le(:,j),re(:,k),mo(:,t));
        % faces=[(1)prob, (2)scale, (3-8)lx,ly,rx,ry,mx,my, (9)theta]
        faces(1,nFace)= p;
        faces(2,nFace)= 0;
        faces(3:4,nFace) = le(3:4,j);
        faces(5:6,nFace) = re(3:4,k);
        faces(7:8,nFace) = mo(3:4,t);
        faces(9,nFace) = geo(3);
     end
  end
end

  noplot = 'n';
  if noplot == 'y'
    return;
  end
 
  plotFaces(faces, le, re, mo, thres, 'b');

  return;