!*********g-2 beam line ********
!*********CD1 version **********
!*********6/2013****************

beam, energy=3.1
parameter[geometry] = open
call, file=parameters_g-2_03112015.lat
call, file = g-2_FODO_04032015.lat
call, file = g-2_horzbend_04032015.lat

!forward = g_2_FFr1 6 quads
!beginning[beta_a] = 6.9758
!beginning[beta_b] = 67.8650
!beginning[alpha_a]  = 0.3759
!beginning[alpha_b] = 8.9266
!beginning[eta_x]  = 0.
!beginning[etap_x] =  0

! forward = FF from doc db 1568 - print element 61 DQ24TOSQ
!beginning[beta_a] = 9.016
!beginning[alpha_a]  = -2.117
!beginning[beta_b]=48.511
!beginning[alpha_b]=-4.184
!beginning[eta_x]=0
!beginning[etap_x]=0

!  CALCULATED BEAM PARAMETERS at:DistStartInjLine
!===================================================================================================
!             emit(m-rad)       4*emit(m-rad)             beta(m)             alpha()          gamma
!---------------------------------------------------------------------------------------------------
!  x:    1.2464487204E-05    4.9857948817E-05    7.7110113584E+00   -1.9978645479E+00    6.473162235
!  y:    1.9282092073E-05    7.7128368293E-05    3.7505192592E+01   -3.8165783977E+00    4.150430804

!beginning[beta_a] = 7.711
!beginning[beta_b]=37.505
!beginning[alpha_a]=-1.99786
!beginning[alpha_b]=-3.816578


!forward = FF   5 quads
beginning[beta_a] = 11.7688
beginning[beta_b] = 10.2033
beginning[alpha_a] = 2.5108
beginning[alpha_b] = 0.8710
beginning[eta_x]  = 0.
beginning[etap_x] =  0

!reverse
!beginning[beta_a] = 57.17
!beginning[beta_b] = 2.539
!beginning[alpha_a]  = -17.47
!beginning[alpha_b] = -0.932
!beginning[eta_x]  = 0.
!beginning[etap_x] =  0


!***************FODO section ***********
PWC021   : line =(B_SWIC)
PWC023   : line =(B_SWIC)
PWC025   : line =(B_SWIC)
IC025   : line =(MIC)

IPW021 : marker
IPW025 : marker
IIC025  : marker

PWC021I   : line =(IPW021,PWC021)
PWC025I   : line =(IPW025,PWC025)
IC025I   : line =(IIC025,IC025)


!IQ021:   marker
!IQ022:   marker
!IQ023:   marker
!IQ024:   marker
!IQ025:   marker



gq21l  =  4.533034E-02
gq22l  = -1.007691E+00
gq23l  =  2.000266E+00
gq24l  = -2.592325E+00
gq25l  =  2.038075E+00



Q021    :    quadrupole, l = lsqa_2,type="sqa",k1 =  gq21l/brho/lsqa_2
Q022    :    quadrupole, l = llqd_2,type="lqd",k1 =  gq22l/brho/llqd_2
Q023    :    quadrupole, l = llqd_2,type="lqd",k1 =  gq23l/brho/llqd_2
Q024    :    quadrupole, l = llqd_2,type="lqd",k1 =  gq24l/brho/llqd_2
Q025    :    quadrupole, l = llqb_2,type="lqb",k1 =  gq25l/brho/llqb_2


IQ021:   marker
IQ022:   marker
IQ023:   marker
IQ024:   marker
IQ025:   marker


Q021I:    line = (dcto24,dflan,dqcoil,IQ021,Q021,Q021,dqcoil,dflan,dcto24)
Q022I:    line = (dflan,dlqcoil,IQ022,Q022,Q022,dlqcoil,dflan)
Q023I:    line = (dflan,dlqcoil,IQ023,Q023,Q023,dlqcoil,dflan)
Q024I:    line = (dflan,dlqcoil,IQ024,Q024,Q024,dlqcoil,dflan)
Q025I:    line = (dflan,dqcoil,IQ025,Q025,Q025,dqcoil,dflan)

dnew1   : drift, l=3.2205214
dnew15  : drift, l=3.3792714
dnew2   : drift, l=0.75862334


dqfftr: drift, l=0.42375

HT024A  =0
VT023A  =0

HT024 : hkicker, l= l_ndb,kick=HT024A
VT023 : vkicker, l= l_ndb,kick=VT023A

IVT023 : marker
IHT024 : marker

HT024I  : line = (dflan1,dndb_coil,IHT024,HT024,dndb_coil,dflan1)
VT023I  : line = (dflan1,dndb_coil,IVT023,VT023,dndb_coil,dflan1)


!************final matrices from Dave ********************

!rbklgmatrix: matrix, RM(1,1)=1.0937,RM(1,2)=2.6811,RM(2,1)=0.29957,&
!             RM(2,2)=1.6497,&
!             RM(3,3)=0.91139,RM(3,4)=2.2901,RM(4,3)=-0.28147,RM(4,4)=0.38998


!rcryomatrix: matrix, RM(1,1)=1.0033,RM(1,2)=0.12004,RM(2,1)=0.058731,&
!             RM(2,2)=1.0036,&
!             RM(3,3)=0.99665,RM(3,4)=0.11975,RM(4,3)=-0.060518,RM(4,4)=0.99609


!rinfAmatrix: matrix, RM(1,1)=1.4502,RM(1,2)=1.6459,RM(2,1)=0.46457,&
!             RM(2,2)=1.2170,&
!             RM(3,3)=0.62617,RM(3,4)=1.3796,RM(4,3)=-0.34479,RM(4,4)=0.83072

!rinfBmatrix: matrix, RM(1,1)=1.0024,RM(1,2)=0.20016,RM(2,1)=0.020742,&
!             RM(2,2)=1.0018,&
!             RM(3,3)=0.99798,RM(3,4)=0.19987,RM(4,3)=-0.016615,RM(4,4)=0.99870


!**************************************************************

MTPT  : marker  ! end of line 30 cm from yoke
ip_l0: marker
end_line: marker

g_2_FFR1: line = (g_2_FODOm,dq24tosqa,Q021I,dq24tosqa,dnew1,PWC021I, Q022I,&    !HT019 - Q020 - Q021 - ring
                 dnew2,Q023I,&
                VT023I,0*PWC023,Q024I,HT024I,Q025I,PWC025I,IC025I,&
                dqfftr,MTPT, end_line)

FF: line = (dq24tosqa,Q021I,dq24tosqa,dnew1,PWC021I,Q022I,&     !Q021 - ring
                dnew2,Q023I,&
                VT023I,0*PWC023,Q024I,HT024I,Q025I,PWC025I,IC025I,&
                dqfftr,MTPT,end_line)

!FF: line = (ip_l0,Q021I,dq24tosqa,dnew1,PWC021I,Q022I,&
!                 dnew2,Q023I,&
!                VT023I,0*PWC023,Q024I,HT024I,Q025I,PWC025I,IC025I,&
!                dqfftr,MTPT,end_line)

rg_2_FFR1 : line = (-g_2_FFR1)

!*************Rubins last matrix*************

g_2_FFRN: line = (g_2_FFR1,rbklgmatrix,rcryomatrix,rinfAmatrix,rinfBmatrix)

rg_2_FFRN: line = (rrinfBmatrix,rrinfAmatrix,rrcryomatrix,rrbklgmatrix,&
                   rg_2_FFR1)


!*********** complete lines***************

g_2_part: line = (ghbend,gFODOmat,g_2_FFRN)


g_2_full:  line = (ECMAG_VMAT,g_2_part)

use,g_2_FFr1
!FFr: line = (-FF)
!use,FF
use, g_2_full

 ! Vars not in original file

!this matches ditkys dist that starts just before Q21 into ring.
!Q021[K1] =      1.0248283650E-01
!Q022[K1] =     -7.9402556528E-02
!Q023[K1] =      4.1493641678E-01
!Q024[K1] =     -6.5308593054E-01
!Q025[K1] =      5.9671824169E-01

!Q022[x_offset] = 0.0005
!Q023[x_offset] = -0.0005
!Q024[x_offset] = -0.0005
!Q025[x_offset] = 0.0005

!forward = FF   5 quads
beginning[beta_a] = 11.7688
beginning[beta_b] = 10.2033
beginning[alpha_a] = 2.5108
beginning[alpha_b] = 0.8710
beginning[eta_x]  = 0.
beginning[etap_x] =  0
