plPolynomBComp

 PLPOLYNOMBCOMP default plotting function for ATPLOT

 Plots polynomB for ring and ring1

0001 function plotdata=plPolynomBComp(lindata,ring,dpp,varargin) %#ok<INUSD>
0002 % PLPOLYNOMBCOMP default plotting function for ATPLOT
0003 %
0004 % Plots polynomB for ring and ring1
0005 
0006 CoD=cat(2,lindata.ClosedOrbit);
0007 PolyBVal=zeros(size(CoD(1,:)));
0008 PolyBVal1=zeros(size(CoD(1,:)));
0009 
0010 PolynomBVal1=zeros(size(CoD(1,:)));
0011 PolynomBVal2=zeros(size(CoD(1,:)));
0012 PolynomBVal3=zeros(size(CoD(1,:)));
0013 PolynomBVal4=zeros(size(CoD(1,:)));
0014 ind=findcells(ring,'PolynomB');
0015 
0016 PolynomBVal1(ind)=getcellstruct(ring,'PolynomB',ind,1,1);
0017 PolynomBVal2(ind)=getcellstruct(ring,'PolynomB',ind,1,2);
0018 PolynomBVal3(ind)=getcellstruct(ring,'PolynomB',ind,1,3);
0019 PolynomBVal4(ind)=getcellstruct(ring,'PolynomB',ind,1,4);
0020 
0021 
0022 plotdata(1).values=[PolynomBVal1' PolynomBVal2'/10 ...];%...
0023                      PolynomBVal3'/100 PolynomBVal4'/5000 ];
0024 plotdata(1).labels={'dipole','quadrupole/10',...
0025                     'sextuople/100','octupole/5000'};
0026 plotdata(1).axislabel='PolynomB';
0027 % dispersion=cat(2,lindata.Dispersion)';
0028 % plotdata(2).values=dispersion(:,3);
0029 % plotdata(2).labels={'\eta_y'};
0030 % plotdata(2).axislabel='vertical dispersion [m]';
0031 
0032 end