!DCO lattice modified with SuperB cell ;  Marica Biagini for LCWS08
!Total lenght 3238.22 m 
! modified : straight sections from  DCO4
! S. Guiducci December 2009

! Work in progress for new version      2 March 2011  

  ASSIGN, PRINT="dsb3_2.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_2.xsif"

  call, filename="DSB3_2.xsif"

! ==============================================================================
! COMMANDS
! ------------------------------------------------------------------------------

  RF_ON  := 0
  SXT_ON := 0


  USE, PHTCR

 savebeta, LABEL=PHTW, PLACE=#E

 twiss,save

 
 USE, PHTR
 twiss,save
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
      STYLE=1, COLOUR=100, SPLINE=.T.,title="PHTR",FILE="dsb3_2"


  USE, PHTCI
 twiss,save
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
      STYLE=1, COLOUR=100, SPLINE=.T.,title="PHTCI",FILE="dsb3_2"


!---------------------------------------------------------------------------
!   Output kvals
!---------------------------------------------------------------------------

 save,filename="kvals.temp.xsif",pattern="K1..*"

   value THETA 
   value THETPI 
   value THETAS


!--------1---------2---------3---------4---------5---------6---------7---------8
! Output the lattice functions
! ------------------------------------------------------------------------------


  USE CX_CELL_HER
   print #e
   TWISS, SAVE

  USE CXpi_CELL_HER
   savebeta, label=midpicell,place=qfbpi[1]
   print #e
   TWISS, SAVE


  USE CXPI_CX  

  savebeta, label=midpicell,place=qfbpi[1]

   print, CLASS=sextupole
   TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
     STYLE=1, COLOUR=100, SPLINE=.T.,title="CXpi and CX CELLs",FILE="dsb3_2"
hinv: string,"((1+(alfx*alfx))/betx)*dx*dx+2*alfx*dx*dpx+betx*dpx*dpx"
!  plot, table=twiss,haxis=s, vaxis= hinv, style=100,FILE="dsb3_2"



   USE SUP0I
   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="SUP0I",FILE="dsb3_2"


   USE SUP2I
   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="SUP2I",FILE="dsb3_2"



  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_2"



   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_2"


  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_2"

! ==============================================================================
!    Beam Parameters
! ==============================================================================

 RF_ON  := 1
  SXT_ON := 1

!----- RINGI begins from the center of INJ/EXTR section ------------------------------
USE RINGI

print, #e
 TWISS, SAVE

  PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY,                &
      COLOUR=100, SPLINE=.T.,title="Ring from INJ",FILE="dsb3_2"
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,                       & 
        COLOUR=100, SPLINE=.T., TITLE="Ring from INJ",FILE="dsb3_2" 

!---- RINGRI begins from dispersion suppressor SUP0R of  STRSECR
  USE RINGRI  


print, #e
 TWISS, SAVE

  PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, COLOUR=100, &
      SPLINE=.T.,title="Ring from Disp. Sup.",FILE="dsb3_2"


    PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX, COLOUR=100,      & 
       SPLINE=.T., TITLE="Ring from Disp. Sup.",FILE="dsb3_2" 

! ==============================================================================

 BEAM, PARTICLE=POSITRON, ENERGY=E0, RADIATE

   BMPM

! ==============================================================================

stop
end