IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Aug 17, 2009, 4:45:20 PM (17 years ago)
Author:
Paul Price
Message:

MOPS needs to know the direction of subtraction. Added a new field in the IPP-MOPS format: POSITIVE, a boolean which indicates whether the subtraction is positive or not.

File:
1 edited

Legend:

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

    r25074 r25100  
    8383my %exp_id;                     # Exposure identifiers
    8484my %exp_name;                   # Exposure names
     85my %direction;                  # Direction of subtraction
    8586{
    8687    my $command;                # Command to run
     
    114115        if ($stage eq 'diff') {
    115116            my $skycell_id = $comp->{skycell_id};
     117
     118            # Positive direction
    116119            $files{"$skycell_id.pos"} = $ipprc->filename( "PPSUB.OUTPUT.SOURCES", $path_base );
    117             $files{"$skycell_id.neg"} = $ipprc->filename( "PPSUB.INVERSE.SOURCES", $path_base ) if
    118                 defined $comp->{bothways} and $comp->{bothways};
    119120            $zp{"$skycell_id.pos"} = $zp;
    120             $zp{"$skycell_id.neg"} = $zp if defined $comp->{bothways} and $comp->{bothways};
    121121            $exp_id{"$skycell_id.pos"} = $comp->{exp_id_1};
    122             $exp_id{"$skycell_id.neg"} = $comp->{exp_id_2} if defined $comp->{bothways} and $comp->{bothways};
    123122            $exp_name{"$skycell_id.pos"} = $comp->{exp_name_1};
    124             $exp_name{"$skycell_id.neg"} = $comp->{exp_name_2} if defined $comp->{bothways} and $comp->{bothways};
     123            $direction{"$skycell_id.pos"} = 1;
     124
     125            # Negative direction
     126            if (defined $comp->{bothways} and $comp->{bothways}) {
     127                $files{"$skycell_id.neg"} = $ipprc->filename( "PPSUB.INVERSE.SOURCES", $path_base );
     128                $zp{"$skycell_id.neg"} = $zp;
     129                $exp_id{"$skycell_id.neg"} = $comp->{exp_id_2};
     130                $direction{"$skycell_id.neg"} = 0;
     131                $exp_name{"$skycell_id.neg"} = $comp->{exp_name_2};
     132            }
     133
    125134        } elsif ($stage eq 'camera') {
    126135            my $cam_id = $comp->{cam_id};
     
    129138            $exp_id{$cam_id} = $comp->{exp_id};
    130139            $exp_name{$cam_id} = $comp->{exp_name};
     140            $direction{$cam_id} = 1;
    131141        }
    132142    }
     
    145155    my $exp_id = $exp_id{$comp};
    146156    my $exp_name = $exp_name{$comp};
     157    my $direction = $direction{$comp};
    147158    if ($product eq "IPP-MOPS") {
    148159        my $outuri = "$outroot.$comp.fits";
     
    150161            &my_die( "Unable to resolve output file $outuri", $pub_id, $PS_EXIT_SYS_ERROR);
    151162
    152         my $command = "$ppMops $file $zp $exp_id $exp_name $out";
     163        my $command = "$ppMops $file $zp $exp_id $exp_name $direction $out";
    153164        my ( $success, $error_code, $full_buf, $stdout_buf, $stderr_buf ) =
    154165            run(command => $command, verbose => $verbose);
Note: See TracChangeset for help on using the changeset viewer.