parameter[lattice]                = g_2
parameter[particle]               = antimuon
parameter[geometry] = open
parameter[p0c]                    = 3.094353005E9
!parameter[absolute_time_tracking] = .true.
bmad_com[aperture_limit_on]       = .true.
 	 
radius = 7.112	 
x0 = 7.112	 

 dThetaQShort = 13.0*pi/180.
 dThetaQGap   =  4.0*pi/180.
 dThetaQLong  = 26.0*pi/180.
radial_field=0.
call, file=dipole.bmad
call, file=quad_curvilinear.dat

!this is the dipole to inherit with multipole quads
FULLRING_B: sbend, rho = 7.112,  tracking_method = runge_kutta, num_steps = 1000, mat6_calc_method = tracking, ds_step=0.001 !, angle = twopi
!call, file = quads/Quad_multipoles.bmad

   r0=0.045
 r04 = r0*r0*r0*r0
 r02=r0*r0

! b1 = 20177.8 * 2/r0!
! b3 = 33.0       * 4 /r0
! b5 =  -45.9     * 6 /(r0)
! b7 =  -5.5      * 8 /(r0)
! b9 =  -391.3   * 10  /(r0)
! b11 = -6.5      * 12 /(r0)
! b13 =  52.3     * 14 /(r0)

! a1=   0.1       * 2 /r0
! a3 = 0.1        * 4 /r0
! a5 =  0.1       * 6 /(r0)
!a7 = 0.2        * 8 /(r0)
!a9 =   0.1     * 10 /(r0)
!a13 = -0.1     * 14 /(r0)


! a_pole_elec

!        a1=  8.60145985E-01
!        a2= -1.24643978E+01
!        a3= -2.60241580E+01
!        a4= -3.81005337E+01
!        a5= -1.31283830E+01
!        a6= -6.72588802E-01
!        a7=  1.98844546E+01
!        a8= -9.09430395E+01
!        a9= -5.92404705E+01
!        a10= -3.62101094E+00
!        a11=  1.21498231E+01
!        a12= -2.26068486E+01
!        a13= -2.04814130E+01
! b_pole_elec

!        b1=  1.00839218E+06
!        b2= -7.04336032E+04
!        b3= -7.04839246E+03
!        b4= -8.56113124E+02
!        b5=  2.51985162E+03
!        b6=  2.60486974E+02
!        b7= -1.44357404E+03
!        b8= -1.82992071E+02
!        b9= -9.85221130E+04
!        b10=  7.39941076E+03
!        b11 = 1.18914893E+03
!        b12 =-4.79468114E+01
!       b13 = 2.55174383E+03


! a_pole_elec
 
         a1 =  -1.70468721E+01
         a2 =   6.81010618E+00
         a3 =  -5.70854043E+01
         a4 =   1.59110986E+01
         a5 =   2.55137629E+01
         a6 =  -3.80217045E+01
         a7 =   8.88223852E+00
         a8 =  -3.21586957E+01
         a9 =  -9.23164800E+01
        a10 =  -6.55720222E+01
        a11 =  -3.95737610E+01
        a12 =  -5.00512342E+01
        a13 =  -1.20906182E+02
        a14 =   1.58963963E+01
        a15 =  -4.74379840E+01
        a16 =  -2.51621736E+01
        a17 =  -2.84584837E+01
        a18 =   1.12377301E+01
        a19 =  -1.05349402E+01
        a20 =  -4.01911746E+00

! b_pole_elec

         b1 =  1.01215278E+06
         b2 = -2.73400385E+03
         b3 = -1.53012977E+04
         b4 = -6.81512724E+02
         b5 = -4.67633366E+03
         b6 =  1.38103377E+02
         b7 =  1.09622838E+03
         b8 = -2.23549991E+02
         b9 = -1.16580419E+05
        b10 =  1.38142966E+02
        b11 =  8.39035218E+02
        b12 =  4.61168344E+01
        b13 = -1.96341601E+04
        b14 = -7.26039029E+01
        b15 = -5.61139445E+03
        b16 =  1.33471239E+01
        b17 = -5.61603384E+03
        b18 = -1.70062655E+01
        b19 = -1.70100175E+03
        b20 = -1.46736280E-01


QUAD1_m:  FULLRING_B, angle = dThetaQLong,  exact_multipoles= vertically_pure,r0_elec=0.045, & 
                             b1_elec = b1, b3_elec = b3, b5_elec=b5, b7_elec=b7, b9_elec = b9, b11_elec=b11, b13_elec=b13, &
                             b2_elec = b2, b4_elec = b4, b6_elec=b6, b8_elec=b8, b10_elec = b10, b12_elec=b12, &
                             b14_elec = b14, b15_elec = b15, b16_elec=b16, b17_elec=b17, b18_elec = b18, b19_elec=b19, b20_elec=b20, &
                             a1_elec = a1, a3_elec = a3, a5_elec=a5, a7_elec=a7, a9_elec = a9, a11_elec = a11,a13_elec=a13, &
                             a2_elec = a2, a4_elec = a4, a6_elec=a6, a8_elec=a8, a10_elec = a10, a12_elec=a12, &
                             a14_elec = a14, a15_elec = a15, a16_elec=a16, a17_elec=a17, a18_elec = a18, a19_elec=a19, a20_elec=a20

test_line: line=(QUAD1, QUAD1_m)
use, test_line