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