clear all;
close all;

nParticles=432;

Pd={'0.01','0.02','0.03'};
cM={'sr','xb','og','vm','*c'};
%cM={'r','b','g','m','c'};

lattice = 'HS1';

%topDir = '/media/georgios/HDD1/Data/runPD';
topDir = './jammingHS1/';
allFile = 'allGaps2.dat';


dim1=500;
dim2=400;
fig=figure('position',[0,0,dim1,dim2]);

for i =1:numel(Pd)
    
    Pd(i)
    dataDir = [ topDir 'N' num2str(nParticles) '_s' strjoin(Pd(i)) filesep ];
    allFilename = [ dataDir allFile ];

    S=[];
    Bs=[];
    Hs=[];
    Hsum=[];
    Xs=[];
    
    S = load(allFilename);
    %S = S - mean(S);
    %S = S(S>0);
    %S = S(S>1.0e-10);

    M = numel(S);
    G = sort(S);
    x = 1:1:M;
    x = x./M;

    loglog( G, x, strjoin(cM(i)), 'DisplayName', ['Pd=' Pd{i}],'LineWidth',2);
    hold on;

end

a3 = 1-0.41269;
b3 = 5.e-1;
PL3 = b3*(G).^a3;
PL3 = PL3(PL3>min(x));
PL3 = PL3(PL3<max(x));
G3 = G(PL3>min(x));
G3 = G3(PL3<max(x));


PL3 = b3*(G3).^a3;
loglog( G3, PL3, 'k-', 'DisplayName', [' slope =' num2str(a3)],'LineWidth',2);

xlim([1e-8,2e+1])
ylim([1e-8,2e+0])
%axis equal
xlabel('gap');
ylabel('CDF(gap)');
title([ lattice '  All Gaps' '   N = ' num2str(nParticles) ]); 
%legend(['all Spectra Sl =' num2str(Sl)], ['slope1 =' num2str(a1)],'Location', 'NorthWest');
legend('-DynamicLegend', 'Location', 'Best');
title([ lattice '  all gaps' '   N = ' num2str(nParticles) ]);