Changeset 12677
- Timestamp:
- Mar 29, 2007, 11:10:42 AM (19 years ago)
- Location:
- trunk/ippScripts/scripts
- Files:
-
- 2 edited
-
detrend_reject_exp.pl (modified) (3 diffs)
-
detrend_reject_imfile.pl (modified) (8 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/ippScripts/scripts/detrend_reject_exp.pl
r12676 r12677 169 169 my $dMean = abs($means[$i] - $meanStats->mean()) ; 170 170 if ($dMean > ($reject_mean * $meanStats->standard_deviation())) { 171 print $logFile "Rejecting $expTag based on ensemble mean value: " . 172 "$means[$i] --> $dMean vs " . $reject_mean * $meanStats->standard_deviation() . "\n"; 171 print $logFile "Rejecting $expTag based on ensemble mean value: "; 173 172 $reject = 1; 174 173 #goto UPDATE; 175 } 174 } else { 175 print $logFile "Ensemble mean meets requirements: "; 176 } 177 print $logFile "$means[$i] --> $dMean vs " . $reject_mean * $meanStats->standard_deviation() . "\n"; 176 178 } else { 177 179 print $logFile "No rejection for ensemble mean\n"; … … 181 183 my $dMean = abs($stdevs[$i] - $stdevStats->mean()); 182 184 if ($dMean > ($reject_stdev * $stdevStats->standard_deviation())) { 183 print $logFile "Rejecting $expTag based on ensemble stdev: " . 184 $stdevs[$i] . " --> $dMean sigma vs " . $reject_stdev . "\n"; 185 print $logFile "Rejecting $expTag based on ensemble stdev: "; 185 186 $reject = 1; 186 187 #goto UPDATE; 187 } 188 } else { 189 print $logFile "Ensemble stdev meets requirements: "; 190 } 191 print $logFile "$stdevs[$i] --> $dMean sigma vs " . $reject_stdev . "\n"; 188 192 } else { 189 193 print $logFile "No rejection for ensemble stdev\n"; … … 193 197 my $dMean = abs($meanStdevs[$i] - $meanStdevStats->mean()); 194 198 if ($dMean > ($reject_meanstdev * $meanStdevStats->standard_deviation())) { 195 print $logFile "Rejecting $expTag based on ensemble mean stdev: " . 196 $meanStdevs[$i] . " --> $dMean sigma vs " . $reject_meanstdev . "\n"; 199 print $logFile "Rejecting $expTag based on ensemble mean stdev: "; 197 200 $reject = 1; 198 201 #goto UPDATE; 199 } 202 } else { 203 print $logFile "Ensemble mean stdev meets requirements: "; 204 } 205 print $logFile "$meanStdevs[$i] --> $dMean sigma vs " . $reject_meanstdev . "\n"; 200 206 } else { 201 207 print $logFile "No rejection for ensemble mean stdev\n"; -
trunk/ippScripts/scripts/detrend_reject_imfile.pl
r12675 r12677 193 193 if ($reject_imfile_mean) { 194 194 if (abs($mean) > $reject_imfile_mean * $stdev) { 195 print $logFile "Rejecting exposure based on bad imfile mean for $name: " . 196 $mean . " vs " . $reject_imfile_mean . "\n"; 197 $reject = 1; 198 } 195 print $logFile "Rejecting exposure based on bad imfile mean for $name: "; 196 $reject = 1; 197 } else { 198 print $logFile "Imfile mean for $name meets requirements: "; 199 } 200 print $logFile "$mean vs $reject_imfile_mean\n"; 199 201 } else { 200 202 print $logFile "No rejection on imfile mean for $name\n"; … … 202 204 if ($reject_imfile_stdev) { 203 205 if ($stdev > $reject_imfile_stdev) { 204 print $logFile "Rejecting exposure based on bad imfile stdev for $name: " . 205 $stdev . " vs " . $reject_imfile_stdev . "\n"; 206 $reject = 1; 207 } 206 print $logFile "Rejecting exposure based on bad imfile stdev for $name: "; 207 $reject = 1; 208 } else { 209 print $logFile "Imfile stdev for $name meets requirements: "; 210 } 211 print "$stdev vs $reject_imfile_stdev\n"; 212 208 213 } else { 209 214 print $logFile "No rejection on imfile stdev for $name\n"; … … 211 216 if ($reject_imfile_meanstdev) { 212 217 if ($meanStdevs[$i] > $reject_imfile_meanstdev) { 213 print $logFile "Rejecting exposure based on bad imfile mean stdev for $name: " . 214 $meanStdevs[$i] . " vs " . $reject_imfile_meanstdev . "\n"; 215 $reject = 1; 216 } 218 print $logFile "Rejecting exposure based on bad imfile mean stdev for $name: "; 219 $reject = 1; 220 } else { 221 print $logFile "Imfile mean stdev for $name meets requirements: "; 222 } 223 print "$meanStdevs[$i] vs $reject_imfile_meanstdev\n"; 217 224 } else { 218 225 print $logFile "No rejection on imfile mean stdev for $name\n"; … … 221 228 if ($mean < $stdev * $reject_imfile_sn) { 222 229 print $logFile "Rejecting exposure based on bad imfile S/N for $name: "; 223 print $logFile "mean: " . $mean . " vs " . "stdev*SNlimit: " . $stdev * $reject_imfile_sn . "\n"; 224 $reject = 1; 225 } 230 $reject = 1; 231 } else { 232 print $logFile "Imfile S/N for $name meets requirements: "; 233 } 234 print $logFile "mean: $mean vs stdev*SNlimit: " . $stdev * $reject_imfile_sn . "\n"; 226 235 } else { 227 236 print $logFile "No rejection on imfile S/N for $name\n"; … … 254 263 if ($reject_exp_mean) { 255 264 if (abs($mean) > $reject_exp_mean * $stdev) { 256 print $logFile "Rejecting exposure based on bad mean: \n"; 257 print $logFile "mean: $mean, stdev: $stdev (limit is: % $reject_exp_mean\n"; 258 $reject = 1; 259 } 265 print $logFile "Rejecting exposure based on bad mean: "; 266 $reject = 1; 267 } else { 268 print $logFile "Exposure mean meets requirements: "; 269 } 270 print $logFile "mean: $mean, stdev: $stdev (limit is: % $reject_exp_mean\n"; 260 271 } else { 261 272 print $logFile "No rejection for exp mean\n"; … … 264 275 if ($reject_exp_stdev) { 265 276 if ($stdev > $reject_exp_stdev) { 266 print $logFile "Rejecting exposure based on bad stdev: " . $stdev . " vs " . 267 $reject_exp_stdev . "\n"; 268 $reject = 1; 269 } 277 print $logFile "Rejecting exposure based on bad stdev: "; 278 $reject = 1; 279 } else { 280 print $logFile "Exposure stdev meets requirements: "; 281 } 282 print $logFile "$stdev vs $reject_exp_stdev\n"; 270 283 } else { 271 284 print $logFile "No rejection for exp stdev\n"; … … 274 287 if ($reject_exp_meanstdev) { 275 288 if ($meanStdev > $reject_exp_meanstdev) { 276 print $logFile "Rejecting exposure based on bad mean stdev: " . $meanStdev . " vs " . 277 $reject_exp_meanstdev . "\n"; 278 $reject = 1; 279 } 289 print $logFile "Rejecting exposure based on bad mean stdev: "; 290 $reject = 1; 291 } else { 292 print $logFile "Exposure mean stdev meets requirements: "; 293 } 294 print $logFile "$meanStdev vs $reject_exp_meanstdev\n"; 280 295 } else { 281 296 print $logFile "No rejection for exp mean stdev\n"; … … 285 300 if (abs($mean) < abs($stdev * $reject_exp_sn)) { 286 301 print $logFile "Rejecting exposure based on poor S/N: \n"; 287 print $logFile "signal: $mean vs noise: $stdev (s/n limit is: $reject_exp_sn)\n"; 288 $reject = 1; 289 } 302 $reject = 1; 303 } else { 304 print $logFile "Exposure S/N meets requirements: \n"; 305 } 306 print $logFile "signal: $mean vs noise: $stdev (s/n limit is: $reject_exp_sn)\n"; 290 307 } else { 291 308 print $logFile "No rejection for exp S/N\n";
Note:
See TracChangeset
for help on using the changeset viewer.
