drawSyntheticSSF.m

function drawSyntheticSSF(A,K,syntheticDataset,plotTitle)

[~,Y] = Cluster.SpectralCluster(A,K);

[~,~,idxTrue] = unique([syntheticDataset.phantomGroup]');
csf = CSF.csf_spatial(Y,idxTrue);

mm = mean(csf);
ss = std(csf);

sym = {'o','x','*'};

hold on
hx = [];

for ff=1:size(Y,1)
    if syntheticDataset(ff).phantomGroup == 1
        color = [1 0 0];
        idxSym = 1;
    elseif syntheticDataset(ff).phantomGroup == 2
        color = [0 1 0];
        idxSym = 2;
    else
        color = [0 0.4470 0.7410];
        idxSym = 3;

    end

    % h1 = plot3(Y(ff,1),Y(ff,2),Y(ff,3),'.','color',color,'UserData',ff,'MarkerSize',8); % for debugging, does not use different markers
    h1 = plot3(Y(ff,1),Y(ff,2),Y(ff,3),sym{idxSym},'color',color,'MarkerSize',6);
    hx(idxSym) = h1;

    % for debugging, lists names for each point
    % phantomID = strcat({' '},syntheticDataset(ff).phantomName);
    % text(Y(ff,1),Y(ff,2),Y(ff,3),phantomID,'Interpreter','none','color',color, 'FontSize',8)

end

axis square

set(gcf,'color','w')
xlabel('k1');ylabel('k2'),zlabel('k3')
lgd = legend(hx,{'v=1','v=3','v=5'});
set(lgd,'color','none');

title([plotTitle,' = [', num2str(round(mm,2,'significant')), ',' ...
    ,num2str(round(ss,2,'significant')),']']);