[430] | 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 | |
---|