% DNNI model around kappa=0.5 multicritical point

close all;
msize = 8;
psize = 20;

widthpix = 250;
heightpix = 200;
linewidth = 1.25;
leftmargin = 95;
bottommargin = 65;
topmargin = 40;
rightmargin = 105;
centermargin = 0;
widthtotal = widthpix*3+leftmargin+rightmargin;
heighttotal = heightpix*2+topmargin+bottommargin;
left = leftmargin/widthtotal;
bottom = bottommargin/heighttotal;
width = widthpix/widthtotal;
height = heightpix/heighttotal;


%% main plot
panelA = axes('Position', [left, bottom, width*2, height*2]);

kappas = [ 0.48, 0.485, 0.49, 0.498, 0.499, 0.5, 0.502];
data = dlmread('../data/2DXDT/heatcap_0.dat');

for rp=1:numel(kappas)
    kappa = kappas(rp);
    kappastr = num2str(kappa);
    vals = data(data(:,1) == kappa, [3, 5]);
    plot(vals(:,1), vals(:,2), '-x');
    if 1==rp, hold on, end
end
xlim([7.8, 36.2]);
ylim([0.5, 1.4]);
xticks([10,20,30]);
yticks([0.5,0.8,1.1,1.4]);
xlabel('$L$', 'interpreter', 'latex');
ylabel('$c^*$', 'interpreter', 'latex');
h=legend(cellfun(@num2str, num2cell(kappas), 'UniformOutput', false), 'location', 'northwest', 'Box','off');
set(h, 'Position', get(h, 'Position')+ [0.07 -0.06 0 0]);
set(gca, 'fontname', 'times new roman', 'fontsize', psize);

%% panel (b)
nowpanel = axes('Position', [left+width*2, bottom+height, width, height]);

kappa = 0.485;
hs = 0;
kappastr = num2str(kappa);

sizesA = 8:4:32;

n=numel(sizesA);
rpp=(1:n)';
corders = [(rpp-1)/(n-1), zeros(n,1), 1 - (rpp-1)/(n-1)];

lgds = {};

for rp=1:n
    asize = sizesA(rp);
    fname = strcat('../data/2DXDT/clength_', kappastr, '_', num2str(hs) , '_', num2str(asize), '.dat');
    if ~exist(fname, 'file'), continue, end
    data = dlmread(fname);
    plot(1./data(:,2), abs(data(:, 6))/asize, 'color', corders(rp,:));
    if 1==rp, hold on, end
end
plot(0.44465, 0.3778, '*', 'color', corders(end,:), 'markersize', msize);
plot([0.433,0.433], [0.12,1], '--k');
xticks([0.1,0.4,7]);
yticks([0, 0.3, 0.6]);
xlim([0.1, 0.7]);
ylim([0, 0.6]);
set(gca, 'xticklabel', []);
%xlabel('$T$', 'interpreter', 'latex');
ylabel('$\xi_1/L$', 'interpreter', 'latex');
xx = gca;
xx.XRuler.TickLength = [0.02, 0];
xx.YRuler.TickLength = [0.03, 0];

set(gca, 'XAxisLocation', 'top', 'YAxisLocation', 'right');
set(gca, 'fontname', 'times new roman', 'fontsize', psize);

%% panel (c)
nowpanel = axes('Position', [left+width*2, bottom, width, height]);

kappa = 0.5;
hs = 0;
kappastr = num2str(kappa);

sizesA = 8:4:32;

n=numel(sizesA);
rpp=(1:n)';
corders = [(rpp-1)/(n-1), zeros(n,1), 1 - (rpp-1)/(n-1)];

lgds = {};

for rp=1:n
    asize = sizesA(rp);
    fname = strcat('../data/2DXDT/clength_', kappastr, '_', num2str(hs) , '_', num2str(asize), '.dat');
    if ~exist(fname, 'file'), continue, end
    data = dlmread(fname);
    xs = [1./data(:,2); 0];
    ys = [ abs(data(:, 6)); 1.44];
    plot(xs, ys/asize, 'color', corders(rp,:));
    if 1==rp, hold on, end
end
%plot([0.972,0.972], [-2,10], '--k');
xticks([0.1,0.4,0.7]);
yticks([0, 0.3]);
xlim([0.1, 0.7]);
ylim([0, 0.6]);
%set(gca, 'yscale', 'log');
xlabel('$T$', 'interpreter', 'latex');
ylabel('$\xi_1/L$', 'interpreter', 'latex');
xx = gca;
xx.XRuler.TickLength = [0.02, 0];
xx.YRuler.TickLength = [0.03, 0];
set(gca, 'YAxisLocation', 'right');
set(gca, 'fontname', 'times new roman', 'fontsize', psize);

%% annotations

annotation('textbox',[0.10, 0.82, 0.1, 0.1], 'String', '(a)', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize);
annotation('textbox',[0.832, 0.82, 0.1, 0.1], 'String', '(b)', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize);
annotation('textbox',[0.832, 0.43, 0.1, 0.1], 'String', '(c)', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize);

annotation('textbox',[0.63, 0.5, 0.1, 0.1], 'String', '$\kappa=0.485$', 'interpreter', 'latex', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize);
annotation('textbox',[0.63, 0.422, 0.1, 0.1], 'String', '$\kappa=0.5$', 'interpreter', 'latex', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize);

print('../figures/artfig_DNNImcritical.eps', '-depsc');