!call,"newnewcolla4.mad"
call,"ff_3tev.mad"
!call,"precolla4.mad"

!shortffnl:line=(totalmis,eff1s)

!shortffnlp:line=(precollima,totalmis)

!shortffnlpp:line=(precollima,totalmis,eff1s)

!beam, particle=electron, energy=1500

srad:= 1.698398e-12


!use, eff1s
!savebeta, label=bqff, place=mff
!twiss, betx=65., bety=18., save
!value, bqff[betx], bqff[bety]

lb4:=lsca*lffb4/2.0/nbfac
lb3:=lsca*lffb3/2.0/nbfac
lb3b:=lsca*lffb3b/2.0/nbfac
lb2:=lsca*lffb2/2.0/nbfac
lb1:=lsca*lffb1/2.0/nbfac
lbs1:=lsca*lffsb1/2.0/nbfac
lbs1a:=lsca*lffsb1a/2.0/nbfac
lbs1b:=lsca*lffsb1b/2.0/nbfac
KFFQ0:=1.0/lsca/lsca*kqd0

comment
!use, eff1
use,eff1s 
print, full
TWISS,SAVE,betx=65.,bety=18.,alfx=0.,alfy=0.
!
setplot ascale=1 lscale=1.5 sscale=1.5 rscale=1.5 xsize=28 ysize=20 &
        post=2 font=-4
!
PLOT,HAXIS=S,VAXIS1=betX,betY, VAXIS2=DX,DY,&
!          hmin=280. hmax=400. vmin=0.,0. vmax=10000.,.1 &
         DELTAP=0.,STYLE=100
endcomment
comment
!==========================================================
!Matrix elements calculation
!==========================================================
use, eff1s
match,betx=65.,bety=18.,&
      alfx=0., alfy=0.
vary,KSF6,step=0.001
vary,KSF5,step=0.001
!vary,KSD4,step=0.001
!vary,KSF1,step=0.001
!vary,KSD0,step=0.001
!-----------------------------------------------------
rmat, #s/#e, rm(1,2)=1., weight(1,2)=0.0
tmatrix, #s/#e, tm(1,1,1)=1., weight(1,1,1)=0.0
!-----------------------------------------------------
Simplex,toler=1.0E-12
lmdif,toler=1.0E-12
endmatch
endcomment

comment
use, eff1s

nbetax:string, "betx/0.008"
nbetay:string, "bety/0.00015"
sizex:string, "sqrt(betx/0.008)"
sizey:string, "sqrt(bety/0.00015)"
delsig:string, &
"1+(1-sqrt(1-1.698398e-12*bety*bety))"

twiss, betx=65.,bety=18.,alfx=0.,alfy=0.,&
       save=IPtwiss, &
       deltap=-0.006:0.006:0.001

table,name=IPtwiss,column=deltap,nbetax,nbetay,sizex,sizey,delsig
archive, table=IPtwiss, filename='betIP.dat'

plot, title='betaIP vs deltap', &
      style=100, symbol=4, &
      haxis=deltap, vaxis1=nbetax,nbetay,range=#e

plot, title='delsig vs deltap', &
      style=100, symbol=4, &
      haxis=deltap, vaxis1=delsig,range=mff

plot, title='alphasIP vs deltap', &
      style=100, symbol=4, &
      haxis=deltap, vaxis1=alfx,alfy, range=#e  
endcomment

value, THETAB2,THETASB1,THETASB1A,THETASB1B
value, lb4,lb3,lb3b,lb2,lb1,lbs1,lbs1a,lbs1b
value,KFFQ0
 
!comment
!use, totalmis
!use, shortffnlpp
use, eff1s
GAMA:= 2.935E6
beam,particle=electron,energy=1500.,&
     ex=.23e-12,ey=.68e-14,radiate
!,radiate
!deltaw:=0.008

track,onepass,damp,quantum
call, filename="madtrack.in"
observe, place=#e, table=colima
run, method=transport,table=colima
tsave, filename="madtrack.out"
!setplot ascale=1 lscale=1.5 sscale=1.5 rscale=1.5 xsize=23 ysize=17 &
!        post=2 font=-4
!plot,table=colima,haxis=x, vaxis=px, multiple
endtrack
!endcomment
!return

stop