no_digested
parameter[lattice] = g_2
parameter[particle] = electron
parameter[geometry]=open
beginning[beta_a]=8.
beginning[beta_b]=22.
beginning[eta_x]=4.
parameter[e_tot] =   3.1e9

gQ18: quadrupole, l=0.35052, k1= -0.374
gQ19: quadrupole, l=0.70104, k1=  0.314
gQ20: quadrupole, l=0.70104, k1= -0.35
gQ21: quadrupole, l=0.70104, k1=  0.314
gQ22: quadrupole, l=0.70104, k1= -0.35
gQ23: quadrupole, l=0.70104, k1=  0.350
gQ24: quadrupole, l=0.70104, k1= -0.374
gQ25: quadrupole, l=0.70104, k1=  0.374
gQ26: quadrupole, l=0.70104, k1= -0.574
gQ27: quadrupole, l=0.70104, k1= 0.774

drM5: drift, l=5.439276
drF1: drift, l=3.959725
drF2: drift, l=0.471276
drF3: drift, l=2.596373
drF4: drift, l=1.964687
drF5: drift, l=0.3

!call, file='GQ18_taylor.dat'
!call, file='GQ19_taylor.dat'
!call, file='GQ20_taylor.dat'
!call, file='GQ21_taylor.dat'
!call, file='GQ22_taylor.dat'
!call, file='GQ23_taylor.dat'
call, file = 'HOLE1_taylor.dat'
call, file = 'INF_taylor.dat'
call, file = 'INF_EXIT_taylor.dat'

m5_line: line=(gQ18,drM5,gQ19,drM5,gQ20,drM5,gQ21,drM5,gQ22,&
               drM5,gQ23,drF1,gQ24,drF2,gQ25,drF3, gQ26,drF4, gQ27, drF5)

!m5_line_t: line=(gQ18_taylor,drM5,gQ19_taylor,drM5,gQ20_taylor,drM5,gQ21_taylor,drM5,gQ22_taylor,&
!               drM5,gQ23_taylor,drF1,gQ24,drF2,gQ25,drF3, gQ26,drF4, gQ27, drF5)


incoming: em_field, descrip = 'injec_fld.dat', field_calc = custom, l = 4.395, &
       tracking_method = runge_kutta, num_steps = 1000 , mat6_calc_method = tracking, ds_step=0.0001

hole1: em_field, descrip = 'injec_fld.dat', field_calc = custom, l = 2.6, & 
       tracking_method = runge_kutta, num_steps = 100 , mat6_calc_method = tracking, ds_step=0.0001

inf: em_field, descrip = 'injec_fld.dat', field_calc = custom, l = 1.7, &
       tracking_method = runge_kutta, num_steps = 1000 , mat6_calc_method = tracking, ds_step=0.001

inf_exit: em_field, descrip = 'injec_fld.dat', field_calc = custom, l = 0.095, &
       tracking_method = runge_kutta, num_steps = 1000 , mat6_calc_method = tracking, ds_step=0.001

inject_3: line=(m5_line,hole1,inf,inf_exit)
inject_t: line=(m5_line,hole1_taylor,inf_taylor,inf_exit_taylor)
inject_t_rev: line = (-inject_t)

!use, m5_line_t
!use, inject_3
!use, inject_t
use, inject_t_rev