!--------1---------2---------3---------4---------5---------6---------7---------8
! WIG (wiggler magnet)
!--------1---------2---------3---------4---------5---------6---------7---------8

! Drifts

  DWIG2:  DRIFT, L=0.1
  DWIG3:  DRIFT, L=0.025


! Bends

!  WIGBMAX := 1.95  ! Peak wiggler field (T)
!  WIGPERD := 0.32  ! Wiggler period (m)
!  WANGLE  := WIGBMAX*(WIGPERD/4)*CLIGHT/E0/1e9  ! Bending angle in one pole

!  BWIGP:  SBEND, L=WIGPERD/4, ANGLE=+WANGLE,   E1=+WANGLE/2, E2=+WANGLE/2
!  BWIGM:  SBEND, L=WIGPERD/4, ANGLE=-WANGLE,   E1=-WANGLE/2, E2=-WANGLE/2
!  BWIGPH: SBEND, L=WIGPERD/8, ANGLE=+WANGLE/2, E1=+WANGLE/4, E2=+WANGLE/4
!  BWIGMH: SBEND, L=WIGPERD/8, ANGLE=-WANGLE/2, E1=-WANGLE/4, E2=-WANGLE/4

!  wang_0 := 0
!  BWIGP_0:  SBEND, L=WIGPERD/4, ANGLE=+WANG_0,   E1=+WANG_0/2, E2=+WANG_0/2
!  BWIGM_0:  SBEND, L=WIGPERD/4, ANGLE=-WANG_0,   E1=-WANG_0/2, E2=-WANG_0/2
!  BWIGPH_0: SBEND, L=WIGPERD/8, ANGLE=+WANG_0/2, E1=+WANG_0/4, E2=+WANG_0/4
!  BWIGMH_0: SBEND, L=WIGPERD/8, ANGLE=-WANG_0/2, E1=-WANG_0/4, E2=-WANG_0/4


! Beamlines

!  WIGPRD: LINE=(BWIGM,DWIG2,BWIGP,DWIG2)

!  WIG:    LINE=(BWIGPH,DWIG3,5*WIGPRD,BWIGM,DWIG2,BWIGP,DWIG3,BWIGMH)
!  WIG_off: LINE=(BWIGPH_0,DWIG3,&
!              5*WIGPRD,BWIGM_0,DWIG2,BWIGP_0,DWIG3,BWIGMH_0)


 WIG: wiggler, n_pole = 12, l_pole=0.16, l=1.92
! wig[b_max]:=1.5
 WIG_off: wiggler, b_max=0., n_pole = 12, l_pole=0.16, l=1.92

  wig1: marker
  wig2: marker
  wig_all: line=(wig1, wig, wig2)

!--------1---------2---------3---------4---------5---------6---------7---------8
! WGC (wiggler cell)
! try to match to
!  beta x      beta y       eta x     alpha x     alpha y
! 14.2899      8.2846      0.0000     -1.2457      0.7401
! 14.2899      8.2846      0.0000     -1.2457      0.7401

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

! Drifts
  lwig_dlr = 0.024 !0.03
  DWGC2A:  DRIFT, L=0.75-LCOR-LDCOR + lwig_dlr -0.5
  DWGC1A:  DRIFT, L=0.75-LDBPM + lwig_dlr      +0.5
  DWGC2B:  DRIFT, L=0.75-LCOR-LDCOR + lwig_dlr +0.5
  DWGC1B:  DRIFT, L=0.75-LDBPM + lwig_dlr      -0.5

! Quadrupoles
k1.qwgc1_off :=  0.58115561
k1.qwgc2_off := -0.55312809

  QWGC1:  QUADRUPOLE, L=QLENGTH !, K1=K1.QWGC1
  QWGC2:  QUADRUPOLE, L=QLENGTH !, K1=K1.QWGC2

QWGC1[K1] =      5.8154065897E-01
QWGC2[K1] =     -5.4978824668E-01

  QWGC1_off:  QUADRUPOLE, L=QLENGTH, K1=K1.QWGC1_off
  QWGC2_off:  QUADRUPOLE, L=QLENGTH, K1=K1.QWGC2_off


! Beamline

  WGC:    LINE=( CQB(QWGC1), DWGC1A, &
                 WIG,        DWGC2A, &
                 CQB(QWGC2), DWGC1B, &
                -WIG,        DWGC2B )

  WGC_off:    LINE=( CQB(QWGC1_off), DWGC1, &
                 WIG_off,        DWGC2, &
                 CQB(QWGC2_off), DWGC1, &
                -WIG_off,        DWGC2 )