parameter[lattice]                = g_2
parameter[particle]               = muon
parameter[geometry] = open
parameter[p0c]                    = 3.094353005E10
parameter[absolute_time_tracking] = .true.
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
call,file= bfield.bmad
call, file=quad_long.bmad
call, file=quad_short.bmad
AFREE: FULLRING, angle = lfree, aperture_type=elliptical ! Gap from injection point to first electric quadrupole
GAP: FULLRING, angle = dThetaQGap
RING: LINE = (AFREE, QUAD_SHORT,GAP,QUAD_LONG,AFREE, QUAD_SHORT,GAP,QUAD_LONG,AFREE, QUAD_SHORT,GAP,QUAD_LONG,AFREE, QUAD_SHORT,GAP,QUAD_LONG)
use, RING