{ TITLE "m4-20101101-412-4W10U" }

define lattice; ringtype = 1;

Energy = 3.0;

{making dP smaller speeds up runtime}
dP = 1e-6; CODeps = 1e-14;

{Meth = 4; Nbend = 4; Nquad = 4; Nsext = 4; Ncorr = 1;}
Meth = 4; Nbend = 10; Nquad = 10; Nsext = 10; Ncorr = 1;

pi = 4.0*arctan(1.0);
c0 = 2.99792458e8;


STR0500:   DRIFT, L= 5.0000000000E-01;
STR0450:   DRIFT, L= 4.5000000000E-01;
STRx403:   DRIFT, L= 4.0311000000E-01;
STR0364:   DRIFT, L= 3.6400000000E-01;
STR0302:   DRIFT, L= 3.0200000000E-01;
STRx203:   DRIFT, L= 2.0311000000E-01;
STR0100:   DRIFT, L= 1.0000000000E-01;
STRx093:   DRIFT, L= 9.2680000000E-02;
STR0075:   DRIFT, L= 7.5000000000E-02;
STR0050:   DRIFT, L= 5.0000000000E-02;
STR0025:   DRIFT, L= 2.5000000000E-02;
STR0020:   DRIFT, L= 2.0000000000E-02;
STR0018:   DRIFT, L= 1.8000000000E-02;
STRx013:   DRIFT, L= 1.2500000000E-02;
STR0010:   DRIFT, L= 1.0000000000E-02;
STRx006:   DRIFT, L= 6.0800000000E-03;


{ Girder Markers }
GS: Multipole, N = 1, Method = Meth;
GE: Multipole, N = 1, Method = Meth;


{ Beam Position Monitors }
BPM_D: Drift, L = 5.0e-2/2.0;
BPM_M: Beam Position Monitor;
BPM:   BPM_D, BPM_M, BPM_D;

{ Horizontal and vertical corrector }
CORR_D:  Drift, L = 10.0e-2/4.0;
CORR_H:  Corrector, Horizontal, Method= Meth;
CORR_V:  Corrector, Vertical, Method= Meth;
CORR:    CORR_D, CORR_V, 2*CORR_D, CORR_H, CORR_D;
CORRs:   CORR_D, CORR_H, 3*CORR_D;


D0:    BENDING, L= 3.6189000000E-01, T= 1.085673,
       K=-0.848412, T1= 0.0, T2= 0.0,
       N = Nbend, Method = Meth;

Df1:   BENDING, L= 5.0000000000E-02, T= 0.149940,
       K=-0.848031, T1= 0.0, T2= 0.0,
       N = Nbend, Method = Meth;

Df2:   BENDING, L= 5.0000000000E-02, T= 0.149685,
       K=-0.847298, T1= 0.0, T2= 0.0,
       N = Nbend, Method = Meth;

Df3:   BENDING, L= 5.0000000000E-02, T= 0.107834,
       K=-0.578736, T1= 0.0, T2= 0.0,
       N = Nbend, Method = Meth;

Df4:   BENDING, L= 5.0000000000E-02, T= 0.005351,
       K=+0.006286, T1= 0.0, T2= 0.0,
       N = Nbend, Method = Meth;

Df5:   BENDING, L= 5.0000000000E-02, T= 0.001543,
       K=+0.000118, T1= 0.0, T2= 0.0,
       N = Nbend, Method = Meth;

Dm6:   BENDING, L= 5.0000000000E-02, T= 0.002217,
       K=-0.000505, T1= 0.0, T2= 0.0,
       N = Nbend, Method = Meth;

Dm5:   BENDING, L= 5.0000000000E-02, T= 0.053563,
       K=-0.286488, T1= 0.0, T2= 0.0,
       N = Nbend, Method = Meth;

Dm4:   BENDING, L= 5.0000000000E-02, T= 0.074264,
       K=-0.413761, T1= 0.0, T2= 0.0,
       N = Nbend, Method = Meth;

Dm3:   BENDING, L= 5.0000000000E-02, T= 0.077218,
       K=-0.420099, T1= 0.0, T2= 0.0,
       N = Nbend, Method = Meth;

Dm2:   BENDING, L= 5.0000000000E-02, T= 0.116500,
       K=-0.594270, T1= 0.0, T2= 0.0,
       N = Nbend, Method = Meth;

Dm1:   BENDING, L= 5.0000000000E-02, T= 0.149165,
       K=-0.836556, T1= 0.0, T2= 0.0,
       N = Nbend, Method = Meth;

Dm0:   BENDING, L= 2.0424000000E-01, T= 0.612721,
       K=-0.846950, T1= 0.0, T2= 0.0,
       N = Nbend, Method = Meth;


DIP:   Df5, Df4, Df3, Df2, Df1, D0;
DIPm:  Dm6, Dm5, Dm4, Dm3, Dm2, Dm1, Dm0, Df1, Df2, Df3, Df4, Df5;


QF:     QUADRUPOLE, L= 1.5000000000E-01, K= 4.018366,
        N = Nquad, Method = Meth;

QFm:    QUADRUPOLE, L= 1.5000000000E-01, K= 3.769247,
        N = Nquad, Method = Meth;

QFend:  QUADRUPOLE, L= 2.5000000000E-01, K= 3.521817,
        N = Nquad, Method = Meth;

QFendW: QUADRUPOLE, L= 2.5000000000E-01, K= 3.649134,
        N = Nquad, Method = Meth;

QFendU: QUADRUPOLE, L= 2.5000000000E-01, K= 3.642797,
        N = Nquad, Method = Meth;

QDend:  QUADRUPOLE, L= 2.5000000000E-01, K=-2.176206,
        N = Nquad, Method = Meth;

QDendW: QUADRUPOLE, L= 2.5000000000E-01, K=-2.513140,
        N = Nquad, Method = Meth;

QDendU: QUADRUPOLE, L= 2.5000000000E-01, K=-2.495706,
        N = Nquad, Method = Meth;


SD:     SEXTUPOLE, L= 1.0000000000E-01, K=-112.549066,
        N = Nsext, Method = Meth;

SDend:  SEXTUPOLE, L= 1.0000000000E-01, K=-162.000000,
        N = Nsext, Method = Meth;

SFm:    SEXTUPOLE, L= 1.0000000000E-01, K=168.000000,
        N = Nsext, Method = Meth;

SFo:    SEXTUPOLE, L= 1.0000000000E-01, K=170.000000,
        N = Nsext, Method = Meth;

SFi:    SEXTUPOLE, L= 1.0000000000E-01, K=209.074683,
        N = Nsext, Method = Meth;


SQFm:   QFm, STR0075, SFm, STRx013, BPM, STRx013, QFm, STR0100, CORR;
SQFo:   QF,  STR0075, SFo, STRx013, BPM, STRx013, QF,  STR0100, CORR;
SQFi:   QF,  STR0075, SFi, STRx013, BPM, STRx013, QF,  STR0100, CORR;


OXX:   MULTIPOLE, L = 1.0e-1, HOM=(4,-58.3029/1.0e-1,0),
       N = Nsext, Method = Meth;

OXY:   MULTIPOLE, L = 1.0e-1, HOM=(4,91.8432/1.0e-1,0),
       N = Nsext, Method = Meth;

OYY:   MULTIPOLE, L = 1.0e-1, HOM=(4,-84.1923/1.0e-1,0),
       N = Nsext, Method = Meth;


h_rf = 176; C =  528.0;
CAV: Cavity, Frequency = c0/C*h_rf, Voltage = 1.8e6, Harnum = h_rf;


lambda_PMDW = 80e-3;
n_PMDW = 25;
L_PMDW = n_PMDW*lambda_PMDW;
B_PMDW=2.22;
Brho = 10.007;
PMDWw: Wiggler, L = L_PMDW, lambda = lambda_PMDW, kxV = 0.0, 
       BoBrhoV = B_PMDW/Brho, N = n_PMDW*4, Method = 2;
{PMDWw: insertion, scaling = 1.0, N = 50, Method = 2,
       file2 = "/home/simon/projects/in/lattice/pmdwkick.dat";}
PMDWd: Drift, L = -L_PMDW/2.0;
PMDW:  PMDWd, PMDWw, PMDWd;


lambda_IVU = 19e-3;
n_IVU = 158;
L_IVU = n_IVU*lambda_IVU;
B_IVU=1.27;
IVUu: Wiggler, L = L_IVU, lambda = lambda_IVU, kxV = 0.0, 
       BoBrhoV = B_IVU/Brho, N = n_IVU*4, Method = 2;
IVUd: Drift, L = -L_IVU/2.0;
IVU:  IVUd, IVUu, IVUd;


LS:    4*STR0500, STR0364;
wLS:   4*STR0500, STR0364;
uLS:   4*STR0500, STR0364;
SS:    STR0450, STR0500, STR0302;
DIPuc: SD, STR0010, DIP, INV(DIP), STR0010, SD;
UC1:   GS, SQFm, STRx203, DIPuc, GE, STRx403;
UC2:   GS, SQFo, STRx203, DIPuc, GE, STRx403;
UC3:   GS, SQFi, STRx203, DIPuc, STRx203, INV(SQFi), GE;
UC4:   STRx403, GS, DIPuc, STRx203, INV(SQFo), GE;
UC5:   STRx403, GS, DIPuc, STRx203, INV(SQFm), GE;
MC1:   BPM, STR0018, CORR, STR0018, OXX, STR0025, QFend, STR0025, OXY, STR0100, QDend,
       STRx006, DIPm, OYY, STRx093, CORRs, BPM, STR0020, SDend, STR0050;
wMC1:  BPM, STR0018, CORR, STR0018, OXX, STR0025, QFendW, STR0025, OXY, STR0100, QDendW,
       STRx006, DIPm, OYY, STRx093, CORRs, BPM, STR0020, SDend, STR0050;
uMC1:  BPM, STR0018, CORR, STR0018, OXX, STR0025, QFendU, STR0025, OXY, STR0100, QDendU,
       STRx006, DIPm, OYY, STRx093, CORRs, BPM, STR0020, SDend, STR0050;
MC2:   BPM, STR0018, CORR, STR0018, OXX, STR0025, QFend, STR0025, OXY, STR0100, QDend,
       STRx006, DIPm, OYY, STRx093, CORR, BPM, STR0020, SDend, STR0050;
wMC2:  BPM, STR0018, CORR, STR0018, OXX, STR0025, QFendW, STR0025, OXY, STR0100, QDendW,
       STRx006, DIPm, OYY, STRx093, CORR, BPM, STR0020, SDend, STR0050;
uMC2:  BPM, STR0018, CORR, STR0018, OXX, STR0025, QFendU, STR0025, OXY, STR0100, QDendU,
       STRx006, DIPm, OYY, STRx093, CORR, BPM, STR0020, SDend, STR0050;
LINE:  LS, GS, MC1, GE, SS, UC1, UC2, UC3, UC4, UC5, INV(SS), GS, INV(MC2), GE, INV(LS);
wLINE1:LS, GS, MC1, GE, SS, UC1, UC2, UC3, UC4, UC5, INV(SS), GS, INV(wMC2), GE, INV(wLS);
wLINE2:wLS, GS, wMC1, GE, SS, UC1, UC2, UC3, UC4, UC5, INV(SS), GS, INV(MC2), GE, INV(LS);
uLINE1:LS, GS, MC1, GE, SS, UC1, UC2, UC3, UC4, UC5, INV(SS), GS, INV(uMC2), GE, INV(uLS);
uLINE2:uLS, GS, uMC1, GE, SS, UC1, UC2, UC3, UC4, UC5, INV(SS), GS, INV(MC2), GE, INV(LS);
wuLINE:wLS, GS, wMC1, GE, SS, UC1, UC2, UC3, UC4, UC5, INV(SS), GS, INV(uMC2), GE, INV(uLS);
uwLINE:uLS, GS, uMC1, GE, SS, UC1, UC2, UC3, UC4, UC5, INV(SS), GS, INV(wMC2), GE, INV(wLS);
uuLINE:uLS, GS, uMC1, GE, SS, UC1, UC2, UC3, UC4, UC5, INV(SS), GS, INV(uMC2), GE, INV(uLS);
wLINES:wLINE1, PMDW, wLINE2;
uLINES:uLINE1, IVU, uLINE2;


{RING: wLINES, CAV;}
{RING: uLINES, CAV;}

RING:  wLINE1, PMDW, wuLINE, IVU, uwLINE,
               PMDW, wuLINE, IVU, uwLINE,
               PMDW, wuLINE, IVU, uwLINE,
               PMDW, wuLINE, IVU, uuLINE,
                             IVU, uuLINE,
                             IVU, uuLINE,
                             IVU, uuLINE,
                             IVU, uuLINE,
                             IVU, uuLINE,
                             IVU, uLINE2,
       5*LINE, CAV;

CELL: RING, SYMMETRY = 1;

END;