no_digested parameter[e_tot] = 1e6 parameter[n_part] = 1e10 parameter[ptc_exact_model] = true parameter[geometry] = open beginning[beta_a] = 10 beginning[beta_b] = 10 particle_start[x] = 0.0001 particle_start[px] = -0.0002 particle_start[y] = 0.0003 particle_start[py] = -0.0004 particle_start[z] = 0.0005 particle_start[pz] = 0.9 ab_multipole1: ab_multipole, a2 = 1e-2, b2 = 1e-2, x_offset = 2e-5, tilt = 0.1 drift1: drift, l = 0.5, num_steps = 10, x_offset = 5e-4, y_pitch = 1e-5, tilt = 0.1 ecollimator1: ecollimator, l = 0.5, num_steps = 10, x_offset = 0.02, y_pitch = 1e-5, tilt = 0.1, hkick = 0.0001, vkick = 0.00002 elseparator1: elseparator, l = 0.5, hkick = -0.00005, gap = 0.11, num_steps = 10, x_offset = 5e-4, y_pitch = 1e-5, tilt = 0.1 elseparator2: elseparator, l = 0.5, vkick = -0.00005, gap = 0.11, num_steps = 10, x_offset = 5e-4, y_pitch = 1e-5, tilt = 0.1 hkicker1: hkicker, l = 0.5, kick = -0.00005, num_steps = 10, x_offset = 5e-4, y_pitch = 1e-5, tilt = 0.1 instrument1: instrument, l = 0.5, num_steps = 10 kicker1: kicker, l = 0.5, hkick = -0.00003, vkick = -0.00003, num_steps = 10, x_offset = 5e-4, y_pitch = 1e-5, tilt = 0.1 lcavity1: lcavity, l = 0.5, rf_frequency = 500e6, voltage = 5e3, num_steps = 10, x_offset = 5e-4, y_pitch = 1e-5, tilt = 0.1, phi0 = 0.15, phi0_err = -0.34 lcavity2: lcavity, l = 0.5, rf_frequency = 500e6, voltage = 5, num_steps = 10, x_offset = 5e-4, y_pitch = 1e-5, tilt = 0.1, phi0 = 0.15, phi0_err = -0.34 monitor1: monitor, l = 0.5, num_steps = 10 multipole1: multipole, k1l = 0.0001, t1, k3l = 0.0045, t3 = 0.0031*pi, x_offset = 5e-4 octupole1: octupole, l = 0.5, k3 = 0.01, tilt, num_steps = 10, x_offset = 0.0001, y_pitch = 1e-5 patch1: patch, x_pitch = 0.0001, t_offset = 1e-8, z_offset = 1, tilt = 1, x_offset = -1.3e-4, y_pitch = 0.0001, y_offset = 5e-4 q0: quadrupole, a0 = 0.00001, a1 = 0.01, vkick = -0.0001, b0_elec = -2e1, a1_elec = 2e1 quadrupole1: q0, l = 0.6, k1 = 0.01, tilt = 0.1, num_steps = 10, x_offset = 5e-4, y_pitch = 1e-5 quadrupole2: q0, l = 0.6, k1 = -0.01, tilt, num_steps = 10, x_offset = 5e-4, y_pitch = 0.0001 quadrupole3: q0, l = 0.6, k1 = 1e-20, tilt, num_steps = 10, x_offset = 5e-4, y_pitch = 0.0001 quadrupole4: q0, l = 0.6, k1 = 0.01, tilt, num_steps = 10, x_offset = 5e-4, y_pitch = 0.0001, fringe_type = full quadrupole5: q0, l = 0.6, k1 = -0.01, tilt, num_steps = 10, x_offset = 5e-4, y_pitch = 0.0001, fringe_type = full rcollimator1: rcollimator, l = 0.5, num_steps = 10, x_offset = 0.02, y_pitch = 1e-5, tilt = 0.1 rfcavity1: rfcavity, l = 0.5, rf_frequency = 500e6, voltage = 1e4, phi0 = 0.1, num_steps = 10, x_offset = 5e-4, y_pitch = 0.0001, tilt = 0.1 rfcavity2: rfcavity, l = 0.5, rf_frequency = 500e6, voltage = 1e4, phi0 = 0.1, num_steps = 10, x_offset = 5e-4, y_pitch = 0.0001, tilt = 0.1, coupler_at = both_ends, coupler_strength = 0.1 sad_mult1: sad_mult, l = 0.1, fringe_type = full, ks = 0.001, a1 = 2e-5, b2 = 4e-3 sad_mult2: sad_mult, fringe_type = full, a1 = 2e-5, b2 = 4e-3 sad_mult3: sad_mult, l = 0.5, ks = 0.001 sb0: sbend, a0 = 1e-5, a1 = 0.1, vkick = -0.00001, b0_elec = -2e1, a1_elec = 2e1 sbend1: sb0, l = 0.6, k1 = 3e-3, num_steps = 10, x_offset = 5e-4, y_pitch = 1e-5, ref_tilt = 0.1, b1 = 0.01, b3 = 0.1 sbend3: sb0, l = 0.6, k1 = 3e-3, k2 = 3e-4, num_steps = 10, x_offset = 5e-4, y_pitch = 1e-5, ref_tilt = 0.1 sbend5: sb0, l = 0.6, g = 0.01, e1 = 0.001, e2 = 0.002, num_steps = 10, fringe_type = full, x_offset = 5e-4, y_pitch = 1e-5, ref_tilt = 0.1, a1 = 0.001, a3 = 0.001, fint = 0.7, hgap = 2 sbend7: sb0, l = 0.6, g = 0.01, e1 = 0.001, e2 = 0.002, num_steps = 10, x_offset = 5e-4, y_pitch = 1e-5, ref_tilt = 0.1, a1 = 0.001, a3 = 0.001 sextupole1: sextupole, l = 0.6, k2 = 0.001, tilt, num_steps = 10, x_offset = 5e-4, y_pitch = 1e-5 solenoid1: solenoid, l = 0.1, ks = 0.0001, num_steps = 10, x_offset = 5e-5, y_pitch = 1e-6, tilt = 0.1 solenoid2: solenoid, l = 0.1, ks = 1e-12, num_steps = 10, x_offset = 5e-5, y_pitch = 1e-5, tilt = 0.1 sol_quad1: sol_quad, l = 0.1, k1 = 0.001, ks = 0.0001, num_steps = 10, x_offset = 5e-4, y_pitch = 1e-5, tilt = 0.1 sol_quad2: sol_quad, l = 0.1, k1 = 1e-5, ks = 0.001, num_steps = 10, x_offset = 5e-4, y_pitch = 1e-5, tilt = 0.1 taylor1: taylor, l = 0.5, {1: 0.1, 1 1 0 0 0 0}, {2: 0.02|}, {3: 0.03|44} vkicker1: vkicker, l = 0.5, kick = 0.00005, num_steps = 10, x_offset = 5e-4, y_pitch = 1e-5, tilt = 0.1 wiggler_map1: wiggler, l = 1.6, num_steps = 10, field_calc = fieldmap, cartesian_map = { term = {3e-4, 3, 4, 5, 0.002, 3e-4, 0.63, x}, term = {3e-4, 3, 4, 5, 0.002, 3e-4, 0.63, y}, term = {3e-4, 5, 3, 4, 0.002, 3e-4, 0.63, x}, term = {3e-4, 5, 3, 4, 0.002, 3e-4, 0.63, y}, term = {3e-4, 4, 5, 3, 0.002, 3e-4, 0.63, x}, term = {3e-4, 4, 5, 3, 0.002, 3e-4, 0.63, y}}, tracking_method = symp_lie_ptc, x_offset = 5e-4, y_pitch = 1e-5, tilt = 0.1 wiggler_periodic1: wiggler, l = 1, b_max = 0.01, n_pole = 20, num_steps = 100, x_offset = 0.01, y_pitch = 1e-5, tilt = 0.1 !----------- ! Put lcavity elements last since they shift reference energy l2: line = (drift1, rfcavity1, rfcavity2, ab_multipole1, drift1, ecollimator1, elseparator1, elseparator2, hkicker1, instrument1, kicker1, monitor1, multipole1, octupole1, patch1, quadrupole1, quadrupole2, quadrupole3, quadrupole4, quadrupole5, rcollimator1, sad_mult1, sad_mult2, sad_mult3, sbend1, sbend3, sbend5, sbend7, sextupole1, solenoid1, solenoid2, sol_quad1, sol_quad2, taylor1, vkicker1, wiggler_map1, wiggler_periodic1, lcavity1, lcavity2) l2[e_tot] = 1e6 l2[geometry] = open use, l2 !particle_start[x] = 0.000 !particle_start[px] = -0.000 !particle_start[y] = 0.000 !particle_start[py] = -0.000 !particle_start[z] = 0.000 !particle_start[pz] = 2 ! !t1: taylor, l = 0.5, {1: 0.1, 1 1 0 0 0 0}, {2: 0.02|}, {3: 0.03|44}, ref_orb = (0,0,0,0,0,0) ! !l3: line = (taylor1) !use, l3