Changeset 34772 for branches/czw_branch/20120906/ippScripts/scripts/regenerate_background.pl

Timestamp:

Dec 6, 2012, 6:32:19 PM (14 years ago)

Author:

watersc1

Message:

merge from trunk. Preliminary versions of stacktool/warptool updates and regenerate_background.pl script. The warp code is finished and tested, and I still need to get the stack version resolved.

File:

• branches/czw_branch/20120906/ippScripts/scripts/regenerate_background.pl (modified) (8 diffs)

branches/czw_branch/20120906/ippScripts/scripts/regenerate_background.pl

@@ -10 +10 @@
 use IPC::Cmd 0.36 qw( can_run run );
 use File::Spec;
+use File::Temp qw( tempfile );
 use PS::IPP::Config 1.01 qw( :standard );
 use Getopt::Long qw (GetOptions :config auto_help auto_version gnu_getopt );
 use Pod::Usage qw( pod2usage );
+use PS::IPP::Metadata::Config;
+use PS::IPP::Metadata::List qw( parse_md_list );
+
+
 
 my $ppConfigDump = can_run('ppConfigDump');
@@ -67 +72 @@
     my $astromSource;               # The astrometry source
     {
-        my $command = "$ppConfigDump -camera $camera -recipe PSWARP DEFAULT -dump-recipe PSWARP -";
+        my $command = "$ppConfigDump -camera $camera -recipe PSWARP BACKGROUND -dump-recipe PSWARP -";
         my ( $success, $error_code, $full_buf, $stdout_buf, $stderr_buf ) =
             run(command => $command, verbose => $verbose);
@@ -97 +102 @@
     }
 
-    my $in_md = $mdcParser->parse_list(join "", @$stdout_buf) or
+    my $in_md = $mdcParser->parse(join "", @$stdout_buf) or
         &my_die("Unable to parse metadata config doc", $stage, $stage_id, $PS_EXIT_PROG_ERROR);
     $imfiles  = parse_md_list($in_md) or
@@ -165 +170 @@
                 &my_die("Unable to generate template skycell",$skycell_id,$warp_id,$PS_EXIT_SYS_ERROR);
         }
-
+        my $bkgOut_log = prepare_output("LOG.EXP",$path_base . ".BKG_REG",$skycell_id,1);
 
         # Construct pswarp command
         my $pswarp_command = "$pswarp ";
         $pswarp_command .= " -list $bkgList_name -astromlist $astromList_name -bkglist $bkgList_name ";
+        $pswarp_command .= " ${path_base} ";
         $pswarp_command .= " $skyCell_name ";
         $pswarp_command .= " -F PSPHOT.PSF.SAVE PSPHOT.PSF.SKY.SAVE -F PSPHOT.OUTPUT PSPHOT.OUT.CMF.MEF ";
@@ -186 +192 @@
         unless ($success) {
             $error_code = (($error_code >> 8)  or $PS_EXIT_PROG_ERROR);
+            &my_die("Unable to perform pswarp: $error_code", $stage, $stage_id, $error_code);
+        }
+        
+        my $update_command = "$warptool ";
+        $update_command .= " -updateskyfile -warp_id $warp_id -skycell_id $skycell_id ";
+        $update_command .= " -set_background_model 1 ";
+        $update_command .= " -dbname $dbname " if defined $dbname;
+
+        ( $success, $error_code, $full_buf, $stdout_buf, $stderrr_buf ) = 
+            run(command => $update_command, verbose => $verbose);
+        unless ($success) {
+            $error_code = (($error_code >> 8)  or $PS_EXIT_PROG_ERROR);
             &my_die("Unable to perform warptool: $error_code", $stage, $stage_id, $error_code);
         }
@@ -223 +241 @@
     
     my $num = 0;
-    
+    my $path_base = '/tmp/x';
+    my $expTime   = 60;
+
     my ($inputMDC_file, $inputMDC_name) = tempfile("/tmp/bkgreg.stack.mdc.${stage_id}.XXXX",
                                                    UNLINK => !$save_temps);
@@ -242 +262 @@
     }
     close($inputMDC_file);
-    
+
+    my $bkgOut_log = prepare_output("LOG.IMFILE",$path_base . "BKG_REG",$skycell_id,1);    
+
     my $ppStack_command = "$ppStackMedian -input $inputMDC_name $path_base ";
     $ppStack_command   .= " -recipe PPSTACK STACK -recipe PPSUB STACK -recipe PSPHOT STACK ";
@@ -253 +275 @@
     $ppStack_command   .= " -Db PHOTOMETRY F -Db VARIANCE F -Db CONVOLVE F ";
     $ppStack_command   .= " -Df DEFAULT.EXPTIME $expTime ";
-    $ppStack_command   .= " -threads $threads -log $bkgOut_log ";
+    $ppStack_command   .= " -threads 1 -log $bkgOut_log ";
     $ppStack_command   .= " -dbname $dbname " if defined $dbname;