pltouschek

plotdata=pltouschek(lindata,ring,dpp)
 plots curly H function and also 1/(sigx sigy) to understand local scattering rate

0001 function plotdata=pltouschek(lindata,ring,dpp) %#ok<INUSD>
0002 %plotdata=pltouschek(lindata,ring,dpp)
0003 % plots curly H function and also 1/(sigx sigy) to understand local scattering rate
0004 beta=cat(1,lindata.beta);
0005 alpha=cat(1,lindata.alpha);
0006 betax=beta(:,1);
0007 betaz=beta(:,2);
0008 alphax=alpha(:,1);
0009 gammax=(1+alphax.*alphax)./betax;
0010 eta=cat(2,lindata.Dispersion)';
0011 etax=eta(:,1);
0012 etaxp=eta(:,2);
0013 Hx=gammax.*etax.*etax+2*alphax.*etax.*etaxp+betax.*etaxp.*etaxp;
0014 
0015 %now compute beamenvelopes
0016 %atxdata=atx(ring,dpp,1:length(ring)+1);
0017 %sig=cat(3,atxdata.beam66);
0018 %sigx=sig(1,1,:);
0019 %sigx=squeeze(sigx);
0020 %sigy=sig(3,3,:);
0021 %sigy=squeeze(sigy);
0022 
0023 %momap_h = varargin(1);
0024 
0025 plotdata(1).values=Hx;
0026 plotdata(1).labels={'H_x'};
0027 plotdata(1).axislabel='H [m]';
0028 
0029 plotdata(2).values=1./betax./betaz;
0030 plotdata(2).labels={'1/(betax*betay)'};
0031 plotdata(2).axislabel={'1/m^2'};
0032 
0033 %plotdata(2).values=1./sigx./sigy;
0034 %plotdata(2).labels={'1/(sigmax*sigmay)'};
0035 %plotdata(2).axislabel={'1/m^2'};
0036 
0037 end