!TITLE,'Test input for PTC_TWISS ';

!=========== RING PARAMETERS ===========================================
option, -echo, -warn, -info;

eg   :=  100;
bg   :=  eg/pmass;
en   := 3.75e-06;
epsx := en/bg;
epsy := en/bg;

beam, particle = proton, energy =   eg        , 
                         sigt=      0.077     , 
                         sige=      1.1e-4    , 
                         npart=     1.05e11   , 
                         exn=4*en,  eyn=4*en  , 
                         kbunch = 10,
                         et = 0.002, bv = 1,
                         ex=epsx,   ey=epsy;

call file="input/fv9.opt";
call file="input/fv9.seq";

option, warn, info;

use,period=fivecell;
// alignment errors

select,flag=error,clear;
select,flag=error,pattern="q.*",range=qf.1;
ealign,dx=0.0001;

acbh1=1e-6;
acbh2=1e-6;
acbh3=1e-6;
acbh4=1e-6;
acbh5=1e-6;

! At the end of the ring the data_block is produced
savebeta, label=twssip, place=#e,sequence=fivecell;

select,flag=twiss,column=name,s,betx,bety,dx,dy;
twiss, file="twiss-table.out"; 

show,twssip;     ! show parameters (see "*.out" file)

! Initialize PTC
ptc_create_universe;
ptc_create_layout,model=2,method=6,nst=10,exact;
!ptc_align;

select,flag=ptc_twiss,column=name,s,beta11,beta21,beta12,beta22,disp1,disp3,x,px,y,py;
ptc_twiss ,closed_orbit,                icase=5,no=3,file="ptc-twiss-table.out";
ptc_normal,closed_orbit,maptable,normal,icase=5,no=3;
ptc_end;
stop;