%% Generate an undecorated version of Figure S8

function []=figureS8()
sigmabar=0.867966783525002;
xi_PTS=2.0322351/sigmabar; % point-to-set length at volume fraction phi=0.608

%% Read in the data
% With parallel tempering
data_name=sprintf('./PT.dat');
fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
NR=round(size(numbers,1)/2); % number of radial data points
R_PT=numbers(2*(1:NR)-1); % cavity radii
tau_PT=numbers(2*(1:NR)-0); % convergence time with parallel tempering

% Without parallel tempering
data_name=sprintf('./NPT.dat');
fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
NR=round(size(numbers,1)/2); % number of radial data points
R_NPT=numbers(2*(1:NR)-1); % cavity radii
tau_NPT=numbers(2*(1:NR)-0); % convergence time without parallel tempering [UNDERESTIMATE for R=R_NPT(1)=2.1]

%% Plot them
figure();
hold on
plot(R_PT,tau_PT);
scatter(R_PT,tau_PT,'s','filled');
plot(R_NPT(2:end),tau_NPT(2:end));
scatter(R_NPT(2:end),tau_NPT(2:end),'o','filled');
scatter(R_NPT(1),tau_NPT(1),'o'); % UNDERESTIMATE for convergence time without parallel tempering R=R_NPT(1)=2.1
plot([xi_PTS,xi_PTS],[10^6,10^10],'Linestyle','--'); % the point-to-set length

xlabel('$R$','Interpreter','latex');
ylabel('Convergence time')
xlim([2,5])
ylim([10^6,10^10])
set(gca,'yscale','log')
axis square
set(gcf, 'PaperPositionMode', 'auto');
print -depsc2 figureS8.eps



end