{ TITLE "m4-20101101-411-4W" }

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.862608, T1= 0.0, T2= 0.0,
       N = Nbend, Method = Meth;

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

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

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

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

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

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

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

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

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

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

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

Dm0:   BENDING, L= 2.0424000000E-01, T= 0.612721,
       K=-0.861122, 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.031557,
        N = Nquad, Method = Meth;

QFm:    QUADRUPOLE, L= 1.5000000000E-01, K= 3.776598,
        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.647797,
        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.513862,
        N = Nquad, Method = Meth;


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

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

SFm:    SEXTUPOLE, L= 1.0000000000E-01, K=156.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=216.891574,
        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,-135.5801/1.0e-1,0),
       N = Nsext, Method = Meth;

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

OYY:   MULTIPOLE, L = 1.0e-1, HOM=(4,-85.5096/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;


LS:    4*STR0500, STR0364;
wLS:   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;
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;
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);
wLINES:wLINE1, PMDW, wLINE2;
wRINGper:wLINES, 3*LINE;


{RING: wLINES, CAV;}
{RING: wRINGper, CAV;}
RING: 4*wRINGper, CAV;

CELL: RING, SYMMETRY = 1;

END;