1 | |
---|
2 | NewFamily_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 | |
---|
16 | Outer_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)); |
---|
35 | B10s = 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)); |
---|
37 | B11s = 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)); |
---|
58 | A10s = 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)); |
---|
60 | A11s = 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 | |
---|
64 | Inner_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)); |
---|
83 | B10s = 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)); |
---|
85 | B11s = 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)); |
---|
106 | A10s = 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)); |
---|
108 | A11s = 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)); |
---|
141 | B10r = 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)); |
---|
143 | B11r = 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)); |
---|
164 | A10r = 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)); |
---|
166 | A11r = 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 | |
---|