comp_luke_coll2_jd

0001 clear
0002 p_opt = -1;
0003 c_opt = true;
0004 
0005 lukedata = load('DKE_RESULTS_Runaway_alpha3_sync_0p6.mat');
0006 lukedata_2 = load('DKE_RESULTS_Runaway_alpha3_sync_0p6_coll2.mat');
0007 
0008 % xlog = 0;
0009 % xlim = [0,max(Ec)];
0010 % xtick = 0:20:200;
0011 ylim = 10.^[-20,0];%max([lukedata.f{1}(:);lukedata_2.f{1}(:)]);
0012 ytick = 10.^(-20:4:0);
0013 
0014 xlog = 1;
0015 xlim = 10.^[-3,2]; 
0016 xtick = 10.^(-3:2);
0017 
0018 
0019 if c_opt,
0020     color1 = 'b';
0021     color2 = 'r';
0022     style1 = '-';
0023     style2 = '--';
0024     size1 = 2;
0025     size2 = 2;
0026     colstr = '_col';
0027 else
0028     color1 = 'k';
0029     color2 = 'k';
0030     style1 = '-';
0031     style2 = '--';
0032     size1 = 0.5;
0033     size2 = 2;
0034     colstr = '';
0035 end
0036 
0037 % PLOT
0038 figure(11),clf
0039 
0040 xlab = '$\mathcal{E}_k$ (MeV)';
0041 ylab = '$f(\mathcal{E}_k,\xi = 1)$';
0042 tit = '';
0043 leg = {'with I_{FP}','without I_{FP}'};
0044 siz = 20+14i;
0045 
0046 graph1D_jd(lukedata_2.Ec,lukedata_2.f{end}(:,end),xlog,1,'','','',NaN,xlim,ylim,style1,'none',color1,size1,siz,gca);
0047 graph1D_jd(lukedata.Ec,lukedata.f{end}(:,end),xlog,1,xlab,ylab,tit,leg,xlim,ylim,style2,'none',color2,size2,siz,gca,0.9,0.7,0.7);
0048 
0049 set(gca,'xtick',xtick)
0050 set(gca,'ytick',ytick)
0051 
0052 set(gca,'YMinorGrid','off')
0053 set(gca,'YMinorTick','on')
0054 if xlog,
0055     set(gca,'XMinorGrid','off')
0056     set(gca,'XMinorTick','on')
0057 end
0058 set(get(gca,'xlabel'),'interpreter','latex')
0059 set(get(gca,'ylabel'),'interpreter','latex')
0060 
0061 print_jd(p_opt,['fig_luke_coll2',colstr],'figures/',11)