!DCO lattice modified with SuperB cell ; Marica Biagini for LCWS08 !modified with straight sections from DCO4 9/12/09 !Total lenght half (3238.22 m) ASSIGN, PRINT="DSB3.print" ASSIGN, ECHO="echo" SETPLOT, XSIZE=25.4, YSIZE=20.32 SETPLOT, LWIDTH=5, LSCALE=1.5, SSCALE=1.5, RSCALE=1.5 ! ============================================================================== ! Load XSIF files ! ------------------------------------------------------------------------------ call, filename="kvals.DSB3.xsif" call, filename="DSB3.xsif" ! ============================================================================== ! COMMANDS ! ------------------------------------------------------------------------------ RF_ON := 0 SXT_ON := 0 COMMENT ! NO MATCH! USE, PHT CELL VARY, K1.QPHT1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QPHT2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QPHT3, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QPHT4, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QPHT5, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QPHT6, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 CONSTRAINT, #E, BETX= 39.087, ALFX=-1.865, MUX=0.557 CONSTRAINT, #E, BETY= 10.671, ALFY= 0.536, MUY=0.563 SIMPLEX, CALLS=1000 MIGRAD, CALLS=1000 ENDMATCH savebeta, LABEL=PHTW, PLACE=#E twiss,save PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & STYLE=1, COLOUR=100, SPLINE=.T.,title="PHT",FILE="DSB3" ! ============================================================================== ! ARC CELLS ! ------------------------------------------------------------------------------ USE CX_CELL_HER ! comment cell vary, KQF, step=.001 vary, KQD, step=.001 constrai, #e,mux=0.718, MUY=0.270 lmdif, tolerance = 1e-10, calls = 2000 endmatch ! endcomment hinv: string,"((1+(alfx*alfx))/betx)*dx*dx+2*alfx*dx*dpx+betx*dpx*dpx" hinvy: string,"((1+(alfy*alfy))/bety)*dy*dy+2*alfy*dy*dpy+bety*dpy*dpy" print #e savebeta, label=arcell,place=#e TWISS, SAVE !PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & ! STYLE=1, COLOUR=100, SPLINE=.T.,title="CELL_1",FILE="DSB" !plot, table=twiss,haxis=s, vaxis= hinv style=100 USE CXpi_CELL_HER !comment MATCH,betx=arcell[betx],bety=arcell[bety],dx=arcell[dx] ! vary, THETpi, step=.001 vary, KQFBpi, step=.001 vary, KQFpi, step=.001 vary, KQDpi, step=.001 RMATRIX, SF2[1]/SF2[3], RM(1,2)=0.0,&!,RM(2,1)=0.0!, & RM(1,1)=-1.0,RM(2,2)=-1.0 constrai, #e,betx=arcell[betx],alfx=arcell[alfx],dx=arcell[dx] constrai, #e,bety=arcell[bety],alfy=arcell[alfy] lmdif, tolerance = 1e-10, calls = 2000 endmatch ! endcomment print #e savebeta, label=arcell,place=#e TWISS, SAVE,betx=arcell[betx],bety=arcell[bety],dx=arcell[dx] !PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & ! STYLE=1, COLOUR=100, SPLINE=.T.,title="CELL_2",FILE="DSB" !plot, table=twiss,haxis=s, vaxis= hinv style=100 USE cXpi_CELL_HER savebeta, label=midpicell,place=qfbpi[1] print,qfbpi[1] TWISS, SAVE,betx=arcell[betx],bety=arcell[bety],dx=arcell[dx] ! ============================================================================== ! DISPERSION SUPPRESSORS ! ------------------------------------------------------------------------------ !Matching arc cells to DSP use DISP_SUP ! comment MATCH,betx=midpicell[betx],bety=midpicell[bety],dx=midpicell[dx] vary, kqfpis, step=.001 vary, kqdpis, step=.001 constrai, MBSUP, dx= 0.,dpx=0. lmdif, tolerance = 1e-10, calls = 5000 migrad,tolerance = 1e-10, calls = 5000 endmatch ! endcomment print #e TWISS, SAVE,betx=midpicell[betx],bety=midpicell[bety],dx=midpicell[dx] ! PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & ! STYLE=1, COLOUR=100, SPLINE=.T., title="DISPSUP",FILE="DSB" use SUPP MATCH,betx=midpicell[betx],bety=midpicell[bety],dx=midpicell[dx] vary, KQFMS2, step=.001 vary, KQDMS2, step=.001 vary, KQFMS1, step=.001 vary, KQDMS1, step=.001 vary, KQDWH0, step=.001 vary, KQFWH0, step=.001 constrai, MRKW, alfx= 0.,alfy=0. constrai,#e, mux=1.71 ,muy=0.94 !!1.75, 0.95 constrai,#e,betx=midpicell[betx],bety=midpicell[bety],dx=midpicell[dx],& alfx=.0,alfy=0.,dpx=0. lmdif, tolerance = 1e-10, calls = 5000 migrad,tolerance = 1e-10, calls = 5000 endmatch savebeta, place=MRKW, label=midsupp savebeta, place=MSUP1, label=midsup1 print #e TWISS, SAVE,betx=midpicell[betx],bety=midpicell[bety],dx=midpicell[dx] !PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & ! STYLE=1, COLOUR=100, SPLINE=.T., title="SUPP",FILE="DSB" !------------------------------------------------------------------------- !Matching arc cells to PHT (DSP per DCO3) use SUP0 MATCH,betx=midpicell[betx],bety=midpicell[bety],dx=midpicell[dx] vary, KQFMS20, step=.001 vary, KQDMS20, step=.001 vary, KQFMS10, step=.001 vary, KQDMS10, step=.001 constrai,#e,betx=phtw[betx],bety= phtw[bety],& alfx= phtw[alfx], alfy=phtw[alfy] !constrai,#e,betx=4.075200613E+01,bety=9.052000311E+00,& ! alfx=-2.141000200E+00,alfy=4.910001143E-01 !dx=0.,dpx=0. lmdif, tolerance = 1e-12, calls = 5000 migrad,tolerance = 1e-12, calls = 5000 endmatch print #e TWISS, SAVE,betx=midpicell[betx],bety=midpicell[bety],dx=midpicell[dx] PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & STYLE=1, COLOUR=100, SPLINE=.T., title="SUP0",FILE="DSB" use DISP_SUP2 ! comment MATCH,betx=midpicell[betx],bety=midpicell[bety],dx=midpicell[dx] vary, kqfpis2, step=.001 vary, kqdpis2, step=.001 constrai, MBSUP, dx= 0.,dpx=0. lmdif, tolerance = 1e-10, calls = 5000 migrad,tolerance = 1e-10, calls = 5000 endmatch ! endcomment print #e TWISS, SAVE,betx=midpicell[betx],bety=midpicell[bety],dx=midpicell[dx] ! PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & ! STYLE=1, COLOUR=100, SPLINE=.T., title="DISPSUP2",FILE="DSB" ! ------------------------------------------------------------------------------ ! USE, WGC print #e savebeta, LABEL=WGCE, PLACE=#E twiss,save ! ! ------------------------------------------------------------------------------ !Matching arc cells to Wiggler (PSD in DCO3) use SUP2 MATCH,betx=midpicell[betx],bety=midpicell[bety],dx=midpicell[dx] vary, KQFMS22, step=.001 vary, KQDMS22, step=.001 vary, KQFMS12, step=.001 vary, KQDMS12, step=.001 constrai,#e,betx=wgce[betx],bety= wgce[bety],& alfx= -wgce[alfx], alfy=-wgce[alfy] lmdif, tolerance = 1e-10, calls = 5000 migrad,tolerance = 1e-10, calls = 5000 endmatch print #e TWISS, SAVE,betx=midpicell[betx],bety=midpicell[bety],dx=midpicell[dx] PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & STYLE=1, COLOUR=100, SPLINE=.T., title="SUP2",FILE="DSB" !--------1---------2---------3---------4---------5---------6---------7---------8 ! Straight Section STRSECR ! ------------------------------------------------------------------------------ USE, MPR MATCH, LINE=PHT VARY, K1.QMPR1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QMPR2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QMPR3, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QMPR4, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 ! VARY, DMPR1[L], STEP=1E-4 ! VARY, DMPR2[L], STEP=1E-4 CONSTRAINT, #E, LINE=RFC SIMPLEX, CALLS=50000 MIGRAD, CALLS=10000 ENDMATCH ! ------------------------------------------------------------------------------ USE, RFC TWISS, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & VMIN=0, VMAX=50 & STYLE=1, COLOUR=100, SPLINE=.T., TITLE="rf CELL", FILE="DSB3" ! ------------------------------------------------------------------------------ USE, MRC MATCH, LINE=RFC VARY, K1.QMRC1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QMRC2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QMRC3, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QMRC4, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 CONSTRAINT, #E, LINE=CCH SIMPLEX, CALLS=1000 MIGRAD, CALLS=1000 ENDMATCH USE, MRC TWISS, LINE=RFC, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, & VMIN=0, VMAX=50 & STYLE=1, COLOUR=100, SPLINE=.T., & TITLE="Matching RF to Chicane", FILE="DSB3" !------------------------------------------------------------------------------ USE, CCH PRINT, #e TWISS, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & VMIN=0,-1, VMAX=50,1 & STYLE=1, COLOUR=100, SPLINE=.T., & TITLE="Circumference Chicane", FILE="DSB3" !------------------------------------------------------------------------------ SCL:= 1.2643 USE, MDW MATCH, LINE=CCH VARY, K1.QMDW1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QMDW2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QMDW3, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QMDW4, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 CONSTRAINT, #E, LINE=WGC SIMPLEX, CALLS=20000 MIGRAD, CALLS=20000 ENDMATCH USE, MDW TWISS, LINE=CCH, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, & VMIN=0, VMAX=50 & STYLE=1, COLOUR=100, SPLINE=.T., & TITLE="Matching Chicane to Wiggler", FILE="DSB3" !--------1---------2---------3---------4---------5---------6---------7---------8 ! Straight Section STRSECI ! ------------------------------------------------------------------------------ USE, MPI MATCH, LINE=PHT VARY, K1.QMPI1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 ! VARY, K1.QMPI2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QMPI3, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 ! VARY, K1.QMPI4, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 CONSTRAINT, #E, LINE=IJ0 SIMPLEX, CALLS=10000 MIGRAD, CALLS=10000 ENDMATCH USE, MPI TWISS, LINE=PHT, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, & VMIN=0, VMAX=100 & STYLE=1, COLOUR=100, SPLINE=.T., & TITLE="Matching PHT to Injection", FILE="DSB3" USE, INJ TWISS, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, & VMIN=0, VMAX=80 & STYLE=1, COLOUR=100, SPLINE=.T., & TITLE="Injection", FILE="DSB3" USE, EXT TWISS, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, & VMIN=0, VMAX=80 & STYLE=1, COLOUR=100, SPLINE=.T., & TITLE="Extraction", FILE="DSB3" ! ------------------------------------------------------------------------------ USE, IFD !comment CELL VARY, K1.QIFD1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QIFD2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 CONSTRAINT, #E, MUX=0.119 CONSTRAINT, #E, MUY=0.23 SIMPLEX, CALLS=1000 MIGRAD, CALLS=1000 ENDMATCH !endcomment TWISS, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, & VMIN=0, VMAX=100 & STYLE=1, COLOUR=100, SPLINE=.T., & TITLE="IFD", FILE="DSB3" ! ------------------------------------------------------------------------------ USE, MFC MATCH, LINE=IFD VARY, K1.QMFC1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QMFC2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QMFC3, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QMFC4, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 CONSTRAINT, #E, LINE=CCH SIMPLEX, CALLS=1000 MIGRAD, CALLS=1000 ENDMATCH USE, MFC TWISS, LINE=IFD, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, & VMIN=0, VMAX=100 & STYLE=1, COLOUR=100, SPLINE=.T., & TITLE="Matching IFD to CCH", FILE="DSB3" USE, MIF MATCH, LINE=IJ0 VARY, K1.QMIF1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QMIF2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QMIF3, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 VARY, K1.QMIF4, STEP=1E-6, LOWER=-0.7, UPPER=+0.7 CONSTRAINT, #E, LINE=IFD SIMPLEX, CALLS=10000 MIGRAD, CALLS=10000 ENDMATCH USE, MIF TWISS, LINE=IJ0, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, & VMIN=0, VMAX=100 & STYLE=1, COLOUR=100, SPLINE=.T., & TITLE="Matching Injection to IFD", FILE="DSB3" ENDCOMMENT ! NO MATCH! !--------------------------------------------------------------------------- ! Output kvals !--------------------------------------------------------------------------- value THETA value THETPI value THETAS value KQF value KQD value KQFB value KQFpi value KQDpi value KQFBpi value kqfpis value kqdpis value KQFMS1 value KQDMS1 value KQFMS2 value KQDMS2 value KQDWH0 value KQFWH0 value KQDW0 value KQFW0 value KQFMS10 value KQDMS10 value KQFMS20 value KQDMS20 value KQDMS22 value KQFMS22 value KQDMS12 value KQFMS12 value kqdpis2 value kqfpis2 VALUE, K1.QPHT1 VALUE, K1.QPHT2 VALUE, K1.QPHT3 VALUE, K1.QPHT4 VALUE, K1.QPHT5 VALUE, K1.QPHT6 VALUE, K1.QMPR1 VALUE, K1.QMPR3 VALUE, K1.QMPR2 VALUE, K1.QMPR4 VALUE, K1.QRFC1 VALUE, K1.QRFC2 VALUE, K1.QCCH1 VALUE, K1.QCCH2 VALUE, K1.QMRC1 VALUE, K1.QMRC2 VALUE, K1.QMRC3 VALUE, K1.QMRC4 VALUE, K1.QMDW1 VALUE, K1.QMDW2 VALUE, K1.QMDW3 VALUE, K1.QMDW4 VALUE, K1.QWGC1 VALUE, K1.QWGC2 VALUE, K1.QMPI1 VALUE, K1.QMPI2 VALUE, K1.QMPI3 VALUE, K1.QMPI4 VALUE, K1.QINJ1 VALUE, K1.QINJ2 VALUE, K1.QMIF1 VALUE, K1.QMIF2 VALUE, K1.QMIF3 VALUE, K1.QMIF4 VALUE, K1.QIFD1 VALUE, K1.QIFD2 VALUE, K1.QMFC1 VALUE, K1.QMFC2 VALUE, K1.QMFC3 VALUE, K1.QMFC4 !--------------------------------------------------------------------------- !--------1---------2---------3---------4---------5---------6---------7---------8 ! Output the lattice functions ! ------------------------------------------------------------------------------ !COMMENT USE, PHT savebeta, LABEL=PHTW, PLACE=#E TWISS, SAVE USE CX_CELL_HER print #e TWISS, SAVE USE CXpi_CELL_HER savebeta, label=midpicell,place=qfbpi[1] print #e TWISS, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & STYLE=1, COLOUR=100, SPLINE=.T.,title="CXpi_CELL",FILE="DSB3" USE SUPP print #e TWISS, SAVE,betx=midpicell[betx],bety=midpicell[bety],dx=midpicell[dx] PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & STYLE=1, COLOUR=100, SPLINE=.T.,title="SUPP",FILE="DSB3" USE SUP0 print #e TWISS, SAVE,betx=midpicell[betx],bety=midpicell[bety],dx=midpicell[dx] PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & STYLE=1, COLOUR=100, SPLINE=.T.,title="SUP0",FILE="DSB3" USE SUP2 print #e TWISS, SAVE,betx=midpicell[betx],bety=midpicell[bety],dx=midpicell[dx] PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & STYLE=1, COLOUR=100, SPLINE=.T.,title="SUP2",FILE="DSB3" USE, WGC savebeta, LABEL=WGCE, PLACE=#E TWISS, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & VMIN=0, VMAX=50 & STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Wiggler", FILE="DSB3" USE STRSECR print, CLASS=MARKER TWISS, SAVE,betx=midpicell[betx],bety=midpicell[bety],dx=midpicell[dx] PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & STYLE=1, COLOUR=100, SPLINE=.T.,title="STRSECR",FILE="DSB3" USE STRSECI print, class=marker TWISS, SAVE,betx=midpicell[betx],bety=midpicell[bety],dx=midpicell[dx] PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & VMIN=0, VMAX=100, & STYLE=1, COLOUR=100, SPLINE=.T.,title="STRSECI",FILE="DSB3" USE HRINGR print #e TWISS, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & STYLE=1, COLOUR=100, SPLINE=.T.,title="HRINGR",FILE="DSB3" USE HRINGI print #e TWISS, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & STYLE=1, COLOUR=100, SPLINE=.T.,title="HRINGI",FILE="DSB3" ! ENDCOMMENT comment USE, RINGI SURVEY, TAPE= "DSB3.survey.txt" USE, EXT_RTML KICK.EXT1= 0.033e-3 KICK.PEXT1= 43.6E-3/2 KICK.PEXT2= 64.5E-3/2 TWISS, SAVE, BETX=18.8627, BETY=72.1899, X=0.0, PX=0.0 ! PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & ! STYLE=1, COLOUR=100,SPLINE=.T.,title="Extraction",FILE="DSB3" PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=X,Y, & STYLE=1, COLOUR=100, title="Extraction orbit",FILE="DSB3" ! SELECT, OPTICS, CLEAR ! SELECT, OPTICS, CLASS=sbend SELECT, OPTICS OPTICS, BETX=18.8627, BETY=72.1899, X=0.0, PX=0.0, & COLUMNS= KEYWORD,NAME, L, S, & BETX, ALFX, MUX, DX, DPX, & BETY, ALFY, MUY, X, PX, & FILENAME= "DSB3.lattice_EXT.txt" KICK.EXT1= 0.0 KICK.PEXT1= 0.0 KICK.PEXT2= 0.0 endcomment ! ============================================================================== BEAM, PARTICLE=POSITRON, ENERGY=E0, RADIATE ! ============================================================================== ! sextupoles for zero chromaticity ! ============================================================================== COMMENT RF_ON := 1 SXT_ON := 1 USE, RINGRI HARMON HVARY, KSF2, STEP=1E-6, LOWER= 0.0, UPPER=+100.0 HVARY, KSD1, STEP=1E-6, LOWER=-100.0, UPPER= 0.0 HTUNE ENDHARM VALUE, KSF2 VALUE, KSD1 TWISS, SAVE ENDCOMMENT ! ============================================================================== ! Beam Parameters ! ============================================================================== USE RINGRI RF_ON := 1 SXT_ON := 0 TWISS, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, & COLOUR=100, SPLINE=.T., TITLE="Ring ALL",FILE="DSB3" EMIT BMPM ! ============================================================================== ! Plot Chromatic functions ! ============================================================================== USE RINGI RF_ON := 0 SXT_ON := 0 TWISS, CHROM, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=WX,WY, & STYLE=1,COLOUR=100,SPLINE=.T., & TITLE="W without sextupoles",FILE="DSB3_W" PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=DDX, & STYLE=1,COLOUR=100,SPLINE=.T., & TITLE="W without sextupoles",FILE="DSB3_W" RF_ON := 0 SXT_ON := 1 TWISS, CHROM, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=WX,WY, & STYLE=1,COLOUR=100,SPLINE=.T., & TITLE="W with sextupoles",FILE="DSB3_W" PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=DDX, & STYLE=1,COLOUR=100,SPLINE=.T., & TITLE="W with sextupoles",FILE="DSB3_W" !------------------------------------------------------------------------- stop end