parameter[lattice]                = g_2
parameter[particle]               = antimuon
parameter[geometry] = open
parameter[p0c]                    = 3.094353005E9
parameter[absolute_time_tracking] = .true.
parameter[custom_attribute4]      = "sbend::quad_voltage"
bmad_com[aperture_limit_on]       = .false.
beginning[beta_a]   = 63.916   ! these values give betax=1.0, betay=9.0 in inflector with Wuzeng inflector map
beginning[beta_b]   = 3.6027
beginning[alpha_a] = 19.567
beginning[alpha_b] = 3.1574

beginning[beta_a]   = 6.3228E+01    ! these values give betax=4.0, betay=16.0 in inflector with Wuzeng inflector map
beginning[beta_b]   = 5.9635E+00  
beginning[alpha_a] = 1.7864E+01  
beginning[alpha_b] = 5.5922E+00 

beginning[eta_x] = 14.262   !these give eta x = 4.0 in inflector
beginning[etap_x] = -3.9241


!==================================
! ELECTRODE SPANNING ANGLES
!==================================
 dThetaQShort = 13.0*pi/180.
 dThetaQGap   =  4.0*pi/180.
 dThetaQLong  = 26.0*pi/180.

lfree = twopi/4 - (dThetaQshort + dThetaQGap + dThetaQLong)

x0=7.112
!FULLRING: sbend, rho=7.112, tracking_method = runge_kutta, num_steps = 1000, mat6_calc_method = tracking, ds_step = 0.001

call, file=dipole.bmad
!call, file=quad_long_a.bmad


kV32=3.44086
!ratio = 22.7/32.
ratio = 1.
quad_inner  = kV32
quad_bottom = kV32 * ratio
quad_outer  = kV32
quad_top    = kV32

quad1_inner  = kV32
quad1_bottom = kV32 * ratio
quad1_outer  = kV32
quad1_top    = kV32

call, file=quad1_short_combined.bmad
call, file=quad1_long_combined.bmad


quad2_inner  = kV32 * ratio
quad2_bottom = kV32 * ratio
quad2_outer  = kV32
quad2_top    = kV32
call, file=quad2_short_combined.bmad
call, file=quad2_long_combined.bmad

quad3_inner  = kV32
quad3_bottom = kV32 * ratio
quad3_outer  = kV32
quad3_top    = kV32
quad4_inner  = kV32
quad4_bottom = kV32 * ratio
quad4_outer  = kV32 * ratio
quad4_top    = kV32



call, file=quad3_short_combined.bmad
call, file=quad4_short_combined.bmad


call, file=quad3_long_combined.bmad
call, file=quad4_long_combined.bmad

quad1_long[quad_voltage]:= KV32 
quad2_long[quad_voltage]:= KV32
quad3_long[quad_voltage]:= kV32
quad4_long[quad_voltage]:= kV32
quad1_short[quad_voltage]:= KV32 
quad2_short[quad_voltage]:= KV32
quad3_short[quad_voltage]:= kV32
quad4_short[quad_voltage]:= kV32
radius = 7.112	 


!QUAD_SHORT: FULLRING,  angle=dThetaQShort
!QUAD_LONG: FULLRING,  angle=dThetaQLong

AFREE: FULLRING,  angle = lfree, aperture_type=elliptical ! Gap from injection point to first electric quadrupole
GAP: FULLRING, angle = dThetaQGap
RING: LINE = (AFREE, QUAD1_SHORT,GAP,QUAD1_LONG,AFREE, QUAD2_SHORT,GAP,QUAD2_LONG,AFREE, QUAD3_SHORT,GAP,QUAD3_LONG,AFREE, QUAD4_SHORT,GAP,QUAD4_LONG)
use, RING