source: PSPA/madxPSPA/tests/share/LHC/LHCRing/SLHCV3.1b/toolkit/rematchCOarc.madx @ 430

Last change on this file since 430 was 430, checked in by touze, 11 years ago

import madx-5.01.00
File size: 4.7 KB
Line 
1!Rematch the orbit in the arcs, S. Fartoukh July 2010
2!Add simplex matching command to ease correction, S. Fartoukh March 2012
3
4kmcb_max=80.e-6;
5
6if(mylhcbeam==1)
7{
8!!IR1-left
9match,sequence=lhcb1,range=e.ds.r8.b1/e.ds.l1.b1,betx=100,bety=100,x=0,px=0,y=0,py=0;
10WEIGHT,     x = 1, px =  10, y = 1, py = 10;
11constraint,sequence=lhcb1,range=e.ds.l1.b1,x=0,y=0,px=0,py=0;
12vary, name=acbh14.l1b1,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
13vary, name=acbh12.l1b1,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
14vary, name=acbv15.l1b1,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
15vary, name=acbv13.l1b1,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
16simplex, calls = 100, tolerance=1.E-30;
17lmdif, calls = 100, tolerance=1.E-30;
18Endmatch;
19!!IR1-right
20match,sequence=lhcb1,range=s.ds.r1.b1/s.ds.r2.b1,betx=100,bety=100,x=0,px=0,y=0,py=0;
21WEIGHT,     x = 1, px =  10, y = 1, py = 10;
22constraint,sequence=lhcb1,range=s.ds.r2.b1,x=0,y=0,px=0,py=0;
23vary, name=acbh13.r1b1,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
24vary, name=acbh15.r1b1,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
25vary, name=acbv12.r1b1,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
26vary, name=acbv14.r1b1,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
27simplex, calls = 100, tolerance=1.E-30;
28lmdif, calls = 100, tolerance=1.E-30;
29Endmatch;
30!!IR5-left
31match,sequence=lhcb1,range=e.ds.r4.b1/e.ds.l5.b1,betx=100,bety=100,x=0,px=0,y=0,py=0;
32WEIGHT,     x = 1, px =  10, y = 1, py = 10;
33constraint,sequence=lhcb1,range=e.ds.l5.b1,x=0,y=0,px=0,py=0;
34vary, name=acbh14.l5b1,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
35vary, name=acbh12.l5b1,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
36vary, name=acbv15.l5b1,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
37vary, name=acbv13.l5b1,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
38simplex, calls = 100, tolerance=1.E-30;
39lmdif, calls = 100, tolerance=1.E-30;
40Endmatch;
41!!IR5-right
42match,sequence=lhcb1,range=s.ds.r5.b1/s.ds.r6.b1,betx=100,bety=100,x=0,px=0,y=0,py=0;
43WEIGHT,     x = 1, px =  10, y = 1, py = 10;
44constraint,sequence=lhcb1,range=s.ds.r6.b1,x=0,y=0,px=0,py=0;
45vary, name=acbh13.r5b1,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
46vary, name=acbh15.r5b1,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
47vary, name=acbv12.r5b1,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
48vary, name=acbv14.r5b1,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
49simplex, calls = 100, tolerance=1.E-30;
50lmdif, calls = 100, tolerance=1.E-30;
51Endmatch;
52};
53
54if(mylhcbeam==2)
55{
56!!IR1-left
57match,sequence=lhcb2,range=e.ds.r8.b2/e.ds.l1.b2,betx=100,bety=100,x=0,px=0,y=0,py=0;
58WEIGHT,     x = 1, px =  10, y = 1, py = 10;
59constraint,sequence=lhcb2,range=e.ds.l1.b2,x=0,y=0,px=0,py=0;
60vary, name=acbh13.l1b2,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
61vary, name=acbh15.l1b2,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
62vary, name=acbv12.l1b2,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
63vary, name=acbv14.l1b2,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
64simplex, calls = 100, tolerance=1.E-30;
65lmdif, calls = 100, tolerance=1.E-30;
66Endmatch;
67!!IR1-right
68match,sequence=lhcb2,range=s.ds.r1.b2/s.ds.r2.b2,betx=100,bety=100,x=0,px=0,y=0,py=0;
69WEIGHT,     x = 1, px =  10, y = 1, py = 10;
70constraint,sequence=lhcb2,range=s.ds.r2.b2,x=0,y=0,px=0,py=0;
71vary, name=acbh12.r1b2,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
72vary, name=acbh14.r1b2,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
73vary, name=acbv13.r1b2,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
74vary, name=acbv15.r1b2,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
75simplex, calls = 100, tolerance=1.E-30;
76lmdif, calls = 100, tolerance=1.E-30;
77Endmatch;
78!!IR5-left
79match,sequence=lhcb2,range=e.ds.r4.b2/e.ds.l5.b2,betx=100,bety=100,x=0,px=0,y=0,py=0;
80WEIGHT,     x = 1, px =  10, y = 1, py = 10;
81constraint,sequence=lhcb2,range=e.ds.l5.b2,x=0,y=0,px=0,py=0;
82vary, name=acbh13.l5b2,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
83vary, name=acbh15.l5b2,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
84vary, name=acbv12.l5b2,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
85vary, name=acbv14.l5b2,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
86simplex, calls = 100, tolerance=1.E-30;
87lmdif, calls = 100, tolerance=1.E-30;
88Endmatch;
89!!IR5-right
90match,sequence=lhcb2,range=s.ds.r5.b2/s.ds.r6.b2,betx=100,bety=100,x=0,px=0,y=0,py=0;
91WEIGHT,     x = 1, px =  10, y = 1, py = 10;
92constraint,sequence=lhcb2,range=s.ds.r6.b2,x=0,y=0,px=0,py=0;
93vary, name=acbh12.r5b2,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
94vary, name=acbh14.r5b2,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
95vary, name=acbv13.r5b2,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
96vary, name=acbv15.r5b2,  step = 0.00001,lower=-kmcb_max,upper=kmcb_max;
97simplex, calls = 100, tolerance=1.E-30;
98lmdif, calls = 100, tolerance=1.E-30;
99Endmatch;
100};
101
102
103return;
104
105
Note: See TracBrowser for help on using the repository browser.