source: PSPA/madxPSPA/tests/test-ptc-twiss-old6/Efcompb1_MB.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: 11.2 KB
Line 
1 
2NewFamily_MB: macro = {
3 b1ss  = 0;a1ss  = 0;
4 b2ss  = 0;a2ss  = 0;
5 b3ss  = 0;a3ss  = 0;
6 b4ss  = 0;a4ss  = 0;
7 b5ss  = 0;a5ss  = 0;
8 b6ss  = 0;a6ss  = 0;
9 b7ss  = 0;a7ss  = 0;
10 b8ss  = 0;a8ss  = 0;
11 b9ss  = 0;a9ss  = 0;
12 b10ss = 0;a10ss = 0;
13 b11ss = 0;a11ss = 0;
14 }
15 
16Outer_MB: macro = {
17 B1s   =  ON_MB*(ON_PERS*( b1M_inj_MB_V1 + ON_DEC * b1M_dec_MB)+
18             (1-ON_PERS)*( b1M_col_MB_V1 + ON_BDOT* b1M_ramp_MB));
19 B2s   =  ON_MB*(ON_PERS*( b2M_inj_MB_V1 + ON_DEC * b2M_dec_MB)+
20             (1-ON_PERS)*( b2M_col_MB_V1 + ON_BDOT* b2M_ramp_MB));
21 B3s   =  ON_MB*(ON_PERS*( b3M_inj_MB_V1 + ON_DEC * b3M_dec_MB)+
22             (1-ON_PERS)*( b3M_col_MB_V1 + ON_BDOT* b3M_ramp_MB));
23 B4s   =  ON_MB*(ON_PERS*( b4M_inj_MB_V1 + ON_DEC * b4M_dec_MB)+
24             (1-ON_PERS)*( b4M_col_MB_V1 + ON_BDOT* b4M_ramp_MB));
25 B5s   =  ON_MB*(ON_PERS*( b5M_inj_MB_V1 + ON_DEC * b5M_dec_MB)+
26             (1-ON_PERS)*( b5M_col_MB_V1 + ON_BDOT* b5M_ramp_MB));
27 B6s   =  ON_MB*(ON_PERS*( b6M_inj_MB_V1 + ON_DEC * b6M_dec_MB)+
28             (1-ON_PERS)*( b6M_col_MB_V1 + ON_BDOT* b6M_ramp_MB));
29 B7s   =  ON_MB*(ON_PERS*( b7M_inj_MB_V1 + ON_DEC * b7M_dec_MB)+
30             (1-ON_PERS)*( b7M_col_MB_V1 + ON_BDOT* b7M_ramp_MB));
31 B8s   =  ON_MB*(ON_PERS*( b8M_inj_MB_V1 + ON_DEC * b8M_dec_MB)+
32             (1-ON_PERS)*( b8M_col_MB_V1 + ON_BDOT* b8M_ramp_MB));
33 B9s   =  ON_MB*(ON_PERS*( b9M_inj_MB_V1 + ON_DEC * b9M_dec_MB)+
34             (1-ON_PERS)*( b9M_col_MB_V1 + ON_BDOT* b9M_ramp_MB));
35B10s   =  ON_MB*(ON_PERS*(b10M_inj_MB_V1 + ON_DEC *b10M_dec_MB)+
36             (1-ON_PERS)*(b10M_col_MB_V1 + ON_BDOT*b10M_ramp_MB));
37B11s   =  ON_MB*(ON_PERS*(b11M_inj_MB_V1 + ON_DEC *b11M_dec_MB)+
38             (1-ON_PERS)*(b11M_col_MB_V1 + ON_BDOT*b11M_ramp_MB));
39 
40 A1s   =  ON_MB*(ON_PERS*( a1M_inj_MB_V1 + ON_DEC * a1M_dec_MB)+
41             (1-ON_PERS)*( a1M_col_MB_V1 + ON_BDOT* a1M_ramp_MB));
42 A2s   =  ON_MB*(ON_PERS*( a2M_inj_MB_V1 + ON_DEC * a2M_dec_MB)+
43             (1-ON_PERS)*( a2M_col_MB_V1 + ON_BDOT* a2M_ramp_MB));
44 A3s   =  ON_MB*(ON_PERS*( a3M_inj_MB_V1 + ON_DEC * a3M_dec_MB)+
45             (1-ON_PERS)*( a3M_col_MB_V1 + ON_BDOT* a3M_ramp_MB));
46 A4s   =  ON_MB*(ON_PERS*( a4M_inj_MB_V1 + ON_DEC * a4M_dec_MB)+
47             (1-ON_PERS)*( a4M_col_MB_V1 + ON_BDOT* a4M_ramp_MB));
48 A5s   =  ON_MB*(ON_PERS*( a5M_inj_MB_V1 + ON_DEC * a5M_dec_MB)+
49             (1-ON_PERS)*( a5M_col_MB_V1 + ON_BDOT* a5M_ramp_MB));
50 A6s   =  ON_MB*(ON_PERS*( a6M_inj_MB_V1 + ON_DEC * a6M_dec_MB)+
51             (1-ON_PERS)*( a6M_col_MB_V1 + ON_BDOT* a6M_ramp_MB));
52 A7s   =  ON_MB*(ON_PERS*( a7M_inj_MB_V1 + ON_DEC * a7M_dec_MB)+
53             (1-ON_PERS)*( a7M_col_MB_V1 + ON_BDOT* a7M_ramp_MB));
54 A8s   =  ON_MB*(ON_PERS*( a8M_inj_MB_V1 + ON_DEC * a8M_dec_MB)+
55             (1-ON_PERS)*( a8M_col_MB_V1 + ON_BDOT* a8M_ramp_MB));
56 A9s   =  ON_MB*(ON_PERS*( a9M_inj_MB_V1 + ON_DEC * a9M_dec_MB)+
57             (1-ON_PERS)*( a9M_col_MB_V1 + ON_BDOT* a9M_ramp_MB));
58A10s   =  ON_MB*(ON_PERS*(a10M_inj_MB_V1 + ON_DEC *a10M_dec_MB)+
59             (1-ON_PERS)*(a10M_col_MB_V1 + ON_BDOT*a10M_ramp_MB));
60A11s   =  ON_MB*(ON_PERS*(a11M_inj_MB_V1 + ON_DEC *a11M_dec_MB)+
61             (1-ON_PERS)*(a11M_col_MB_V1 + ON_BDOT*a11M_ramp_MB));
62}
63 
64Inner_MB: macro = {
65 B1s   =  ON_MB*(ON_PERS*( b1M_inj_MB_V2 + ON_DEC * b1M_dec_MB)+
66             (1-ON_PERS)*( b1M_col_MB_V2 + ON_BDOT* b1M_ramp_MB));
67 B2s   =  ON_MB*(ON_PERS*( b2M_inj_MB_V2 + ON_DEC * b2M_dec_MB)+
68             (1-ON_PERS)*( b2M_col_MB_V2 + ON_BDOT* b2M_ramp_MB));
69 B3s   =  ON_MB*(ON_PERS*( b3M_inj_MB_V2 + ON_DEC * b3M_dec_MB)+
70             (1-ON_PERS)*( b3M_col_MB_V2 + ON_BDOT* b3M_ramp_MB));
71 B4s   =  ON_MB*(ON_PERS*( b4M_inj_MB_V2 + ON_DEC * b4M_dec_MB)+
72             (1-ON_PERS)*( b4M_col_MB_V2 + ON_BDOT* b4M_ramp_MB));
73 B5s   =  ON_MB*(ON_PERS*( b5M_inj_MB_V2 + ON_DEC * b5M_dec_MB)+
74             (1-ON_PERS)*( b5M_col_MB_V2 + ON_BDOT* b5M_ramp_MB));
75 B6s   =  ON_MB*(ON_PERS*( b6M_inj_MB_V2 + ON_DEC * b6M_dec_MB)+
76             (1-ON_PERS)*( b6M_col_MB_V2 + ON_BDOT* b6M_ramp_MB));
77 B7s   =  ON_MB*(ON_PERS*( b7M_inj_MB_V2 + ON_DEC * b7M_dec_MB)+
78             (1-ON_PERS)*( b7M_col_MB_V2 + ON_BDOT* b7M_ramp_MB));
79 B8s   =  ON_MB*(ON_PERS*( b8M_inj_MB_V2 + ON_DEC * b8M_dec_MB)+
80             (1-ON_PERS)*( b8M_col_MB_V2 + ON_BDOT* b8M_ramp_MB));
81 B9s   =  ON_MB*(ON_PERS*( b9M_inj_MB_V2 + ON_DEC * b9M_dec_MB)+
82             (1-ON_PERS)*( b9M_col_MB_V2 + ON_BDOT* b9M_ramp_MB));
83B10s   =  ON_MB*(ON_PERS*(b10M_inj_MB_V2 + ON_DEC *b10M_dec_MB)+
84             (1-ON_PERS)*(b10M_col_MB_V2 + ON_BDOT*b10M_ramp_MB));
85B11s   =  ON_MB*(ON_PERS*(b11M_inj_MB_V2 + ON_DEC *b11M_dec_MB)+
86             (1-ON_PERS)*(b11M_col_MB_V2 + ON_BDOT*b11M_ramp_MB));
87 
88 A1s   =  ON_MB*(ON_PERS*( a1M_inj_MB_V2 + ON_DEC * a1M_dec_MB)+
89             (1-ON_PERS)*( a1M_col_MB_V2 + ON_BDOT* a1M_ramp_MB));
90 A2s   =  ON_MB*(ON_PERS*( a2M_inj_MB_V2 + ON_DEC * a2M_dec_MB)+
91             (1-ON_PERS)*( a2M_col_MB_V2 + ON_BDOT* a2M_ramp_MB));
92 A3s   =  ON_MB*(ON_PERS*( a3M_inj_MB_V2 + ON_DEC * a3M_dec_MB)+
93             (1-ON_PERS)*( a3M_col_MB_V2 + ON_BDOT* a3M_ramp_MB));
94 A4s   =  ON_MB*(ON_PERS*( a4M_inj_MB_V2 + ON_DEC * a4M_dec_MB)+
95             (1-ON_PERS)*( a4M_col_MB_V2 + ON_BDOT* a4M_ramp_MB));
96 A5s   =  ON_MB*(ON_PERS*( a5M_inj_MB_V2 + ON_DEC * a5M_dec_MB)+
97             (1-ON_PERS)*( a5M_col_MB_V2 + ON_BDOT* a5M_ramp_MB));
98 A6s   =  ON_MB*(ON_PERS*( a6M_inj_MB_V2 + ON_DEC * a6M_dec_MB)+
99             (1-ON_PERS)*( a6M_col_MB_V2 + ON_BDOT* a6M_ramp_MB));
100 A7s   =  ON_MB*(ON_PERS*( a7M_inj_MB_V2 + ON_DEC * a7M_dec_MB)+
101             (1-ON_PERS)*( a7M_col_MB_V2 + ON_BDOT* a7M_ramp_MB));
102 A8s   =  ON_MB*(ON_PERS*( a8M_inj_MB_V2 + ON_DEC * a8M_dec_MB)+
103             (1-ON_PERS)*( a8M_col_MB_V2 + ON_BDOT* a8M_ramp_MB));
104 A9s   =  ON_MB*(ON_PERS*( a9M_inj_MB_V2 + ON_DEC * a9M_dec_MB)+
105             (1-ON_PERS)*( a9M_col_MB_V2 + ON_BDOT* a9M_ramp_MB));
106A10s   =  ON_MB*(ON_PERS*(a10M_inj_MB_V2 + ON_DEC *a10M_dec_MB)+
107             (1-ON_PERS)*(a10M_col_MB_V2 + ON_BDOT*a10M_ramp_MB));
108A11s   =  ON_MB*(ON_PERS*(a11M_inj_MB_V2 + ON_DEC *a11M_dec_MB)+
109             (1-ON_PERS)*(a11M_col_MB_V2 + ON_BDOT*a11M_ramp_MB));
110}
111 
112 
113 R_ALL =  3.0;    !default cut of the distribution (RANDOM);
114 GCUTR =  R_ALL;
115 GCUTS =  1.5;    !default cut of the distribution (SYSTEMATIC);
116 Rr    =  0.017;
117 
118 ON_MB:=1;
119 
120!****************Main and Disp. Suppr. Dipoles (MB)*******************
121!  System,"echo ' Efcomp : Main and Disp. Suppr. Dipoles (MB)'";
122! Define sigma for random componentSet, bscale, 10E-4;
123 B1r  =  ON_MB*(ON_PERS*sqrt( b1R_inj_MB^2+ ON_DEC^2* b1R_dec_MB^2)+
124            (1-ON_PERS)*sqrt( b1R_col_MB^2+ON_BDOT^2* b1R_ramp_MB^2));
125 B2r  =  ON_MB*(ON_PERS*sqrt( b2R_inj_MB^2+ ON_DEC^2* b2R_dec_MB^2)+
126            (1-ON_PERS)*sqrt( b2R_col_MB^2+ON_BDOT^2* b2R_ramp_MB^2));
127 B3r  =  ON_MB*(ON_PERS*sqrt( b3R_inj_MB^2+ ON_DEC^2* b3R_dec_MB^2)+
128            (1-ON_PERS)*sqrt( b3R_col_MB^2+ON_BDOT^2* b3R_ramp_MB^2));
129 B4r  =  ON_MB*(ON_PERS*sqrt( b4R_inj_MB^2+ ON_DEC^2* b4R_dec_MB^2)+
130            (1-ON_PERS)*sqrt( b4R_col_MB^2+ON_BDOT^2* b4R_ramp_MB^2));
131 B5r  =  ON_MB*(ON_PERS*sqrt( b5R_inj_MB^2+ ON_DEC^2* b5R_dec_MB^2)+
132            (1-ON_PERS)*sqrt( b5R_col_MB^2+ON_BDOT^2* b5R_ramp_MB^2));
133 B6r  =  ON_MB*(ON_PERS*sqrt( b6R_inj_MB^2+ ON_DEC^2* b6R_dec_MB^2)+
134            (1-ON_PERS)*sqrt( b6R_col_MB^2+ON_BDOT^2* b6R_ramp_MB^2));
135 B7r  =  ON_MB*(ON_PERS*sqrt( b7R_inj_MB^2+ ON_DEC^2* b7R_dec_MB^2)+
136            (1-ON_PERS)*sqrt( b7R_col_MB^2+ON_BDOT^2* b7R_ramp_MB^2));
137 B8r  =  ON_MB*(ON_PERS*sqrt( b8R_inj_MB^2+ ON_DEC^2* b8R_dec_MB^2)+
138            (1-ON_PERS)*sqrt( b8R_col_MB^2+ON_BDOT^2* b8R_ramp_MB^2));
139 B9r  =  ON_MB*(ON_PERS*sqrt( b9R_inj_MB^2+ ON_DEC^2* b9R_dec_MB^2)+
140            (1-ON_PERS)*sqrt( b9R_col_MB^2+ON_BDOT^2* b9R_ramp_MB^2));
141B10r  =  ON_MB*(ON_PERS*sqrt(b10R_inj_MB^2+ ON_DEC^2*b10R_dec_MB^2)+
142            (1-ON_PERS)*sqrt(b10R_col_MB^2+ON_BDOT^2*b10R_ramp_MB^2));
143B11r  =  ON_MB*(ON_PERS*sqrt(b11R_inj_MB^2+ ON_DEC^2*b11R_dec_MB^2)+
144            (1-ON_PERS)*sqrt(b11R_col_MB^2+ON_BDOT^2*b11R_ramp_MB^2));
145 
146 A1r  =  ON_MB*(ON_PERS*sqrt( a1R_inj_MB^2+ ON_DEC^2* a1R_dec_MB^2)+
147            (1-ON_PERS)*sqrt( a1R_col_MB^2+ON_BDOT^2* a1R_ramp_MB^2));
148 A2r  =  ON_MB*(ON_PERS*sqrt( a2R_inj_MB^2+ ON_DEC^2* a2R_dec_MB^2)+
149            (1-ON_PERS)*sqrt( a2R_col_MB^2+ON_BDOT^2* a2R_ramp_MB^2));
150 A3r  =  ON_MB*(ON_PERS*sqrt( a3R_inj_MB^2+ ON_DEC^2* a3R_dec_MB^2)+
151            (1-ON_PERS)*sqrt( a3R_col_MB^2+ON_BDOT^2* a3R_ramp_MB^2));
152 A4r  =  ON_MB*(ON_PERS*sqrt( a4R_inj_MB^2+ ON_DEC^2* a4R_dec_MB^2)+
153            (1-ON_PERS)*sqrt( a4R_col_MB^2+ON_BDOT^2* a4R_ramp_MB^2));
154 A5r  =  ON_MB*(ON_PERS*sqrt( a5R_inj_MB^2+ ON_DEC^2* a5R_dec_MB^2)+
155            (1-ON_PERS)*sqrt( a5R_col_MB^2+ON_BDOT^2* a5R_ramp_MB^2));
156 A6r  =  ON_MB*(ON_PERS*sqrt( a6R_inj_MB^2+ ON_DEC^2* a6R_dec_MB^2)+
157            (1-ON_PERS)*sqrt( a6R_col_MB^2+ON_BDOT^2* a6R_ramp_MB^2));
158 A7r  =  ON_MB*(ON_PERS*sqrt( a7R_inj_MB^2+ ON_DEC^2* a7R_dec_MB^2)+
159            (1-ON_PERS)*sqrt( a7R_col_MB^2+ON_BDOT^2* a7R_ramp_MB^2));
160 A8r  =  ON_MB*(ON_PERS*sqrt( a8R_inj_MB^2+ ON_DEC^2* a8R_dec_MB^2)+
161            (1-ON_PERS)*sqrt( a8R_col_MB^2+ON_BDOT^2* a8R_ramp_MB^2));
162 A9r  =  ON_MB*(ON_PERS*sqrt( a9R_inj_MB^2+ ON_DEC^2* a9R_dec_MB^2)+
163            (1-ON_PERS)*sqrt( a9R_col_MB^2+ON_BDOT^2* a9R_ramp_MB^2));
164A10r  =  ON_MB*(ON_PERS*sqrt(a10R_inj_MB^2+ ON_DEC^2*a10R_dec_MB^2)+
165            (1-ON_PERS)*sqrt(a10R_col_MB^2+ON_BDOT^2*a10R_ramp_MB^2));
166A11r  =  ON_MB*(ON_PERS*sqrt(a11R_inj_MB^2+ ON_DEC^2*a11R_dec_MB^2)+
167            (1-ON_PERS)*sqrt(a11R_col_MB^2+ON_BDOT^2*a11R_ramp_MB^2));
168 
169! Define systematic component per arc and set errors for MB's
170   select, flag=error, clear; select, flag=error,  range= IP1/IP2,        class= MB;
171      eoption, seed= SEEDSYS    ; exec NewFamily_MB; exec Outer_MB;
172      eoption, seed= SEEDRAN    ; exec SetEfcomp_B;
173 
174   select, flag=error, clear; select, flag=error,  range= IP2/IP3.B1,     class= MB;
175      eoption, seed= SEEDSYS+100; exec NewFamily_MB; exec Inner_MB;
176      eoption, seed= SEEDRAN+1100; exec SetEfcomp_B;
177 
178   select, flag=error, clear; select, flag=error,  range= IP3.B1/IP4.B1,  class= MB;
179      eoption, seed= SEEDSYS+200; exec NewFamily_MB; exec Inner_MB;
180      eoption, seed= SEEDRAN+1200; exec SetEfcomp_B;
181 
182   select, flag=error, clear; select, flag=error,  range= IP4.B1/IP5,     class= MB;
183      eoption, seed= SEEDSYS+300; exec NewFamily_MB; exec Inner_MB;
184      eoption, seed= SEEDRAN+1300; exec SetEfcomp_B;
185 
186   select, flag=error, clear; select, flag=error,  range= IP5/IP6.B1,     class= MB;
187      eoption, seed= SEEDSYS+400; exec NewFamily_MB; exec Outer_MB;
188      eoption, seed= SEEDRAN+1400; exec SetEfcomp_B;
189 
190   select, flag=error, clear; select, flag=error,  range= IP6.B1/IP7.B1,  class= MB;
191      eoption, seed= SEEDSYS+500; exec NewFamily_MB; exec Outer_MB;
192      eoption, seed= SEEDRAN+1500; exec SetEfcomp_B;
193 
194   select, flag=error, clear; select, flag=error,  range= IP7.B1/IP8,     class= MB;
195      eoption, seed= SEEDSYS+600; exec NewFamily_MB; exec Outer_MB;
196      eoption, seed= SEEDRAN+1600; exec SetEfcomp_B;
197 
198   select, flag=error, clear; select, flag=error,  range= IP8/IP1.L1,     class= MB;
199      eoption, seed= SEEDSYS+700; exec NewFamily_MB; exec Inner_MB;
200      eoption, seed= SEEDRAN+1700; exec SetEfcomp_B;
201 
202 
203 
Note: See TracBrowser for help on using the repository browser.