Commit 413dbc0a authored by Simone Liuzzo's avatar Simone Liuzzo

corrected cell naming, fixed figures for plot

parent b4814f88
......@@ -14,7 +14,7 @@ dkk=(k-k0l)./k0l;
if largefig
figure;
figure(41);
subplot(2,1,1);
extfigax = gca;
bar(extfigax,[dkk]);
......@@ -49,7 +49,9 @@ bar(ax,[cor_per_cell_quad],'EdgeColor','none');
ylabel(ax,'\DeltaK/K');
grid(ax,'on');
set(ax,'Xlim',[0 32+1],'FontSize',6);
ax.XTick=[1:32];
ax.XTickLabel=arrayfun(@(a)num2str(a),[4:32,1:3],'un',0);
ax.XTickLabelRotation=90;
......
......@@ -82,8 +82,8 @@ end
if isfield(handles,'axes4') % displays corr. strengths
if get(handles.largefigures,'Value')
figure; extfigax = gca;
bar(extfigax,qemb.cor);
figure(51); extfigax = gca;
plot(extfigax,qemb.cor);
legend(extfigax,'initial','corrected');
title(extfigax,'Corrector strengths');
grid(extfigax,'on');
......@@ -93,6 +93,9 @@ if isfield(handles,'axes4') % displays corr. strengths
cor_per_cell_std = squeeze(std(reshape(qemb.cor((1+extracor/2):(Nc-extracor/2),:),32,[],2),1,2));
% cor_per_cell_mean = squeeze(mean(reshape(qemb.cor((1+extracor/2):(Nc-extracor/2),:),32,[],2),2));
bar(handles.axes4,cor_per_cell_std,'EdgeColor','none');
handles.axes4.XTick=[1:32];
handles.axes4.XTickLabel=arrayfun(@(a)num2str(a),[4:32,1:3],'un',0);
handles.axes4.XTickLabelRotation=90;
title(handles.axes4,'std Corrector strengths (per cell)');
legend(handles.axes4,'initial','corrected');
grid(handles.axes4,'on');
......
......@@ -46,8 +46,9 @@ end
if isfield(handles,'axes4') % displays corr. strengths
if get(handles.largefigure,'Value')
figure; extfigax = gca;
bar(extfigax,qemb.skewcor);
figure(21);
extfigax = gca;
plot(extfigax,qemb.skewcor);
legend(extfigax,'initial','corrected');
title(extfigax,'Corrector strengths');
grid(extfigax,'on');
......@@ -60,6 +61,10 @@ if isfield(handles,'axes4') % displays corr. strengths
32,[],size(qemb.skewcor,2)),1,2));
% cor_per_cell_mean = squeeze(mean(reshape(qemb.cor((1+extracor/2):(Nc-extracor/2),:),32,[],2),2));
bar(handles.axes4,cor_per_cell_std,'EdgeColor','none');
handles.axes4.XTick=[1:32];
handles.axes4.XTickLabel=arrayfun(@(a)num2str(a),[4:32,1:3],'un',0);
handles.axes4.XTickLabelRotation=90;
title(handles.axes4,'std Corrector strengths (per cell)');
legend(handles.axes4,'initial','corrected');
grid(handles.axes4,'on');
......
......@@ -33,7 +33,7 @@ if nargin >= 4 % displays quad rotations
if largefig
figure;
figure(31);
subplot(2,1,1); extfigax = gca;
plot(extfigax,tilts);
......@@ -68,6 +68,10 @@ if nargin >= 4 % displays quad rotations
ylabel(ax,'\Theta');
yyaxis(ax,'right');
bar(ax,[cor_per_cell_dip],'EdgeColor','none');
ax.XTick=[1:32];
ax.XTickLabel=arrayfun(@(a)num2str(a),[4:32,1:3],'un',0);
ax.XTickLabelRotation=90;
legend(ax,'std quad','std dip');
grid(ax,'on');
set(ax,'Xlim',[0 32+1],'FontSize',6);
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment