function show_paper_figure(dataset)
% dataset is 1, 12, 123, or 1234 for 1x1, 2x2, 3x3, 4x4

dataOrigEps2=load_xgrow_data('coptdata_orig0.2N');              % original tile set
dataOrigGse65=load_xgrow_data('renatFeb12/orig_low.data');   % these are 100 runs, 10000 tiles
dataOrigGse85=load_xgrow_data('renatFeb12/orig_high.data');  % these are 100 runs, 10000 tiles
dataRSqEps1=load_xgrow_data('2x2RSqEps0.1All');                 % full no-slippage tile set, 2x2RSq, variable number of trials
dataRSqGse65=load_xgrow_data('2x2RSqGse6.5All');                % full no-slippage tile set, 2x2RSq, variable number of trials
dataRSqGse65b=load_xgrow_data('2x2GseRSq6.5BIGeps');            % full no-slippage tile set, 2x2RSq, variable number of trials
dataRSqGse75=load_xgrow_data('2x2RSqGse7.5All');                % full no-slippage tile set, 2x2RSq, variable number of trials
dataRSqGse75b=load_xgrow_data('2x2GseRSq7.5BIGeps');            % full no-slippage tile set, 2x2RSq, variable number of trials
dataRSqGse75b=dataRSqGse75b(162:end,:);        %% let's avoid cluttering with high-error, high eps data
%data3x3Eps2=load_xgrow_data('coptdata_3x3proof0.2All');         % full no-slippage tile set, 3x3, 30 trials
data3x3Eps15=load_xgrow_data('3x3Eps0.15All');                  % full no-slippage tile set, 3x3, in progress
data3x3Gse55=load_xgrow_data('3x3Gse5.5All');                   % full no-slippage tile set, 3x3, variable number of trials
data3x3Gse65=load_xgrow_data('3x3Gse6.5All');                   % full no-slippage tile set, 3x3, variable number of trials

%data4x4Eps15=load_xgrow_data('4x4Eps0.15All');                  % full no-slippage tile set, 4x4, in progress
data4x4Gse55=load_xgrow_data('4x4Gse5.5All');                   % full no-slippage tile set, 4x4, in prog
data4x4Gse60=load_xgrow_data('4x4Gse6.0All');                   % full no-slippage tile set, 4x4, in prog

clf
loglog([1e-4 1e-4 10 10], [1e-7 5e-2 1e-7 5e-2], 'k.'); 

%% lines identical to what SAkinplot would do, but underneath other data
h=line([1e-4 10], [.98e-4 3.1e-2]); set(h,'Color',[0 0 0]);
if dataset>1
  h=line([1.8e-3 10], [10^(-7) 5.5*10^(-4)]); set(h,'Color',[0 0 0]);
  if dataset>12
    h=line([.03 10], [10^(-7) 7*10^(-5)]); set(h,'Color',[0 0 0]);
    %line([.1 10], [10^(-7) 10^(-4)])
  end
end

%%%%% draw it pretty
%hold off
hold on; legendH=[];

[h,so]=plot_speed_vs_errors(dataOrigEps2,1e-7,'g ',1-1); legendH=[legendH h];     % linefit=1?
[h,soea]=plot_speed_vs_errors(dataOrigGse65,1e-7,'b ',0); legendH=[legendH h];
[h,soeb]=plot_speed_vs_errors(dataOrigGse85,1e-7,'r ',0); legendH=[legendH h];
%text(1,3e-2,['orig slope=' num2str(so)])
text(2e-3,3e-3,['1x1 slope=' num2str(so)])

if dataset>1
 [h,sp4]=plot_speed_vs_errors(dataRSqEps1,1e-7,'g ',1-1); legendH=[legendH h];     % linefit=1?
 [h,spe4a]=plot_speed_vs_errors(dataRSqGse65,1e-7,'b ',0); legendH=[legendH h];
 [h,spe4aa]=plot_speed_vs_errors(dataRSqGse65b,1e-7,'b ',0); legendH=[legendH h];
 [h,spe4b]=plot_speed_vs_errors(dataRSqGse75,1e-7,'r ',0); legendH=[legendH h];
 [h,spe4bb]=plot_speed_vs_errors(dataRSqGse75b,1e-7,'r ',0); legendH=[legendH h];
 %text(1,4e-4,['RSq Mar24 slope=' num2str(sp4)])
 text(.9,8e-4,['2x2 slope=' num2str(sp4)])

 if dataset > 12
  [h,sp3x3]=plot_speed_vs_errors(data3x3Eps15,1e-7,'g ',1); legendH=[legendH h];
  [h,spe3x3a]=plot_speed_vs_errors(data3x3Gse55,1e-7,'b ',0); legendH=[legendH h];
  [h,spe3x3b]=plot_speed_vs_errors(data3x3Gse65,1e-7,'r ',0); legendH=[legendH h];
  %text(1,1e-6,['3x3 June slope=' num2str(sp3x3)])
  sp3x3=(log(7e-5)-log(1e-7))  / (log(10)-log(.03)) ;
  text(1,3e-6,['3x3 slope=' num2str(sp3x3)])
 end

 if dataset > 123
  %[h,sp4x4]=plot_speed_vs_errors(data4x4Eps15,1e-7,'g ',1); legendH=[legendH h];
  [h,spe4x4a]=plot_speed_vs_errors(data4x4Gse55,1e-7,'b ',0); legendH=[legendH h];
  [h,spe4x4b]=plot_speed_vs_errors(data4x4Gse60,1e-7,'r ',0); legendH=[legendH h];
  %text(1,1e-6,['4x4 June slope=' num2str(sp4x4)])
  sp4x4=(log(7e-5)-log(1e-7))  / (log(10)-log(.03)) ;
  text(1,3e-6,['4x4 slope=' num2str(sp4x4)])
 end

end

if 0
legend(legendH,...
       '1x1, eps=0.2 (size=100000, 20 flakes)', ...
       '1x1, Gse=6.5 (size=10000, 100 flakes)',  ...
       '1x1, Gse=8.5 (size=10000, 100 flakes)',  ...
       '2x2, eps=0.1 (size=100000, ~200 flakes)',  ...
       '2x2, Gse=6.5 (size=100000, ~100 flakes)',  ...
       '2x2, Gse=7.5 (size=100000, ~250 flakes)',  ...
       '3x3, eps=.15 (size=100000, ~20 flakes?)',  ...
       '3x3, Gse=5.5 (size=100000, ~20 flakes)',  ...
       '3x3, Gse=6.5 (size=100000, ~60 flakes)',  ...
       '4x4, eps=.15 (size=100000, ~20 flakes?)',  ...
       '4x4, Gse=5.5 (size=100000, ~20 flakes?)',  ...
       '4x4, Gse=6.0 (size=100000, ~60 flakes?)',  ...
       2 );
end

axis([1e-4 10 1e-7 5e-2])


if dataset>1
 SAkinplotAllMp([6.5 8.5; 6.5 7.5], [0.05:0.1:3.95; 0.03:0.03:1.20], ...
                [1/3 0], [], [2 12], []);
%                [1/3 0], [], [2 12], [1e-4 10]);

 h=title('{\bf Simulation and Kinetic Trap Predictions for Proofreading Tiles}');
 set(h,'FontSize',15);
else
 SAkinplotAllMp([6.5 8.5], [0.05:0.1:.95], ...
                [1/3 0], [], [], []);
%                [1/3 0], [], [], [1e-4 10]);

 h=title('{\bf Simulation and Kinetic Trap Predictions}');
 set(h,'FontSize',15);
end

h=xlabel('\bf speed (layers/sec)'); 
set(h,'FontSize',15);
h=ylabel('\bf error rate (mismatches/tile)');
set(h,'FontSize',15);