TITLE "ILC Damping Ring (FODO4)"

  ASSIGN, PRINT="ILCDR.print.txt"
  ASSIGN, ECHO="ILCDR.echo.txt"
  OPTION, -ECHO, DOUBLE, VERIFY

! ==============================================================================
! load LATTICE AND MAGNETS file
! ------------------------------------------------------------------------------

  CALL, FILENAME="FODO4.xsif"
  CALL, FILENAME="k1FODO4-72degree.txt"

!  call "k1FODO4-72degree.txt" for alpha=4*10^-4
!  call "k1FODO4-90degree.txt" for alpha=2*10^-4
!  call "k1FODO4-60degree.txt" for alpha=6*10^-4

! ==============================================================================
! COMMANDs
! ------------------------------------------------------------------------------

  OPTION, ECHO

  BEAM, PARTICLE=ELECTRON, ENERGY=E0


  RF_ON := 0
  USE, RING

  SELECT, OPTICS, CLEAR
  SELECT, OPTICS, #S/#E
  OPTICS, COLUMNS= NAME, KEYWORD, S, L, K0L, E1, E2, &
                   K1L, KS1L, K2L, K3L, TILT, TYPE, &
          FILENAME= "FODO4.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= "FODO4.lfunctions.txt"

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

  USE, ACELL
  RF_ON  := 0
  SXT_ON := 1

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

!  STOP

  USE, RING
  RF_ON  := 1
  SXT_ON := 1

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


  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