source: MML/trunk/at/simulator/element/user/sextcorr.m @ 4

function z=sextupole_corr(fname,L,S,KickAngle,method)

%SEXTUPOLE('FAMILYNAME',Length [m],S,Kickangle,'METHOD')
%       creates a new family in the FAMLIST - a structure with fields%          FamName    
%       FamName                 family name
%       Length                  length[m]
%       S                                       S-strngth of the sextupole
%       NumIntSteps             Number of integration steps
%       MaxOrder
%       R1                                      6 x 6 rotation matrix at the entrance
%       R2                      6 x 6 rotation matrix at the entrance
%       T1                                      6 x 1 translation at entrance 
%       T2                                      6 x 1 translation at exit4
%       ElemData.PolynomA= [0 0 0 0];    
%       ElemData.PolynomB= [0 0 S 0]; 
%   KickAngle      same as corrector however in this case it is a thick
%                  element kick, ie it adds PolynomB(1) = sin(angle)/L.
%       PassMethod     name of the function to use for tracking
% returns assigned address in the FAMLIST that is uniquely identifies
% the family


ElemData.FamName = fname;  % add check for identical family names
ElemData.Length = L;
ElemData.MaxOrder = 3;
ElemData.NumIntSteps = 10;
ElemData.R1 = diag(ones(6,1));
ElemData.R2 = diag(ones(6,1));
ElemData.T1 = zeros(1,6);
ElemData.T2 = zeros(1,6);
ElemData.KickAngle = KickAngle;
ElemData.PolynomA= [0 0 0 0];    
ElemData.PolynomB= [0 0 S 0]; 
ElemData.PassMethod=method;

global FAMLIST
z = length(FAMLIST)+1; % number of declare families including this one
FAMLIST{z}.FamName = fname;
FAMLIST{z}.NumKids = 0;
FAMLIST{z}.KidsList= [];
FAMLIST{z}.ElemData= ElemData;