source: MML/trunk/at/simulator/element/sbend.m

function z=sbend(fname,L,A,A1,A2,K,method)
%BEND('FAMILYNAME',  Length[m], BendingAngle[rad], EntranceAngle[rad],
%       ExitAngle[rad], K, 'METHOD')
%       creates a new family in the FAMLIST - a structure with fields
%               FamName         family name
%               Length          length of the arc for an on-energy particle [m]
%               BendingAngle    total bending angle [rad]
%               EntranceAngle   [rad] (0 - for sector bends)
%               ExitAngle               [rad] (0 - for sector bends)
%               ByError                 error in the dipole field relative to the design value 
%               K                               quadrupole K-value for combined funtion bends
%               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.BendingAngle   = A;
ElemData.EntranceAngle  = A1;
ElemData.ExitAngle      = A2;
ElemData.ByError        = 0;
ElemData.K                  = K;

ElemData.R1 = diag(ones(6,1));
ElemData.R2 = diag(ones(6,1));
ElemData.T1 = zeros(1,6);
ElemData.T2 = zeros(1,6);

ElemData.PolynomA                       = [0 0 0 0];     
ElemData.PolynomB                       = [0 K 0 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;