function plotHemoSpeedHist(summit, nonsummit, expDir, var, expName)
    figure;
    hold on;
    histSum = summit(find(summit));
    histSum = histSum(~isnan(histSum));
    a = histogram(log(histSum),100,'Normalization','Probability');

    nonHistSum = nonsummit(find(nonsummit));
    nonHistSum = nonHistSum(~isnan(nonHistSum));
    b = histogram(log(nonHistSum),100,'Normalization','Probability');

    legend({'Summit','Non-summit'},'Location','Best');
    xlabel('Log(speed)');
    ylabel('Probability (normalized)');
    title([expName,"Speed distribution"]);
    
    ylim([0 0.04]);
    
    set(gcf,'PaperUnits','inches','PaperPosition',[0 0 7.5 6]);
    outName = strcat(expDir,var,'-histogam');
    print(outName,'-dpng');  
    close;

end