TITLE, "LHC Version 6.503 -  July 2009" ;

option, -echo, -warn, -info;

real const l.tanal = 0;
real const l.tanar = 0;
real const l.tanc = 0;
real const l.tcdqa = 0;
real const l.tcp   = 0.0;

real const l.mbxwt = 0.0;
real const l.mbls2 = 0.0;
real const l.tasb  = 0.0;
real const l.btvss075 = 0.0;
real const l.mblw = 0.0;

real const l.tctvb = 0;
real const l.x5zdc002 = 0;
real const l.tcddm = 0;

call, file="../share/LHC/LHCRing/V6.503/V6.5.seq";
call, file="../share/LHC/LHCRing/V6.503/V6.5.inj.str";

beam, sequence=lhcb1, particle=proton, pc=450;
beam, sequence=lhcb2, particle=proton, pc=450, bv=-1;

on_x1:=1;
on_sep1:=1;
on_atlas:=0;
on_x2:=1;
on_sep2:=1;
on_alice:=0;
on_x5:=1;
on_sep5:=1;
on_cms  :=0;
on_x8:=1;
on_sep8:=1;
on_lhcb :=0;

option, warn, info;

use, sequence=lhcb1;
select, flag=twiss, clear;
select,flag=twiss,column=name,s,beta11,beta21,beta12,beta22,disp1,disp3,x,px,y,py;
twiss, sequence=lhcb1, file="twiss-table.out";

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

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

ptc_end;