{ 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;
Brho = 1e9*Energy/c0;


STR0250:   DRIFT, L= 2.5000000000E-01;
STR0190:   DRIFT, L= 1.9000000000E-01;
STR0140:   DRIFT, L= 1.4000000000E-01;
STR0090:   DRIFT, L= 9.0000000000E-02;
STR0080:   DRIFT, L= 8.0000000000E-02;
STR0025:   DRIFT, L= 2.5000000000E-02;
STR0015:   DRIFT, L= 1.5000000000E-02;


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


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

{ Horizontal and vertical corrector }
CORRH: Corrector, Horizontal, Method= Meth;
CORRV: Corrector, Vertical, Method= Meth;


{T3 and OPA quad/sext strengths are normalized}
{m->m/2 for MAD8->T3}
{m->m/2 for MAD8->OPA}


QFmaster=1.0;QDmaster=1.0;NUXmaster=1.0;NUYmaster=1.0;NUXmaster2=1.0;NUYmaster2=1.0;
{to get tunes 11.22/3.13}
QFmaster=1.000002803;QDmaster=1.000021214;
{from 501 lattice -> to get chroma 2.0/2.0}
NUXmaster=1.038248585;NUYmaster=1.030848314;
{to get chroma 1.0/1.0 in 503 lattice}
NUXmaster2=1.006159465;NUYmaster2=1.001797131;


DIP: BENDING, L= 1.00000E+00, T= 15.0,
     K=-1.347963*QDmaster, T1= 0.0, T2= 0.0,
     N = Nbend, Method = Meth;


SQFo: BENDING, L= 2.00E-01, K= 5.736667*QFmaster,
      HOM = (3, 36.676825, 0.0),
      N = Nquad, Method = Meth;

SQFi: BENDING, L= 2.00E-01, K= 4.998490*QFmaster,
      HOM = (3, 28.040076*NUXmaster, 0.0),
      N = Nquad, Method = Meth;


SDi: SEXTUPOLE, L= 1.0000000000E-01, K=-68.175522*NUYmaster*NUYmaster2,
     N = Nsext, Method = Meth;

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


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

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


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


lambda_SCW = 61e-3;
n_SCW = 25;
L_SCW = n_SCW*lambda_SCW;
B_SCW = 3.0;
SCWw: Wiggler, L = L_SCW, lambda = lambda_SCW, kxV = 0.0, 
       BoBrhoV = B_SCW/Brho, N = n_SCW*4, Method = 2;
SCWd: Drift, L = -L_SCW/2.0;
SCW:  SCWd, SCWw, SCWd;


SCi1: STR0015, SCi, Sci, STR0015;
SCi2: HVPU, STR0015, SCi, CORRV, CORRH, SCi, STR0015;
SCo1: STR0015, SCo, CORRH, CORRV, SCo, STR0015;

STR:  7*STR0250;

hACHR1: STR, GS, HVPU, SQFo, SCo1, STR0080, SDo, STR0090, DIP,
        STR0190, SDi, STR0140, SCi1, SQFi;
hACHR2: STR, GE, HVPU, SQFo, SCo1, STR0080, SDo, STR0090, DIP,
        STR0190, SDi, STR0090, SCi2, SQFi;

ACHR:  hACHR1, INV(hACHR2);

{RING: ACHR, CAV;}
RING: 12*ACHR, CAV;

CELL: RING, SYMMETRY = 1;

END;