function []=figure1c()
close all
figure()
hold on
%% Loading data
tau0=2479;
% ceiling_T_VFT=0.0907;
% ceiling_T_A=0.0738;
% T_g=mean([ceiling_T_VFT,ceiling_T_A]);
% T_MCT=0.123458;

% local MC
data_name=sprintf('./tau_psi_1000_local.dat');
fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Ndata=round(size(numbers,1)/2);
Ts_local=numbers(2*(1:Ndata)-1);
taus_local=numbers(2*(1:Ndata));
T_inverses_local=1./Ts_local;

% Fit
T_fit_inverse=0:0.01:40;
T_fit=1./T_fit_inverse;
lnA_local=-7.25558;
B_local=3.14764;
taufit_local=(exp(lnA_local))*exp(B_local./(T_fit));

% swap MC
data_name=sprintf('./tau_psi_1000_swap.dat');
fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Ndata=round(size(numbers,1)/2);
Ts_swap=numbers(2*(1:Ndata)-1);
taus_swap=numbers(2*(1:Ndata));
T_inverses_swap=1./Ts_swap;
%% Plotting
scatter(T_inverses_local,taus_local/tau0,100,[0,0,0],'filled')
scatter(T_inverses_swap,taus_swap/tau0,100,[0.9,0,0],'filled','marker','s')
plot(T_fit_inverse,taufit_local/tau0,'Color',[0,0,0])
xlim([0,40])
ylim([10^(-1),10^(6)])
xlabel('$1/T$','Interpreter','latex');
ylabel('$\tau_{\alpha}/\tau_{0}$','Interpreter','latex')
set(gca,'Yscale','Log');
print(gcf, '-dpdf', 'figure1c.pdf');
end