1 | !S. Fartoukh Novembre 2009 |
---|
2 | !Rematch the orbit at all IP after the error abd bb implementation |
---|
3 | !S. Fartoukh March 2012 |
---|
4 | ! Use acbh5.l8b1 for IR8 H crossing |
---|
5 | ! Add a SIMPLEX matching command to ease the rematching |
---|
6 | |
---|
7 | limitMCBY := 96.3 *1.e-6*1; ! Existing MCBY in IR2/8 |
---|
8 | limitMCBYY := 96.3 *1.e-6*1; ! 2.25Tm for MCBYY (2.5T *0.9 m) @4.5K in IR1/5 |
---|
9 | |
---|
10 | limitMCBC := 89.87*1.e-6*1; |
---|
11 | limitMCBW := 80.14*1.e-6*1; |
---|
12 | |
---|
13 | if(mylhcbeam==1) |
---|
14 | { |
---|
15 | |
---|
16 | !!IR1 |
---|
17 | match,sequence=lhcb1,range=e.ds.l1.b1/s.ds.r1.b1,betx=100,bety=100,x=0,px=0,y=0,py=0; |
---|
18 | WEIGHT, x = 1, px = 10, y = 1, py = 10; |
---|
19 | constraint,sequence=lhcb1,range=IP1,x=xnom1,y=ynom1,px=pxnom1,py=pynom1; |
---|
20 | constraint,sequence=lhcb1,range=s.ds.r1.b1,x=0.,y=0.,px=0.,py=0.; |
---|
21 | |
---|
22 | vary, name=ACBCH6.L1B1, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
23 | vary, name=ACBCV5.L1B1, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
24 | vary, name=ACBCH5.R1B1, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
25 | vary, name=ACBCV6.R1B1, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
26 | vary, name=ACBYHS4.L1B1, step = 0.00001,lower=-limitMCBYY,upper=limitMCBYY; |
---|
27 | vary, name=ACBYHS4.R1B1, step = 0.00001,lower=-limitMCBYY,upper=limitMCBYY; |
---|
28 | vary, name=ACBYVS4.L1B1, step = 0.00001,lower=-limitMCBYY,upper=limitMCBYY; |
---|
29 | vary, name=ACBYVS4.R1B1, step = 0.00001,lower=-limitMCBYY,upper=limitMCBYY; |
---|
30 | |
---|
31 | simplex, calls = 100, tolerance=1.E-30; |
---|
32 | lmdif, calls = 100, tolerance=1.E-30; |
---|
33 | Endmatch; |
---|
34 | |
---|
35 | !!IR2 |
---|
36 | match,sequence=lhcb1,range=e.ds.l2.b1/s.ds.r2.b1,betx=100,bety=100,x=0,px=0,y=0,py=0; |
---|
37 | WEIGHT, x = 1, px = 10, y = 1, py = 10; |
---|
38 | constraint,sequence=lhcb1,range=IP2,x=xnom2,y=ynom2,px=pxnom2,py=pynom2; |
---|
39 | constraint,sequence=lhcb1,range=s.ds.r2.b1,x=0.,y=0.,px=0.,py=0.; |
---|
40 | |
---|
41 | vary, name=acbyhs5.l2b1, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
42 | vary, name=acbchs5.r2b1, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
43 | vary, name=acbyvs5.l2b1, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
44 | vary, name=acbcvs5.r2b1, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
45 | vary, name=acbyhs4.l2b1, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
46 | vary, name=acbyhs4.r2b1, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
47 | vary, name=acbyvs4.l2b1, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
48 | vary, name=acbyvs4.r2b1, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
49 | |
---|
50 | simplex, calls = 100, tolerance=1.E-30; |
---|
51 | lmdif, calls = 100, tolerance=1.E-30; |
---|
52 | Endmatch; |
---|
53 | |
---|
54 | !!IR5 |
---|
55 | match,sequence=lhcb1,range=e.ds.l5.b1/s.ds.r5.b1,betx=100,bety=100,x=0,px=0,y=0,py=0; |
---|
56 | WEIGHT, x = 1, px = 10, y = 1, py = 10; |
---|
57 | constraint,sequence=lhcb1,range=IP5,x=xnom5,y=ynom5,px=pxnom5,py=pynom5; |
---|
58 | constraint,sequence=lhcb1,range=s.ds.r5.b1,x=0.,y=0.,px=0.,py=0.; |
---|
59 | |
---|
60 | vary, name=ACBCH6.L5B1, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
61 | vary, name=ACBCV5.L5B1, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
62 | vary, name=ACBCH5.R5B1, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
63 | vary, name=ACBCV6.R5B1, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
64 | vary, name=ACBYHS4.L5B1, step = 0.00001,lower=-limitMCBYY,upper=limitMCBYY; |
---|
65 | vary, name=ACBYHS4.R5B1, step = 0.00001,lower=-limitMCBYY,upper=limitMCBYY; |
---|
66 | vary, name=ACBYVS4.L5B1, step = 0.00001,lower=-limitMCBYY,upper=limitMCBYY; |
---|
67 | vary, name=ACBYVS4.R5B1, step = 0.00001,lower=-limitMCBYY,upper=limitMCBYY; |
---|
68 | |
---|
69 | simplex, calls = 100, tolerance=1.E-30; |
---|
70 | lmdif, calls = 100, tolerance=1.E-30; |
---|
71 | Endmatch; |
---|
72 | |
---|
73 | !!IR8 |
---|
74 | match,sequence=lhcb1,range=e.ds.l8.b1/s.ds.r8.b1,betx=100,bety=100,x=0,px=0,y=0,py=0; |
---|
75 | WEIGHT, x = 1, px = 10, y = 1, py = 10; |
---|
76 | constraint,sequence=lhcb1,range=IP8,x=xnom8,y=ynom8,px=pxnom8,py=pynom8; |
---|
77 | constraint,sequence=lhcb1,range=s.ds.r8.b1,x=0.,y=0.,px=0.,py=0.; |
---|
78 | |
---|
79 | !vary, name=acbchs5.l8b1, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; ! temporarily unavailable |
---|
80 | !vary, name=acbwh5.l8b1, step = 0.00001,lower=-limitMCBW,upper=limitMCBW; ! temporary replacement |
---|
81 | vary, name=acbch5.l8b1, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
82 | vary, name=acbyhs4.l8b1, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
83 | vary, name=acbyhs4.r8b1, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
84 | vary, name=acbyhs5.r8b1, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
85 | vary, name=acbcvs5.l8b1, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
86 | vary, name=acbyvs4.l8b1, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
87 | vary, name=acbyvs4.r8b1, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
88 | vary, name=acbyvs5.r8b1, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
89 | |
---|
90 | simplex, calls = 100, tolerance=1.E-30; |
---|
91 | lmdif, calls = 100, tolerance=1.E-30; |
---|
92 | Endmatch; |
---|
93 | |
---|
94 | }; |
---|
95 | |
---|
96 | if(mylhcbeam==2) |
---|
97 | { |
---|
98 | |
---|
99 | !!IR1 |
---|
100 | match,sequence=lhcb2,range=e.ds.l1.b2/s.ds.r1.b2,betx=100,bety=100,x=0,px=0,y=0,py=0; |
---|
101 | WEIGHT, x = 1, px = 10, y = 1, py = 10; |
---|
102 | constraint,sequence=lhcb2,range=IP1,x=xnom1,y=ynom1,px=pxnom1,py=pynom1; |
---|
103 | constraint,sequence=lhcb2,range=s.ds.r1.b2,x=0.,y=0.,px=0.,py=0.; |
---|
104 | |
---|
105 | vary, name=ACBCH6.R1b2, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
106 | vary, name=ACBCV5.R1b2, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
107 | vary, name=ACBCH5.L1b2, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
108 | vary, name=ACBCV6.L1b2, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
109 | vary, name=ACBYHS4.L1b2, step = 0.00001,lower=-limitMCBYY,upper=limitMCBYY; |
---|
110 | vary, name=ACBYHS4.R1b2, step = 0.00001,lower=-limitMCBYY,upper=limitMCBYY; |
---|
111 | vary, name=ACBYVS4.L1b2, step = 0.00001,lower=-limitMCBYY,upper=limitMCBYY; |
---|
112 | vary, name=ACBYVS4.R1b2, step = 0.00001,lower=-limitMCBYY,upper=limitMCBYY; |
---|
113 | |
---|
114 | simplex, calls = 100, tolerance=1.E-30; |
---|
115 | lmdif, calls = 2000, tolerance=1.E-30; |
---|
116 | Endmatch; |
---|
117 | |
---|
118 | !!IR2 |
---|
119 | match,sequence=lhcb2,range=e.ds.l2.b2/s.ds.r2.b2,betx=100,bety=100,x=0,px=0,y=0,py=0; |
---|
120 | WEIGHT, x = 1, px = 10, y = 1, py = 10; |
---|
121 | constraint,sequence=lhcb2,range=IP2,x=xnom2,y=ynom2,px=pxnom2,py=pynom2; |
---|
122 | constraint,sequence=lhcb2,range=s.ds.r2.b2,x=0.,y=0.,px=0.,py=0.; |
---|
123 | |
---|
124 | vary, name=acbyhs5.l2b2, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
125 | vary, name=acbchs5.r2b2, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
126 | vary, name=acbyvs5.l2b2, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
127 | vary, name=acbcvs5.r2b2, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
128 | vary, name=acbyhs4.l2b2, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
129 | vary, name=acbyhs4.r2b2, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
130 | vary, name=acbyvs4.l2b2, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
131 | vary, name=acbyvs4.r2b2, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
132 | |
---|
133 | simplex, calls = 100, tolerance=1.E-30; |
---|
134 | lmdif, calls = 2000, tolerance=1.E-30; |
---|
135 | Endmatch; |
---|
136 | |
---|
137 | !!IR5 |
---|
138 | match,sequence=lhcb2,range=e.ds.l5.b2/s.ds.r5.b2,betx=100,bety=100,x=0,px=0,y=0,py=0; |
---|
139 | WEIGHT, x = 1, px = 10, y = 1, py = 10; |
---|
140 | constraint,sequence=lhcb2,range=IP5,x=xnom5,y=ynom5,px=pxnom5,py=pynom5; |
---|
141 | constraint,sequence=lhcb2,range=s.ds.r5.b2,x=0.,y=0.,px=0.,py=0.; |
---|
142 | |
---|
143 | vary, name=ACBCH6.R5b2, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
144 | vary, name=ACBCV5.R5b2, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
145 | vary, name=ACBCH5.L5b2, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
146 | vary, name=ACBCV6.L5b2, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
147 | vary, name=ACBYHS4.L5b2, step = 0.00001,lower=-limitMCBYY,upper=limitMCBYY; |
---|
148 | vary, name=ACBYHS4.R5b2, step = 0.00001,lower=-limitMCBYY,upper=limitMCBYY; |
---|
149 | vary, name=ACBYVS4.L5b2, step = 0.00001,lower=-limitMCBYY,upper=limitMCBYY; |
---|
150 | vary, name=ACBYVS4.R5b2, step = 0.00001,lower=-limitMCBYY,upper=limitMCBYY; |
---|
151 | |
---|
152 | simplex, calls = 100, tolerance=1.E-30; |
---|
153 | lmdif, calls = 2000, tolerance=1.E-30; |
---|
154 | Endmatch; |
---|
155 | |
---|
156 | !!IR8 |
---|
157 | match,sequence=lhcb2,betx=100,bety=100,x=0,px=0,y=0,py=0; |
---|
158 | WEIGHT, x = 1, px = 10, y = 1, py = 10; |
---|
159 | constraint,sequence=lhcb2,range=IP8,x=xnom8,y=ynom8,px=pxnom8,py=pynom8; |
---|
160 | constraint,sequence=lhcb2,range=s.ds.r8.b2,x=0.,y=0.,px=0.,py=0.; |
---|
161 | |
---|
162 | vary, name=acbchs5.l8b2, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
163 | vary, name=acbyhs5.r8b2, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
164 | vary, name=acbcvs5.l8b2, step = 0.00001,lower=-limitMCBC,upper=limitMCBC; |
---|
165 | vary, name=acbyvs5.r8b2, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
166 | vary, name=acbyhs4.l8b2, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
167 | vary, name=acbyhs4.r8b2, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
168 | vary, name=acbyvs4.l8b2, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
169 | vary, name=acbyvs4.r8b2, step = 0.00001,lower=-limitMCBY,upper=limitMCBY; |
---|
170 | |
---|
171 | simplex, calls = 100, tolerance=1.E-30; |
---|
172 | lmdif, calls = 2000, tolerance=1.E-30; |
---|
173 | Endmatch; |
---|
174 | |
---|
175 | }; |
---|
176 | |
---|
177 | return; |
---|
178 | |
---|
179 | |
---|