!Quad multipoles
 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)

QUAD1_LONG:  FULLRING_B, angle = dThetaQLong,  field_index=0.185, 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, &
                             a1_elec = a1, a3_elec = a3, a5_elec=a5, a7_elec=a7, a9_elec = a9, a13_elec=a13  !, &

!                             wall = {superimpose = .true., thickness = 0.0, & 
!                             section={type=wall_start, s =  0, v(1) = {0.07,-0.025}, v(2) = {0.07,0.025}}, &
!                             section={type=wall_end, s =  dThetaQLong*radius, v(1) = {0.07,-0.025}, v(2) = {0.07,0.025}}}

QUAD1_SHORT: FULLRING_B, angle = dThetaQShort, field_index=0.185, 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, &
                             a1_elec = a1, a3_elec = a3, a5_elec=a5, a7_elec=a7, a9_elec = a9, a13_elec=a13 !, &
!                 wall = { superimpose = .true., thickness = 0.0,       &
!                 section = { type = wall_start, s = 0.,                  v(1)={0.05025,-0.0235}, v(2)={0.05025, 0.0235} }, &
!                 section = { type = wall_end,   s = dThetaQShort*radius, v(1)={0.05025,-0.0235}, v(2)={0.05025, 0.0235} }}


QUAD1:       LINE = (QUAD1_SHORT, QFREE, QUAD1_LONG)


QUAD2_SHORT:  FULLRING_B, angle = dThetaQShort, field_index = 0.185, 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, &
                             a1_elec = a1, a3_elec = a3, a5_elec=a5, a7_elec=a7, a9_elec = a9, a13_elec=a13
QUAD2_LONG:   FULLRING_B, angle = dThetaQLong, field_index=0.185,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, &
                             a1_elec = a1, a3_elec = a3, a5_elec=a5, a7_elec=a7, a9_elec = a9, a13_elec=a13
QUAD3_SHORT:  FULLRING_B, angle = dThetaQShort, field_index = 0.185, 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, &
                             a1_elec = a1, a3_elec = a3, a5_elec=a5, a7_elec=a7, a9_elec = a9, a13_elec=a13
QUAD3_LONG:   FULLRING_B, angle = dThetaQLong, field_index=0.185,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, &
                             a1_elec = a1, a3_elec = a3, a5_elec=a5, a7_elec=a7, a9_elec = a9, a13_elec=a13
QUAD4_SHORT:  FULLRING_B, angle = dThetaQShort, field_index = 0.185, 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, &
                             a1_elec = a1, a3_elec = a3, a5_elec=a5, a7_elec=a7, a9_elec = a9, a13_elec=a13
QUAD4_LONG:   FULLRING_B, angle = dThetaQLong, field_index=0.185,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, &
                             a1_elec = a1, a3_elec = a3, a5_elec=a5, a7_elec=a7, a9_elec = a9, a13_elec=a13
!QUAD:        LINE = (QUAD_SHORT, QFREE, QUAD_LONG)