no_digested

particle_start[x] = 1e-4
particle_start[px] = 0e-4

! Branch 0

rfc: rfcavity, l = 1, phi0 = 0.1, voltage = 1e6, rf_frequency = 500e6
rfc2: rfc
rfc3: rfc

mm: marker, superimpose, ref = rfc, offset = 0.1
mm3: marker, superimpose, ref = rfc3

branch0: line = (rfc, rfc2, rfc3)
branch0[geometry] = open
branch0[e_tot] = 2e6
branch0[beta_a]  = 10  
branch0[beta_b]  = 12
branch0[eta_x]   =  1.3
branch0[etap_x]  =  0.1

! Branch 1

bend: sbend, rho = 1, angle = pi/10, y_pitch = 1e-4, ref_tilt = 0.7, roll = 2e-4,
          x_offset = 0.01, y_offset = 0.02, e1 = 0.3, e2 = 0.5
bend2: bend
bend3: bend, rho = 1, angle = 1, b1 = 1

m: marker, superimpose, ref = bend, offset = 0.1
m3: marker, superimpose, ref = bend3

branch1: line = (bend, bend2, bend3)
branch1[e_tot]  = 7e6
branch1[geometry] = open
branch1[beta_a]  = 10  
branch1[beta_b]  = 12

! branch 2

sol: solenoid, ks = 0.3, l = 1
s: sextupole, k2 = 0.7, l = 1, num_steps = 6, integrator_order = 8
b: sbend, angle = 0.3, l = 1, e1 = 0.1

branch2: line = (sol, s, b)
branch2[e_tot] = 7e6
branch2[geometry] = open
branch2[beta_a]  = 10  
branch2[beta_b]  = 12
branch2[eta_x]   =  1.3
branch2[etap_x]  =  0.1

! branch 3

e_tot = 5e6
lcc: lcavity, l = 2, voltage = e_tot, phi0 = 0.1, rf_frequency = 500e6

branch3: line = (lcc)
branch3[e_tot] = e_tot
branch3[geometry] = open
branch3[beta_a]  = 20  
branch3[beta_b]  = 12
branch3[eta_x]   =  1.3
branch3[etap_x]  =  0.1

! branch 4

eg: e_gun, l = 1
d: drift, l = 1
ss: solenoid, l = 1, superimpose, ref = eg, offset = 0.5

branch4: line = (eg, d)
branch4[e_tot] = 1e6
branch4[geometry] = open
branch4[beta_a]  = 20  
branch4[beta_b]  = 12

use, branch0, branch1, branch2, branch3, branch4