%% Generate an undecorated version of Figure 4

function []=Figure4()
%% Read in data for Kob-Andersen liquid
% Relaxation time data
data_name=sprintf('./tau_alphas_KA.dat'); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Ts=numbers(3*(1:5)-2); % temperatures
tau_alphas_KA=numbers(3*(1:5)-1); % tau_alphas, i.e., structural relaxation time
tau_alphas_error=numbers(3*(1:5)); % error bars
% Structural PTS data
data_name=sprintf('./point_to_set_KA.dat'); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Ts=numbers(3*(1:5)-2); % temperatures
xi_PTSs_KA=numbers(3*(1:5)-1); % xi_PTSs, i.e., structural point-to-set length
xi_PTSs_error=numbers(3*(1:5)); % error bars

%% Read in data for monodisperse hard disks
% Relaxation time data
data_name=sprintf('./HexTime0.dat'); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
phis=numbers(3*(1:5)-2); % packing fractions
tau_alphas_0=numbers(3*(1:5)-1); % tau_alphas, i.e., hexatic relaxation time
tau_alphas_error=numbers(3*(1:5)); % error bars
% Hexatic PTS data
data_name=sprintf('./HexLength0.dat'); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
phis=numbers(3*(1:5)-2); % packing fractions
xi_PTSs_0=numbers(3*(1:5)-1); % xi_PTSs, i.e., hexatic point-to-set length
xi_PTSs_error=numbers(3*(1:5)); % error bars

%% Read in data for 11% polydisperse hard disks
% Relaxation time data
data_name=sprintf('./HexTime11.dat'); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
phis=numbers(3*(1:5)-2); % packing fractions
tau_alphas_11=numbers(3*(1:5)-1); % tau_alphas, i.e., hexatic relaxation time
tau_alphas_error=numbers(3*(1:5)); % error bars
% Hexatic PTS data
data_name=sprintf('./HexLength11.dat'); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
phis=numbers(3*(1:5)-2); % packing fractions
xi_PTSs_11=numbers(3*(1:5)-1); % xi_PTSs, i.e., hexatic point-to-set length
xi_PTSs_error=numbers(3*(1:5)); % error bars

%% Normalize them at high temperature/low density
tau_alphas_KA=tau_alphas_KA/tau_alphas_KA(2);
xi_PTSs_KA=xi_PTSs_KA/xi_PTSs_KA(2);
tau_alphas_0=tau_alphas_0/tau_alphas_0(2);
xi_PTSs_0=xi_PTSs_0/xi_PTSs_0(2);
tau_alphas_11=tau_alphas_11/tau_alphas_11(2);
xi_PTSs_11=xi_PTSs_11/xi_PTSs_11(2);

%% Plot them
figure();
hold on
plot(xi_PTSs_KA,tau_alphas_KA);
plot(xi_PTSs_0,tau_alphas_0);
plot(xi_PTSs_11,tau_alphas_11);
xlabel('$\xi/\xi_0$','Interpreter','latex')
ylabel('$\tau_{\alpha}(\xi)/\tau_{\alpha}(\xi_0)$','Interpreter','latex')
xlim([0,8])
ylim([0.1,1000])
set(gca,'yscale','log')
axis square
set(gcf, 'PaperPositionMode', 'auto');
print -depsc2 Figure4.eps

end