%% Generate an undecorated version of Figure 3f

function []=Figure3f()
ylimhigh=0.0003; % y-axis range
angular_index=8; % angular momentum

figure();
hold on
for temperature_index=1:5
    % temperature_index=1 corresponds to temperature T1=1.000;
    % temperature_index=2 corresponds to temperature T2=0.800;
    % temperature_index=3 corresponds to temperature T3=0.600;
    % temperature_index=4 corresponds to temperature T4=0.510;
    % temperature_index=5 corresponds to temperature T5=0.450;
    
    %% Read in the data
    data_name=sprintf('./l=%dChi_T%d.dat',angular_index,temperature_index); fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
    Nr=round(size(numbers,1)/3);% number of radial data points
    numbers=reshape(numbers,3,Nr);
    R=numbers(1,1:Nr).'; % cavity radii where positional point-to-set correlation functions are measured
    gPTS=numbers(2,1:Nr).'; % positional point-to-set correlation functions
    gPTS_error=numbers(3,1:Nr).'; % error bars
      

    %% Plot them
    errorbar(R,gPTS,gPTS_error,'Linestyle','-','Marker','.')

end
xlabel('$R$','Interpreter','latex')
ylabel('$\chi^{\rm PTS}_{8}$','Interpreter','latex');
xlim([1,4])
ylim([0,ylimhigh])
axis square
set(gcf, 'PaperPositionMode', 'auto');
print -depsc2 Figure3f.eps
    
end