!*********g-2 beam line ********
!*********CD1 version **********
!*********4/2015****************

!***************FODO transport section ***********

gqFl  = 1.241606E+00 
gqDl  =-1.157250E+00 

gq20l  = 1.708738E+00

IQ016:   marker
IQ017:   marker
IQ018:   marker
IQ019:   marker
IQ020:   marker


Q016 : quadrupole, l= l4q24_2,type="4q24",k1 = gqFl/brho/l4q24_2    
Q017 : quadrupole, l= l4q24_2,type="4q24",k1 = gqDl/brho/l4q24_2     
Q018 : quadrupole, l= l4q24_2,type="4q24",k1 = gqFl/brho/l4q24_2    
Q019 : quadrupole, l= l4q24_2,type="4q24",k1 = gqDl/brho/l4q24_2     
Q020 : quadrupole, l= l4q24_2,type="4q24",k1 = gq20l/brho/l4q24_2    

Q016I:    line = (dflan,dqcoil,IQ016,Q016,Q016,dflan,dqcoil)   
Q017I:    line = (dflan,dqcoil,IQ017,Q017,Q017,dflan,dqcoil)   
Q018I:    line = (dflan,dqcoil,IQ018,Q018,Q018,dflan,dqcoil)   
Q019I:    line = (dflan,dqcoil,IQ019,Q019,Q019,dflan,dqcoil)   
Q020I:    line = (dflan,dqcoil,IQ020,Q020,Q020,dflan,dqcoil)   


HT019A  =0
VT017A  =0

HT019 : hkicker, l= l_nda,kick=HT019A
VT017 : hkicker, l= l_nda,kick=VT017A

IVT017 :  marker
IHT019 :  marker

HT019I  : line = (dflan1,dnda_coil,IHT019,HT019,dnda_coil,dflan1)
VT017I  : line = (dflan1,dnda_coil,IVT017,VT017,dnda_coil,dflan1)


VT018A  =0
HT020A  =0

VT018 : vkicker, l= l_nda,kick=VT018A
HT020 : hkicker, l= l_nda,kick=HT020A

IVT018 : marker
IHT020 : marker

VT018I  : line = (dflan1,dnda_coil,IVT018,VT018,dnda_coil,dflan1)
HT020I  : line = (dflan1,dnda_coil,IHT020,HT020,dnda_coil,dflan1)

PWC014   : line =(B_SWIC)
PWC019   : line =(B_SWIC)

IPW019  :  marker
IPW014  :  marker

PWC014I   : line =(IPW014,PWC014)
PWC019I   : line =(IPW019,PWC019)

    

gFOdrf_1m : drift, l=3.621751349
gFOdrf_2m : drift, l=4.062386965
gFOdrf_3m : drift, l=2.045822985

dr4q24  : drift, l=-d4q24

dcto24    : drift, l=dc4q24

MOM:   marker
MOM2:   marker

MAT: marker


g_2_FODOm: line = (dm1,HT019I,gFOdrf_1m,PWC019I,Q020I,&
      gFOdrf_2m,HT020I,gFOdrf_3m)


!********matching section***********

gq12l  =  1.425462E+00 
gq13l  = -1.600554E+00  
gq14l  =  1.065328E+00  
gq15l  = -1.070378E+00 

IQ012:   marker
IQ013:   marker
IQ014:   marker
IQ015:   marker


Q012 : quadrupole, l = l4q24_2,type="4q24",k1 =  gq12l/brho/l4q24_2   
Q013 : quadrupole, l = l4q24_2,type="4q24",k1 =  gq13l/brho/l4q24_2   
Q014 : quadrupole, l = l4q24_2,type="4q24",k1 =  gq14l/brho/l4q24_2   
Q015 : quadrupole, l = l4q24_2,type="4q24",k1 =  gq15l/brho/l4q24_2 

Q012I:    line = (dflan,dqcoil,IQ012,Q012,Q012,dflan,dqcoil)   
Q013I:    line = (dflan,dqcoil,IQ013,Q013,Q013,dflan,dqcoil)   
Q014I:    line = (dflan,dqcoil,IQ014,Q014,Q014,dflan,dqcoil)   
Q015I:    line = (dflan,dqcoil,IQ015,Q015,Q015,dflan,dqcoil)   



HT011A  =0
VT011A  =0

HT011 : hkicker, l= l_nda,kick=HT011A
VT011 : hkicker, l= l_nda,kick=VT011A

IHT011: marker
IVT011: marker

HT011I  : line = (dflan1,dnda_coil,IHT011,HT011,dnda_coil,dflan1)
VT011I  : line = (0*dflan1,dnda_coil,IVT011,VT011,dnda_coil,dflan1)


dpdctoa   : drift, l=-pdatoc


drFF =0.

drF1 =5.864880E-01 
drF2 =1.354272E-01 
drF3 =-8.178969E-01 
drF4 =2.964839E-01


drFOD1 : drift, l=-0.132712
drFOD2 : drift, l=1.0936-drF1+drF2+d4q24
drFOD3 : drift, l=6.0936-drF2+drF3
!drFOD4 : drift, l=5.739728-drF3/5+drF4/5
drFOD5 : drift, l=0.2968-drF4


drFOD4 : drift, l=5.608288-drF3/5+drF4/5-0.274664158
drFOD4a : drift, l=5.608288-drF3/5+drF4/5-0.579500158
drFOD4aa : drift, l=5.608288-drF3/5+drF4/5-0.782664158
drFOD4aaa : drift, l=5.608288-drF3/5+drF4/5-0.433414158



gFODOmat:  line = (HT011I,VT011I,drFOD1,Q012I,drFOD2,&
                  Q013I,drFOD3,Q014I,drFOD4aaa,PWC014I,Q015I,&
                  drFOD4,Q016I,drFOD4,Q017I,VT017I,drFOD4aa,&
                  Q018I,drFOD4a,VT018I,Q019I,drFOD5)

rgFODOmat: line = (-gFODOmat)