function []=figure3a()
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
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;

%% Fitting
x_grain=0.02*(T_g);
Ts_fit=0:x_grain:T_g;
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;
Sconf_fits4=Sconf_fits;

%% Plotting
scatter(Sconfs4,Ts4,'Marker','s')
line4=line(Sconf_fits4,Ts_fit);
set(line4,'Color',[0,0,0.7],'LineStyle',':');
xlim([0,0.5])
ylim([0,0.15])
xlabel('$1/R$','Interpreter','latex');
ylabel('$T$','Interpreter','latex')
print(gcf, '-dpdf', 'figure3a.pdf');
end