function makeEpsSet22(filename)

load(filename);
close all

figure(1)
set(gcf,'paperunits','centimeters','papersize',[21 6],'paperposition',[0 0 21 6]);

hoffset = .07;
voffset = .16;
% voffset2 = .05;
fw = .28;
fh = .73;

steps = length(interaction{1}(100:end));

part1 = 1:(ceil(steps/4));
part2 = ((ceil(steps/4))+1):(ceil(2*steps/4));
part3 = (ceil(2*steps/4)+1):(ceil(3*steps/4));
part4 = (ceil(3*steps/4)+1):steps;

axes('position',[0+hoffset voffset fw fh]);
[hits1,bins] = hist(log2(interaction{1}(part1)),0:10); %#ok<*USENS>
[hits2,bins] = hist(log2(interaction{1}(part2)),0:10); %#ok<*USENS>
[hits3,bins] = hist(log2(interaction{1}(part3)),0:10); %#ok<*USENS>
[hits4,bins] = hist(log2(interaction{1}(part4)),0:10); %#ok<*USENS>
ymax = max(max(hits1),max(hits2))*1.1;
% mu1 = log2(mean(interaction{1}(100:end)));

dw = .2;
h = bar(bins-1.5*dw,[hits1(:)],1);
hold on
set(h,'facecolor',[.9 .9 .9],'barwidth',dw);
h = bar(bins-dw/2,[hits2(:)],1);
set(h,'facecolor',[.5 .5 .5],'barwidth',dw);
h = bar(bins+dw/2,[hits3(:)],1);
set(h,'facecolor',[.2 .2 .2],'barwidth',dw);
h = bar(bins+1.5*dw,[hits4(:)],1);
set(h,'facecolor',[0 0 0],'barwidth',dw);
grid on
hold off

title('ARPF');
xlabel('$K_n$','interpreter','latex');
ylabel('Count');
set(gca,'ylim',[0 ymax],'xlim',[-.5 10.5],'box','on','fontsize',9);
set(gca,'box','off');

axes('position',[.31+hoffset voffset fw fh]);
[hits1,bins] = hist(log2(interaction{2}(part1)),0:10); %#ok<*USENS>
[hits2,bins] = hist(log2(interaction{2}(part2)),0:10); %#ok<*USENS>
[hits3,bins] = hist(log2(interaction{2}(part3)),0:10); %#ok<*USENS>
[hits4,bins] = hist(log2(interaction{2}(part4)),0:10); %#ok<*USENS>

h = bar(bins-1.5*dw,[hits1(:)],1);
hold on
set(h,'facecolor',[.9 .9 .9],'barwidth',dw);
h = bar(bins-dw/2,[hits2(:)],1);
set(h,'facecolor',[.5 .5 .5],'barwidth',dw);
h = bar(bins+dw/2,[hits3(:)],1);
set(h,'facecolor',[.2 .2 .2],'barwidth',dw);
h = bar(bins+1.5*dw,[hits4(:)],1);
set(h,'facecolor',[0 0 0],'barwidth',dw);
grid on
hold off

title('Simple');
xlabel('$K_n$','interpreter','latex');
set(gca,'ylim',[0 ymax],'yticklabel',{},'xlim',[-.5 10.5],'fontsize',9);
set(gca,'box','off');

axes('position',[.62+hoffset voffset fw fh]);
[hits1,bins] = hist(log2(interaction{3}(part1)),0:10); %#ok<*USENS>
[hits2,bins] = hist(log2(interaction{3}(part2)),0:10); %#ok<*USENS>
[hits3,bins] = hist(log2(interaction{3}(part3)),0:10); %#ok<*USENS>
[hits4,bins] = hist(log2(interaction{3}(part4)),0:10); %#ok<*USENS>

h = bar(bins-1.5*dw,[hits1(:)],1);
hold on
set(h,'facecolor',[.9 .9 .9],'barwidth',dw);
h = bar(bins-dw/2,[hits2(:)],1);
set(h,'facecolor',[.5 .5 .5],'barwidth',dw);
h = bar(bins+dw/2,[hits3(:)],1);
set(h,'facecolor',[.2 .2 .2],'barwidth',dw);
h = bar(bins+1.5*dw,[hits4(:)],1);
set(h,'facecolor',[0 0 0],'barwidth',dw);
grid on
hold off

title('Greedy');
xlabel('$K_n$','interpreter','latex');
set(gca,'ylim',[0 ymax],'yticklabel',{},'xlim',[-.5 10.5],'fontsize',9);
set(gca,'box','off');

print -depsc -loose figure2

return
