% Plot comparison between NN and NNN transfer matrix results

function Figure2()
    %% compare xi1 for NN and NNN
    close all;
    addpath('../helpfunctions');
    widthpix = 375;
    heightpix = 300;
    leftmargin = 50;
    bottommargin = 50;
    topmargin = 10;
    rightmargin = 30;
    centermargin = 80;
    widthtotal = widthpix+leftmargin+rightmargin;
    heighttotal = heightpix*2+topmargin+bottommargin;
    left = leftmargin/widthtotal;
    bottom = bottommargin/heighttotal;
    cmwidth = centermargin/widthtotal;
    cmheight = centermargin/heighttotal;
    width = widthpix/widthtotal;
    height = heightpix/heighttotal;
    psize=2;
    fig = figure('rend','painters','pos',[100 100 widthtotal heighttotal]);
    axbg = axes('Position',[0,0,1,1], 'Fontname', 'Times New Roman', 'visible', 'off');
    text(axbg, 0.01, 0.99, '(a)', 'Fontname', 'Times New Roman','fontsize',18);
    text(axbg, 0.01, 0.52, '(b)', 'Fontname', 'Times New Roman','fontsize',18);
    
    %% 3D
    ax2 = axes('Position', [left,bottom+(1-1)*(height+cmheight),width,height], 'Fontname', 'Times New Roman','fontsize',16);
    tdNN = dlmread('../data3D/xi3DD1.5_NN.dat','',1,0);
    tdNNN = dlmread('../data3D/xi3DD1.5.dat','',1,0);
    plot(tdNN(:,1), tdNN(:,3), '-k', tdNNN(:,1), tdNNN(:,3), 'ok');
    xlabel('$P_{x}^*$', 'Interpreter', 'Latex');
    ylabel('$\xi_r$', 'Interpreter', 'Latex');
    xticks(0:10:50);yticks(0:0.1:0.4);
    unifyfigureend;
    
    % inset draw NNN/NN 
    ax2small = axes('Position', [left+width*0.45,bottom+(1-1)*(height+cmheight)+height*0.2,width*0.45,height*0.45], 'Fontname', 'Times New Roman','fontsize',16);
    plot(tdNN(:,1), tdNNN(:,3) ./ tdNN(:,3)-1, '-k');
    xlabel('$P_{x}^*$', 'Interpreter', 'Latex');
    ylabel('$\xi_r^{\mathrm{NNN}} / \xi_r^{\mathrm{NNN}}-1$', 'Interpreter', 'Latex');
    xlim([0 50]);ylim([-0.01 0.01]);
    xticks(0:25:50);yticks(-0.01:0.01:0.01);
    unifyfigureend;
    set(ax2small, 'FontSize', 14);
    
    %% 2D
    ax1 = axes('Position', [left,bottom+(2-1)*(height),width,height], 'Fontname', 'Times New Roman','fontsize',16);
    sdNN = dlmread('../data2D/xiD1.5_NN.dat','',1,0);
    sdNNN = dlmread('../data2D/xiD1.5.dat','',1,0);
    plot(sdNN(:,1), sdNN(:,2), '-k', sdNNN(:,1), sdNNN(:,3), 'ok');
    xlim([0 50]);
    %xlabel('$P_{x}^*$', 'Interpreter', 'Latex');
    ylabel('$\xi_y$', 'Interpreter', 'Latex');
    set(gca, 'XTickLabel', []);
    yticks(0:20:80);
    unifyfigureend;
    
    % inset draw NNN/NN 
    ax1small = axes('Position', [left+width*0.2,bottom+(2-1)*(height)+height*0.5,width*0.45,height*0.45], 'Fontname', 'Times New Roman','fontsize',16);
    plot(sdNN(:,1), sdNNN(:,3) ./ sdNN(:,2)-1, '-k');
    xlim([0 50]);ylim([0 0.04]);
    xticks(0:25:50);yticks(0:0.02:0.04);
    xlabel('$P_{x}^*$', 'Interpreter', 'Latex');
    ylabel('$\xi_y^{\mathrm{NNN}} / \xi_y^{\mathrm{NNN}}-1$', 'Interpreter', 'Latex');
    unifyfigureend;
    set(ax1small, 'FontSize', 14);
    
    print('Figure2.eps','-depsc');
end