
%tumorSetSize=105;

figure(99);clf;mode=1;indOfInterest=[1:tumorSetSize];
for te=1:10
    subplot(5,2,te)
    TumorFail=projMatOnTumor(indOfInterest,mode,te)<0;
    stem(TumorFail);hold on;
    temp=flipud(squeeze(projMatOnNormal(:,mode,te)));
    NormalFail=(temp(indOfInterest)>0);
    stem(NormalFail*2,'r');
    BothFail=TumorFail.*NormalFail;
    stem(BothFail*3,'k');
    CombFail=((projMatOnTumor(indOfInterest,mode,te)-temp(indOfInterest))<0);
    stem(CombFail*4,'m');
    hold off;grid on
    legend(['Tumor ' num2str(sum(TumorFail))],...
        ['Normal ' num2str(sum(NormalFail))],...
        ['Both ' num2str(sum(BothFail))],...
        ['Combination ' num2str(sum(CombFail))])
end

figure(999);clf;mode=1;indOfInterest=tumorSetSize+[1:normalSetSize];
for te=1:10
    subplot(5,2,te)
    TumorFail=projMatOnTumor(indOfInterest,mode,te)>0;
    stem(TumorFail);hold on;
    temp=flipud(squeeze(projMatOnNormal(:,mode,te)));
    NormalFail=(temp(indOfInterest)<0);
    stem(NormalFail*2,'r');
    BothFail=TumorFail.*NormalFail;
    stem(BothFail*3,'k');
    CombFail=((projMatOnTumor(indOfInterest,mode,te)-temp(indOfInterest))>0);
    stem(CombFail*4,'m');
    hold off;grid on
    legend(['Tumor ' num2str(sum(TumorFail))],...
        ['Normal ' num2str(sum(NormalFail))],...
        ['Both ' num2str(sum(BothFail))],...
        ['Combination ' num2str(sum(CombFail))])
end

