Changeset 9097 for trunk/ippScripts/scripts/detrend_reject_exp.pl
- Timestamp:
- Oct 2, 2006, 12:16:03 PM (20 years ago)
- File:
-
- 1 edited
-
trunk/ippScripts/scripts/detrend_reject_exp.pl (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/ippScripts/scripts/detrend_reject_exp.pl
r9091 r9097 73 73 } 74 74 75 my @exp Ids; # Array of exposure IDs75 my @expTags; # Array of exposure IDs 76 76 my @means; # Array of means 77 77 my @stdevs; # Array of stdevs 78 78 my @meanStdevs; # Array of mean stdevs 79 79 foreach my $exposure (@$exposures) { 80 die "Unable to find exposure id.\n" if not defined $exposure->{exp_ id};80 die "Unable to find exposure id.\n" if not defined $exposure->{exp_tag}; 81 81 die "Unable to find mean.\n" if not defined $exposure->{bg}; 82 82 die "Unable to find stdev.\n" if not defined $exposure->{bg_stdev}; 83 83 die "Unable to find mean stdev.\n" if not defined $exposure->{bg_mean_stdev}; 84 push @exp Ids, $exposure->{exp_id};84 push @expTags, $exposure->{exp_tag}; 85 85 push @means, $exposure->{bg}; 86 86 push @stdevs, $exposure->{bg_stdev}; … … 105 105 my $numRejected = 0; # Number of exposures rejected 106 106 for (my $i = 0; $i < scalar @means; $i++) { 107 my $exp Id = $expIds[$i]; # Exposure ID108 my $command = "$dettool -updateresidexp -det_id $det_id -iteration $iter -exp_ id $expId"; # Command to run107 my $expTag = $expTags[$i]; # Exposure ID 108 my $command = "$dettool -updateresidexp -det_id $det_id -iteration $iter -exp_tag $expTag"; # Command to run 109 109 my $reject = 0; # Reject this exposure? 110 110 if (defined REJECT_MEAN->{$det_type} and … … 112 112 $meanStats->standard_deviation() > 0 and 113 113 ($means[$i] - $meanStats->mean()) / $meanStats->standard_deviation() > REJECT_MEAN->{$det_type}) { 114 print "Rejecting $exp Idbased on bad mean: " .114 print "Rejecting $expTag based on bad mean: " . 115 115 (($means[$i] - $meanStats->mean()) / $meanStats->standard_deviation()) . 116 116 " vs " . REJECT_MEAN->{$det_type} . "\n"; … … 121 121 ($stdevs[$i] - $stdevStats->mean()) / $stdevStats->standard_deviation() > 122 122 REJECT_STDEV->{$det_type}) { 123 print "Rejecting $exp Idbased on bad stdev: " .123 print "Rejecting $expTag based on bad stdev: " . 124 124 (($stdevs[$i] - $stdevStats->mean()) / $stdevStats->standard_deviation()) . 125 125 " vs " . REJECT_STDEV->{$det_type} . "\n"; … … 130 130 ($meanStdevs[$i] - $meanStdevStats->mean()) / $meanStdevStats->standard_deviation() > 131 131 REJECT_MEAN_STDEV->{$det_type}) { 132 print "Rejecting $exp Idbased on bad stdev: " .132 print "Rejecting $expTag based on bad stdev: " . 133 133 (($meanStdevs[$i] - $meanStdevStats->mean()) / $meanStdevStats->standard_deviation()) . 134 134 " vs " . REJECT_MEAN_STDEV->{$det_type} . "\n";
Note:
See TracChangeset
for help on using the changeset viewer.
