IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
May 10, 2011, 4:43:58 PM (15 years ago)
Author:
watersc1
Message:

missed cases for quality flagging.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/ippScripts/scripts/lap_science.pl

    r31511 r31515  
    475475        }
    476476        # Do quality checks here
    477         if ((($exposure->{chipRun_state} eq 'full')&&($exposure->{chip_bad_quality} / $exposure->{chip_component_count} > 0.05)) ||
    478             (($exposure->{camRun_state} eq 'full')&&($exposure->{cam_bad_quality} / $exposure->{cam_component_count} > 0)) ||
    479             (($exposure->{warpRun_state} eq 'full')&&($exposure->{warp_bad_quality} / $exposure->{warp_component_count} > 0.2)) ||
    480             (($exposure->{diffRun_state} eq 'full')&&($exposure->{diff_bad_quality} / $exposure->{diff_component_count} > 0.5)))  {
    481             printf( "QUALITY: $exposure->{exp_id} has bad quality fractions: chip: %d/%d cam: %d/%d warp: %d/%d diff: %d/%d\n",
    482                     $exposure->{chip_bad_quality} , $exposure->{chip_component_count},
    483                     $exposure->{cam_bad_quality} , $exposure->{cam_component_count},
    484                     $exposure->{warp_bad_quality} , $exposure->{warp_component_count},
     477        my $is_bad_quality = 0;
     478        if ((defined($exposure->{chipRun_state}))&&($exposure->{chipRun_state} eq 'full')&&
     479            ($exposure->{chip_component_count} > 0)&&($exposure->{chip_bad_quality} / $exposure->{chip_component_count} > 0.05)) {
     480            printf("QUALITY: $exposure->{exp_id} has bad chip quality: %d / %d\n",
     481                   $exposure->{chip_bad_quality} , $exposure->{chip_component_count});
     482            $is_bad_quality = 1;
     483        }
     484        elsif ((defined($exposure->{camRun_state}))&&($exposure->{camRun_state} eq 'full')&&
     485               ($exposure->{cam_bad_quality} / $exposure->{cam_component_count} > 0)) {
     486            printf("QUALITY: $exposure->{exp_id} has bad cam quality: %d / %d\n",
     487                   $exposure->{cam_bad_quality} , $exposure->{cam_component_count});
     488            $is_bad_quality = 1;
     489        }
     490        elsif ((defined($exposure->{warpRun_state}))&&($exposure->{warpRun_state} eq 'full')&&
     491               ($exposure->{warp_bad_quality} / $exposure->{warp_component_count} > 0.2)) {
     492            printf("QUALITY: $exposure->{exp_id} has bad warp quality: %d / %d\n",
     493                   $exposure->{warp_bad_quality} , $exposure->{warp_component_count});
     494            $is_bad_quality = 1;
     495        }
     496        elsif ((defined($exposure->{diffRun_state}))&&($exposure->{diffRun_state} eq 'full')&&
     497               ($exposure->{diff_bad_quality} / $exposure->{diff_component_count} > 0.5)) {
     498            printf("QUALITY: $exposure->{exp_id} has bad diff quality: %d / %d\n",
    485499                    $exposure->{diff_bad_quality} , $exposure->{diff_component_count});
    486 
    487             unless ($exposure->{diffRun_state} eq 'full') {
     500            $is_bad_quality = 1;
     501        }
     502        if ($is_bad_quality) {
     503            unless ((defined($exposure->{diffRun_state}))&&
     504                    ($exposure->{diffRun_state} eq 'full')) {
    488505                $needs_qstack = 1;
    489506                $needs_something_private = 1;
Note: See TracChangeset for help on using the changeset viewer.