option,-echo

!use,eff1s
use,BDS

option,rbarc=.false.

GAMA:= 2.935E6
beam,radiate,energy= 1500.0,ex=.2e-12,ey=.33e-14
!
! historical initial conditions
!
betx0 := 65
alfx0 := 0.
bety0 := 18
alfy0 := 0.
!
!  new initial conditions from Daniel (PLACET)
!  note that the file madtrack.in has not been
!  matched to these new parameters 
!
!betx0 := 8.44404
!alfx0 := 0.62103
!bety0 := 28.4579 
!alfy0 :=  -1.91574
BETX0:=64.171
BETY0:=18.0
ALFX0:=-1.95133
ALFY0:=0.605865

savebeta, label = ipip, place = #e
twiss, betx = betx0, bety = bety0, alfx=alfx0, alfy=alfy0
value, ipip[betx], ipip[alfx], ipip[bety], ipip[alfy]


    MATCH, orbit, betx=betx0,bety=bety0,alfx=alfx0,alfy=alfy0
      weight,betx=1000,bety=1000000,alfx=1000,alfy=1000
!        VARY, qd0[k1], STEP=1.E-5
        VARY, bqf8, STEP=1.E-5
        VARY, bqmd11, STEP=1.E-5
        VARY, bqmd12, STEP=1.E-5
        VARY, bqmd13, STEP=1.E-5
        VARY, bqmd14,   STEP=1.E-5
        WEIGHT, BETX=10, ALFX=0.1, BETY=1000, ALFY=0.1
        CONSTR, #E, ALFX=0, ALFY=0, BETX=0.008, BETY=0.00015
        simplex, TOL=1.0E-15,calls=150000
!       LMDIF, TOL=1.0E-12,calls=150000
        MIGRAD, TOL=1.0E-14,calls=150000
    ENDMATCH


savebeta, label = ipip, place = #e
twiss, betx = betx0, bety = bety0, alfx=alfx0, alfy=alfy0
value, ipip[betx], ipip[alfx], ipip[bety], ipip[alfy]


!
! GAMA:=0.
!
! beam,energy= 250.0,ex=.2e-12,ey=.33e-14
deltaw:=0.008
!
! TRACKING  
!
track,onepass,damp,quantum
!
call, filename="madtrack.in"
RUN, METHOD=TRANSPORT,TABLE=TRMP
TSAVE, FILENAME="mad.trackout"
!ENDTRACK
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

return