option echo = f

BEAM, PARTICLE=ELECTRON, ENERGY=2.000, RADIATE, &
NPART=1.0E8, KBUNCH=244, BUNCHED
   SETPLOT, FONT=2, LSCALE=1.5, RSCALE=1.2, SSCALE=1.0,&
   ASCALE=1.2,  LWIDTH=5

!***********************!
!  General definitions  !
!***********************!
   D2R  :=  PI/180

!**********!
!  Markers !
!**********!


!**********!
!  Drifts  !
!**********!

   D11   : drift, L= 2.6225
   D12   : drift, L= 0.075
   D13A  : drift, L= 0.075
   D13B  : drift, L= 0.075
   D13C  : drift, L= 0.075
   D14A  : drift, L= 0.075
   D14B  : drift, L= 0.075
   D15  : drift, L= 0.25
   D33   : drift, L= 0.075
   DQF2  : drift, L= 0.075

!**************!
!  Sextupoles  !
!**************!

   !ksf :=   947.8157892742132*2
   !ksd :=  -757.4570494303927*2
   !ksf :=   934.9901601929981*2
   !ksd :=  -701.4701059669954*2
   !ksf :=   953.0604543385057*2
   !ksd :=  -762.2039821059361*2
   !ksf :=   951.0211513669213*2
   !ksd :=  -761.5381496911989*2
   ksf :=   948.4701611340747*2
   ksd :=  -760.0489843237999*2
   !ksf :=   836.4099889519214*2
   !ksd :=  -670.7368638458836*2
   !ksf :=   932.3190552042591*2
   !ksd :=  -784.0000868009072*2

   kshh := 0.01
   kshh2 := 0.01

   SF : sextupole, L= 0.28/2, K2= ksf
   SD : sextupole, L= 0.28/2, K2= ksd
   SHH : sextupole, L= 0.025, K2= kshh
   SHH2 : sextupole, L= 0.025, K2= kshh2

!***************!
!  Quadrupoles  !
!***************!  

   !v2
   !fac :=  0.9760586927306200
   !fac2 :=   1.002864674306957
   !kqf1 :=   13.85586313759599
   !kqd1 :=  -11.91316631639049
   !kqf2 :=   11.33706156351902
   !kqd2 :=   9.070176720528069

   !kqf1 :=   13.78757384393482
   !kqd1 :=  -11.83539140296732

   !fac :=  0.9979179530614163
   !fac2 :=   1.007500000000000
   !kqf1 :=   13.98953608558538
   !kqd1 :=  -11.37897211383264
   !kqf2 :=   11.34727374406785
   !kqd2 :=   9.117623537882565

   !kqf1 :=   13.66558790241613
   !kqd1 :=  -11.27106388494587

   !fac :=  0.9569391294489886
   !fac2 :=  0.9988000000000000
   !kqf1 :=   13.86641899834866
   !kqd1 :=  -11.69753739933391
   !kqf2 :=   11.39746713308766
   !kqd2 :=   8.778403017773806

   !kqf1 :=   13.75433759960392
   !kqd1 :=  -11.69019292603233

   !fac :=  0.9625603175522259
   !fac2 :=   1.000000000000000
   !kqf1 :=   13.85953450871226
   !kqd1 :=  -11.78967812295987
   !kqf2 :=   11.37674497530826
   !kqd2 :=   8.875288413848681

   !kqf1 :=   13.77220998186793
   !kqd1 :=  -11.76856637709610

   !fac :=  0.9859816025334327
   !fac2 :=   1.005000000000000
   !kqf1 :=   13.89098836380313
   !kqd1 :=  -11.90386600025346
   !kqf2 :=   11.31973355516637
   !kqd2 :=   9.170443401543441

   !kqf1 :=   13.78330613851694
   !kqd1 :=  -11.82050400628746

   !fac :=  0.9625603175522259
   !fac2 :=   1.000000000000000
   !kqf1 :=   13.71316968063578
   !kqd1 :=  -11.83796429858221
   !kqf2 :=   11.37946950495304
   !kqd2 :=   8.865447328995909

   !kqf1 :=   13.81395120403858
   !kqd1 :=  -11.93018932763247

   !fac :=  0.9685812430711895
   !fac2 :=   1.004694249298262
   !kqf1 :=   13.75748503114882
   !kqd1 :=  -11.95634394156803
   !kqf2 :=   11.40686648799266
   !kqd2 :=   8.808747104466098

   !v4
   !fac :=   1.044713469188013
   !fac2 :=   1.020925579957719
   !kqf1 :=   13.76124249473754
   !kqd1 :=  -11.95634394156803
   !kqf2 :=   11.28140137045003
   !kqd2 :=   9.548672864810342

   !fac_l := 5 ; fac_r := 5
   !fac :=   1.042424364729266
   !fac2 :=   1.018967070561279
   !kqf1 :=   13.77628051320097
   !kqd1 :=  -11.95634394156803
   !kqf2 :=   11.38207668826223
   !kqd2 :=   9.105455768236210

   !fac_l := 5 ; fac_r := 1.1
   !fac :=   1.044640289713829
   !fac2 :=   1.020867282248639
   !kqf1 :=   13.76176543664228
   !kqd1 :=  -11.95634394156803
   !kqf2 :=   11.28479677466247
   !kqd2 :=   9.534009698457746

   !fac_l := 5 ; fac_r := 1.25
   !fac :=   1.044532442283790
   !fac2 :=   1.020780793557143
   !kqf1 :=   13.76253045151756
   !kqd1 :=  -11.95634394156803
   !kqf2 :=   11.28977737131393
   !kqd2 :=   9.512465191345191

   !fac :=   1.050000000000000
   !fac2 :=   1.050000000000000
   !kqf1 :=   13.73052334803463
   !kqd1 :=  -11.95634394156803
   !kqf2 :=   11.63149920976716
   !kqd2 :=   8.423366425009345

   !fac :=  0.9590104122466128
   !fac2 :=  0.9953568101806221
   !kqf1 :=   15.66028347713443
   !kqd1 :=  -9.384202079748363
   !kqf2 :=   11.57615436745611
   !kqd2 :=   8.104253958504543

   !fac_l := 5 ; fac_r := 0.75
   !fac :=   1.044901055712621
   !fac2 :=   1.021073555900503
   !kqf1 :=   13.75988750701056
   !kqd1 :=  -11.95634394156803
   !kqf2 :=   11.27263816106433
   !kqd2 :=   9.586426226152060
   !kq31 :=  -2.295754413817030
   !kq32 :=  -6.858696000000000

   !fac_l := 1 ; fac_r := 0.75
   fac :=   1.045370130208120
   fac2 :=   1.021444443552734
   kqf1 :=   13.75646581259283
   kqd1 :=  -11.95634394156803
   kqf2 :=   11.25063041240804
   kqd2 :=   9.680862787452575
   kq31 :=  -2.293027930708193
   kq32 :=  -6.858696000000000

   !fac_l := 1 ; fac_r := 1.25
   !fac :=   1.044149118887584
   !fac2 :=   1.020475374919297
   !kqf1 :=   13.76523792989202
   !kqd1 :=  -11.95634394156803
   !kqf2 :=   11.30746407663828
   !kqd2 :=   9.435760711243207
   !kq31 :=  -2.300054358039655
   !kq32 :=  -6.858696000000000

   !fac_l := 2 ; fac_r := 1.25
   !fac :=   1.044337246767068
   !fac2 :=   1.020625081452635
   !kqf1 :=   13.76391261638910
   !kqd1 :=  -11.95634394156803
   !kqf2 :=   11.29879258693808
   !kqd2 :=   9.473410717491975
   !kq31 :=  -2.298985042582590
   !kq32 :=  -6.858696000000000

   !fac_l := 2 ; fac_r := 0.75
   !fac :=   1.045124100181271
   !fac2 :=   1.021249814617294
   !kqf1 :=   13.75826681602572
   !kqd1 :=  -11.95634394156803
   !kqf2 :=   11.26219211072885
   !kqd2 :=   9.631319582868953
   !kq31 :=  -2.294460877359949
   !kq32 :=  -6.858696000000000

   !fac_l := 3 ; fac_r := 0.75
   !fac :=   1.045008429786731
   !fac2 :=   1.021158372803517
   !kqf1 :=   13.75910865338619
   !kqd1 :=  -11.95634394156803
   !kqf2 :=   11.26761323594375
   !kqd2 :=   9.608036677041911
   !kq31 :=  -2.295132301247568
   !kq32 :=  -6.858696000000000

   !fac_l := 3 ; fac_r := 1.25
   !fac :=   1.044435767339106
   !fac2 :=   1.020703623617029
   !kqf1 :=   13.76321587641105
   !kqd1 :=  -11.95634394156803
   !kqf2 :=   11.29424456861910
   !kqd2 :=   9.493124000137684
   !kq31 :=  -2.298423984170595
   !kq32 :=  -6.858696000000000

   !fac_l := 4 ; fac_r := 1.25
   !fac :=   1.044494137240050
   !fac2 :=   1.020750204827775
   !kqf1 :=   13.76280223732737
   !kqd1 :=  -11.95634394156803
   !kqf2 :=   11.29154789512622
   !kqd2 :=   9.504802086432433
   !kq31 :=  -2.298091243926975
   !kq32 :=  -6.858696000000000

   !fac_l := 4 ; fac_r := 0.75
   !fac :=   1.044942941678994
   !fac2 :=   1.021106633753415
   !kqf1 :=   13.75958397219493
   !kqd1 :=  -11.95634394156803
   !kqf2 :=   11.27067877926536
   !kqd2 :=   9.594856151624587
   !kq31 :=  -2.295511858534988
   !kq32 :=  -6.858696000000000

   QF1  : quadrupole,   L= 0.19,   K1= kqf1
   QF2  : quadrupole,   L= 0.19,   K1= kqf2
   QF3  : quadrupole,   L= 0.305,   K1= 15*fac2
   QF4  : quadrupole,   L= 0.1525,   K1= 15*fac2
   QF5  : quadrupole,   L= 0.1525,   K1= 15*fac2
   QD1  : quadrupole,   L= 0.18,   K1= kqd1
   QD2  : quadrupole,   L= 0.09,   K1= kqd2

!***********!
!  Dipoles  !
!***********!

   !kq31 :=  -2.247900620775051
   !kq32 :=  -6.858696000000000

   !kq31 :=  -2.582922394689871
   !kq32 :=  -6.858696000000000

   !kq31 :=  -2.455286567150831
   !kq32 :=  -6.858696000000000

   !kq31 :=  -2.372348353751067
   !kq32 :=  -6.858696000000000

   !kq31 :=  -2.238633903134690
   !kq32 :=  -6.858696000000000

   !kq31 :=  -2.392157746473020
   !kq32 :=  -6.858696000000000

   !kq31 :=  -2.334443393637003
   !kq32 :=  -6.858696000000000

   !v4
   !kq31 :=  -2.296838833044971
   !kq32 :=  -6.858696000000000

   !fac_l := 5 ; fac_r := 5
   !kq31 :=  -2.309106162788844
   !kq32 :=  -6.858696000000000

   !fac_l := 5 ; fac_r := 1.1
   !kq31 :=  -2.297258259025512
   !kq32 :=  -6.858696000000000

   !fac_l := 5 ; fac_r := 1.25
   !kq31 :=  -2.297872752973702
   !kq32 :=  -6.858696000000000

   !kq31 :=  -2.219496475675777
   !kq32 :=  -6.858696000000000

   !kq31 :=  -4.180856017430765
   !kq32 :=  -6.858696000000000

   lbend := 0.5
   theta := 0.05817764173

   fac_l := 1
   fac_r := 0.75
  
   l1 := lbend/(1+2*fac_l)
   l2 := fac_l*l1

   theta1 := theta/(1+2*fac_l/fac_r)
   theta2 := (theta-theta1)/2

   B3_1  : sbend,    L= l1, angle= theta1, &
                     K1= -6.858696*fac, E1=0, E2=0 
   B3_2  : sbend,    L= l2, angle= theta2, &
                     K1= -6.858696*fac, E1=0, E2=0 
   B3_1H  : sbend,    L= l1/2, angle= theta1/2, &
                     K1= -6.858696*fac, E1=0, E2=0 
   B31_1 : sbend,    L= 0.34*l1/lbend, angle= theta1, &
                     K1= kq31, E1=0, E2=0
   B31_2 : sbend,    L= 0.34*l2/lbend, angle= theta2, &
                     K1= kq31, E1=0, E2=0
   B32_1 : sbend,    L= l1, angle= theta1, &
                     K1= kq32*fac, E1=0, E2=0
   B32_2 : sbend,    L= l2, angle= theta2, &
                     K1= kq32*fac, E1=0, E2=0

!*************!
!  Beamlines  !
!*************!

B3: LINE=(B3_2, B3_1, B3_2)
B3U: LINE=(B3_2, B3_1H)
B3D: LINE=(B3_1H, B3_2)
B31: LINE=(B31_2, B31_1, B31_2)
B32: LINE=(B32_2, B32_1, B32_2)

CELL1: LINE=(D11,SHH, &
   D12,QF1,D13A,QD1,D13B,SHH2,D13C,B31, &
   D14A,SD,SD,D14B,QF2,DQF2, &
   SF,SF,DQF2,QD2,D15,B32,D33, &
   QF3,D33,B3,D33,QF4, &
   QF4,D33,B3,D33,QF5, &
   QF5,D33,B3,D33,QF5, &
   QF5,D33,B3,D33,QF4, &
   QF4,D33,B3,D33,QF3, &
   D33,B32,D15,QD2,DQF2,SF,SF, &
   DQF2,QF2,D14B,SD,SD,D14A, &
   B31,D13C,SHH2,D13B,QD1,D13A,QF1,D12, &
   SHH,D11)

UNITCELL: LINE=(QF5,D33,B3U,B3D,D33,QF5)

!***************!
!  RF Cavities  !
!***************!

 RFC: RFCAV, L=1e-9, VOLT=1, LAG=161.2379/360, HARMON=244

!*********!
!  Rings  !
!*********!

   Ring1 : line=(12*cell1,rfc)

use, Ring1

print, FULL
SAVEBETA, LABEL=twissparam,place=#E
TWISS,SAVE
EMIT, DELTAP=0.0
BMPM
PLOT,table=twiss,haxis=s,vaxis1=betx,bety,vaxis2=dx,colour=100,spline

!match
!rmatrix, #s/#e
!endmatch

STOP