function []=figure2a()
close all
figure()
hold on
%% Loading data
ceiling_T_VFT=0.0907;
ceiling_T_A=0.0738;
T_g=mean([ceiling_T_VFT,ceiling_T_A]);

% PEL 1
data_name=sprintf('./sconf_PEL1.dat');
fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Ndata=round(size(numbers,1)/2);
Ts1=numbers(2*(1:Ndata)-1);
Sconfs1=numbers(2*(1:Ndata));

% PEL 2
data_name=sprintf('./sconf_PEL2.dat');
fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Ndata=round(size(numbers,1)/2);
Ts2=numbers(2*(1:Ndata)-1);
Sconfs2=numbers(2*(1:Ndata));

% Frankel-Ladd
data_name=sprintf('./sconf_FL.dat');
fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Ndata=round(size(numbers,1)/4);
Ts3=numbers(4*(1:Ndata)-3);
Sconfs3=numbers(4*(1:Ndata)-2);
% [value-L3, value+U3]
ErrorU3=numbers(4*(1:Ndata)-1);
ErrorL3=numbers(4*(1:Ndata));

% Point-to-set correlation length
data_name=sprintf('./cavityPTSthres.dat');
fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Ndata=round(size(numbers,1)/2);
Ts4=numbers(2*(1:Ndata)-1);
xis=numbers(2*(1:Ndata));

xi_g=xis(5)+(xis(4)-xis(5))*((T_g-Ts4(5))/(Ts4(4)-Ts4(5)));
Sconfs4=xi_g./xis;

%% Plotting
scatter(Ts1,Sconfs1,'Marker','o')
scatter(Ts2,Sconfs2,'Marker','d')
errorbar(Ts3,Sconfs3,ErrorL3,ErrorU3, 'LineStyle','none','Marker','^')
scatter(Ts4,Sconfs4,'Marker','s')
xlim([0,0.15])
ylim([0,2])
xlabel('$T$','Interpreter','latex');
ylabel('$s_{\mathrm{conf}}$','Interpreter','latex')
print(gcf, '-dpdf', 'figure2a.pdf');
end