hsr.bmad

call, file='rhic-lat.bmad'

rhic_rip = 590.581639053215
dx0_eic_ip6 = +0.81
dtheta0_eic_ip6 = -0.017

lwdl = 6.0
lwds = 4.4
lwc = 2

lwd2 = 2.0
lww2 = 0.4

lc2w_wd = 1.36941
lgv = 0.11
lgvfl = 0.20
lgvwd = 0.42

valve_wd: instrument, l=lgv
ogvfl: drift, l=lgvfl
ogvwd: drift, l=lgvwd
oc2w_wd: drift, l=lc2w_wd-lbq1-lflq1x
oww2: drift, l=lww2

call, file='ir2-lat.bmad'
call, file='ir4-lat.bmad'
call, file='ir6-lat.bmad'
call, file='ir8c-lat.bmad'
call, file='ir8n-lat.bmad'
call, file='ir10-lat.bmad'
call, file='ir12-lat.bmad'
call, file='rhic-ps.bmad'
call, file='ir2-ps.bmad'
call, file='ir10-ps.bmad'

qmain_ps: overlay={ y_qmain_ps[i] : +i, b_qmain_ps[i] : -i }, var={i}
qtrim_ps: overlay={ y_qtrim_ps[i] : +i, b_qtrim_ps[i] : -i }, var={i}

ip6d: marker
ip6w: marker

parameter[geometry]        = closed
parameter[particle]        = antiproton
parameter[ptc_exact_model] = t

beginning[x_position] = -rhic_rip + dx0_eic_ip6
beginning[theta_position] = dtheta0_eic_ip6 + pi

hsr_hi_n: line=(ip6d,hsr6d,
 arc07y,ins07y,ir08n,ins08y,arc09y,
 ins09y,ir10,ins10b,
 arc11b,
 mat11b,bo11_snk7,ins11b,bo11_du3_hsr,trp11b,ir12,trp12b,bi12_du3_hsr,mat12b,
 arc01b,
 ir2,
 arc03y,
 mat03y,yi3_du3_hsr,trp03y,ir4,trp04y,yo4_du3_hsr,mat04y,
 arc05y,
 hsr6w,ip6w)

hsr_lo_n: line=(ip6d,hsr6d,
 arc07y,ins07y,ir08n,ins08y,arc09y,
 ins09y,ir10_41,ins10y,
 arc11y,
 ins11y,ir12_41,trp12b,bi12_du3_hsr,mat12b,
 arc01b,
 ir2,
 arc03y,
 mat03y,yi3_du3_hsr,trp03y,ir4,trp04y,yo4_du3_hsr,mat04y,
 arc05y,
 hsr6w,ip6w)

hsr: line=(ip6d,hsr6d,
 arc07y,ir8,arc09y,
 ins09y,ir10,ins10b,
 arc11b,
 mat11b,bo11_snk7,ins11b,bo11_du3_hsr,trp11b,ir12,trp12b,bi12_du3_hsr,mat12b,
 arc01b,
 ir2,
 arc03y,
 mat03y,yi3_du3_hsr,trp03y,ir4,trp04y,yo4_du3_hsr,mat04y,
 arc05y,
 hsr6w,ip6w)

hsr41: line=(ip6d,hsr6d,
 arc07y,ir8,arc09y,
 ins09y,ir10_41,ins10y,
 arc11y,
 ins11y,ir12_41,trp12b,bi12_du3_hsr,mat12b,
 arc01b,
 ir2,
 arc03y,
 mat03y,yi3_du3_hsr,trp03y,ir4,trp04y,yo4_du3_hsr,mat04y,
 arc05y,
 hsr6w,ip6w)