—
—
—
—
—
—
—
—
—
—
—
—
—

              
# -*-perl-*-
package Astro::FITS::HdrTrans::LCOFLOYDS;
HideShow 17 lines of Pod
=head1 NAME
Astro::FITS::HdrTrans::LCOFLOYDS - LCO 2.0m FLOYDS translations
=head1 SYNOPSIS
  use Astro::FITS::HdrTrans::LCOFLOYDS;
  %gen = Astro::FITS::HdrTrans::LCOFLOYDS->translate_from_FITS( %hdr );
=head1 DESCRIPTION
This class provides a generic set of translations that are specific to
2.0m FLOYDSs at LCO.
=cut
use 5.006;
use warnings;
use strict;
use Carp;
# Inherit from LCO base class.
use base qw/ Astro::FITS::HdrTrans::LCO /;
our $VERSION = "1.66";
# for a constant mapping, there is no FITS header, just a generic
# header that is constant
# NULL mappings used to override base-class implementations.
my @NULL_MAP = qw/ /;
my %CONST_MAP = ( OBSERVATION_MODE => 'spectroscopy',
                 GRATING_NAME => 'FLOYDS GRATING',
                 GRATING_ORDER       => 1,
                 GRATING_DISPERSION  => 1.73/10000.0,
                 GRATING_WAVELENGTH  => 0.556,
                 NSCAN_POSITIONS     => 1,
                 NUMBER_OF_READS     => 1,
                 SCAN_INCREMENT      => 1,
                );
my %UNIT_MAP = (
                 SLIT_NAME  => "APERTURE",
                 SLIT_WIDTH => "APERWID",
                 X_DIM      => "NAXIS1",
                 Y_DIM      => "NAXIS2",
             );
# Create the translation methods
__PACKAGE__->_generate_lookup_methods( \%CONST_MAP, \%UNIT_MAP, \@NULL_MAP );
HideShow 17 lines of Pod
=head1 METHODS
=over 4
=item B<this_instrument>
The name of the instrument required to match (case insensitively)
against the INSTRUME/INSTRUMENT keyword to allow this class to
translate the specified headers. Called by the default
C<can_translate> method.
  $inst = $class->this_instrument();
Returns "LCOFLOYDS".
=cut
sub this_instrument {
   return qr/(^en05)|(^en06)/i
}
HideShow 23 lines of Pod
=back
=head1 COMPLEX CONVERSIONS
These methods are more complicated than a simple mapping. We have to
provide both from- and to-FITS conversions All these routines are
methods and the to_ routines all take a reference to a hash and return
the translated value (a many-to-one mapping) The from_ methods take a
reference to a generic hash and return a translated hash (sometimes
these are many-to-many)
=over 4
=cut
=item B<to_DEC_SCALE>
Sets the declination scale in arcseconds per pixel.  The C<PIXSCALE>
is used when it's defined.  Otherwise it returns a default value of 0.2320
arcsec/pixel, multiplied by C<YBINNING> assuming this is defined
=cut
sub to_DEC_SCALE {
   my $self = shift;
   my $FITS_headers = shift;
   my $decscale = 0.3860;
# Assumes either x-y scales the same or the y corresponds to
# declination.
   my $ccdscale = $self->via_subheader( $FITS_headers, "PIXSCALE" );
   if ( defined $ccdscale ) {
      $decscale = $ccdscale;
   } else {
      my $ybinning = $self->via_subheader( $FITS_headers, "YBINNING" );
      if ( defined $ybinning ) {
        $decscale = $decscale * $ybinning;
      }
   }
   return $decscale;
}
HideShow 8 lines of Pod
=item B<to_DEC_TELESCOPE_OFFSET>
Sets the declination telescope offset in arcseconds.   It uses the
C<CAT-DEC> and C<DEC> keywords to derive the offset, and if either
does not exist, it returns a default of 0.0.
=cut
sub to_DEC_TELESCOPE_OFFSET {
   my $self = shift;
   my $FITS_headers = shift;
   my $decoffset = 0.0;
   if ( exists $FITS_headers->{"CAT-DEC"} && exists $FITS_headers->{DEC} ) {
# Obtain the reference and telescope declinations positions measured in degrees.
      my $refdec = $self->dms_to_degrees( $FITS_headers->{"CAT-DEC"} );
      my $dec = $self->dms_to_degrees( $FITS_headers->{DEC} );
# Find the offsets between the positions in arcseconds on the sky.
      $decoffset = 3600.0 * ( $dec - $refdec );
   }
# The sense is reversed compared with UKIRT, as these measure the
# places on the sky, not the motion of the telescope.
   return -1.0 * $decoffset;
}
HideShow 9 lines of Pod
=item B<to_DR_RECIPE>
Returns the data-reduction recipe name.  The selection depends on the
values of the C<OBJECT> and C<OBSTYPE> keywords.  The default is
"QUICK_LOOK".  A dark returns "REDUCE_DARK", and an object's recipe is
"JITTER_SELF_FLAT".
=cut
sub to_DR_RECIPE {
  my $self = shift;
  my $FITS_headers = shift;
  my $recipe = "QUICK_LOOK";
  if ( exists $FITS_headers->{OBSTYPE} ) {
    if ( $FITS_headers->{OBSTYPE} =~ /ARC/i ) {
      $recipe = "REDUCE_ARC";
    } elsif ( $FITS_headers->{OBSTYPE} =~ /BIAS/i ) {
      $recipe = "REDUCE_BIAS";
    } elsif ( $FITS_headers->{OBSTYPE} =~ /DARK/i ) {
      $recipe = "REDUCE_DARK";
    } elsif ( $FITS_headers->{OBSTYPE} =~ /LAMPFLAT/i ) {
      $recipe = "REDUCE_FLAT";
    } elsif ( $FITS_headers->{OBSTYPE} =~ /EXPOSE/i ) {
      #     $recipe = "JITTER_SELF_FLAT";
      $recipe = "POINT_SOURCE_NOSTD";
    } elsif ( $FITS_headers->{OBSTYPE} =~ /STANDARD/i ) {
      $recipe = "STANDARD_STAR";
    }
  }
  return $recipe;
}
HideShow 8 lines of Pod
=item B<to_RA_SCALE>
Sets the RA scale in arcseconds per pixel.  The C<PIXSCALE>
is used when it's defined.  Otherwise it returns a default value of 0.2320
arcsec/pixel, multiplied by C<XBINNING> assuming this is defined (1.0 otherwise)
=cut
sub to_RA_SCALE {
   my $self = shift;
   my $FITS_headers = shift;
   my $rascale = 0.3860;
# Assumes either x-y scales the same or the x corresponds to
# ra.
   my $ccdscale = $self->via_subheader( $FITS_headers, "PIXSCALE" );
   if ( defined $ccdscale ) {
      $rascale = $ccdscale;
   } else {
      my $xbinning = $self->via_subheader( $FITS_headers, "XBINNING" );
      if ( defined $xbinning ) {
        $rascale = $rascale * $xbinning;
      }
   }
   return $rascale;
}
HideShow 8 lines of Pod
=item B<to_RA_TELESCOPE_OFFSET>
Sets the right-ascension telescope offset in arcseconds.   It uses the
C<CAT-RA>, C<RA>, C<CAT-DEC> keywords to derive the offset, and if any
of these keywords does not exist, it returns a default of 0.0.
=cut
sub to_RA_TELESCOPE_OFFSET {
   my $self = shift;
   my $FITS_headers = shift;
   my $raoffset = 0.0;
   if ( exists $FITS_headers->{"CAT-DEC"} &&
        exists $FITS_headers->{"CAT-RA"} && exists $FITS_headers->{RA} ) {
# Obtain the reference and telescope sky positions measured in degrees.
      my $refra = $self->hms_to_degrees( $FITS_headers->{"CAT-RA"} );
      my $ra = $self->hms_to_degrees( $FITS_headers->{RA} );
      my $refdec = $self->dms_to_degrees( $FITS_headers->{"CAT-DEC"} );
# Find the offset between the positions in arcseconds on the sky.
      $raoffset = 3600.0 * ( $ra - $refra ) * $self->cosdeg( $refdec );
   }
# The sense is reversed compared with UKIRT, as these measure the
# place son the sky, not the motion of the telescope.
   return -1.0 * $raoffset;
}
HideShow 7 lines of Pod
=item B<to_SLIT_ANGLE>
Returns the slit PA, trapping UNKNOWN values and setting them to -999 (could
in principle be recalculated from HA, Dec and Latitude via some SLALIBing...
=cut
sub to_SLIT_ANGLE {
   my $self = shift;
   my $FITS_headers = shift;
   my $slit_angle = -999.0;
   if ( exists $FITS_headers->{APERPA} ) {
      my $slit_pa = $FITS_headers->{APERPA};
      if ( $slit_pa !~ /UNKNOWN/i ) {
         $slit_angle = $slit_pa;
      }
   }
   return $slit_angle;
}
HideShow 7 lines of Pod
=item B<to_X_LOWER_BOUND>
Returns the lower bound along the X-axis of the area of the detector
as a pixel index.
=cut
sub to_X_LOWER_BOUND {
   my $self = shift;
   my $FITS_headers = shift;
   my @bounds = $self->getbounds( $FITS_headers );
   return $bounds[ 0 ];
}
HideShow 7 lines of Pod
=item B<to_X_UPPER_BOUND>
Returns the upper bound along the X-axis of the area of the detector
as a pixel index.
=cut
sub to_X_UPPER_BOUND {
   my $self = shift;
   my $FITS_headers = shift;
   my @bounds = $self->getbounds( $FITS_headers );
   return $bounds[ 1 ];
}
HideShow 7 lines of Pod
=item B<to_Y_LOWER_BOUND>
Returns the lower bound along the Y-axis of the area of the detector
as a pixel index.
=cut
sub to_Y_LOWER_BOUND {
   my $self = shift;
   my $FITS_headers = shift;
   my @bounds = $self->getbounds( $FITS_headers );
   return $bounds[ 2 ];
}
HideShow 7 lines of Pod
=item B<to_Y_UPPER_BOUND>
Returns the upper bound along the Y-axis of the area of the detector
as a pixel index.
=cut
sub to_Y_UPPER_BOUND {
   my $self = shift;
   my $FITS_headers = shift;
   my @bounds = $self->getbounds( $FITS_headers );
   return $bounds[ 3 ];
}
# Supplementary methods for the translations
# ------------------------------------------
# Obtain the detector bounds from a section in [xl:xu,yl:yu] syntax.
# If the TRIMSEC header is absent, use a default which corresponds
# to the useful part of the array (minus bias strips).
sub getbounds{
   my $self = shift;
   my $FITS_headers = shift;
   my @bounds = ( 1, 2048, 1,  512 );
   if ( exists $FITS_headers->{CCDSUM} ) {
      my $binning = $FITS_headers->{CCDSUM};
      if ( $binning eq '2 2' ) {
         @bounds = ( 1, 1024,  1,  256 );
      }
   }
   if ( exists $FITS_headers->{TRIMSEC} ) {
      my $section = $FITS_headers->{TRIMSEC};
      if ( $section !~ /UNKNOWN/i ) {
        $section =~ s/\[//;
        $section =~ s/\]//;
        $section =~ s/,/:/g;
        my @newbounds = split( /:/, $section );
        if (@newbounds == grep { $_ == 0 } @newbounds) {
          print "ERR: TRIMSEC all 0\n";
        } else {
          if ( $FITS_headers->{INSTRUME} !~ /^en0/i ) {
# Unless this is any data (which has a bad TRIMSEC), update bounds array
            @bounds = @newbounds;
          }
        }
      }
   }
#   print("DBG: Bounds=@bounds\n");
   return @bounds;
}
HideShow 14 lines of Pod
=back
=head1 SEE ALSO
C<Astro::FITS::HdrTrans>, C<Astro::FITS::HdrTrans::LCO>.
=head1 AUTHOR
Tim Lister E<lt>tlister@lcogt.netE<gt>
=head1 COPYRIGHT
=cut
1;