no_digested

particle_start[x]  = 1e-4
particle_start[px] = 2e-4
particle_start[y]  = 3e-4
particle_start[py] = 4e-4
particle_start[z]  = 5e-3
particle_start[pz] = 0.3
particle_start[spin_x] = 0.3
particle_start[spin_y] = 0.4
particle_start[spin_z] = 0.5

beginning[beta_a] =  10
beginning[beta_b] =  10

e_tot = 1e10
parameter[geometry] = open
beginning[e_tot] = e_tot
parameter[absolute_time_tracking] = T

ab_multipole1:     ab_multipole, a2 = 20, b2 = 10, x_offset = 0.02, tilt = 0.1
beambeam1:         beambeam, sig_x = 3e-3, sig_y = 3e-4, x_offset = 0.05, y_offset = 0.05
drift1:            drift, l = 4.5, num_steps = 10, x_offset = 0.05, y_pitch = 0.001, tilt = 0.1
ecollimator1:      ecollimator, l = 4.5, x_limit = 0.09/2, y_limit = 0.05/2, num_steps = 10, x_offset = 0.02, y_pitch = 0.001, tilt = 0.1
elseparator1:      elseparator, l = 4.5, hkick = 0.005, gap = 0.11, num_steps = 10, x_offset = 0.05, y_pitch = 0.001, tilt = 0.1
elseparator2:      elseparator, l = 4.5, vkick = 0.005, gap = 0.11, num_steps = 10, x_offset = 0.05, y_pitch = 0.001, tilt = 0.1
hkicker1:          hkicker, l = 4.5, kick = 0.005, num_steps = 10, x_offset = 0.05, y_pitch = 0.001, tilt = 0.1
instrument1:       instrument, l = 4.5, num_steps = 10
kicker1:           kicker, l = 4.5, hkick = 0.003, vkick = 0.003, num_steps = 10, x_offset = 0.05, y_pitch = 0.001, tilt = 0.1
monitor1:          monitor, l = 4.5, num_steps = 10
multipole1:        multipole, k1l = 1, t1, k3l = 4.5, t3 = 0.31*pi, x_offset = 0.05
octupole1:         octupole, l = 4.5, k3 = 1e3, tilt, num_steps = 10, x_offset = 0.01, y_pitch = 0.001
patch1:            patch, x_offset = 1e-5, y_offset = 2e-5, tilt = 0.3, z_offset = 1, x_pitch = 3e-5, y_pitch = 4e-5
quadrupole1:       quadrupole, l = 0.6, k1 = 1, tilt, num_steps = 10, x_offset = 0.05, y_pitch = 0.001
quadrupole2:       quadrupole, l = 0.6, k1 = -1, tilt, num_steps = 10, x_offset = 0.05, y_pitch = 0.001
quadrupole3:       quadrupole, l = 0.6, k1 = 1e-20, tilt, num_steps = 10, x_offset = 0.05, y_pitch = 0.001
quadrupole4:       quadrupole, l = 0.6, k1 = 1, tilt, num_steps = 10, x_offset = 0.05, y_pitch = 0.001, fringe_type = full
quadrupole5:       quadrupole, l = 0.6, k1 = -1, tilt, num_steps = 10, x_offset = 0.05, y_pitch = 0.001, fringe_type = full
rcollimator1:      rcollimator, l = 4.5, x_limit = 0.09/2, y_limit = 0.05/2, num_steps = 10, x_offset = 0.02, y_pitch = 0.001, tilt = 0.1
rfcavity1:         rfcavity, l = 4.5, rf_frequency = 500e6, voltage = 5e5, phi0 = 0.1, num_steps = 10, x_offset = 0.05, 
                             y_pitch = 0.001, tilt = 0.1
rfcavity2:         rfcavity, l = 4.5, rf_frequency = 500e6, voltage = 5e5, phi0 = 0.1, num_steps = 10, x_offset = 0.05, 
                             y_pitch = 0.001, tilt = 0.1, coupler_at = both_ends, coupler_strength = 0.1
sad_mult1:         sad_mult, l = 0.1, fringe_type = full, ks = 0.1, a1 = 2e-5, b2 = 4e-3, fringe_at = exit_end
sad_mult2:         sad_mult, fringe_type = full, a1 = 2e-5, b2 = 4e-3
sbend1:            sbend, l = 0.6, k1 = 0.3, num_steps = 10, x_offset = 0.05, y_pitch = 0.001, ref_tilt = 0.1
sbend2:            sbend, l = 0.6, k2 = 3, num_steps = 10, x_offset = 0.05, y_pitch = 0.001, ref_tilt = 0.1
sbend3:            sbend, l = 0.6, k1 = 0.3, k2 = 3, num_steps = 10, x_offset = 0.05, y_pitch = 0.001, ref_tilt = 0.1
sbend4:            sbend, l = 0.6, g = 0.01, dg = -0.0099999, e1 = 0.1, e2 = 0.2, num_steps = 10, 
                          fringe_type = full, x_offset = 0.05, y_pitch = 0.001, ref_tilt = 0.1
sbend5:            sbend, l = 0.6, g = 0.01, e1 = 0.1, e2 = 0.2, num_steps = 10, fringe_type = full, 
                          x_offset = 0.05, y_pitch = 0.001, ref_tilt = 0.1
sbend6:            sbend, l = 0.6, g = 0.01, dg = -0.0099999, e1 = 0.1, e2 = 0.2, num_steps = 10, 
                          x_offset = 0.05, y_pitch = 0.001, ref_tilt = 0.1
sbend7:            sbend, l = 0.6, g = 0.01, e1 = 0.1, e2 = 0.2, num_steps = 10, x_offset = 0.05, y_pitch = 0.001, ref_tilt = 0.1
sextupole1:        sextupole, l = 0.6, k2 = 1, tilt, num_steps = 10, x_offset = 0.05, y_pitch = 0.001
solenoid1:         solenoid, l = 2.6, ks = 1, num_steps = 10, x_offset = 0.05, y_pitch = 0.001, tilt = 0.1
solenoid2:         solenoid, l = 2.6, ks = 1e-12, num_steps = 10, x_offset = 0.05, y_pitch = 0.001, tilt = 0.1
sol_quad1:         sol_quad, l = 2.6, k1 = 1, ks = 1, num_steps = 10, x_offset = 0.05, y_pitch = 0.001, tilt = 0.1
sol_quad2:         sol_quad, l = 2.6, k1 = 1e-5, ks = 1, num_steps = 10, x_offset = 0.05, y_pitch = 0.001, tilt = 0.1
vkicker1:          vkicker, l = 4.5, kick = 0.005, num_steps = 10, x_offset = 0.05, y_pitch = 0.001, tilt = 0.1
wiggler_map1:      wiggler, l = 0.1, term(1) = {0.03, 3.00, 4.00, 5.00, 0.63}, num_steps = 10, tracking_method = taylor
wiggler_periodic1: wiggler, l = 1, b_max = 0.2, n_pole = 20, num_steps = 10, x_offset = 0.01, y_pitch = 0.001, tilt = 0.1

WIG_01E: WIGGLER, 
   cartesian_map = call::wig1_18w_1.cartesian_map, L = 1.6, TILT = -2.61431456E-4, n_pole = 1,
   B_MAX = 1.9515008405509, POLARITY = 0.9765258, X_PITCH = -3.308429488E-4,
   Y_PITCH = 1.504456327E-4, X_OFFSET = 5.642326478E-4, Y_OFFSET = 5.754310708E-6,
   Z_OFFSET = -4.53921036E-4, DS_STEP = 0.16, mat6_calc_method = Symp_Lie_Bmad,
   tracking_method = Symp_Lie_Bmad, A1 = -1.522837748E-5

*[num_steps] = 3

l1: line = (ab_multipole1,beambeam1,drift1,ecollimator1,elseparator1,elseparator2,hkicker1,instrument1,kicker1,monitor1,
             multipole1,octupole1,patch1, quadrupole1,quadrupole2,quadrupole3)
l2: line = (quadrupole4,quadrupole5,rcollimator1,rfcavity1,rfcavity2,sbend1,sbend2,sbend3,sbend4,sbend5,sbend6,sbend7,
!l2: line = (quadrupole4,quadrupole5,rcollimator1,sbend1,sbend2,sbend3,sbend4,sbend5,sbend6,sbend7,
             sextupole1,solenoid1,solenoid2,sol_quad1,sol_quad2,vkicker1, sad_mult1, sad_mult2, wiggler_map1, wiggler_periodic1)
ll: line = (l1,l2)

use, ll