!**********triplet lattice*****************
call "mta_diag_10_7_4.lat"
BEAM, PARTICLE="PROTON",CHARGE=1

use,diagnos
print,range=diagnos
!twiss functions at middle of Q74

!twiss,save, betx=3.773, alfx=1.445, bety=9.575, alfy=-3.445
!twiss,save, betx=4.914, alfx=1.738, bety=7.278, alfy=-2.963
twiss,save, betx=3.691, alfx=0.403, bety=9.767, alfy=-0.822

plot,table=twiss,colour=100,haxis=s,vaxis1=betx,bety,&
vaxis2=dx
plot,table=twiss,colour=100,haxis=s,vaxis1=betx,bety,&
vaxis2=dy

print, range=diagnos

!survey

!survey,x0=30331.68916,z0=30368.79956,y0=225.552451,theta0=3.141680967
!!survey,z0=30331.68915,x0=30368.76781,y0=225.552451,theta0=4.712477293
!survey,x0=30331.68917,z0=30368.86306,y0=225.552451,theta0=3.141680967
!!!survey,z0=30331.688936,x0=30368.766650,y0=225.546093,theta0=4.712477293

!stop

!!match,betx=4,914,alfx=1.738,bety=7.278,alfy=-2.963
!match,betx=3.773,alfx=1.445,bety=9.575,alfy=-3.445
match, betx=3.691, alfx=0.403, bety=9.767, alfy=-0.822

vary uq01[K1], step=.01, lower =-2, upper=2
vary uq02[K1], step=.01, lower =-2, upper=2
vary uq03[K1], step=.01, lower =-2, upper=2
vary uq04[K1], step=.01, lower =-2, upper=2
vary uq05[K1], step=.01, lower =-3, upper=4
vary uq06[K1], step=.01, lower =-3, upper=4
vary uq07[K1], step=.01, lower =-3, upper=4
vary uq08[K1], step=.01, lower =-3, upper=4
vary uq09[K1], step=.01, lower =-3, upper=4
!vary QM11[K1], step=.01, lower =-3, upper=4
!!vary T1, step=.01,lower=0,upper=0.18

constraint, diagnos, betx<180,bety<180
constraint, #E, dy=0,dpy=0
constraint, #E, dx=0,dpx=0
!!constraint, #E,betx=2,alfx=0,bety=2,alfy=0
!constraint, #E,betx>2,bety>1,alfx=0,alfy=0
!constraint, diagnos, betx<30,bety<30
constraint, UMW5_C,betx=0.5,alfx=0,bety=0.5,alfy=0

LMDIF, calls=1000, tole=1E-08
endmatch
!twiss,save, betx=3.773, alfx=1.445, bety=9.575, alfy=-3.445
!twiss,save, betx=4.914, alfx=1.738, bety=7.278, alfy=-2.963
twiss,save, betx=3.691, alfx=0.403, bety=9.767, alfy=-0.822

plot,table=twiss,colour=100,haxis=s,vaxis1=betx,bety,&
vaxis2=dx
plot,table=twiss,colour=100,haxis=s,vaxis1=betx,bety,&
vaxis2=dy

stop