Changeset 20483 for trunk/ippScripts/scripts/magic_tree.pl

Timestamp:

Oct 30, 2008, 7:36:29 PM (18 years ago)

Author:

bills

Message:

require --warp_id so we can look up the skycell files

File:

• trunk/ippScripts/scripts/magic_tree.pl (modified) (9 diffs)

trunk/ippScripts/scripts/magic_tree.pl

@@ -33 +33 @@
 
 # Parse the command-line arguments
-my ($magic_id, $tess_id, $camera, $ra0, $dec0, $dbname, $outroot,
+my ($magic_id, $warp_id, $tess_id, $camera, $ra0, $dec0, $dbname, $outroot,
     $save_temps, $verbose, $no_update, $no_op, $redirect);
 GetOptions(
            'magic_id=s'      => \$magic_id,   # Magic identifier
+           'warp_id=s'       => \$warp_id,    # id for corresponding warps
            'tess_id=s'       => \$tess_id,    # Tessellation identifier
            'camera=s'        => \$camera,     # Camera name
@@ -51 +52 @@
 
 pod2usage( -msg => "Unknown option: @ARGV", -exitval => 2 ) if @ARGV;
-pod2usage( -msg => "Required options: --magic_id --camera --tess_id --ra --dec --outroot",
+pod2usage( -msg => "Required options: --magic_id --warp_id --camera --tess_id --ra --dec --outroot",
            -exitval => 3) unless
     defined $magic_id and
+    defined $warp_id and
     defined $tess_id and
     defined $ra0 and
@@ -69 +71 @@
 my $missing_tools;
 my $magictool = can_run('magictool') or (warn "Can't find magictool" and $missing_tools = 1);
+my $warptool = can_run('warptool') or (warn "Can't find warptool" and $missing_tools = 1);
 if ($missing_tools) {
     warn("Can't find required tools.");
@@ -76 +79 @@
 my $mdcParser = PS::IPP::Metadata::Config->new; # Parser for metadata config files
 
+### Get a list of warpSkyfiles
+my %warpSkyfiles;                   # hash of warps
+{
+    my $command = "$warptool -warped -warp_id $warp_id"; # Command to run
+    $command .= " -dbname $dbname" if defined $dbname;
+    my ( $success, $error_code, $full_buf, $stdout_buf, $stderr_buf ) =
+        run(command => $command, verbose => $verbose);
+    unless ($success) {
+        $error_code = (($error_code >> 8) or $PS_EXIT_PROG_ERROR);
+        &my_die("Unable to perform warptool -warped: $error_code", $magic_id, $error_code);
+    }
+
+    my $metadata = $mdcParser->parse(join "", @$stdout_buf) or
+        &my_die("Unable to parse metadata config doc", $magic_id, $PS_EXIT_PROG_ERROR);
+
+    my $warps = parse_md_list($metadata) or
+        &my_die("Unable to parse metadata list", $magic_id, $PS_EXIT_PROG_ERROR);
+
+    # make a hash indexed by skycell_id 
+    foreach my $warp ( @$warps ) {
+        my $skycell_id = $warp->{skycell_id};
+        die "failed to lookup skycell_id from warp" if !$skycell_id;
+        $warpSkyfiles{$skycell_id} = $warp;
+    }
+}
 ### Get a list of skycells
 my @skycells;                   # List of skycells
@@ -102 +130 @@
 my @fields;
 foreach my $skycell_id ( @skycells ) {
-    my $skyfile = $ipprc->filename("SKYCELL.TEMPLATE", $outroot, $skycell_id );
-    $ipprc->skycell_file($tess_id, $skycell_id, $skyfile, $verbose) or &my_die("Unable to generate skycells $skycell_id", $magic_id, $PS_EXIT_PROG_ERROR);
+    # the filename method doesn't add the $skycell_id
+
+    # my $skyfile = $ipprc->filename("SKYCELL.TEMPLATE", $outroot, $skycell_id );
+    #my $skyfile = $ipprc->filename("SKYCELL.TEMPLATE", $outroot, $skycell_id );
+
+#    my $skycellBase = "$outroot/$skycell_id";
+#    my $skyfile = $ipprc->filename("SKYCELL.TEMPLATE", $skycellBase);
+#    $ipprc->skycell_file($tess_id, $skycell_id, $skyfile, $verbose) or &my_die("Unable to generate skycells $skycell_id", $magic_id, $PS_EXIT_PROG_ERROR);
+
+    my $warp = $warpSkyfiles{$skycell_id};
+    die "warpSkyfile for $skycell_id not found" if !$warp;
+    my $skyfileBase = $warp->{path_base};
+    my $skyfile = $ipprc->filename("SKYCELL.TEMPLATE", $skyfileBase);
     my $skyfileResolved = $ipprc->file_resolve( $skyfile );
     my ($header, $status) = Astro::FITS::CFITSIO::fits_read_header( $skyfileResolved );
@@ -181 +220 @@
 };
 my @tasks = ( $root );
-while (scalar @tasks > 1) {
+while (scalar @tasks) {
     my $node = shift @tasks;
     divide_node($node, \@tasks);
-}
+} 
 
 ### Format tree for magictool
@@ -193 +232 @@
 
 ### Input tree into database
-{
+if (!$no_update) {
     my $command = "$magictool -inputtree";
     $command   .= " -magic_id $magic_id";
@@ -241 +280 @@
     my @upper = splice(@sorted, $median); # Upper part of the sorted list
 
+#    print " median $median\n";
+
     return (\@sorted, \@upper);
 }
@@ -257 +298 @@
         children => {},
     };
+
+#    my $n = scalar @$contents;
+#    print "new node: $node->{position}: $n\n";
 
     $parent->{children}->{$position} = $node;