! BMAD input file

lattice := FODO

beam, energy = 1

beginning[beta_x] =11.27161797
beginning[beta_y] = 4.31731041
beginning[alpha_x]=-1.61152116
beginning[alpha_y]= 0.67030841
beginning[eta_x]  = 2.65296216
beginning[etap_x] = 0.38439728

PARAMETER[LATTICE_TYPE] = LINEAR_LATTICE

ip:        marker

d1 : drift, l=0.25
d2 : drift, l=1.0

bend: sbend, l=2.5, angle = pi*2/40
db: drift, l=2.5 ! drift for the missing bend scheme

Quad_x: quadrupole, l = 0.50
Quad_y: quadrupole, l = 0.50
Quad_x[k1]= 5.7264354E-01
Quad_y[k1]=-5.9079362E-01

Quad_x0: quadrupole, l = 0.50
Quad_y0: quadrupole, l = 0.50
Quad_x1: quadrupole, l = 0.50
Quad_y1: quadrupole, l = 0.50
Quad_x0[k1]= 5.7264354E-01
Quad_y0[k1]=-3.1166922E-01
Quad_x1[k1]= 4.0289801E-01
Quad_y1[k1]=-5.4692642E-01

Quad_1: quadrupole, l = 0.50
Quad_2: quadrupole, l = 0.50
quad_3: quadrupole, l = 0.50
Quad_4: quadrupole, l = 0.50
Quad_1[k1]= 2.1406135E-01
Quad_2[k1]=-5.3874299E-01
Quad_3[k1]= 1.5367782E+00
Quad_4[k1]=-2.3965273E+00

det_y: MARKER, ALIAS = "det_y"
det_x: MARKER, ALIAS = "det_x"

fodo0: line = (Quad_x0, d1, db, d1, Quad_y0, d1, bend, d1)

fodo1: line = (Quad_x1, d1, db, d1, det_y, Quad_y1, d1, bend, d1, det_x)

insert: line = (Quad_1, d2, Quad_2, d2, Quad_3, d2, Quad_4, d2, ip)

module: line = (fodo0, fodo1, insert)

use, module