%% Generate an undecorated version of Figure S7A

function []=figureS7A()
volume_fraction_index=1; % corresponds to volume fraction phi=0.568
cavity_radius_indices=[2,3,4,107]; % corresponds to cavity radii R=[1.38,1.73,2.07,3.11]

for ii=1:4
    %% Read in the data
    cavity_radius_index=cavity_radius_indices(ii);
    data_name=sprintf('./PDF%d_%d.dat',volume_fraction_index,cavity_radius_index);
    fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
    Ng=round(size(numbers,1)/2); % number of grid points for the histogram
    q_c=numbers(2*(1:Ng)-1); % core overlap
    P=numbers(2*(1:Ng)-0); % probability distribution function
         
    %% Plot them
    plot(q_c,P)
    hold on
end
xlim([0,1])
ylim([0,6])
xlabel('$q_{\rm c}$','Interpreter','latex')
ylabel('$P$','Interpreter','latex')
axis square
set(gcf, 'PaperPositionMode', 'auto');
print -depsc2 FigureS7A.eps
end