1 | !---|----1----|----2----|----3----|----4----|----5----|----6----|----7----|---8; |
---|
2 | Option, warn,info; |
---|
3 | |
---|
4 | Title, "Matching chromaticity for beam1 and beam2 of the LHC 6.4 lattice with MADX"; |
---|
5 | |
---|
6 | |
---|
7 | !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++; |
---|
8 | ! TARGET CHROMATICITY ; |
---|
9 | !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++; |
---|
10 | QPRIME = 2.0; |
---|
11 | |
---|
12 | !---------------------------------------------------------------; |
---|
13 | ! CALL IN SELECTION, SEQUENCE and OPTICS FILES; |
---|
14 | !---------------------------------------------------------------; |
---|
15 | |
---|
16 | option, -echo,-warn; |
---|
17 | call, file = "../share/LHC/LHCRing/V6.503/V6.5.seq"; |
---|
18 | call, file = "../share/LHC/LHCRing/V6.503/V6.5.coll.str"; |
---|
19 | option, echo; |
---|
20 | |
---|
21 | Beam, particle = proton, sequence=lhcb1, energy = 450.0, |
---|
22 | NPART=1.05E11, sige= 4.5e-4 ; |
---|
23 | Beam, particle = proton, sequence=lhcb2, energy = 450.0, |
---|
24 | bv = -1, NPART=1.05E11, sige= 4.5e-4; !<<<< |
---|
25 | ! the 'bv' parameter indicates that beam2 rotates counter-clockwise |
---|
26 | |
---|
27 | ! Select beam1: |
---|
28 | USE, period=lhcb1; |
---|
29 | |
---|
30 | ! Generate a TWISS table: |
---|
31 | select, flag=twiss, clear; |
---|
32 | select, flag=twiss, pattern="IP\.*",column=name,s,betx,bety; |
---|
33 | TWISS,sequence=lhcb1,file=twiss.b1.data; |
---|
34 | |
---|
35 | ! Select beam2: |
---|
36 | USE, period=lhcb2; |
---|
37 | |
---|
38 | ! Generate a TWISS table: |
---|
39 | select, flag=twiss, clear; |
---|
40 | select, flag=twiss, pattern="IP\.*",column=name,s,betx,bety; |
---|
41 | TWISS,sequence=lhcb2,file=twiss.b2.data; |
---|
42 | |
---|
43 | ! Match the chromaticity: |
---|
44 | match,sequence=lhcb1,lhcb2; |
---|
45 | vary, name=KSF.B1, step=0.00001; |
---|
46 | vary, name=KSD.B1, step=0.00001; |
---|
47 | vary, name=KSF.B2, step=0.00001; |
---|
48 | vary, name=KSD.B2, step=0.00001; |
---|
49 | global,sequence=lhcb1,DQ1=QPRIME; |
---|
50 | global,sequence=lhcb1,DQ2=QPRIME; |
---|
51 | global,sequence=lhcb2,DQ1=QPRIME; |
---|
52 | global,sequence=lhcb2,DQ2=QPRIME; |
---|
53 | Lmdif, calls=10, tolerance=1.0e-21; |
---|
54 | endmatch; |
---|
55 | |
---|
56 | ! Generate a TWISS table: |
---|
57 | select, flag=twiss, clear; |
---|
58 | select, flag=twiss, pattern="IP\.*",column=name,s,betx,bety; |
---|
59 | TWISS,sequence=lhcb1,file=twiss.b1.data; |
---|
60 | |
---|
61 | ! Generate a TWISS table: |
---|
62 | select, flag=twiss, clear; |
---|
63 | select, flag=twiss, pattern="IP\.*",column=name,s,betx,bety; |
---|
64 | TWISS,sequence=lhcb2,file=twiss.b2.data; |
---|
65 | |
---|
66 | stop; |
---|
67 | |
---|