package Subs::XrtTdrss; ############################################################################## # # DESCRIPTION: This subroutine runs the xrttdrss script to process xrt tdrss # data. # # HISTORY: # HISTORY: $Log: XrtTdrss.pm,v $ # HISTORY: Revision 1.4 2006/01/31 16:49:45 apsop # HISTORY: Fix bugs in renaming and combining postage stamp images. # HISTORY: # HISTORY: Revision 1.3 2006/01/29 19:39:47 apsop # HISTORY: Code for combining multiple tdrss images into one file. # HISTORY: # HISTORY: Revision 1.2 2005/11/08 20:10:42 apsop # HISTORY: <Previous comment bogus>New module for processing xrt tdrss data. # HISTORY: # HISTORY: Revision 1.1 2005/11/08 19:22:28 apsop # HISTORY: Populate the TIMELIST and DATALIST hashes. Used to be an SWCheckInput. # HISTORY: # HISTORY: # # VERSION: 0.0 # ############################################################################## use Subs::Sub; use Util::SwiftTags; @ISA = ("Subs::Sub"); use strict; sub new { my $proto=shift; my $self=$proto->SUPER::new(); $self->{DESCRIPTION}="Process XRT tdrss messages"; return $self; } ################## # METHODS: ################## sub body { my $self=shift; my $log =$self->log(); my $filename=$self->filename(); my $procpar =$self->procpar(); my $jobpar =$self->jobpar(); my $tdrss = Util::HEAdas->new('xrttdrss') ->params({spec1file => 'NONE', spec2file => 'NONE', posfile => 'CALDB', imbiasfile => 'CALDB', pdbiasfile => 'CALDB'}) ->is_script(1); my @images = $filename->get('tdimage', 'x', 'im', '*'); foreach my $image (@images){ my $index = ( $filename->parse($image, 'tdimage') )[2]; my ($sk_image, $rw_image) = ($image) x 2; unless($index){ $sk_image =~ s/\.img/_sk\.img/; $rw_image =~ s/\.img/_rw\.img/; }else{ $sk_image =~ s/(_\d{2}\.img)/_sk$1/; $rw_image =~ s/(_\d{2}\.img)/_rw$1/; } rename $image, $rw_image; $tdrss->params({imagefile => $rw_image, outimagefile => $sk_image}) ->run(); } ############################################### # Check it we need to reformat the image files ############################################### my $tempim = 'xrttdrss_fits.tmp'; my $create = Util::Ftool->new('fimgcreate') ->params({bitpix => 8, naxes => '0', datafile => 'none', outfile => $tempim}); my $append = Util::HEAdas->new('ftappend'); my $cphead = Util::Ftool->new('cphead'); foreach my $type ('tdrawimage', 'tdskyimage'){ @images = sort( $filename->get($type, 'x', 'im', '*') ); if(@images){ my $image = $images[0]; my $imfits = Util::FITSfile->new($image); if( $imfits->nhdus() <= 1 || $imfits->keyword('EXTNAME') !~ /(LNG|SHT)\d{9}I\d{2}/ ){ $create->run(); foreach $image (@images){ $append->params({infile => $image .'[0]', outfile => $tempim}) ->run(); } unlink @images; rename $tempim, $image; $imfits = Util::FITSfile->new($image); my ($tstart, $tstop) = (1E10, 0); for(my $ext=1; $ext < $imfits->nhdus(); $ext++){ $imfits->ext($ext); my ($start, $stop) = ($imfits->keyword('TSTART'), $imfits->keyword('TSTOP')); $tstart = $start if $start < $tstart; $tstop = $stop if $stop > $tstop; my $int = int($start); my $frac = int( ($start-$int)/2E-2 ); my $mode = $imfits->keyword('DATAMODE'); my $pneu = 'UNK'; $pneu = 'LNG' if $mode =~ /LONGIMA/; $pneu = 'SHT' if $mode =~ /SHORTIMA/; $imfits->keyword('EXTNAME', $pneu.$int.'I'.$frac); } my $start_date = Util::Date->new($tstart); my $stop_date = Util::Date->new($tstop); $imfits->ext(0); $imfits->begin_many_keywords(); $imfits->keyword('INSTRUME', 'XRT '); $imfits->keyword('TSTART', $tstart); $imfits->keyword('TSTOP', $tstop); $imfits->keyword('DATE-OBS', $start_date->date().'T'.$start_date->time() ); $imfits->keyword('DATE-END', $stop_date->date().'T'.$stop_date->time() ); $imfits->end_many_keywords(); } } } } 1;