1 | title, "V6.5 Thin Lens" ; // USE ONLY WITH MAD-X version > June 2005 |
---|
2 | |
---|
3 | option,-echo,-warn; |
---|
4 | |
---|
5 | set, format="12.6f"; |
---|
6 | select, flag=twiss, clear; |
---|
7 | select, flag=twiss, pattern="IP\.*", column=name,s,betx,bety,alfx,alfy,dx,dpx,mux,muy,x,px,y,py; |
---|
8 | select, flag=twiss, pattern="s.ds.l\.*", column=name,s,betx,bety,alfx,alfy,dx,dpx,mux,muy,x,px,y,py; |
---|
9 | select, flag=twiss, pattern="e.ds.r\.*", column=name,s,betx,bety,alfx,alfy,dx,dpx,mux,muy,x,px,y,py; |
---|
10 | |
---|
11 | call, file="../share/LHC/LHCRing/V6.503/V6.5.seq"; |
---|
12 | call, file="../share/LHC/LHCRing/V6.503/V6.5.inj.str"; |
---|
13 | option,echo,warn,info; |
---|
14 | |
---|
15 | beam, particle = proton, sequence=lhcb1, energy = 450.0; |
---|
16 | |
---|
17 | ! without this we get poor tune (and more?) in twiss table after makethin |
---|
18 | use, sequence=lhcb1; |
---|
19 | twiss, sequence=lhcb1, file=lhc_twiss_b1.before.dat, betx=beta.ip1, bety=beta.ip1, x=+x.ip1, y=+y.ip1, py=+py.ip1; |
---|
20 | |
---|
21 | ! convert to THIN |
---|
22 | option, info, warn, -echo; |
---|
23 | call, file="../share/LHC/LHCRing/V6.503/toolkit/slice.madx"; |
---|
24 | option, echo, warn, info; |
---|
25 | makethin, sequence=lhcb1; |
---|
26 | |
---|
27 | ! set new strength file |
---|
28 | option,-echo,-warn,-info; |
---|
29 | call, file="../share/LHC/LHCRing/V6.503/V6.5.thin.inj.str"; |
---|
30 | option,echo,warn,info; |
---|
31 | |
---|
32 | use, sequence=lhcb1; |
---|
33 | twiss, sequence=lhcb1, file=lhc_twiss_b1.dat, betx=beta.ip1, bety=beta.ip1, x=+x.ip1, y=+y.ip1, py=+py.ip1; |
---|
34 | save,sequence=lhcb1,file="V6.5.thin.b1.seq"; |
---|
35 | |
---|
36 | stop; |
---|