source: MML/trunk/machine/SOLEIL/LT1/Lattices/maille_LT1.m @ 4

function maille_LT1
% Lattice definition file
% Lattice for LT1
%
% Laurent S. Nadolski, SOLEIL, 03/04

global FAMLIST THERING GLOBVAL
GLOBVAL.E0 = 100e6;
GLOBVAL.LatticeFile = mfilename;
FAMLIST = cell(0);
disp(['** Loading lattice in ' mfilename])

%% markers           
%% Ecrans
ECR1 =  marker('BPM', 'IdentityPass');
ECR2 =  marker('BPM', 'IdentityPass');
ECR3 =  marker('BPM', 'IdentityPass');
ECR4 =  marker('BPM', 'IdentityPass');
%% Collimateurs 
COLL1 =  marker('COLL1', 'IdentityPass');
COLL2 =  marker('COLL2', 'IdentityPass');
%% perturbateurs
KICKQD =  marker('KICKQD', 'IdentityPass');
PERTM =  marker('PERTM', 'IdentityPass');
%% divers
FIN =  marker('FIN', 'IdentityPass');
DEBUT =  marker('DEBUT', 'IdentityPass');

% DRIFT SPACES

SDL1  = drift('SDL1', 1.119000, 'DriftPass');
SDL1A = drift('SDL1A',  0.550000, 'DriftPass');
SDL2  = drift('SDL2',  0.180000, 'DriftPass');
SDL2B = drift('SDL2B',  0.674150, 'DriftPass');
SDL3  = drift('SDL3',  0.180000, 'DriftPass');
SDL3B = drift('SDL3B',  2.220000, 'DriftPass');
SDL4  = drift('SDL4',  0.180000, 'DriftPass');
SDL4B = drift('SDL4B', 0.770000, 'DriftPass');
SDL5  = drift('SDL5', 0.180000, 'DriftPass');
SDL5B = drift('SDL5B', 0.570000, 'DriftPass');
SDL5H1= drift('SDL5H1',  0.180000, 'DriftPass');
SDL5H2= drift('SDL5H2',  0.570000, 'DriftPass');
SDL5C = drift('SDL5C',  0.300000, 'DriftPass');
SDL6  = drift('SDL6', 1.983680, 'DriftPass');
SDL6B = drift('SDL6B',  0.380000, 'DriftPass');
SDL6T = drift('SDL6T',  0.350000, 'DriftPass');
SDL7  = drift('SDL7',  0.130000, 'DriftPass');
SDL7B = drift('SDL7B',  0.971910, 'DriftPass');
SDL8  = drift('SDL8',  0.180000, 'DriftPass');
SDL8B = drift('SDL8B', 0.300000, 'DriftPass');
SDL9  = drift('SDL9',  0.300000, 'DriftPass');
SDL9B = drift('SDL9B',  0.180000, 'DriftPass');
SDL0  = drift('SDL0',  0.243340, 'DriftPass');
SDL0B = drift('SDL0B',  0.500000, 'DriftPass');
SDL0B1= drift('SDL0B1', 0.500000, 'DriftPass');
SDL0C = drift('SDL0C',  2.393090, 'DriftPass');
SBOO1 = drift('SBOO1',  1.294100, 'DriftPass');
SBOO2 = drift('SBOO2',  0.764434, 'DriftPass');
SBOO3 = drift('SBOO3', 0.150000, 'DriftPass');
SBOO4 = drift('SBOO4',  2.000000, 'DriftPass');
PERTEBR = drift('PERTEBR',  0.130000, 'DriftPass');
PERTELM = drift('PERTELM',  0.30000, 'DriftPass');
PERTSLM = drift('PERTSLM',  0.30000, 'DriftPass');
PERTSBR = drift('PERTSBR',  0.130000, 'DriftPass');

% QUADRUPOLES 
L= 0.15/2;
QP1  =  quadrupole('QP1', L,  4.17582 , 'QuadLinearPass');
QP2  =  quadrupole('QP2', L, -2.78267 , 'QuadLinearPass');
QP3  =  quadrupole('QP3', L, -2.61001 , 'QuadLinearPass');
QP4  =  quadrupole('QP4', L,  4.29826   , 'QuadLinearPass');
QP5  =  quadrupole('QP5', L,  4.94250   , 'QuadLinearPass');
QP6  =  quadrupole('QP6', L, -4.94777  , 'QuadLinearPass');
QP7  =  quadrupole('QP7', L,  4.92238  , 'QuadLinearPass');
QDBOO=  quadrupole('QDBOO', 0.2,  -0.891375 , 'QuadLinearPass');

% Q1 = [QP1 QP1];
% Q2 = [QP2 QP2];
% Q3 = [QP3 QP3];
% Q4 = [QP4 QP4];
% Q5 = [QP5 QP5];
% Q6 = [QP6 QP6];
% Q7 = [QP7 QP7];
%% EN attente d'éléments coupés en deux
L= 0.15; 
Q1  =  quadrupole('QP', L,  4.17582 , 'QuadLinearPass');
Q2  =  quadrupole('QP', L, -2.78267 , 'QuadLinearPass');
Q3  =  quadrupole('QP', L, -2.61001 , 'QuadLinearPass');
Q4  =  quadrupole('QP', L,  4.29826   , 'QuadLinearPass');
Q5  =  quadrupole('QP', L,  4.94250   , 'QuadLinearPass');
Q6  =  quadrupole('QP', L, -4.94777  , 'QuadLinearPass');
Q7  =  quadrupole('QP', L,  4.92238  , 'QuadLinearPass');


%% Vertical Correctors 
CV  =  corrector('COR',0.0,[0 0],'CorrectorPass');
CV1  =  corrector('CV1',0.0,[0 0],'CorrectorPass');
CV2  =  corrector('CV2',0.0,[0 0],'CorrectorPass');
CV3  =  corrector('CV3',0.0,[0 0],'CorrectorPass');

%% Horizontal Correctors 
CH  =  corrector('CH',0.0,[0 0],'CorrectorPass');
CH1  =  corrector('CH1',0.0,[0 0],'CorrectorPass');
CH2  =  corrector('CH2',0.0,[0 0],'CorrectorPass');
CH3  =  corrector('CH3',0.0,[0 0],'CorrectorPass');
CH4  =  corrector('CH4',0.0,[0 0],'CorrectorPass');

% DIPOLES
% {** 1.3815 factor to fit with BETA ??? strange **}
theta = -15*pi/180;
fullgap = 0.105*0.724*2/6*1.3815;
BEND  =  rbend2('BEND', 0.30, theta, theta/2, theta/2, 0.0, ...
                fullgap,'BendLinearFringeTiltPass');

% {** septa booster **}
% beta_gap=0.03/6;
% tracy_gap=beta_gap*0.724*2;
theta = -7.5*pi/180;
fullgap= 0.03/6*0.724*2; %*1.3815;
L = 0.30;
SEP     = rbend2('SEP', L,  theta, theta/2, theta/2, 0.0, fullgap, 'BendLinearFringeTiltPass');
SEPup   = rbend2('SEPup', L/2, theta/2, theta/2, 0.0, 0.0,  fullgap, 'BendLinearFringeTiltPass');
SEPdown = rbend2('SEPdown', L/2, theta/2, 0.0, theta/2, 0.0,  fullgap, 'BendLinearFringeTiltPass');

% Begin Lattice
%%  Superperiods  
LINE2 = [ ... 
    SDL1       SDL1A     QP1       QP1       SDL2       CH1        SDL2B     ...
    QP2        QP2       SDL3      CV1       SDL3B      QP3        QP3    ...
    SDL4       SDL4B     QP4       QP4       SDL5H1     CH2        SDL5H2     ...
    ECR1       SDL5C     BEND      SDL6      COLL1      SDL6B      COLL2      ...
    SDL6T      QP5       QP5       SDL7      CH3        SDL7B      CV2        ...
    SDL8       QP6       QP6       SDL8B     BEND       SDL9       ECR2    ...
    SDL0       QP7       QP7       SDL9B     SDL0C      CV3        SDL0B1     ...
    ECR3       SDL0B     SEPup     CH4       SEPdown    SBOO1      QDBOO      ...
    KICKQD     QDBOO  ...    
    SBOO2      SBOO3     PERTEBR   PERTELM   PERTM      PERTSLM    PERTSBR    ...
    SBOO3      ECR4      SBOO4
];

LINE = [ ... 
    SDL1       SDL1A     Q1       SDL2       CH        SDL2B     ...
    Q2         SDL3      CV       SDL3B      Q3        ...
    SDL4       SDL4B     Q4       SDL5H1     CH        SDL5H2     ...
    ECR1       SDL5C     BEND      SDL6      COLL1      SDL6B      COLL2      ...
    SDL6T      Q5       SDL7      CH        SDL7B      CV        ...
    SDL8       Q6       SDL8B     BEND       SDL9       ECR2    ...
    SDL0       Q7       SDL9B     SDL0C      CV        SDL0B1     ...
    ECR3       SDL0B     SEPup     CH       SEPdown    SBOO1      QDBOO      ...
    KICKQD     QDBOO  ...    
    SBOO2      SBOO3     PERTEBR   PERTELM   PERTM      PERTSLM    PERTSBR    ...
    SBOO3      ECR4      SBOO4
];

ELIST = [DEBUT, LINE, FIN];

THERING=cell(size(ELIST));
for i=1:length(THERING)
    THERING{i} = FAMLIST{ELIST(i)}.ElemData;
    FAMLIST{ELIST(i)}.NumKids=FAMLIST{ELIST(i)}.NumKids+1;
    FAMLIST{ELIST(i)}.KidsList = [FAMLIST{ELIST(i)}.KidsList i];
end

%% compute total length

L0_tot=0;
for i=1:length(THERING)
    L0_tot=L0_tot+THERING{i}.Length;
end

disp('** Done **');