MomAperture_allRing

 all Ring momentum aperture
points=1:10:length(THERING);

0001 function [map_l,map_h]=MomAperture_allRing(THERING,points,varargin)
0002 % all Ring momentum aperture
0003 %points=1:10:length(THERING);
0004 map_h=zeros(length(points),1);
0005 map_l=zeros(length(points),1);
0006 
0007 if numel(varargin)>0
0008     nturns=varargin{1};
0009 else
0010     nturns=100;
0011 end
0012 
0013 for i=1:length(points)
0014    % disp([i length(points) i/length(points)*100])
0015     %cycle ring
0016      THERING_cycl=[THERING(points(i):end); THERING(1:points(i)-1)]';
0017         try
0018             map_h(i)=momentum_aperture_at(THERING_cycl,+0.1,[10^-6 10^-6],0,0,+0.01,2,10,nturns);
0019         catch
0020             map_h(i)=0;
0021         end
0022         
0023         try
0024             map_l(i)=momentum_aperture_at(THERING_cycl,-0.1,[10^-6 10^-6],0,0,-0.01,2,10,nturns);
0025         catch
0026             map_l(i)=0;
0027         end
0028        
0029 
0030 end
0031 
0032 
0033 return