%% Generate an undecorated version of Figure 2d

function []=Figure2d()
polydispersity_index=0;  % polydispersity, in unit of percentage.

%% Read in the data
% Positional point-to-set correlation length
data_name=sprintf('./xi_posPTS_poly%d.dat',polydispersity_index); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Nd=round(size(numbers,1)/3);% number of data points
numbers=reshape(numbers,3,Nd);
phis_posPTS=numbers(1,1:Nd).'; % packing fractions
xiposPTS=numbers(2,1:Nd).'; % positional point-to-set correlation length
xiposPTS_error=numbers(3,1:Nd).'; % error bars

% Positional pair correlation length
data_name=sprintf('./xi_pospair_poly%d.dat',polydispersity_index); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Nd=round(size(numbers,1)/3);% number of data points
numbers=reshape(numbers,3,Nd);
phis_pospair=numbers(1,1:Nd).'; % packing fractions
xipospair=numbers(2,1:Nd).'; % positional pair correlation length
xipospair_error=numbers(3,1:Nd).'; % error bars

% Hexatic point-to-set correlation length
data_name=sprintf('./xi_hexPTS_poly%d.dat',polydispersity_index); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Nd=round(size(numbers,1)/3);% number of data points
numbers=reshape(numbers,3,Nd);
phis_hexPTS=numbers(1,1:Nd).'; % packing fractions
xihexPTS=numbers(2,1:Nd).'; % hexatic point-to-set correlation length
xihexPTS_error=numbers(3,1:Nd).'; % error bars

% Hexatic pair correlation length
data_name=sprintf('./xi_hexpair_poly%d.dat',polydispersity_index); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Nd=round(size(numbers,1)/3);% number of data points
numbers=reshape(numbers,3,Nd);
phis_hexpair=numbers(1,1:Nd).'; % packing fractions
xihexpair=numbers(2,1:Nd).'; % hexatic pair correlation length
xihexpair_error=numbers(3,1:Nd).'; % error bars

%% Plot them
figure();
hold on
errorbar(phis_posPTS,xiposPTS/xiposPTS(2),xiposPTS_error/xiposPTS(2),'r','Linestyle','none','Marker','s')
plot(phis_pospair,xipospair/xipospair(2),'r','Linestyle','--')
errorbar(phis_hexPTS,xihexPTS/xihexPTS(2),xihexPTS_error/xihexPTS(2),'b','Linestyle','none','Marker','d')
plot(phis_hexpair,xihexpair/xihexpair(2),'b','Linestyle','--')
xlabel('$\phi$','Interpreter','latex')
ylabel('$\xi/\xi_0$','Interpreter','latex');
axis square
set(gcf, 'PaperPositionMode', 'auto');
print -depsc2 Figure2d.eps


end