figure(1) %%%%I want to have a plot with fixed dimension for icons in legend hplot = plot(yH(1,:),xH(1,:),'^', yC(:,1),xC(:,1),'*',yC(:,2),xC(:,2),'*',... yC(:,3),xC(:,3),'*',yC(:,4),xC(:,4),'*',yC(:,5),xC(:,5),'*',... yC(:,6),xC(:,6),'*','MarkerSize',s); % Markersize: specifys the size of the marker in points (s in questo caso) hold on plot(Ymcporigine,Xmcporigine,'k-','MarkerEdgeColor','k','MarkerSize',1); %Plot contorno MCP hold on plot(Yh, Xh, 'b-', 'MarkerSize', s); %Plot alone circolare hold off %Label main axes - It is necessary to give the label instructions after plot in order to avoid overlap xlabel(gca, 'Deflessione magnetica [m]'); % label lower x axis ylabel(gca,'Deflessione elettrica [m]'); %label left y axis %particles outside MCP radius won't be appear in figure xlim([0, Rmcp]) ylim([0, Rmcp]) %%%% Legend l=legend(hplot, 'H^+','C^+','C^{+2}','C^{+3}','C^{+4}','C^{+5}','C^{+6}', 'Location','BestOutside'); a=get(l,'children'); set(a(1:3:end),'MarkerSize',10); %%%% Secindary x axis %xlabel(gca, 'Deflessione magnetica [m]'); % label asse x principale %set secondary x limit as the momentum of a H+ at distance equal to the MCP radius % Secondo Harres y=(q*B*LB*L)/sqrt(2mEkin) ==> mv=q*B*LB*L/y mv_max = (q*B*LB*L)/Rmcp; %mv_max = 1; %Layout instruction set(gca,'Box','off'); % Turn off the box surrounding the whole axes axesPosition = get(gca,'Position'); %# Get the current axes position hNewAxes = axes('Position',axesPosition,... %# Place a new axes on top... 'Color','none',... %# ... with no background color 'XAxisLocation','top',... %# ... located on the top 'Ytick', [],... %# ... with no y tick marks 'Xlim', [0, mv_max],... %# ... should define x axis scale (need to set xmax = mv_max) 'Box','off'); %# ... and no surrounding box xlabel(hNewAxes,'Momentum (H^+)'); %# Add a label to the top axis set(gca, 'XTickLabel', num2str(get(gca,'XTick')','%d')) %%%%%Plot title - It is necessary to give the title instruction after secondary x axis in order to avoid overlap title(['Calcolo approssimato interazione ioni campo magnetico B=', num2str(B), 'Tesla']); %%%% Secondary y axis %ylabel(gca,'Deflessione elettrica [m]'); %label asse y principale %set secondary y limit as the energy of a H+ at distance equal to the MCP radius % Secondo Harres x=(q*E*Le*L)/(2mEkin) ==> Ekin=q*E*Le*L/2mx Le = 0.07; %Estensione del C.E. : 70 mm E = 100000; %campo TP.m Ekin_max = (q*E*Le*L)/(2*m_H*Rmcp); %mv_max = 1; set(gca,'Box','off'); % Turn off the box surrounding the whole axes axesPosition = get(gca,'Position'); %# Get the current axes position hNewAxes = axes('Position',axesPosition,... %# Place a new axes on top... 'Color','none',... %# ... with no background color 'YAxisLocation','right',... %# ... located on the right 'Ylim', [0, Ekin_max],... %# ... should define y axis scale (need to set ymax=Ekin_max) 'Xtick', [],... %# ... with no y tick marks 'Box','off'); %# ... and no surrounding box ylabel(hNewAxes,'Energy (H^+)'); %# Add a label to the top axis set(gca, 'YTickLabel', num2str(get(gca,'YTick')','%d'))