! ILCDR FODO LATTICE, version FODO4 ! Created by Sun Yi-Peng (PKU/IHEP), Aug. 01, 2007 ! Decrease the number of wiggler sections from 8 to 2 ! Use the new dispersion suppressor, injection\extraction, and Chicane ! Momentum compaction can be freely tuned between 2 and 6*10^-4 ! Adjust the arc cell length and number. ! Half quads and sextupoles cleaned by Dr. Weishi Wan of LBNL. ! Optimization work is done together with Dr. Weishi Wan of LBNL. ! Wiggler and RF design from OCS lattice. ! Baseline arc cell Phase advance 72/72 ! Beam Energy E0 := 5.0 !-------------------------------------Magnets strength and RF voltage--------- !--common-- K1QDINJ1 := -0.417487446272 K1QFINJ1 := 0.294579826789 !---------- !--90 degree arc cell, for alpha=2e-4-- !------------------------------- comment K1QFA := -0.4820134000000003 K1QDA := 0.4882812000000000 K1QFA1 := -0.5056879815423139 K1QDA1 := 0.4905187924276428 K1QFA2 := -0.5056879815423139 K1QDA2 := 0.4905187924276428 K1QFAs := -0.5056879815423139 K1QDMS1 := 0.2347994569919224 K1QFMS1 := -0.1715394783876347 K1QDMS2 := 0.1895012119349973 K1QFMS2 := -0.2261768333441385 K1QDMW1 := 0.5408258475354180 K1QFMW1 := -0.3176534008589770 K1QDMW2 := 0.1703463156704844 K1QFMW2 := -0.3904019959250564 K1QDMW3 := 0.1530383772463411 K2.SF = 0.4590075143*2 K2.SD = -1.1155334423*2 !--RF voltage-- VRF := 15 endcomment !------------------------------- !--60 degree arc cell, for alpha=6e-4-- !------------------------------- comment K1QFA := -0.3346122060699204 K1QDA := 0.3393936767839698 K1QFA1 := -0.2374856816866342 K1QDA1 := 0.3560718711285560 K1QFA2 := -0.3560373627441347 K1QDA2 := 0.1179405107729569E-01 K1QFAs := 0.8301031302299172 K1QDMS1 := 0.6128539796453565E-01 K1QFMS1 := -0.3672718181809480 K1QDMS2 := 0.1408553491155269 K1QFMS2 := -0.1690566864432421 K1QDMW1 := 0.3312941181105442 K1QFMW1 := -0.5272714713269512 K1QDMW2 := 0.1070898647670208 K1QFMW2 := -0.3018323164123268 K1QDMW3 := 0.4617437525864292E-01 K2.SF = 0.1967426185728997*2 K2.SD = -0.3456509856057895*2 !--RF voltage-- VRF := 31 endcomment !------------------------------- !--72 degree arc cell, for alpha=4e-4-- !------------------------------- !comment K1QFA := -0.3985402996900860 K1QDA := 0.3988894940400053 K1QFA1 := -0.3600204591990463 K1QDA1 := 0.4173115393767193 K1QFA2 := -0.3789683877324434 K1QDA2 := 0.2035430894222265 K1QFAs := 0.4138387120708208 K1QDMS1 := 0.1522762649374542 K1QFMS1 := -0.3139571798303095 K1QDMS2 := 0.1412389008290027 K1QFMS2 := -0.1651688175183891 K1QDMW1 := 0.4191381700233475 K1QFMW1 := -0.4683848561066384 K1QDMW2 := 0.1119024327168110 K1QFMW2 := -0.2902800023702720 K1QDMW3 := 0.8484382914563864E-02 K2.SF = 0.2765984551082348*2 K2.SD = -0.5541814175344777*2 !--RF voltage-- VRF := 22 !endcomment !------------------------------- !---------------------------------------------------- ! BPMs, KICKERs, and MARKERs !---------------------------------------------------- MARKSUN1: MARKER MARKSUN2: MARKER MARKSUN3: MARKER MARKSUN4: MARKER MARKSUN5: MARKER MARKSUN6: MARKER M0: MARKER MRK: MARKER MRKL: MARKER MRKA: MARKER MRKAD: MARKER MRKW: MARKER MRKT: MARKER MRKRF: MARKER MARKINJ1: MARKER MARKINJ2: MARKER !-------------------------------------------------------------------- !-----------start of lattice----------------------------------------- !---------------------------------------------------- ! RF Cavities !---------------------------------------------------- RF_ON := 0 RFVOLT := RF_ON*VRF/32 ! 650MHz RF, 12 RF Cavities RFC : RFCAVITY, L=0.23, VOLT=RFVOLT, LAG=0.5, HARMON=14042 RFC1: DRIFT, L=0.23 ! --------Arc section !---------------------------------------------------- BANGLE := TWOPI/184 HBANGLE := BANGLE/2 QBANGLE := BANGLE/4 BA2: SBEND, L=2, ANGLE=HBANGLE, E1=QBANGLE, E2=QBANGLE BA2F2: SBEND, L=1, ANGLE=QBANGLE, E1=0.0, E2=QBANGLE BA2B2: SBEND, L=1, ANGLE=QBANGLE, E1=QBANGLE, E2=0.0 QFA: QUADRUPOLE, L=0.2, K1= K1QFA QDA: QUADRUPOLE, L=0.2, K1= K1QDA SXT_ON := 1 SF : SEXTUPOLE, L=0.25, K2= K2.SF*SXT_ON SD : SEXTUPOLE, L=0.25, K2= K2.SD*SXT_ON SF2 : SEXTUPOLE, L=0.25, K2= K2.SF*SXT_ON SD2 : SEXTUPOLE, L=0.25, K2= K2.SD*SXT_ON DSF: DRIFT, L=0.2 DSD: DRIFT, L=0.2 DLA1: DRIFT, L=6.65 DLA2: DRIFT, L=3.6 DLA3: DRIFT, L=4.7 DFODO1: DRIFT,L=9.55 DFODO2: DRIFT,L=2.05 DFODO3: DRIFT,L=10 DFODO4: DRIFT,L=2.5 ACELL: Line= ( QFA, DSD, SD, DFODO2, BA2, DFODO3, & QDA, DSF, SF, DFODO1, BA2, DFODO4 ) ACELL1st: Line= ( DSD, SD, DFODO2, BA2, DFODO3, & QDA, DSF, SF, DFODO1, BA2, DFODO4 ) ACELL2: Line= ( QFA, DSD, SD2, DFODO2, BA2, DFODO3, & QDA, DSF, SF2, DFODO1, BA2, DFODO4 ) ! --------Dispersion Suppressor section !---------------------------------------------------- BANGLEs1 := TWOPI/184/2 HBANGLEs1 := BANGLEs1/2 QBANGLEs1 := BANGLEs1/4 BA2s1: SBEND, L=2, ANGLE=HBANGLEs1, E1=QBANGLEs1, E2=QBANGLEs1 BA2F2s1: SBEND, L=1, ANGLE=QBANGLEs1, E1=0.0, E2=QBANGLEs1 BA2B2s1: SBEND, L=1, ANGLE=QBANGLEs1, E1=QBANGLEs1, E2=0.0 QFA1: QUADRUPOLE, L=0.2, K1= K1QFA1 QDA1: QUADRUPOLE, L=0.2, K1= K1QDA1 SF1 : SEXTUPOLE, L=0.25, K2= 0 SD1 : SEXTUPOLE, L=0.25, K2= 0 ACELLs1: Line= ( QFA1, DSD, SD1, DFODO2, BA2s1, DFODO3, & QDA1, DSF, SF1, DFODO1, BA2s1, DFODO4 ) ACELLs1R: Line= ( DSD, SD1, DFODO2, BA2s1, DFODO3, & QDA1, DSF, SF1, DFODO1, BA2s1, DFODO4, QFA1 ) !---------------------------------------------------- BANGLEs2 := TWOPI/184/2 HBANGLEs2 := BANGLEs2/2 QBANGLEs2 := BANGLEs2/4 BA2s2: SBEND, L=2, ANGLE=HBANGLEs2, E1=QBANGLEs2, E2=QBANGLEs2 BA2F2s2: SBEND, L=1, ANGLE=QBANGLEs2, E1=0.0, E2=QBANGLEs2 BA2B2s2: SBEND, L=1, ANGLE=QBANGLEs2, E1=QBANGLEs2, E2=0.0 QFA1H: QUADRUPOLE, L=0.1, K1= K1QFAs QFA2: QUADRUPOLE, L=0.2, K1= K1QFA2 QDA2: QUADRUPOLE, L=0.2, K1= K1QDA2 ACELLs2: Line= ( QFA2, DSD, SD1, DFODO2, BA2s2, DFODO3, & QDA2, DSF, SF1, DFODO1, BA2s2, DFODO4, QFA1H ) ACELLs2R: Line= ( QFA1H,DSD, SD1, DFODO2, BA2s2, DFODO3, & QDA2, DSF, SF1, DFODO1, BA2s2, DFODO4, QFA2 ) ACELLs0: LINE= (ACELLs1,ACELLs2) ACELLs0R: Line= (ACELLs2R,ACELLs1R) !---------------------------------------------------- ! Wiggler Section !---------------------------------------------------- D1W : DRIFT, L=0.5 D2W : DRIFT, L=0.1 D3W : DRIFT, L=0.025 BWP : SBEND, L=0.05*2, ANGLE=0.473418870904E-2*2, & E1=0.473418870904E-2, E2=0.473418870904E-2 BWPU : SBEND, L=0.05, ANGLE=0.473418870904E-2, & E1=0.473418870904E-2, E2=0.0 BWPD : SBEND, L=0.05, ANGLE=0.473418870904E-2, & E1=0.0, E2=0.473418870904E-2 BWM : SBEND, L=0.05*2, ANGLE=-0.473418870904E-2*2, & E1=-0.473418870904E-2, E2=-0.473418870904E-2 BWMU : SBEND, L=0.05, ANGLE=-0.473418870904E-2,& E1=-0.473418870904E-2, E2=0.0 BWMD : SBEND, L=0.05, ANGLE=-0.473418870904E-2,& E1=0.0, E2=-0.473418870904E-2 BWPH : SBEND, L=0.05, ANGLE=0.473418870904E-2, & E1=0.236709435452E-2, E2=0.236709435452E-2 BWMH : SBEND, L=0.05, ANGLE=-0.473418870904E-2,& E1=-0.236709435452E-2, E2=-0.236709435452E-2 !BWP : LINE=(BWPU,BWPD) !BWM : LINE=(BWMU,BWMD) WIGPRD: Line= ( BWM, D2W, BWP, D2W ) WIGPRD5: Line= ( WIGPRD, WIGPRD, WIGPRD, WIGPRD, WIGPRD ) WIGSECU: Line= ( BWPH, D3W, WIGPRD5, BWM, D2W, BWP, D3W, BWMH ) WIGSECD: Line= ( BWMH, D3W, BWP, D2W, WIGPRD5, BWM, D3W, BWPH ) !QDWH : QUADRUPOLE, L=0.15, K1= -5.956951E-01 QFWH : QUADRUPOLE, L=0.15, K1= 6.565377E-01 QDW : QUADRUPOLE, L=0.3, K1= -5.956951E-01 QFW : QUADRUPOLE, L=0.3, K1= 6.565377E-01 CWIG: Line= ( MRKW, QFW, D1W, WIGSECU, D1W, & QDW, D1W, WIGSECD, D1W ) !---------------------------------------------------- ! RF Section !---------------------------------------------------- DRF : DRIFT, L = 0.75 !QDRF : QUADRUPOLE, L=0.15, K1= -6.624226E-01 !QFRF : QUADRUPOLE, L=0.15, K1= 7.015743E-01 QDRF : QUADRUPOLE, L=0.3, K1= -6.624226E-01 QFRF : QUADRUPOLE, L=0.3, K1= 7.015743E-01 RFCELLU: Line= ( MRKRF, & QFWH, DRF, RFC, DRF, DRF, RFC, DRF, & QDRF, DRF, RFC, DRF, DRF, RFC, DRF, & QFRF, DRF, RFC, DRF, DRF, RFC, DRF, & QDRF, DRF, RFC, DRF, DRF, RFC, DRF ) RFCELLD: Line= ( MRKRF, & QFW, DRF, RFC, DRF, DRF, RFC, DRF, & QDRF, DRF, RFC, DRF, DRF, RFC, DRF, & QFRF, DRF, RFC, DRF, DRF, RFC, DRF, & QDRF, DRF, RFC, DRF, DRF, RFC, DRF ) RFCELL1U: Line= ( MRKRF, & QFWH, DRF, RFC1, DRF, DRF, RFC1, DRF, & QDRF, DRF, RFC1, DRF, DRF, RFC1, DRF, & QFRF, DRF, RFC1, DRF, DRF, RFC1, DRF, & QDRF, DRF, RFC1, DRF, DRF, RFC1, DRF ) RFCELL1D: Line= ( MRKRF, & QFW, DRF, RFC1, DRF, DRF, RFC1, DRF, & QDRF, DRF, RFC1, DRF, DRF, RFC1, DRF, & QFRF, DRF, RFC1, DRF, DRF, RFC1, DRF, & QDRF, DRF, RFC1, DRF, DRF, RFC1, DRF ) !---------------------------------------------------- ! Matching to Wiggler/RF !---------------------------------------------------- QDMW1: QUADRUPOLE, L=0.3, K1= K1QDMW1 QFMW1: QUADRUPOLE, L=0.3, K1= K1QFMW1 QDMW2: QUADRUPOLE, L=0.3, K1= K1QDMW2 QFMW2: QUADRUPOLE, L=0.3, K1= K1QFMW2 QDMW3: QUADRUPOLE, L=0.3, K1= K1QDMW3 DLMW1: DRIFT, L = 7-0.886775+1.6589+0.525 DLMW2: DRIFT, L = 4+0.1 DLMW3: DRIFT, L = 12-1.2 DLMW4: DRIFT, L = 5 DLMW5: DRIFT, L = 1.5+0.6589 SUPW : LINE = (DLMW1, QFMW1, DLMW2, QDMW1, DLMW3, QFMW2, & DLMW4, QDMW2, DLMW5, QDMW3) !---------------------------------------------------- ! Matching to Injection/extraction !---------------------------------------------------- QDMS1: QUADRUPOLE, L=0.3, K1= K1QDMS1 QFMS1: QUADRUPOLE, L=0.3, K1= K1QFMS1 QDMS2: QUADRUPOLE, L=0.3, K1= K1QDMS2 QFMS2: QUADRUPOLE, L=0.3, K1= K1QFMS2 DLMS1: DRIFT, L = 7 DLMS2: DRIFT, L = 5 DLMS3: DRIFT, L = 14 DLMS4: DRIFT, L = 5 DLMS5: DRIFT, L = 7 SUPS : LINE = (DLMS1, QFMS1, DLMS2, QDMS1, DLMS3, QFMS2, & DLMS4, QDMS2) !---------------------------------------------------- ! Injection/extraction section !---------------------------------------------------- ! Fast strip line kicker ! Theta = 2*10E3(Voltage)*.3(Length)/5E9(E0)/.03(Gap)= 4E-5(rad) KICK_ON := 0 D1KA : DRIFT, L=0.05 BKM : HKICK, L=0.3, KICK=-4E-5*KICK_ON BKP : HKICK, L=0.3, KICK= 4E-5*KICK_ON KINJ1 : LINE = (D1KA, BKM, D1KA) KINJ2 : LINE = (D1KA, BKP, D1KA) KEXT2 : LINE = (D1KA, BKP, D1KA) ! Septum data from APS storage ring ! Thin septum: 2mm thickness. Peak field 0.73 T ! Thick septum: 30mm thickness. Peak field 1.08 T BSPI1 : HKICK, L=1, KICK=0.0436*KICK_ON BSPI2 : HKICK, L=1, KICK=0.0645*KICK_ON BSPE1 : HKICK, L=1, KICK=0.0436*KICK_ON BSPE2 : HKICK, L=1, KICK=0.0645*KICK_ON QDINJ1: QUADRUPOLE, L=0.3, K1= K1QDINJ1 QFINJ1: QUADRUPOLE, L=0.3, K1= K1QFINJ1 DLINJ1: DRIFT, L = 5.6 DLINJ1a: DRIFT, L = 1.1 DLINJ1b: DRIFT, L = 2 DLINJ1c: DRIFT, L = 1 DLINJ2: DRIFT, L = 3+1.5 DLINJ3: DRIFT, L = 17.706875-1.5 DINJ7: DRIFT, L = 0.5 INJ: LINE =(DLINJ1a,BSPI2,DINJ7,BSPI1, DLINJ1b, & QFINJ1, DLINJ2, QDINJ1, DLINJ3,DLINJ3, QDINJ1, & DLINJ2, QFINJ1, DLINJ1c,5*KINJ2, 15*KINJ1,3*KINJ2, & DLINJ1c,QFINJ1,DLINJ2,QDINJ1,DLINJ3,DLINJ3,QDINJ1, & DLINJ2,QFINJ1,DLINJ1b,BSPE2,DINJ7,BSPE1,DLINJ1a) !---------------------------------------------------- !ring !---------------------------------------------------- ACELL22_1st: Line= ( ACELL1st, ACELL2, ACELL, ACELL2, & ACELL, ACELL2, ACELL, ACELL2, & ACELL, ACELL2, ACELL, ACELL2, & ACELL, ACELL2, ACELL, ACELL2, & ACELL, ACELL2, ACELL, ACELL2, & ACELL, ACELL2 ) ACELL22: Line= ( ACELL, ACELL2, ACELL, ACELL2, & ACELL, ACELL2, ACELL, ACELL2, & ACELL, ACELL2, ACELL, ACELL2, & ACELL, ACELL2, ACELL, ACELL2, & ACELL, ACELL2, ACELL, ACELL2, & ACELL, ACELL2 ) CWIG10: Line= ( CWIG, CWIG, CWIG, CWIG, CWIG, & CWIG, CWIG, CWIG, CWIG, CWIG ) SECTORW: Line= ( ACELL22, ACELLs0, SUPW, RFCELLU, RFCELLD, & CWIG10, CWIG10, QFWH, -SUPW, & ACELLs0R, ACELL22_1st ) SECTORINJ: Line= ( ACELL22, ACELLs0, SUPS, INJ, & -SUPS, ACELLs0R, ACELL22_1st ) SECTORSF: Line= ( -SUPS, & ACELLs0R, ACELL22_1st ) SECTORSB: Line= ( ACELL22, ACELLs0, SUPS, INJ ) SECTORINJSHORT: Line= ( ACELL, ACELL2, ACELLs0, SUPS, INJ, -SUPS, & ACELLs0R, ACELL1st, ACELL2 ) SECTORWSHORT: Line= ( ACELL, ACELL2, ACELLs0, SUPW, RFCELLU, RFCELLD, & CWIG10, CWIG10, QFWH, -SUPW, & ACELLs0R, ACELL1st, ACELL2 ) RING : LINE = (SECTORSF, SECTORW, SECTORINJ, SECTORW, SECTORSB) !-----------end of lattice------------------------------------------- !-------------------------------------------------------------------- BEAM, PARTICLE=ELECTRON, ENERGY=E0 SETPLOT, FONT=2, LSCALE=1.5, RSCALE=1.2, SSCALE=1.0,& ASCALE=1.2, LWIDTH=5 use, ACELL print, full TWISS, SAVE PLOT,table=twiss,haxis=s,vaxis1=betx,bety,vaxis2=dx,colour=100,spline, & FILE="ACELL" use, RING print, full TWISS, SAVE PLOT,table=twiss,haxis=s,vaxis1=betx,bety,vaxis2=dx,colour=100,spline, & FILE="FODO4" use, SECTORINJSHORT print, full TWISS, SAVE PLOT,table=twiss,haxis=s,vaxis1=betx,bety,vaxis2=dx,colour=100,spline, & FILE="Injection" use, SECTORWSHORT print, full TWISS, SAVE PLOT,table=twiss,haxis=s,vaxis1=betx,bety,vaxis2=dx,colour=100,spline, & FILE="Wiggler" ASSIGN, PRINT="ILCDR.bmpm.txt" USE, RING RF_ON := 0 SXT_ON := 0 PRINT, #E TWISS, TUNES, SAVE RF_ON := 1 SXT_ON := 1 PRINT, #E TWISS, COUPLE, SAVE BMPM BEAM, PARTICLE=electron, ENERGY=E0, NPART=2.0E10, RADIATE EMIT STOP