! 9 First cryomodule 
% 0.2563      sto cell_l1
% 0.1154      sto cell_l2
% 0.2193      sto cell_l3

% 0.0507      sto drilen
% 0.2982      sto drilen2
% 15.813e6    sto vpeak
% 93.270      sto phase1  
% 149.27      sto phase2 

% 87.270      sto phase11  
% 143.27      sto phase12 

% 1300e6      sto freq1

linA1: rfca, l = "cell_l1", phase = "phase1 ", volt = " vpeak cell_l1 * 0.54 *", phase_reference =1,&
            freq = "freq1", change_p0 =1, END1_FOCUS = 1, END2_FOCUS = 1, BODY_FOCUS_MODEL = SRS

linA2: rfca, l = "cell_l2", phase = "phase1 180.0 + ", volt = " vpeak cell_l2 *",phase_reference =2,&
            freq = "freq1", change_p0 =1, END1_FOCUS = 1, END2_FOCUS = 1, BODY_FOCUS_MODEL = SRS

linA3: rfca, l = "cell_l2", phase = "phase1 ", volt = " vpeak cell_l2 *",phase_reference =3,&
            freq = "freq1", change_p0 =1, END1_FOCUS = 1, END2_FOCUS = 1, BODY_FOCUS_MODEL = SRS

linA4: rfca, l = "cell_l3", phase ="phase1 ", volt = " vpeak cell_l3 * 0.62 *",phase_reference =4,&
            freq = "freq1", change_p0 =1, END1_FOCUS = 1, END2_FOCUS = 1, BODY_FOCUS_MODEL = SRS

linB1: rfca, l = "cell_l3", phase ="phase2 ", volt = " vpeak cell_l3 * 0.62 *",phase_reference =5,&
            freq = "freq1", change_p0 =1, END1_FOCUS = 1, END2_FOCUS = 1, BODY_FOCUS_MODEL = SRS

linB2: rfca, l = "cell_l2" phase = "phase2 180.0 + ", volt = " vpeak cell_l2 *",phase_reference =6,&
            freq = "freq1", change_p0 =1, END1_FOCUS = 1, END2_FOCUS = 1, BODY_FOCUS_MODEL = SRS

linB3: rfca, l = "cell_l2", phase = "phase2 ", volt = " vpeak cell_l2 *",phase_reference =7,&
            freq = "freq1", change_p0 =1, END1_FOCUS = 1, END2_FOCUS = 1, BODY_FOCUS_MODEL = SRS

linB4: rfca, l = "cell_l1", phase = "phase2", volt = " vpeak cell_l1 * 0.52 *",phase_reference =8,&
            freq = "freq1", change_p0 =1, END1_FOCUS = 1, END2_FOCUS = 1, BODY_FOCUS_MODEL = SRS

l1: drift, l = "drilen"
l2: drift, l = "drilen2"
LinacA: line = (linA1, linA2, linA3, linA2, linA3, linA3, linA3,linA2,linA4)

LinacB: line = (linB1, linB2, linB3, linB2, linB3, linB3, linB3,linB2,linB4)

cryo1:line= (LinacA, l1,l1,LinacB)

linA11: rfca, l = "cell_l1", phase = "phase11 ", volt = " vpeak cell_l1 * 0.54 *", phase_reference =9,&
            freq = "freq1", change_p0 =1, END1_FOCUS = 1, END2_FOCUS = 1, BODY_FOCUS_MODEL = SRS

linA12: rfca, l = "cell_l2", phase = "phase11 180.0 + ", volt = " vpeak cell_l2 *",phase_reference =10,&
            freq = "freq1", change_p0 =1, END1_FOCUS = 1, END2_FOCUS = 1, BODY_FOCUS_MODEL = SRS

linA13: rfca, l = "cell_l2", phase = "phase11 ", volt = " vpeak cell_l2 *",phase_reference =11,&
            freq = "freq1", change_p0 =1, END1_FOCUS = 1, END2_FOCUS = 1, BODY_FOCUS_MODEL = SRS

linA14: rfca, l = "cell_l3", phase ="phase11 ", volt = " vpeak cell_l3 * 0.62 *",phase_reference =12,&
            freq = "freq1", change_p0 =1, END1_FOCUS = 1, END2_FOCUS = 1, BODY_FOCUS_MODEL = SRS

linB11: rfca, l = "cell_l3", phase ="phase12 ", volt = " vpeak cell_l3 * 0.62 *",phase_reference =13,&
            freq = "freq1", change_p0 =1, END1_FOCUS = 1, END2_FOCUS = 1, BODY_FOCUS_MODEL = SRS

linB12: rfca, l = "cell_l2" phase = "phase12 180.0 + ", volt = " vpeak cell_l2 *",phase_reference =14,&
            freq = "freq1", change_p0 =1, END1_FOCUS = 1, END2_FOCUS = 1, BODY_FOCUS_MODEL = SRS

linB13: rfca, l = "cell_l2", phase = "phase12 ", volt = " vpeak cell_l2 *",phase_reference =15,&
            freq = "freq1", change_p0 =1, END1_FOCUS = 1, END2_FOCUS = 1, BODY_FOCUS_MODEL = SRS

linB14: rfca, l = "cell_l1", phase = "phase12", volt = " vpeak cell_l1 * 0.52 *",phase_reference =16,&
            freq = "freq1", change_p0 =1, END1_FOCUS = 1, END2_FOCUS = 1, BODY_FOCUS_MODEL = SRS

LinacA1: line = (linA11, linA12, linA13, linA12, linA13, linA13, linA13,linA12,linA4)

LinacB1: line = (linB11, linB12, linB13, linB12, linB13, linB13, linB13,linB12,linB14)

cryo2:line = (LinacA1,l1,l1,LinacB1)

! Drift for beam sepration
T1Sd1: lscdrift, l = 0.5,lsc = 0, bins = 1024
T1Sd2: lscdrift, l = 0.5, lsc = 0, bins = 1024
T1Sd3: lscdrift, l = 0.15, lsc = 0, bins = 1024
T1Sd4: lscdrift, l = 0.5, lsc = 0, bins = 1024

! Matching
T1Md1: lscdrift, l = 0.51231,lsc=0,bins = 1024
T1Md2: lscdrift, l = 0.05,lsc = 0, bins = 1024
T1Md3: lscdrift, l = 0.5,lsc = 0, bins = 1024

!Bending
T1Bd1: lscdrift, l = 0.2,lsc=0,bins=1024
T1Bd2: lscdrift, l = 0.3,lsc=0,bins=1024 
T1Bd3: lscdrift, l = 0.2,lsc=0,bins=1024
T1Bd4: lscdrift, l = 0.5,lsc=0,bins=1024

! Beam Sepration 
T1Squad1: quadrupole, l =0.2, k1 =-8.428991514
T1Squad2: quadrupole, l =0.1, k1 =14.87149813
T1Squad2b: quadrupole, l = 0.1, k1 = 14.87149813
T1Squad3: quadrupole, l = 0.2, k1 = -8.428991514
T1Squad4: quadrupole, l = 0.2, k1 = -8.428991514
T1Squad5: quadrupole, l = 0.1, k1 = 14.87149813
T1Squad5b: quadrupole, l = 0.1, k1 = 14.87149813
T1Squad6: quadrupole, l = 0.2, k1 = -8.428991514

! Matching
T1Mquad1: quadrupole, l = 0.2, k1 = 8.927990563
T1Mquad2: quadrupole, l = 0.2, k1 = -12.97499714
T1Mquad3: quadrupole, l = 0.2, k1 = -12.97499714
T1Mquad4: quadrupole, l = 0.2, k1 = 8.927990563

!Bending 
T1Bquad1: quadrupole, l = 0.2, k1 = -6.117992651
T1Bquad2: quadrupole, l = 0.2, k1 = 7.039998651
T1Bquad3: quadrupole, l = 0.2, k1 = -6.117992651
T1Bquad4: quadrupole, l = 0.1, k1 = 2.402997054
T1Bquad4b: quadrupole, l= 0.1, k1 = 2.402997054
T1Bquad5: quadrupole, l = 0.2, k1 = -6.117992651
T1Bquad6: quadrupole, l = 0.2, k1 = 7.039998651
T1Bquad7: quadrupole, l = 0.2, k1 = -6.117992651

T1SDip1: sbend, l= 0.761543549467, angle = 0.872664625997, tilt=-1.57079632679,&
 hgap = 0.025, fint = 0.4, e1 = -0.174532925199, e2 = 0.349065850399,  k2 = 4.4
T1SDip2: sbend, l = 0.697651735253, angle = 0.872664625997, tilt = 1.57079632679,&
 hgap = 0.025, fint = 0.4, e1 = 0.349065850399, k2 = 4.4
T1SDip3: sbend, l = 0.697651735253, angle = 0.872664625997, tilt= -1.57079632679, &
  hgap = 0.025, fint = 0.4, e2 = 0.349065850399, k2 = 4.4
T1SDip4: sbend, l = 0.697651735253, angle = 0.872664625997, tilt = 1.57079632679, &
 hgap = 0.025, fint = 0.4, e1 = 0.349065850399 , k2 = 4.4
T1Bdip1: sbend, l = 0.943106114608, angle = 0.785398163397, hgap = 0.025, fint = 0.4, k2 = 4.4
T1Bdip2: sbend, l = 0.943106114608, angle = 0.785398163397, hgap = 0.025, fint = 0.4, k2 = 4.4
T1SDip5: sbend, l = 0.697651735253, angle = 0.872664625997, tilt = 1.57079632679, &
 hgap = 0.025, fint = 0.4, e2 = 0.349065850399, k2 = 4.4
T1SDip6: sbend, l = 0.697651735253, angle = 0.872664625997, tilt = -1.57079632679, &
 hgap = 0.025, fint = 0.4, e1 =0.349065850399, k2 = 4.4
T1SDip7: sbend, l = 0.697651735253, angle = 0.872664625997, tilt = 1.57079632679,&
 hgap = 0.025, fint = 0.4, e2 = 0.349065850399, k2 = 4.4
T1SDip8: sbend, l = 0.761543549467 , angle = 0.872664625997, tilt = -1.57079632679, &
 hgap = 0.025, fint = 0.4, e1 = 0.349065850399, e2 = -0.174532925199 , k2 = 4.4

T1Start: marker
T1Mstart: marker
T1Bstart: marker
T1mitte: marker
T1Bend: marker
T1Mend: marker
T1end: marker

! drift I've used for compensation
T1_comp: lscdrift, l= 0.04788, lsc=0,bins=1024  

T1_Beam_Sep_Up: line = (T1Start, T1Sd1, T1_comp, &
 T1SDip1, T1Sd1, T1SDip2, &
 T1Sd2, T1Squad1, T1Sd3, T1Squad2, T1Squad2b, T1Sd3, T1Squad3, T1Sd2,&
 T1SDip3, T1Sd4, T1SDip4)

T1_Matching1: line  = (T1Mstart, T1Md1, T1Mquad1, T1Md2, T1Mquad2, T1Md3)

T1_Bend: line = (T1Bstart, &
 T1Bdip1, T1Bd1, T1Bquad1, T1Bd2, T1Bquad2, T1Bd2, T1Bquad3, T1Bd3,& 
 T1Bdip2, T1Bd4, T1Bquad4, T1mitte, T1Bquad4b, T1Bd4, T1Bdip2, T1Bd3,&
 T1Bquad5, T1Bd2, T1Bquad6, T1Bd2, T1Bquad7, T1Bd1, T1Bdip1, T1Bend)

T1_Matching2: line = ( T1Md3, T1Mquad3, T1Md2, T1Mquad4, T1Md1, T1Mend)

T1_Beam_Sep_Down: line = (T1SDip5, T1Sd4, T1SDip6,&
 T1Sd2, T1Squad4, T1Sd3, T1Squad5, T1Squad5b, T1Sd3, T1Squad6, T1Sd2,&
 T1SDip7, T1Sd1, T1SDip8, &
 T1Sd1, T1_comp, T1end)

T1: line = (T1_Beam_Sep_Up, T1_Matching1, T1_Bend, T1_Matching2, T1_Beam_Sep_Down)   

!==============================================!
!  T2 arc for 55 MeV
!==============================================!

T2dri1: lscdrift, l=0.50, lsc=0,bins=1024
T2dri2: lscdrift, l=0.05, lsc=0,bins=1024
T2dri3: lscdrift, l=0.0947, lsc=0,bins=1024
T2dri4: lscdrift, l=0.89, lsc=0,bins=1024
T2dri5: lscdrift, l=0.15, lsc=0,bins=1024
T2dri6: lscdrift, l=0.20, lsc=0,bins=1024
T2dri7: lscdrift, l=0.30, lsc=0,bins=1024
T2dri8: lscdrift, l=0.00203, lsc=0,bins=1024
T2dri9: lscdrift, l=1.6323, lsc=0,bins=1024

T2_comp: lscdrift, l=0.0771930769, lsc=0, bins=1024

T2quad1: quadrupole, l = 0.2, k1 = 1.06289890909
T2quad2: quadrupole, l = 0.2, k1 = -1.11904177455
T2quad3: quadrupole, l = 0.2, k1 = -8.40048388593
T2quad4: quadrupole, l = 0.2, k1 = 12.7499764587
T2quad5: quadrupole, l = 0.2, k1 = -6.75098791634
T2quad6: quadrupole, l = 0.2, k1 = 6.60648980332
T2quad7: quadrupole, l = 0.2, k1 = -6.10448901472
T2quad8: quadrupole, l = 0.2, k1 = 4.9144912119
T2quad9: quadrupole, l = 0.2, k1 = 2.40299130623

T2dip1: sbend, l= 0.351036, angle = 0.102101761242,hgap = 0.025, &
 e2 =0.102101761242, fint = 0.4, k2 = 4.4
T2dip2: sbend, l= 0.351036, angle = 0.204203522483, tilt=-3.14159265359,&
 hgap = 0.025, fint = 0.4, e1 = 0.102101761242, e2 = 0.102101761242,  k2 = 4.4
T2dip3: sbend, l= 0.351036, angle = 0.102101761242,hgap = 0.025, &
 e1=0.102101761242, fint = 0.4, k2 = 4.4
T2dip4: sbend, l= 0.697695, angle = 0.872664625997, tilt=-1.57079632679,&
 hgap = 0.025, fint = 0.4, e1 = -0.174532925199, e2 = 0.349065850399,  k2 = 4.4
T2dip5: sbend, l= 0.697695, angle = 0.872664625997, tilt=1.57079632679,&
 hgap = 0.025, fint = 0.4, e1 = 0.349065850399,k2 = 4.4
T2dip6: sbend, l= 0.697695, angle = 0.872664625997, tilt=-1.57079632679,&
 hgap = 0.025, fint = 0.4, e2 = 0.349065850399,k2 = 4.4
T2dip7: sbend, l= 0.943106,angle =0.785398163397,hgap = 0.025, fint = 0.4,k2 = 4.4
T2dip8: sbend, l= 0.697695, angle = 0.872664625997, tilt=1.57079632679,&
 hgap = 0.025, fint = 0.4, e2 = 0.349065850399,k2 = 4.4
T2dip9: sbend, l= 0.697695, angle = 0.872664625997, tilt=-1.57079632679,&
 hgap = 0.025, fint = 0.4, e1 = 0.349065850399,k2 = 4.4
T2dip10: sbend, l= 0.697695, angle = 0.872664625997, tilt=-1.57079632679,&
 hgap = 0.025, fint = 0.4, e1 = 0.349065850399, e2 = -0.174532925199,  k2 = 4.4

T2ECStart: marker
T2EcEnd: marker
T2BSStart: marker
T2MStart: marker
T2BStart: marker
T2BSEnd: marker
T2Mend: marker
T2ICStart: marker
T2ICEnd: marker

T2_Extraction_Chicane: line= (T2ECStart,T2_comp, T2dri1, T2quad1, T2dri2,T2quad2,T2dri1,&
 T2dip1, T2dri3, T2dip2, T2dri3,T2dip3, T2dri1, T2ECend)

T2_Beam_Sepration_Up: line = ( T2BSStart, T2dip4, T2dri4, T2dip5,T2dri1, T2quad3,&
  T2dri5, T2quad4, T2dri5,T2quad3,T2dri1,T2dip6, T2dri4,T2dip5,T2dri1)

T2_Matching1: line = (T2MStart, T2dri8, T2quad5, T2dri2, T2quad6, T2dri1)

T2_Bend: line = ( T2BStart, T2dip7,T2dri6, &
  T2quad7, T2dri7, T2quad8, T2dri7, T2quad7,&
  T2dri6, T2dip7, T2dri1,&
  T2quad9, &
  T2dri1, T2dip7, T2dri6, &
  T2quad7, T2dri7, T2quad8, T2dri7, T2quad7, &
  T2dri6, T2dip7)

T2_Matching2: line = ( T2dri1, T2quad6, T2dri2, T2quad5, T2dri8, T2dri1, T2Mend)

T2_Beam_Sepration_Down: line = ( T2dip8, T2dri4, T2dip9, T2dri1, T2quad3, &
  T2dri5, T2quad4, T2dri5, T2quad3, T2dri1, T2dip8, T2dri4, T2dip10, T2dri1, T2BSEnd)

T2_Injection_Chicane: line= (T2ICStart,T2dip1,T2dri3,T2dip2,T2dri3,T2dip3,T2dri1,&
  T2quad1,T2dri2,T2quad2,T2dri1,T2_comp, T2ICEnd)

T2: line= (T2_Extraction_Chicane, T2_Beam_Sepration_Up, T2_Matching1, &
  T2_Bend, T2_Matching2, T2_Beam_Sepration_Down, T2_Injection_Chicane)

!=====================================!
!T3 arc 80 MeV
!=====================================!

T3Sd1: lscdrift, l= 0.724497,LSC=0,bins=1024
T3Sd2: lscdrift, l=0.5,LSC=0,bins=1024
T3Sd3: lscdrift,l=0.15,LSC=0,bins=1024
T3Sd4: lscdrift,l= 0.724497,LSC=0,bins=1024
T3Md1: lscdrift,l= 0.5336,LSC=0,bins=1024
T3Md2: lscdrift,l=0.05,LSC=0,bins=1024
T3Md3: lscdrift,l=0.5,LSC=0,bins=1024
T3Bd1: lscdrift,l=0.2,LSC=0,bins=1024
T3Bd2: lscdrift,l=0.3,LSC=0,bins=1024
T3Bd3: lscdrift,l=0.2,LSC=0,bins=1024
T3Bd4: lscdrift,l=0.5,LSC=0,bins=1024
T3d1: lscdrift,l=0.5,LSC=0,bins=1024
T3d2: lscdrift,l=1.6323,LSC=0,bins=1024

T3Sdip1: sbend, l = 0.660443745972, angle = 0.309795942229, tilt = -1.57079632679,&
 hgap=0.025, fint=0.4, e1= -0.174532925199, e2= 0.309795942229
T3Sdip2: sbend, l= 0.386073899125, angle = 0.309795942229,&
 tilt = 1.57079632679, hgap=0.025, fint=0.4, e1= 0.309795942229
T3Sdip3: sbend, l= 0.386073899125, angle = 0.309795942229, &
 tilt= -1.57079632679, hgap=0.025, fint=0.4, e2=0.309795942229
T3Sdip4: sbend, l= 0.660443745972, angle = 0.309795942229, &
 tilt= 1.57079632679 , hgap=0.025, fint=0.4, e1= 0.309795942229
T3Bdip1: sbend, l= 0.943106114608, angle= 0.785398163397, hgap=0.025, fint=0.4
T3Bdip2: sbend, l= 0.943106114608, angle= 0.785398163397, hgap=0.025, fint=0.4
T3Sdip5: sbend, l= 0.660443745972, angle= 0.309795942229, &
 tilt= 1.57079632679, hgap=0.025, fint=0.4, e2= 0.309795942229
T3Sdip6: sbend, l= 0.386073899125,angle= 0.309795942229, &
 tilt = -1.57079632679, hgap=0.025, fint=0.4, e1= 0.309795942229
T3Sdip7: sbend, l=0.386073899125, angle= 0.309795942229,&
 tilt= 1.57079632679, hgap=0.025, fint=0.4, e2= 0.309795942229
T3Sdip8: sbend, l= 0.660443745972, angle= 0.309795942229, &
 tilt= -1.57079632679, hgap=0.025, fint=0.4, e1= 0.309795942229, e2 = -0.174532925199

T3Squad1: quadrupole, l =0.2, k1= -15.9108382
T3Squad2: quadrupole, l= 0.1, k1= 18.20915278
T3Squad2b: quadrupole, l= 0.1, k1=18.20915278
T3Squad3: quadrupole, l=0.2, k1= -15.9108382
T3Squad4: quadrupole, l=0.2, k1= -15.9108382
T3Squad5: quadrupole, l=0.1, k1= 18.20915278
T3Squad5b: quadrupole, l= 0.1 , k1=18.20915278
T3Squad6: quadrupole, l=0.2, k1=-15.9108382
T3Mquad1: quadrupole, l=0.2 , k1=-5.967123856
T3Mquad2: quadrupole, l= 0.2, k1= 9.008396855
T3Mquad3: quadrupole, l=0.2, k1=9.008396855
T3Mquad4: quadrupole, l=0.2, k1= -5.967123856
T3Bquad1: quadrupole, l= 0.2, k1=-6.113518723
T3Bquad2: quadrupole, l=0.2, k1= 11.44111386
T3Bquad3: quadrupole, l=0.2, k1=-6.113518723
T3Bquad4: quadrupole, l= 0.1, k1= 2.3982392
T3Bquad4b: quadrupole, l=0.1, k1= 2.3982392
T3Bquad5: quadrupole, l=0.2, k1=-6.113518723
T3Bquad6: quadrupole, l=0.2, k1=11.44111386
T3Bquad7: quadrupole, l=0.2, k1= -6.113518723

T3start: marker
T3Mstart: marker
T3Bstart: marker
T3mitte: marker
T3Bend: marker
T3Mend: marker
T3end: marker

T3_comp: lscdrift, l= 0.0745753846, LSc=0, bins=1024

T3: line = (T3start,&
T3d1, T3_comp, &
T3Sdip1, T3Sd1, T3Sdip2,& 
T3Sd2, T3Squad1, T3Sd3, T3Squad2, T3Squad2b, T3Sd3, T3Squad3, T3Sd2,&
T3Sdip3, T3Sd4, T3Sdip4,& 
T3Mstart,& 
T3Md1, T3Mquad1, T3Md2, T3Mquad2, T3Md3,&
T3Bstart,&
T3Bdip1, T3Bd1,& 
T3Bquad1, T3Bd2, T3Bquad2, T3Bd2, T3Bquad3,& 
T3Bd3, T3Bdip2, T3Bd4,&
T3Bquad4, T3mitte, T3Bquad4b, &
T3Bd4, T3Bdip2, T3Bd3,& 
T3Bquad5, T3Bd2, T3Bquad6, T3Bd2, T3Bquad7, &
T3Bd1, T3Bdip1, &
T3Bend,&
T3Md3, T3Mquad3, T3Md2, T3Mquad4, T3Md1, &
T3Mend, &
T3Sdip5, T3Sd4, T3Sdip6,& 
T3Sd2, T3Squad4, T3Sd3, T3Squad5, T3Squad5b, T3Sd3, T3Squad6, T3Sd2,&
T3Sdip7, T3Sd1, T3Sdip8, T3_comp, T3d1)

!================================================!
! Internal experiment set-up at 105 MeV
! It will give-->180 degree phase
!================================================!

intdri1: lscdrift,l=0.107517,LSC=0,bins=1024
intdri2: lscdrift,l=0.05,LSC=0,bins=1024
intdri3: lscdrift,l=0.50,LSC=0,bins=1024
intdri4: lscdrift,l=0.70 ,LSC=0,bins=1024
intdri5: lscdrift,l=0.85,LSC=0,bins=1024
intdri6: lscdrift,l=1.35,LSC=0,bins=1024
intdri7: lscdrift,l=0.15,LSC=0,bins=1024
intdri8: lscdrift,l=0.20,LSC=0,bins=1024
intdri9: lscdrift,l=1.90,LSC=0,bins=1024
intdri10: lscdrift,l=1.20,LSC=0,bins=1024
intdri11: lscdrift,l=0.30,LSC=0,bins=1024
intdri12: lscdrift,l=1.4375,LSC=0,bins=1024
intdri13: lscdrift,l=2.875,LSC=0,bins=1024
intdri14: lscdrift,l=0.598854,LSC=0,bins=1024
intdri15: lscdrift,l= 0.0780764,LSC=0,bins=1024
intdri16:  lscdrift,l=0.77367,LSC=0,bins=1024
! 0.09128 intdri15
int_comp: lscdrift, l=0.031285 , LSC=0, bins=1024
! 0.031285  0.0097384615
intquad1: Quadrupole, l= 0.2, k1= -4.17198038095
intquad2: Quadrupole, l= 0.2, k1= 5.03186025158
intquad3: Quadrupole, l=0.2, k1=-6.07943343025
intquad4:Quadrupole, l= 0.2, k1= 7.9886773408
intquad5:Quadrupole, l= 0.2, k1= 7.16468618057
intquad6: Quadrupole, l= 0.2, k1=-13.5907990796
intquad7: Quadrupole, l= 0.2, k1=-1.33747489783
intquad8: Quadrupole, l= 0.2, k1= 1.1128878624
intquad9: Quadrupole, l= 0.2, k1= -5.90408366187
intquad10: Quadrupole, l= 0.2, k1=6.02918829105
intquad11: Quadrupole, l= 0.2, k1=11.7117028236
intquad12: Quadrupole, l= 0.2, k1=-11.7117028236
intquad13: Quadrupole, l= 0.2, k1= -4.93991547276
intquad14: Quadrupole, l= 0.2, k1= 4.89720796099
intquad15: Quadrupole, l= 0.2, k1= -8.79910077166
intquad16: Quadrupole, l= 0.2, k1= 12.6884308699
intquad17: Quadrupole, l= 0.2, k1= 9.6115342365
intquad18: Quadrupole, l= 0.2, k1= -9.6115342365
intquad19:Quadrupole, l= 0.2, k1= -15.9833819284
intquad20: Quadrupole, l=0.2, k1= 17.1816848143
intquad21:Quadrupole, l= 0.2, k1= -5.00925022446
intquad22: Quadrupole, l= 0.2, k1= 5.7779740157
intquad23: Quadrupole, l= 0.2, k1= 0.55675657142
intquad24: Quadrupole, l= 0.2, k1= -0.58530819047

intdip1: sbend, l =0.628737, angle = 0.523598775598 , tilt = 3.14159265359, &
 hgap = 0.025, fint = 0.4, k2 = 4.4
intdip2: sbend, l = 0.838317 ,angle = 0.698131700798 ,&
 hgap = 0.025, fint = 0.4, k2 = 4.4
intdip3: sbend, l = 0.943106 ,angle= 0.785398163397 , &
 hgap = 0.025, fint = 0.4, k2 = 4.4
intdip4: sbend, l =0.350553, angle = 0.053407075111 , e2 = 0.053407075111 ,&
 hgap = 0.025, fint = 0.4, k2 = 4.4
intdip5: sbend,l=0.350553,angle= 0.106814150222 ,tilt = 3.14159265359 ,&
 hgap = 0.025,e1 = 0.053407075111,e2 =0.053407075111,fint = 0.4, k2 = 4.4
intdip6: sbend, l =0.350553,angle = 0.053407075111 , e1 =0.053407075111,&
 hgap = 0.025, fint = 0.4, k2 = 4.4
intdip7: sbend, l = 0.671524,angle = 0.439997504428 , tilt = -1.57079632679, &
 hgap = 0.025, e1 = 0.439997504428 ,e2= -0.174532925199 , fint = 0.4, k2 = 4.4
intdip8: sbend, l =0.39261,angle = 0.439997504428 , tilt = 1.57079632679 ,&
 hgap = 0.025, e2 = 0.439997504428 , fint = 0.4, k2 = 4.4
intdip9: sbend, l =0.39261,angle = 0.439997504428 ,tilt = -1.57079632679 ,&
 hgap = 0.025, e1 =0.439997504428, fint = 0.4, k2 = 4.4
intdip10: sbend, l =0.671524 ,angle =0.439997504428 , tilt= 1.57079632679, &
 hgap = 0.025, e2= 0.439997504428 , fint = 0.4, k2 = 4.4
intdip11: sbend,l= 0.671524,angle = 0.439997504428,tilt = -1.57079632679, &
 hgap=0.025, e1 = -0.174532925199 ,e2= 0.439997504428 ,fint=0.4, k2 =4.4
intdip12: sbend, l =0.39261, angle = 0.439997504428 , tilt= 1.57079632679, &
 hgap = 0.025, e1 = 0.439997504428, fint = 0.4, k2 = 4.4
intdip13: sbend,l=0.39261,angle= 0.439997504428,tilt= -1.57079632679, &
 hgap = 0.025, e2 = 0.439997504428, fint = 0.4, k2 = 4.4
intdip14: sbend, l =  0.671524  , angle = 0.439997504428 ,tilt=1.57079632679, &
 hgap = 0.025, e1= 0.439997504428, fint = 0.4, k2 = 4.4

IntStart: marker
IntChiStart: marker
IntChiend: marker
IntBS1Start:marker
IntGarc1Start: marker
IntGarc2Start: marker
Inttel1Start: marker
IntGarc3Start: marker
Inttel2Start: marker
IntBMStart: marker
IntBS2Start:marker
intEnd : marker

Exp_guidance1: line = (IntStart, &
   intdri3, intquad23,intdri2, intquad24,intdri3)

Extraction_Chicane: line = ( IntChiStart, &
   intdip4, intdri15, intdip5, intdri15,intdip6,intdri3,IntChiend)

Beam_Sepration1: line = (IntBS1Start, &
   intdip11, intdri16, intdip12,intdri3,&
   intquad19, intdri7, intquad20, intdri7,intquad19, intdri3,&
   intdip13, intdri16,intdip14, intdri14, intquad21, intdri2,&
   intquad22,intdri5, intdri3, intdri3)

Exp_guidance_arc1: line = (IntGarc1Start, &
   intdri1, intquad1, intdri2, intquad2, intdri3,intdip1,&
   intdri3,intquad3,intdri4,intquad4,&
   intdri4, intquad3, intdri3, intdip1, intdri3,&
   intquad2, intdri2, intquad1, intdri5,&
   intquad5, intdri2, intquad6, intdri2,intquad5,& 
   intdri6,intdri7, intdri7, intdri6,& 
   intquad5, intdri2, intquad6, intdri2, intquad5)

Exp_guidance_arc2: line = (IntGarc2Start, &
   intdri5, intquad7, intdri2, intquad8, intdri3, intdip2,& 
   intdri3, intquad9, intdri8, intquad10,&
   intdri3, intdip2, intdri3, intquad10, intdri8,& 
   intquad9, intdri3,intdip2,intdri3,&
   intquad8, intdri2, intquad7, intdri3)

Telescope1: line = (Inttel1Start, &
   intquad11, intdri2, intquad12,intdri9,&
   intquad11,intdri2, intquad12, intdri10)

Exp_guidance_arc3: line = ( IntGarc3Start, &
   intquad13,intdri2,intquad14, intdri3,&
    intdip3, intdri3,intquad15,&
   intdri11, intquad16, intdri11, intquad15, intdri3,& 
   intdip3,intdri3,intquad14,intdri2,intquad13)

Telescope2: line = (Inttel2Start, &
   intdri12,intquad17,intdri2, intquad18, intdri13, &
   intquad17, intdri2, intquad18,intdri13,intquad17,&
   intdri2, intquad18, intdri13, intquad17, intdri2, intquad18, intdri12)

Beam_matching: line = (IntBMstart, &
   intdri1, intquad22, intdri2, intquad21,intdri14)

Beam_Sepration2: line = (IntBS2Start, &
   intdip10, intdri16, intdip9,intdri3,intquad19,intdri7,intquad20,&
   intdri7,intquad19,intdri3,&
   intdip8,intdri16, intdip7,intdri3)

Injection_Chicane: line = ( IntChiStart, &
   intdip4, intdri15,&
   intdip5, intdri15,intdip6)

Exp_guidance2: line = (intdri3,intquad24,intdri2,intquad23,intdri3, intEnd) 

intexp: line = ( Exp_guidance1, Extraction_Chicane, Beam_Sepration1, &
    Exp_guidance_arc1, Exp_guidance_arc2, Telescope1, Exp_guidance_arc3, &
    Telescope2, Beam_matching, Beam_Sepration2, Injection_Chicane, Exp_guidance2)

! beamint: line = (T1, l2, cryo1, l2, T2, l2, cryo1, l2, T3,l2, cryo1, l2, intexp) !intexp for 180 phase  
! beamreturn: line = (l2, cryo1, l2, T3)   ! T3 is creating problem 


! just to test

l3: drift ,l =2.307692308    !  multiple of wavelength 
l4: drift ,l =1.153846154

beamline: line = (cryo1,l3,cryo1,l3,l4,cryo1)