!! TPS lattice 'd24p79h2' is used for CSR simulation
!! Revolution period= 1.72919627e-6 sec.
!! Nominal total beam current of 400 mA and 600 bunches are assumed
!! Single bunch current(ave.)= 2/3 mA for nominal operation
!! 2/3 mA/bunch== 1.1528e-9 C
!! 1 mA/bunch== 1.72919627e-9 C, 2 mA/bunch== 3.458392539e-9 C
!! 4 mA/bunch== 6.916785078e-9 C, 5 mA/bunch== 8.64598135e-9 C
!! PJ. Chou, March 13, 2008
!! modified by M. Borland, April 14, 2008
!! modified by PJ. Chou to keep the same period length(MB added an rf cavity w/o
!!   keeping the same period length). 2008-4-15

q: charge, total= 1.72919627e-9    !bunch charge[C]--> 1 mA
DDD      : csrdrift, l = 0.125000, dz=0.01, use_stupakov=1
DL1      : csrdrift, l = 6.000000, dz=0.01, use_stupakov=1
DL1rf    : csrdrift, l = 5.800000, dz=0.01, use_stupakov=1
DL2      : csrdrift, l = 0.730000, dz=0.01, use_stupakov=1
DL2A     : csrdrift, l = 0.355000, dz=0.01, use_stupakov=1
DL2B     : csrdrift, l = 0.125000, dz=0.01, use_stupakov=1
DL3      : csrdrift, l = 0.700000, dz=0.01, use_stupakov=1 
DL4      : csrdrift, l = 0.750000, dz=0.01, use_stupakov=1
DL4A     : csrdrift, l = 0.125000, dz=0.01, use_stupakov=1
DL4B     : csrdrift, l = 0.375000, dz=0.01, use_stupakov=1 
DC5      : csrdrift, l = 0.550000, dz=0.01, use_stupakov=1 
DC6      : csrdrift, l = 0.950000, dz=0.01, use_stupakov=1
DC6A     : csrdrift, l = 0.275000, dz=0.01, use_stupakov=1
DC6B     : csrdrift, l = 0.425000, dz=0.01, use_stupakov=1
DC7      : csrdrift, l = 0.255000, dz=0.01, use_stupakov=1 
DS1      : csrdrift, l = 3.500000, dz=0.01, use_stupakov=1 
DS2      : csrdrift, l = 0.680000, dz=0.01, use_stupakov=1
DS2A     : csrdrift, l = 0.305000, dz=0.01, use_stupakov=1
DS2B     : csrdrift, l = 0.125000, dz=0.01, use_stupakov=1 
DS3      : csrdrift, l = 0.800000, dz=0.01, use_stupakov=1
DS3A     : csrdrift, l = 0.425000, dz=0.01, use_stupakov=1
DS3B     : csrdrift, l = 0.125000, dz=0.01, use_stupakov=1
DS4      : csrdrift, l = 0.320000, dz=0.01, use_stupakov=1

!!{** Quadrupole **}      
QS5      : kquad, l = 0.300000, k1 =   1.601058, n_kicks=20
QS4      : kquad, l = 0.300000, k1 =  -0.975030, n_kicks=20
QS3      : kquad, l = 0.300000, k1 =  -1.353903, n_kicks=20
QS2      : kquad, l = 0.600000, k1 =   1.545859, n_kicks=20
QS1      : kquad, l = 0.300000, k1 =  -1.696006, n_kicks=20
QL3      : kquad, l = 0.300000, k1 =  -1.237476, n_kicks=20
QL2      : kquad, l = 0.600000, k1 =   1.297832, n_kicks=20
QL1      : kquad, l = 0.300000, k1 =  -1.078493, n_kicks=20

!!{** Dipole **}
!!MB: reduce number of kicks from 200 to 30
BEND: csrcsbend, l=1.1, angle=0.1308996939, k1=0.0, e1=0.06544984695, e2=&
0.06544984695, integration_order=4, n_kicks=30,nonlinear=1,&
bins=600,high_frequency_cutoff0=0.2,high_frequency_cutoff1=0.25,&
derbenev_criterion_mode="evaluate"
bend1: csrcsbend, l=1.1, angle=0.1308996939, k1=0.0, e1=0.06544984695, e2=&
0.06544984695, integration_order=4, n_kicks=30,nonlinear=1,&
bins=600,high_frequency_cutoff0=0.2,high_frequency_cutoff1=0.25,&
derbenev_criterion_mode="evaluate"
bend8: csrcsbend, l=1.1, angle=0.1308996939, k1=0.0, e1=0.06544984695, e2=&
0.06544984695, integration_order=4, n_kicks=30,nonlinear=1,&
bins=600,high_frequency_cutoff0=0.2,high_frequency_cutoff1=0.25,&
derbenev_criterion_mode="evaluate"
! original bins=600

!!{** Sextupole **}
!SS1       : ksext,l=0.025, k2 = 12.351504, n_kicks=10  
!SS2       : ksext,l=0.025, k2 = -28.418568, n_kicks=10 
!SS3       : ksext,l=0.025, k2 = 38.783568, n_kicks=10  
!SS4       : ksext,l=0.025, k2 = -42.464976, n_kicks=10 
!SS5       : ksext,l=0.025, k2 = 32.451576, n_kicks=10  
!SS6       : ksext,l=0.025, k2 = -40.431864, n_kicks=10 
!SSD       : ksext,l=0.025, k2 = -28.2858, n_kicks=10   
!SSF       : ksext,l=0.025, k2 = 38.693352, n_kicks=10  
!S1  : LINE=(SS1,SS1,SS1,SS1,SS1)
!S2  : LINE=(SS2,SS2,SS2,SS2,SS2) 
!S3  : LINE=(SS3,SS3,SS3,SS3,SS3) 
!S4  : LINE=(SS4,SS4,SS4,SS4,SS4) 
!S5  : LINE=(SS5,SS5,SS5,SS5,SS5)
!S6  : LINE=(SS6,SS6,SS6,SS6,SS6)
!SD  : LINE=(SSD,SSD,SSD,SSD,SSD)
!SF  : LINE=(SSF,SSF,SSF,SSF,SSF)

S1       : ksext,l=0.125, k2 = 12.351504, n_kicks=10  
S2       : ksext,l=0.125, k2 = -28.418568, n_kicks=10 
S3       : ksext,l=0.125, k2 = 38.783568, n_kicks=10  
S4       : ksext,l=0.125, k2 = -42.464976, n_kicks=10 
S5       : ksext,l=0.125, k2 = 32.451576, n_kicks=10  
S6       : ksext,l=0.125, k2 = -40.431864, n_kicks=10 
SD       : ksext,l=0.125, k2 = -28.2858, n_kicks=10   
SF       : ksext,l=0.125, k2 = 38.693352, n_kicks=10  

BPM      : MONI, L=0.0
KHV      : KICKER, HKICK=0.000,VKICK = 0.000

!!MB: Need rf cavity for synchrotron motion
rf       : rfca,l=0.0,volt=3.5e6,phase="180 8.526886888778189e-01 3.5 / dasin -",freq="864 1.72919627e-6  /",change_t=1

LLONG    :LINE=(DL1,BPM, QL1, DL2A, S1,KHV,S1, DL2B, QL2, DL3,&
               BPM,QL3, DL4A, S2,KHV,S2, DL4B)
RLONG    :LINE=(DS1,BPM, QS1, DS2A, S3,KHV,S3, DS2B, QS2, DS3A,&
               BPM,S4,KHV,S4, DS3B, QS3,DS4)
HLACH    :LINE=(DC5, QS4,BPM, DC6A, SD,KHV,SD, DC6B, QS5, DC7, &
               SF,KHV,SF,BPM)
SHORT    :LINE=(DS1,BPM, QS1, DS2A, S5,KHV,S5, DS2B, QS2, DS3A,&
               BPM,S6,KHV,S6, DS3B, QS3, DS4)
HSACH    :LINE=(DC5, QS4,BPM, DC6A, SD,KHV,SD, DC6B, QS5, DC7, &
               SF,KHV,SF,BPM)
HSACH1   :LINE=(DC5,QS4,BPM, DC6A, SD,KHV,SD, DC6B, QS5, DC7)
LCELL    :LINE=(LLONG, BEND1, HLACH, -HSACH1, BEND, -RLONG)
SCELL1    :LINE=(SHORT, BEND, HSACH, -HSACH1, BEND, -SHORT)
SCELL2    :LINE=(SHORT, BEND, HSACH, -HSACH1, BEND, -SHORT)
ILCELL   :LINE=(RLONG, BEND, HSACH1,-HLACH, BEND8, -LLONG)

w1: watch,filename="%s.w1",mode="coordinates"
w2: watch,filename="%s.w2",mode="parameters",flush_interval=1

!!MB: Include lumped classical/quantum radiation effects, with defined emittance ratio ("coupling")
sr       : sreffects,coupling=0.01
SUPERA   :LINE=(LCELL, SCELL1, SCELL2, ILCELL)
ring: line=(q, sr, rf, w1, w2, 6*SUPERA)

return