function []=figure2b()
close all
figure()
hold on
%% Data parameters
ceiling_T_VFT=0.0907;
ceiling_T_A=0.0738;
T_g=mean([ceiling_T_VFT,ceiling_T_A]);

%% Loading data
x_grain=0.001*(T_g);
x_grain_plot=0.03*(T_g);
Ts_fit=0:x_grain:T_g;
Ts_fit_plot=0:x_grain_plot:T_g;

% 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));
where=logical(Ts1<T_g);
Ts=Ts1(where);
Sconfs=Sconfs1(where);
p=polyfit(Ts,Sconfs,2);
c=p(1); b=p(2); a=p(3);
Sconf_fits=a+b*Ts_fit+c*Ts_fit.^2;
[~,i_where]=min(abs(Ts_fit-T_g));
Sconf_collapse=Sconf_fits(i_where);
Sconfs1=Sconfs1/Sconf_collapse;

% 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));
where=logical(Ts2<T_g);
Ts=Ts2(where);
Sconfs=Sconfs2(where);
p=polyfit(Ts,Sconfs,2);
c=p(1); b=p(2); a=p(3);
Sconf_fits=a+b*Ts_fit+c*Ts_fit.^2;
[~,i_where]=min(abs(Ts_fit-T_g));
Sconf_collapse=Sconf_fits(i_where);
Sconfs2=Sconfs2/Sconf_collapse;

% 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));
where=logical(Ts3<T_g);
Ts=Ts3(where);
Sconfs=Sconfs3(where);
p=polyfit(Ts,Sconfs,2);
c=p(1); b=p(2); a=p(3);
Sconf_fits=a+b*Ts_fit+c*Ts_fit.^2;
[~,i_where]=min(abs(Ts_fit-T_g));
Sconf_collapse=Sconf_fits(i_where);
Sconfs3=Sconfs3/Sconf_collapse;
ErrorU3=ErrorU3/Sconf_collapse;
ErrorL3=ErrorL3/Sconf_collapse;

% 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));
Sconfs4=1./xis;
where=logical(Ts4<T_g);
Ts=Ts4(where);
Sconfs=Sconfs4(where);
p=polyfit(Ts,Sconfs,2);
c=p(1); b=p(2); a=p(3);
Sconf_fits=a+b*Ts_fit+c*Ts_fit.^2;
[~,i_where]=min(abs(Ts_fit-T_g));
Sconf_collapse=Sconf_fits(i_where);
Sconfs4=Sconfs4/Sconf_collapse;

Sconf_fits_plot=a+b*Ts_fit_plot+c*Ts_fit_plot.^2;
Sconf_fits_plot=Sconf_fits_plot/Sconf_collapse;

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

end