{ TITLE "m5-ring" }

define lattice; ringtype = 1;

Energy = 1.5;

{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;


STR0250:   DRIFT, L= 2.5000000000E-01;
STRa250:   DRIFT, L= 2.5000000000E-01;
STR0248:   DRIFT, L= 2.4800000000E-01;
STR0200:   DRIFT, L= 2.0000000000E-01;
STR0150:   DRIFT, L= 1.5000000000E-01;
STR0130:   DRIFT, L= 1.3000000000E-01;
STR0100:   DRIFT, L= 1.0000000000E-01;
STR0098:   DRIFT, L= 9.8000000000E-02;
STR0095:   DRIFT, L= 9.5000000000E-02;
STR0080:   DRIFT, L= 8.0000000000E-02;
STR0050:   DRIFT, L= 5.0000000000E-02;
STR0025:   DRIFT, L= 2.5000000000E-02;
STR0015:   DRIFT, L= 1.5000000000E-02;
STRa015:   DRIFT, L= 1.5000000000E-02;
STRa010:   DRIFT, L= 1.0000000000E-02;
STR0002:   DRIFT, L= 2.0000000000E-03;
STRa002:   DRIFT, L= 2.0000000000E-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_Da: Drift, L = 5.0e-2/2.0;
BPM:    BPM_D, BPM_M, BPM_D;
BPMa:   BPM_Da, BPM_M, BPM_Da;

{ Horizontal and vertical corrector }
CORR_H:  Corrector, Horizontal, Method= Meth;
CORR_V:  Corrector, Vertical, Method= Meth;
CORR_Ha: Corrector, Horizontal, Method= Meth;
CORR_Va: Corrector, Vertical, Method= Meth;


bm14: BENDING, L= 2.00000E-02, T= 0.000684,
      K=0.009083, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
bm13: BENDING, L= 2.50000E-02, T= 0.001035,
      K=0.010144, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
bm12: BENDING, L= 2.50000E-02, T= 0.014773,
      K=-0.013719, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
bm11: BENDING, L= 2.50000E-02, T= 0.125664,
      K=-0.424484, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
bm10: BENDING, L= 2.50000E-02, T= 0.248086,
      K=-0.894212, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
bm9:  BENDING, L= 2.50000E-02, T= 0.355434,
      K=-1.370585, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
bm8:  BENDING, L= 2.50000E-02, T= 0.375289,
      K=-1.350285, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
bm7:  BENDING, L= 2.50000E-02, T= 0.375477,
      K=-1.348728, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
bm6:  BENDING, L= 2.50000E-02, T= 0.375486,
      K=-1.349533, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
bm5:  BENDING, L= 2.50000E-02, T= 0.375472,
      K=-1.349813, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
bm4:  BENDING, L= 2.50000E-02, T= 0.375446,
      K=-1.350099, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
bm3:  BENDING, L= 2.50000E-02, T= 0.375406,
      K=-1.350067, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
bm2:  BENDING, L= 1.00000E-01, T= 1.501131,
      K=-1.349876, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
bm1:  BENDING, L= 1.00000E-01, T= 1.500354,
      K=-1.349340, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
bm0:  BENDING, L= 1.00000E-01, T= 1.499919,
      K=-1.349033, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
b0:   BENDING, L= 1.00000E-01, T= 1.499942,
      K=-1.349122, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
b1:   BENDING, L= 1.00000E-01, T= 1.500434,
      K=-1.349474, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
b2:   BENDING, L= 1.00000E-01, T= 1.501257,
      K=-1.349976, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
b3:   BENDING, L= 2.50000E-02, T= 0.375446,
      K=-1.350157, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
b4:   BENDING, L= 2.50000E-02, T= 0.375486,
      K=-1.350161, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
b5:   BENDING, L= 2.50000E-02, T= 0.375512,
      K=-1.349995, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
b6:   BENDING, L= 2.50000E-02, T= 0.375526,
      K=-1.349559, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
b7:   BENDING, L= 2.50000E-02, T= 0.375520,
      K=-1.348866, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
b8:   BENDING, L= 2.50000E-02, T= 0.375331,
      K=-1.350210, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
b9:   BENDING, L= 2.50000E-02, T= 0.355459,
      K=-1.369168, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
b10:  BENDING, L= 2.50000E-02, T= 0.248126,
      K=-0.894396, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
b11:  BENDING, L= 2.50000E-02, T= 0.125713,
      K=-0.427264, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
b12:  BENDING, L= 2.50000E-02, T= 0.014833,
      K=-0.023614, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
b13:  BENDING, L= 2.50000E-02, T= 0.001044,
      K=0.010325, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;
b14:  BENDING, L= 2.00000E-02, T= 0.000716,
      K=0.009315, T1= 0.0, T2= 0.0,
      N = Nbend, Method = Meth;

DIP:  bm14, bm13, bm12, bm11, bm10, bm9, bm8, bm7, bm6,
      bm5, bm4, bm3, bm2, bm1, bm0, b0, b1, b2, b3, b4,
      b5, b6, b7, b8, b9, b10, b11, b12, b13, b14;


SQFi_en: QUADRUPOLE, L= 5.00E-02, K= 5.004739,
         HOM = (3, 26.403945, 0.0),
         N = Nquad, Method = Meth;
SQFi_m:  QUADRUPOLE, L= 1.50E-01, K= 4.999557,
         HOM = (3, 29.112205, 0.0),
         N = Nquad, Method = Meth;
SQFi:    SQFi_en, SQFi_m;


SQFo_en: QUADRUPOLE, L= 5.00E-02, K= 5.745075,
         HOM = (3, 41.824688, 0.0),
         N = Nquad, Method = Meth;
SQFo_m:  QUADRUPOLE, L= 1.00E-01, K= 5.734628,
         HOM = (3, 45.280000, 0.0),
         N = Nquad, Method = Meth;
SQFo_ex: QUADRUPOLE, L= 5.00E-02, K= 5.744196,
         HOM = (3, 41.757489, 0.0),
         N = Nquad, Method = Meth;
SQFo:    SQFo_en, SQFo_m, SQFo_ex;


SDi: SEXTUPOLE, L= 1.0000000000E-01, K=-68.730946,
     N = Nsext, Method = Meth;

SDo: SEXTUPOLE, L= 1.0000000000E-01, K=-98.320000,
     N = Nsext, Method = Meth;

SCi: SEXTUPOLE, L= 1.0000000000E-02, K=-21.389322,
     N = Nsext, Method = Meth;

SCo: SEXTUPOLE, L= 1.0000000000E-02, K=-22.000000,
     N = Nsext, Method = Meth;


h_rf = 32; C =  96.0;
CAV: Cavity, Frequency = c0/C*h_rf, Voltage = 5.6e5, Harnum = h_rf;


IP:   MARKER;
PSMc: SEXTUPOLE, L= 0.0, K=-0*46.780,
      N = Nsext, Method = Meth;
KIc:  MARKER;
PVc:  MARKER;

PSM:  STR0200, PSMc, STR0200;
KI:   STR0100, KIc, STR0100;
PV:   STR0100, PVc, STR0100;
SEPT: 5*STRa250, STRa002;

SCi1: STRa015, SCi, SCi, STRa015;
SCi2: BPMa, STRa015, SCi, CORR_Va, CORR_Ha, SCi, STRa015;
SCo1: STR0015, SCo, CORR_H, CORR_V, SCo, STR0015;

hSTR:    7*STR0250;
hSTR1a:  SEPT, IP, STR0248, STR0250;
hSTR2b:  4*STR0250, STR0002, PSM, STR0050, KI, STR0098;
hSTR12b: 4*STR0250, STR0150, PV, 4*STR0100;

hBlock1: GS, BPM, SQFo, SCo1, STR0050, SDo, STR0025, DIP,
         STR0095, SDi, STR0130, SCi1, STRa010, SQFi;
hBlock2: GE, BPM, SQFo, SCo1, STR0050, SDo, STR0025, DIP,
         STR0095, SDi, STR0080, SCi2, STRa010, SQFi;

ACHR:   hSTR, hBlock1, INV(hBlock2), INV(hSTR);
ACHR1:  hSTR1a, hBlock1, INV(hBlock2), INV(hSTR);
ACHR2:  hSTR, hBlock1, INV(hBlock2), INV(hSTR2b);
ACHR12: hSTR, hBlock1, INV(hBlock2), INV(hSTR12b);


{for injection studies}
{starts at IP}
{RING: IP, STR0248, STR0250, hBlock1, INV(hBlock2), INV(hSTR), ACHR2, 9*ACHR, CAV, ACHR12, SEPT;}
{starts at KI}
{RING: KIc, STR0100, STR0050, PSM, STR0002, 4*STR0250, 9*ACHR, CAV, ACHR12, ACHR1, hSTR, hBlock1, INV(hBlock2),
      STR0098, STR0100;}
{starts at PSM}
{RING: PSMc, STR0200, STR0002, 4*STR0250, 9*ACHR, CAV, ACHR12, ACHR1, hSTR, hBlock1, INV(hBlock2), STR0098,
      KI, STR0050, STR0200;}


{RING: ACHR, CAV;}
{RING: 12*ACHR, CAV;}
RING: ACHR1, ACHR2, 9*ACHR, CAV, ACHR12;

CELL: RING, SYMMETRY = 1;

END;