!--------1---------2---------3---------4---------5---------6---------7--------8
  TITLE "ILC Damping Ring (OCS8) 20 June 2006"
!--------1---------2---------3---------4---------5---------6---------7--------8

  ASSIGN, PRINT="OCS8.print.txt"
  ASSIGN, ECHO="OCS8.echo.txt"

  SETPLOT, XSIZE=25.4, YSIZE=20.32
  SETPLOT, LWIDTH=5, LSCALE=1.5,SSCALE=1.5,RSCALE=1.5 

!--------1---------2---------3---------4---------5---------6---------7--------8

  OPTION, -ECHO, DOUBLE, VERIFY

  CALL, FILENAME="OCS8.xsif"

  OPTION, ECHO

!--------1---------2---------3---------4---------5---------6---------7--------8

COMMENT
  USE, RING
  RF_ON  := 0
  SXT_ON := 1

  SELECT, ERROR, CLEAR
  SELECT, ERROR, CLASS=QUADRUPOLE

  EOPT, SEED=123456789
  EALIGN, DY=1e-6*TGAUSS(5)

  SELECT, OPTICS, CLEAR
  SELECT, OPTICS, #S/#E
  OPTICS, COLUMNS= S, L, X, PX, Y, PY, &
          FILENAME= "OCS8.closedorbit.txt"

ENDCOMMENT

!--------1---------2---------3---------4---------5---------6---------7--------8

  BEAM, PARTICLE=ELECTRON, ENERGY=5.0

  RF_ON := 0

  USE, RING

  SELECT, OPTICS, CLEAR
  SELECT, OPTICS, CLASS=DRIFT
  OPTICS, COLUMNS= NAME, KEYWORD, S, L, K0L, &
          FILENAME= "OCS8.layout.txt"

  SELECT, OPTICS, CLEAR
  SELECT, OPTICS, #S/#E
  OPTICS, COLUMNS= NAME, KEYWORD, S, L, K0L, E1, E2, &
                   K1L, KS1L, K2L, K3L, TILT, TYPE, &
          FILENAME= "OCS8.lattice.txt"

  SELECT, OPTICS, CLEAR
  SELECT, OPTICS, #S/#E
  OPTICS, COLUMNS= NAME, S,    L,            &
                   BETX, ALFX, MUX, DX, DPX, &
                   BETY, ALFY, MUY,          &
                   K0L,  E1, E2, K1L,  K2L   &
          FILENAME= "OCS8.lfunctions.txt"

!--------1---------2---------3---------4---------5---------6---------7--------8

  RF_ON  := 0
  SXT_ON := 1

  USE, ACELL
  PRINT, #S/#E
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        STYLE=1, COLOUR=100, SPLINE=.T., FILE="OCS8.ACELL.ps"

  HX: STRING,"((1+(alfx*alfx))/betx)*dx*dx+2*alfx*dx*dpx+betx*dpx*dpx"

  plot, table=twiss, haxis=s, vaxis1=hx, &
!        vmin=0, vmax=0.01, &
     spline=.F., colour=100, style=1, range=#s/#e

  USE, CWIG
  PRINT, #S/#E
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        STYLE=1, COLOUR=100, SPLINE=.T., FILE="OCS8.CWIG.ps"

  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=HINV, &
        STYLE=1, COLOUR=100, SPLINE=.T., FILE="OCS8.CWIG.Hfn.ps"

  USE, SECTORWIG
  PRINT, #S/#E
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        STYLE=1, COLOUR=100, SPLINE=.T., FILE="OCS8.SECTORW.ps"

  USE, RFCELL
  PRINT, #S/#E
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        STYLE=1, COLOUR=100, SPLINE=.T., FILE="OCS8.RFCELL.ps"

  SXT_ON := 0

  USE, RING
  PRINT, #E
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        STYLE=1, COLOUR=100, SPLINE=.T., FILE="OCS8.RING.ps"

stop

!--------1---------2---------3---------4---------5---------6---------7--------8

!  TWISS, TUNES, DELTAP=-0.01:0.01:0.001, SAVE
!  PLOT, HAXIS=DELTAP, VAXIS1=QX, VAXIS2=QY, &
!                      SPLINE, COLOUR=100, STYLE=1

!--------1---------2---------3---------4---------5---------6---------7--------8

  ASSIGN, PRINT="OCS8.bmpm.txt"

  USE, RING

  RF_ON  := 0
  SXT_ON := 0
  PRINT, #E
  TWISS, TUNES, SAVE

  RF_ON  := 1
  SXT_ON := 1
  PRINT, #E
  TWISS, TUNES, SAVE

! Note: BMPM does not work when skew quadrupoles are present in lattice

  BMPM

  BEAM, PARTICLE=electron, ENERGY=E0, NPART=2.0E10, RADIATE
  EMIT

  STOP