%% Generate an undecorated version of Figure 3g

function []=Figure3g()

%% Read in the data
% Positional point-to-set correlation length
data_name=sprintf('./xi_posPTS.dat'); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Nd=round(size(numbers,1)/3);% number of data points
numbers=reshape(numbers,3,Nd);
T_posPTS=numbers(1,1:Nd).'; % temperature
xiposPTS=numbers(2,1:Nd).'; % positional point-to-set correlation length
xiposPTS_error=numbers(3,1:Nd).'; % error bars

% Positional pair correlation length
data_name=sprintf('./xi_pospair.dat'); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Nd=round(size(numbers,1)/3);% number of data points
numbers=reshape(numbers,3,Nd);
T_pospair=numbers(1,1:Nd).'; % temperature
xipospair=numbers(2,1:Nd).'; % positional pair correlation length
xipospair_error=numbers(3,1:Nd).'; % error bars

% l=6 point-to-set correlation length
data_name=sprintf('./xi_l=6PTS.dat'); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Nd=round(size(numbers,1)/3);% number of data points
numbers=reshape(numbers,3,Nd);
T_6PTS=numbers(1,1:Nd).'; % temperature
xi6PTS=numbers(2,1:Nd).'; % l=6 point-to-set correlation length
xi6PTS_error=numbers(3,1:Nd).'; % error bars

% l=6 pair correlation length
data_name=sprintf('./xi_l=6pair.dat'); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Nd=round(size(numbers,1)/3);% number of data points
numbers=reshape(numbers,3,Nd);
T_6pair=numbers(1,1:Nd).'; % temperature
xi6pair=numbers(2,1:Nd).'; % l=6 pair length
xi6pair_error=numbers(3,1:Nd).'; % error bars

% l=7 point-to-set correlation length
data_name=sprintf('./xi_l=7PTS.dat'); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Nd=round(size(numbers,1)/3);% number of data points
numbers=reshape(numbers,3,Nd);
T_7PTS=numbers(1,1:Nd).'; % temperature
xi7PTS=numbers(2,1:Nd).'; % l=7 point-to-set correlation length
xi7PTS_error=numbers(3,1:Nd).'; % error bars

% l=7 pair correlation length
data_name=sprintf('./xi_l=7pair.dat'); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Nd=round(size(numbers,1)/3);% number of data points
numbers=reshape(numbers,3,Nd);
T_7pair=numbers(1,1:Nd).'; % temperature
xi7pair=numbers(2,1:Nd).'; % l=7 pair length
xi7pair_error=numbers(3,1:Nd).'; % error bars

% l=8 point-to-set correlation length
data_name=sprintf('./xi_l=8PTS.dat'); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Nd=round(size(numbers,1)/3);% number of data points
numbers=reshape(numbers,3,Nd);
T_8PTS=numbers(1,1:Nd).'; % temperature
xi8PTS=numbers(2,1:Nd).'; % l=8 point-to-set correlation length
xi8PTS_error=numbers(3,1:Nd).'; % error bars

% l=8 pair correlation length
data_name=sprintf('./xi_l=8pair.dat'); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Nd=round(size(numbers,1)/3);% number of data points
numbers=reshape(numbers,3,Nd);
T_8pair=numbers(1,1:Nd).'; % temperature
xi8pair=numbers(2,1:Nd).'; % l=8 pair length
xi8pair_error=numbers(3,1:Nd).'; % error bars

% l=12 point-to-set correlation length
data_name=sprintf('./xi_l=12PTS.dat'); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Nd=round(size(numbers,1)/3);% number of data points
numbers=reshape(numbers,3,Nd);
T_12PTS=numbers(1,1:Nd).'; % temperature
xi12PTS=numbers(2,1:Nd).'; % l=12 point-to-set correlation length
xi12PTS_error=numbers(3,1:Nd).'; % error bars

% l=12 pair correlation length
data_name=sprintf('./xi_l=12pair.dat'); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Nd=round(size(numbers,1)/3);% number of data points
numbers=reshape(numbers,3,Nd);
T_12pair=numbers(1,1:Nd).'; % temperature
xi12pair=numbers(2,1:Nd).'; % l=12 pair length
xi12pair_error=numbers(3,1:Nd).'; % error bars

%% Plot them
figure();
hold on
errorbar(1./T_posPTS,xiposPTS/xiposPTS(2),xiposPTS_error/xiposPTS(2),'r','Linestyle','-')
plot(1./T_pospair,xipospair/xipospair(2),'r','Linestyle','--')
errorbar(1./T_6PTS,xi6PTS/xi6PTS(2),xi6PTS_error/xi6PTS(2),'b','Linestyle','-')
plot(1./T_6pair,xi6pair/xi6pair(2),'b','Linestyle','--')
errorbar(1./T_7PTS,xi7PTS/xi7PTS(2),xi7PTS_error/xi7PTS(2),'b','Linestyle','-')
plot(1./T_7pair,xi7pair/xi7pair(2),'b','Linestyle','--')
errorbar(1./T_8PTS,xi8PTS/xi8PTS(2),xi8PTS_error/xi8PTS(2),'b','Linestyle','-')
plot(1./T_8pair,xi8pair/xi8pair(2),'b','Linestyle','--')
errorbar(1./T_12PTS,xi12PTS/xi12PTS(2),xi12PTS_error/xi12PTS(2),'b','Linestyle','-')
plot(1./T_12pair,xi12pair/xi12pair(2),'b','Linestyle','--')


xlim([1,2.5])
ylim([0.5,2.2])
xlabel('$1/T$','Interpreter','latex')
ylabel('$\xi/\xi_0$','Interpreter','latex');
axis square
set(gcf, 'PaperPositionMode', 'auto');
print -depsc2 Figure3g.eps


end