MOS Spectroscopy of Extragalactic Field#

Use case: emission-line measurements and template matching on 1D spectra.
Data: LEGA-C spectra and galaxy template spectra; optical rest-frame.
Tools: specutils, astropy, matplotlib, pandas.
Cross-intrument: NIRSpec
Documentation: This notebook is part of a STScI’s larger post-pipeline Data Analysis Tools Ecosystem.

Introduction#

This notebook will perform a seris of spectroscopic analyses on multiple spectra, including smoothing, continuum fitting and subtraction, line identification, centroiding and flux measurements, gaussian fitting, equivalent widths, and template fitting.

Note: This notebook is intended to ultimately be compatible with the final data products (1D and 2D spectra) from the JWST pipeline. These data products are not available yet, so the notebook uses LEGA-C data (van der Wel et al. 2016, Straatmann et al. 2018) for now.

LEGA-C is a galaxy survey of about 3000 galaxies at z~0.8 and M* > 10^10 M_sun in the COSMOS field. The spectra sample the rest-frame optical between ~3000A and 5000A at high resolution and very high signal-to-noise ratio. More information about the survey can be found here: http://www.mpia.de/home/legac/

First, set the environment with astroconda and specutils.

conda create -n astroutils stsci

source activate astroutils

Imports#

#general os
import os
import zipfile
import urllib.request

#general plotting
from matplotlib import pyplot as plt

params={'legend.fontsize':'18','axes.labelsize':'18',
        'axes.titlesize':'18','xtick.labelsize':'18',
        'ytick.labelsize':'18','lines.linewidth':2,'axes.linewidth':2,'animation.html': 'html5'}
plt.rcParams.update(params)
plt.rcParams.update({'figure.max_open_warning': 0})

#table/math handling
import pandas as pd
import numpy as np
np.seterr(all='ignore')  # hides irrelevant warnings about divide-by-zero, etc

#astropy
import astropy
import astropy.units as u
from astropy.table import QTable
from astropy.io import fits,ascii
from astropy.nddata import StdDevUncertainty
from astropy.modeling import models
from astropy.visualization import quantity_support
from astropy import constants as const

#specutils
import specutils
from specutils import Spectrum1D, SpectralRegion
from specutils.manipulation import gaussian_smooth
from specutils.fitting import fit_generic_continuum
from specutils.fitting import find_lines_derivative
from specutils.fitting import find_lines_threshold
from specutils.fitting import fit_lines
from specutils.manipulation import noise_region_uncertainty
from specutils.analysis import centroid
from specutils.analysis import line_flux
from specutils.analysis import equivalent_width
from specutils.analysis import template_comparison

quantity_support();  # auto-recognizes units on matplotlib plots

Check versions. Should be:#

Pandas: 1.0.1

Numpy: 1.18.1

Astropy: 4.0

Specutils: 0.7

print("Pandas: ",pd.__version__)
print("Numpy: ",np.__version__)
print("Astropy: ",astropy.__version__)
print("Specutils: ",specutils.__version__)
Pandas:  2.1.4
Numpy:  1.26.2
Astropy:  6.0.0
Specutils:  1.12.0

Choose one galaxy#

file1d = observedfiles + 'legac_M1_v3.7_spec1d_130902.fits'
file1dwht = observedfiles + 'legac_M1_v3.7_wht1d_130902.fits'
file2d = observedfiles + 'legac_M1_v3.7_spec2d_130902.fits'

Inspect its 2D spectrum#

hdu2d = fits.open(file2d)
hdu2d.info()
Filename: ./mos_spectroscopy/observed/legac_M1_v3.7_spec2d_130902.fits
No.    Name      Ver    Type      Cards   Dimensions   Format
  0  PRIMARY       1 PrimaryHDU     716   (6166, 81)   float32   
hdu2d[0].header
SIMPLE  =                    T / file does conform to FITS standard             
BITPIX  =                  -32 / number of bits per data pixel                  
NAXIS   =                    2 / number of data axes                            
NAXIS1  =                 6166 / length of data axis 1                          
NAXIS2  =                   81 / length of data axis 2                          
EXTEND  =                    T / FITS dataset may contain extensions            
COMMENT   FITS (Flexible Image Transport System) format is defined in 'Astronomy
COMMENT   and Astrophysics', volume 376, page 359; bibcode: 2001A&A...376..359H 
DATE    = '2015-01-26T12:24:07' / file creation date (YYYY-MM-DDThh:mm:ss UT)   
ORIGIN  = 'ESO     '           / European Southern Observatory                  
MJD-OBS =       57040.16930682 / Obs start 2015-01-18T04:03:48.109              
DATE-OBS= '2015-01-18T04:03:48.108' / Date of observation                       
EXPTIME =             1199.997 / Total integration time. 00:19:59.997           
TELESCOP= 'ESO-VLT-U3'         / ESO <TEL>                                      
RA      =           150.220555 / 10:00:52.9 RA (J2000) pointing                 
DEC     =              2.41104 / 02:24:39.7  DEC (J2000) pointing               
EQUINOX =                2000. / Standard FK5                                   
RADECSYS= 'FK5     '           / Coordinate reference frame                     
LST     =            25866.379 / 07:11:06.379 LST at start                      
UTC     =               14623. / 04:03:43.000 UT at start                       
OBSERVER= 'VIMOS   '           / Name of observer                               
INSTRUME= 'VIMOS   '           / Instrument used                                
PI-COI  = '555555555'          / Name of PI and COI                             
OBJECT  = 'COSMOS_M01'         / Original target                                
CRVAL1  =               5800.3 / 10:00:52.9, RA at ref pixel                    
CRPIX1  =                   1. / Reference pixel in X                           
CTYPE1  = 'LINEAR  '           / pixel coordinate system                        
CRVAL2  =                   1. / 02:24:39.7, DEC at ref pixel                   
CRPIX2  =                   1. / Reference pixel in Y                           
CTYPE2  = 'PIXEL   '           / pixel coordinate system                        
INHERIT =                    T / denotes the INHERIT keyword convention         
CD1_1   =                  0.6 / Translation matrix element                     
CD1_2   =                   0. / Translation matrix element                     
CD2_1   =                   0. / Translation matrix element                     
CD2_2   =                   1. / Translation matrix element                     
DATAMD5 = '77a1ef86719040e2fea8d1d96fef35e0' / MD5 checksum                     
PIPEFILE= 'mos_science_sky_extracted.fits' / Filename of data product           
ARCFILE = 'VIMOS.2015-01-18T04:03:48.110.fits' / Archive File Name              
HIERARCH ESO OBS DID = 'ESO-VLT-DIC.OBS-1.12' / OBS Dictionary                  
HIERARCH ESO OBS EXECTIME = 6835 / Expected execution time                      
HIERARCH ESO OBS GRP = '0       ' / linked blocks                               
HIERARCH ESO OBS ID =  1163772 / Observation block ID                           
HIERARCH ESO OBS NAME = 'legac_M01_2' / OB name                                 
HIERARCH ESO OBS OBSERVER = 'UNKNOWN ' / Observer Name                          
HIERARCH ESO OBS PI-COI ID = 73498 / ESO internal PI-COI ID                     
HIERARCH ESO OBS PI-COI NAME = 'UNKNOWN ' / PI-COI name                         
HIERARCH ESO OBS PROG ID = '194.A-2005(B)' / ESO program identification         
HIERARCH ESO OBS START = '2015-01-18T03:38:14' / OB start time                  
HIERARCH ESO OBS NTPL =      3 / Number of templates within OB                  
HIERARCH ESO OBS TPLNO =     2 / Template number within OB                      
HIERARCH ESO OBS AIRM =     2. / Req. max. airmass                              
HIERARCH ESO OBS AMBI FWHM = 1.3 / Req. max. seeing                             
HIERARCH ESO OBS AMBI TRANS = '2CLR    ' / Req. sky transparency                
HIERARCH ESO OBS ATM = 'NATM    ' / Req. Atmospheric Turbulence Model           
HIERARCH ESO OBS CONTAINER ID = -999 / Scheduling container ID                  
HIERARCH ESO OBS CONTAINER TYPE = 'U       ' / Scheduling container type        
HIERARCH ESO OBS MOON DIST = 120 / Req. min. angular dist. from moon            
HIERARCH ESO OBS MOON FLI = 0.3 / Req. max. fractional lunar illum.             
HIERARCH ESO OBS STREHLRATIO = 0. / Req. strehl ratio                           
HIERARCH ESO OBS TARG NAME = 'COSMOS_M01' / OB target name                      
HIERARCH ESO OBS TWILIGHT = -15 / Req. twilight                                 
HIERARCH ESO OBS WATERVAPOUR = 0. / Req. water vapour                           
HIERARCH ESO TPL FILE DIRNAME = '$INS_ROOT/$INS_USER/ADF' / Storage filename    
HIERARCH ESO TPL FILE KEEP = T / paramfile keep-flag                            
HIERARCH ESO TPL DID = 'ESO-VLT-DIC.TPL-1.9' / Data dictionary for TPL          
HIERARCH ESO TPL ID = 'VIMOS_mos_obs_Offset' / Template signature ID            
HIERARCH ESO TPL NAME = 'MOS Observation' / Template name                       
HIERARCH ESO TPL PRESEQ = 'VIMOS_mos_obs_Offset' / Sequencer script             
HIERARCH ESO TPL START = '2015-01-18T04:01:12' / TPL start time                 
HIERARCH ESO TPL VERSION = '@(#) $Revision: 77996 $' / Version of the templa    
HIERARCH ESO TPL NEXP =      4 / Number of exposures within templat             
HIERARCH ESO TPL EXPNO =     1 / Exposure number within template                
HIERARCH ESO TEL DID = 'ESO-VLT-DIC.TCS' / Data dictionary for TEL              
HIERARCH ESO TEL ID = 'v 261159' / TCS version number                           
HIERARCH ESO TEL DATE = '2000-01-01T00:00:00' / TCS installation date           
HIERARCH ESO TEL ALT =  40.643 / Alt angle at start                             
HIERARCH ESO TEL AZ =  243.107 / Az angle at start S=0,W=90                     
HIERARCH ESO TEL GEOELEV = 2648. / Elevation above sea level (m)                
HIERARCH ESO TEL GEOLAT = -24.6268 / Tel geo latitute (+=North)                 
HIERARCH ESO TEL GEOLON = -70.4045 / Tel geo longitude (+=East)                 
HIERARCH ESO TEL OPER = 'I, Condor' / Telescope Operator                        
HIERARCH ESO TEL FOCU ID = 'NB      ' / Telescope focus station ID              
HIERARCH ESO TEL FOCU LEN = 120. / Focal length                                 
HIERARCH ESO TEL FOCU SCALE = 1.718 / Focal scale                               
HIERARCH ESO TEL FOCU VALUE = -31.348 / M2 setting                              
HIERARCH ESO TEL PARANG START = -125.768 / Parallactic angle at start           
HIERARCH ESO TEL AIRM START = 1.533 / Airmass at start                          
HIERARCH ESO TEL AMBI FWHM START = 1.08 / Observatory Seeing queried from AS    
HIERARCH ESO TEL AMBI PRES START = 744.4 / Observatory ambient air pressure q   
HIERARCH ESO TEL AMBI WINDSP = 0.82 / Observatory ambient wind speed que        
HIERARCH ESO TEL AMBI WINDDIR = 325.5 / Observatory ambient wind direction      
HIERARCH ESO TEL AMBI RHUM = 3. / Observatory ambient relative humid            
HIERARCH ESO TEL AMBI TEMP = 15.77 / Observatory ambient temperature qu         
HIERARCH ESO TEL MOON RA = 265.331996 / 17:41:19.6 RA (J2000)                   
HIERARCH ESO TEL MOON DEC = -17.93387 / -17:56:01.9 DEC (J2000)                 
HIERARCH ESO TEL TH M1 TEMP = 13.63 / M1 superficial temperature                
HIERARCH ESO TEL TRAK STATUS = 'NORMAL  ' / Tracking status                     
HIERARCH ESO TEL DOME STATUS = 'FULLY-OPEN' / Dome status                       
HIERARCH ESO TEL CHOP ST =   F / True when chopping is active                   
HIERARCH ESO TEL TARG ALPHA = 100052.8 / Alpha coordinate for the target        
HIERARCH ESO TEL TARG DELTA = 22436. / Delta coordinate for the target          
HIERARCH ESO TEL TARG EPOCH = 2000. / Epoch                                     
HIERARCH ESO TEL TARG EPOCHSYSTEM = 'J       ' / Epoch system (default J=Julian)
HIERARCH ESO TEL TARG EQUINOX = 2000. / Equinox                                 
HIERARCH ESO TEL TARG PMA = 0. / Proper Motion Alpha                            
HIERARCH ESO TEL TARG PMD = 0. / Proper motion Delta                            
HIERARCH ESO TEL TARG RADVEL = 0. / Radial velocity                             
HIERARCH ESO TEL TARG PARALLAX = 0. / Parallax                                  
HIERARCH ESO TEL TARG COORDTYPE = 'M       ' / Coordinate type (M=mean A=apparen
HIERARCH ESO TEL PARANG END = -128.789 / Parallactic angle at end               
HIERARCH ESO TEL AIRM END = 1.421 / Airmass at end                              
HIERARCH ESO TEL AMBI FWHM END = 1.03 / Observatory Seeing queried from AS      
HIERARCH ESO TEL IA FWHM =  1. / Delivered seeing corrected by airm             
HIERARCH ESO TEL IA FWHMLIN = 1.01 / Delivered seeing on IA detector (l         
HIERARCH ESO TEL IA FWHMLINOBS = 1.25 / Delivered seeing on IA detector (l      
HIERARCH ESO TEL AMBI PRES END = 744.3 / Observatory ambient air pressure q     
HIERARCH ESO TEL AMBI TAU0 = 0.0039 / Average coherence time                    
HIERARCH ESO TEL GS1 ALPHA = 100057.969 / Guide Star alpha                      
HIERARCH ESO TEL GS1 DELTA = 23539.05 / Guide Star delta                        
HIERARCH ESO TEL GS1 PPOS = 'NEG     ' / Selected probe position (POS or NE     
HIERARCH ESO TEL TARG NAME = '        ' / Name of special target position       
HIERARCH ESO INS SWSIM = 'NORMAL  ' / Software simulation                       
HIERARCH ESO INS ID = 'VIMOS   ' / Instrument identification                    
HIERARCH ESO INS DID = 'ESO-VLT-DIC.VIMOS_ICS-242857' / Data dictionary         
HIERARCH ESO INS SOFT DATE = 'Mar 20 2014' / Instrument SW Release date         
HIERARCH ESO INS SOFT ID = '186976  ' / Version of ICS                          
HIERARCH ESO INS FOCU1 TEMP = 14.61792 / Temperature(C) for the beam            
HIERARCH ESO INS ADF TYPE = 'MOS     ' / Type of ADF                            
HIERARCH ESO INS ADF VERSION = 3.28 / Version of ADF                            
HIERARCH ESO INS ADF UNIT = 'MILLIMETER' / Unit of slit/object coordinates      
HIERARCH ESO INS ADF ID = 'legac_M1Q1_HR_Red_M1Q1.adp' / Name for the ADF       
HIERARCH ESO INS ADF ADMID = '2015-01-18T03:38:06-101-1.adm' / Name for the *.  
HIERARCH ESO INS ADF COMMENT = '        ' / Optional user comment for ADF       
HIERARCH ESO INS MASK1 ID = 103063 / Mask identification number                 
HIERARCH ESO INS ADF GRISM NAME = 'HR_Red  ' / Grism name                       
HIERARCH ESO INS ADF GRISM SPECTLEN = 4096. / Spectrum length in pixels         
HIERARCH ESO INS ADF SKYREG = 1.8 / Sky region in arcsec                        
HIERARCH ESO INS REF NO =    1 / Number of reference apertures                  
HIERARCH ESO INS ARC NO =    0 / Number of curved apertures                     
HIERARCH ESO INS SLIT NO =  33 / Number of all apertures but refere             
HIERARCH ESO INS SHU NO =    1 / Number of shutter positions                    
HIERARCH ESO INS SHU1 POSL = 0. / Lower shutter positions                       
HIERARCH ESO INS SHU1 POSH = 2440. / Upper shutter positions                    
HIERARCH ESO INS REF1 TYPE = 'SQUARE  ' / Type of REF i                         
HIERARCH ESO INS REF1 ID = 132308 / ID of REF i                                 
HIERARCH ESO INS REF1 OBJ RA = 150.345742 / RA (J2000) of object                
HIERARCH ESO INS REF1 OBJ DEC = 2.347825 / DEC (J2000) of object                
HIERARCH ESO INS REF1 X = -27.353 / x co-ordinate of REF i                      
HIERARCH ESO INS REF1 Y = 105.937 / y co-ordinate of REF i                      
HIERARCH ESO INS REF1 DIMX = 3.576 / Size in x direction of REF i               
HIERARCH ESO INS REF1 DIMY = 3.579 / Size in y direction of REF i               
HIERARCH ESO INS SLIT1 TYPE = 'RECTANGLE' / Type of SLIT i                      
HIERARCH ESO INS SLIT1 ID = 135357 / ID of SLIT i                               
HIERARCH ESO INS SLIT1 OBJ RA = 150.295375 / RA (J2000) of object               
HIERARCH ESO INS SLIT1 OBJ DEC = 2.387456 / DEC (J2000) of object               
HIERARCH ESO INS SLIT1 X = -114.091 / x co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT1 Y = 0.598 / y co-ordinate of SLIT i                      
HIERARCH ESO INS SLIT1 DIMX = 11.926 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT1 DIMY = 0.582 / Size in y direction of SLIT i             
HIERARCH ESO INS SLIT2 TYPE = 'RECTANGLE' / Type of SLIT i                      
HIERARCH ESO INS SLIT2 ID = 135418 / ID of SLIT i                               
HIERARCH ESO INS SLIT2 OBJ RA = 150.294267 / RA (J2000) of object               
HIERARCH ESO INS SLIT2 OBJ DEC = 2.384864 / DEC (J2000) of object               
HIERARCH ESO INS SLIT2 X = -105.266 / x co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT2 Y = -1.718 / y co-ordinate of SLIT i                     
HIERARCH ESO INS SLIT2 DIMX = 5.247 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT2 DIMY = 0.582 / Size in y direction of SLIT i             
HIERARCH ESO INS SLIT3 TYPE = 'RECTANGLE' / Type of SLIT i                      
HIERARCH ESO INS SLIT3 ID = 134929 / ID of SLIT i                               
HIERARCH ESO INS SLIT3 OBJ RA = 150.314113 / RA (J2000) of object               
HIERARCH ESO INS SLIT3 OBJ DEC = 2.381319 / DEC (J2000) of object               
HIERARCH ESO INS SLIT3 X = -98.833 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT3 Y = 39.782 / y co-ordinate of SLIT i                     
HIERARCH ESO INS SLIT3 DIMX = 7.155 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT3 DIMY = 0.582 / Size in y direction of SLIT i             
HIERARCH ESO INS SLIT4 TYPE = 'RECTANGLE' / Type of SLIT i                      
HIERARCH ESO INS SLIT4 ID = 134657 / ID of SLIT i                               
HIERARCH ESO INS SLIT4 OBJ RA = 150.333579 / RA (J2000) of object               
HIERARCH ESO INS SLIT4 OBJ DEC = 2.377983 / DEC (J2000) of object               
HIERARCH ESO INS SLIT4 X = -91.636 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT4 Y = 80.503 / y co-ordinate of SLIT i                     
HIERARCH ESO INS SLIT4 DIMX = 6.799 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT4 DIMY = 0.582 / Size in y direction of SLIT i             
HIERARCH ESO INS SLIT5 TYPE = 'RECTANGLE' / Type of SLIT i                      
HIERARCH ESO INS SLIT5 ID = 134371 / ID of SLIT i                               
HIERARCH ESO INS SLIT5 OBJ RA = 150.297429 / RA (J2000) of object               
HIERARCH ESO INS SLIT5 OBJ DEC = 2.374828 / DEC (J2000) of object               
HIERARCH ESO INS SLIT5 X = -84.818 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT5 Y = 4.892 / y co-ordinate of SLIT i                      
HIERARCH ESO INS SLIT5 DIMX = 6.318 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT5 DIMY = 0.581 / Size in y direction of SLIT i             
HIERARCH ESO INS SLIT6 TYPE = 'RECTANGLE' / Type of SLIT i                      
HIERARCH ESO INS SLIT6 ID = 134169 / ID of SLIT i                               
HIERARCH ESO INS SLIT6 OBJ RA = 150.29875 / RA (J2000) of object                
HIERARCH ESO INS SLIT6 OBJ DEC = 2.372181 / DEC (J2000) of object               
HIERARCH ESO INS SLIT6 X = -78.739 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT6 Y = 7.657 / y co-ordinate of SLIT i                      
HIERARCH ESO INS SLIT6 DIMX = 5.364 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT6 DIMY = 0.581 / Size in y direction of SLIT i             
HIERARCH ESO INS SLIT7 TYPE = 'RECTANGLE' / Type of SLIT i                      
HIERARCH ESO INS SLIT7 ID = 133783 / ID of SLIT i                               
HIERARCH ESO INS SLIT7 OBJ RA = 150.279029 / RA (J2000) of object               
HIERARCH ESO INS SLIT7 OBJ DEC = 2.368242 / DEC (J2000) of object               
HIERARCH ESO INS SLIT7 X = -71.828 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT7 Y = -33.576 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT7 DIMX = 7.987 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT7 DIMY = 0.582 / Size in y direction of SLIT i             
HIERARCH ESO INS SLIT8 TYPE = 'RECTANGLE' / Type of SLIT i                      
HIERARCH ESO INS SLIT8 ID = 133480 / ID of SLIT i                               
HIERARCH ESO INS SLIT8 OBJ RA = 150.338367 / RA (J2000) of object               
HIERARCH ESO INS SLIT8 OBJ DEC = 2.364314 / DEC (J2000) of object               
HIERARCH ESO INS SLIT8 X = -63.618 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT8 Y = 90.513 / y co-ordinate of SLIT i                     
HIERARCH ESO INS SLIT8 DIMX = 7.988 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT8 DIMY = 0.582 / Size in y direction of SLIT i             
HIERARCH ESO INS SLIT9 TYPE = 'RECTANGLE' / Type of SLIT i                      
HIERARCH ESO INS SLIT9 ID = 133163 / ID of SLIT i                               
HIERARCH ESO INS SLIT9 OBJ RA = 150.336546 / RA (J2000) of object               
HIERARCH ESO INS SLIT9 OBJ DEC = 2.36105 / DEC (J2000) of object                
HIERARCH ESO INS SLIT9 X = -56.106 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT9 Y = 86.692 / y co-ordinate of SLIT i                     
HIERARCH ESO INS SLIT9 DIMX = 6.557 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT9 DIMY = 0.582 / Size in y direction of SLIT i             
HIERARCH ESO INS SLIT10 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT10 ID = 133098 / ID of SLIT i                              
HIERARCH ESO INS SLIT10 OBJ RA = 150.292 / RA (J2000) of object                 
HIERARCH ESO INS SLIT10 OBJ DEC = 2.358644 / DEC (J2000) of object              
HIERARCH ESO INS SLIT10 X = -50.193 / x co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT10 Y = -6.454 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT10 DIMX = 4.767 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT10 DIMY = 0.581 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT11 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT11 ID = 132633 / ID of SLIT i                              
HIERARCH ESO INS SLIT11 OBJ RA = 150.276625 / RA (J2000) of object              
HIERARCH ESO INS SLIT11 OBJ DEC = 2.355419 / DEC (J2000) of object              
HIERARCH ESO INS SLIT11 X = -44.294 / x co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT11 Y = -38.594 / y co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT11 DIMX = 6.555 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT11 DIMY = 0.581 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT12 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT12 ID = 132419 / ID of SLIT i                              
HIERARCH ESO INS SLIT12 OBJ RA = 150.280325 / RA (J2000) of object              
HIERARCH ESO INS SLIT12 OBJ DEC = 2.351606 / DEC (J2000) of object              
HIERARCH ESO INS SLIT12 X = -35.117 / x co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT12 Y = -30.854 / y co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT12 DIMX = 11.322 / Size in x direction of SLIT i           
HIERARCH ESO INS SLIT12 DIMY = 0.581 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT13 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT13 ID = 132048 / ID of SLIT i                              
HIERARCH ESO INS SLIT13 OBJ RA = 150.254533 / RA (J2000) of object              
HIERARCH ESO INS SLIT13 OBJ DEC = 2.343506 / DEC (J2000) of object              
HIERARCH ESO INS SLIT13 X = -20.753 / x co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT13 Y = -84.791 / y co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT13 DIMX = 9.297 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT13 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT14 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT14 ID = 131393 / ID of SLIT i                              
HIERARCH ESO INS SLIT14 OBJ RA = 150.300879 / RA (J2000) of object              
HIERARCH ESO INS SLIT14 OBJ DEC = 2.340411 / DEC (J2000) of object              
HIERARCH ESO INS SLIT14 X = -12.715 / x co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT14 Y = 12.107 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT14 DIMX = 6.315 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT14 DIMY = 0.581 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT15 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT15 ID = 131127 / ID of SLIT i                              
HIERARCH ESO INS SLIT15 OBJ RA = 150.265312 / RA (J2000) of object              
HIERARCH ESO INS SLIT15 OBJ DEC = 2.337044 / DEC (J2000) of object              
HIERARCH ESO INS SLIT15 X = -5.918 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT15 Y = -62.244 / y co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT15 DIMX = 6.793 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT15 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT16 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT16 ID = 131013 / ID of SLIT i                              
HIERARCH ESO INS SLIT16 OBJ RA = 150.283083 / RA (J2000) of object              
HIERARCH ESO INS SLIT16 OBJ DEC = 2.334078 / DEC (J2000) of object              
HIERARCH ESO INS SLIT16 X = 0.693 / x co-ordinate of SLIT i                     
HIERARCH ESO INS SLIT16 Y = -25.086 / y co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT16 DIMX = 5.958 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT16 DIMY = 0.581 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT17 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT17 ID = 130497 / ID of SLIT i                              
HIERARCH ESO INS SLIT17 OBJ RA = 150.326658 / RA (J2000) of object              
HIERARCH ESO INS SLIT17 OBJ DEC = 2.330731 / DEC (J2000) of object              
HIERARCH ESO INS SLIT17 X = 7.242 / x co-ordinate of SLIT i                     
HIERARCH ESO INS SLIT17 Y = 66.006 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT17 DIMX = 6.673 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT17 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT18 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT18 ID = 130276 / ID of SLIT i                              
HIERARCH ESO INS SLIT18 OBJ RA = 150.348 / RA (J2000) of object                 
HIERARCH ESO INS SLIT18 OBJ DEC = 2.326914 / DEC (J2000) of object              
HIERARCH ESO INS SLIT18 X = 14.692 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT18 Y = 110.66 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT18 DIMX = 7.748 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT18 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT19 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT19 ID = 129957 / ID of SLIT i                              
HIERARCH ESO INS SLIT19 OBJ RA = 150.3202 / RA (J2000) of object                
HIERARCH ESO INS SLIT19 OBJ DEC = 2.323739 / DEC (J2000) of object              
HIERARCH ESO INS SLIT19 X = 22.018 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT19 Y = 52.509 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT19 DIMX = 6.435 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT19 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT20 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT20 ID = 129553 / ID of SLIT i                              
HIERARCH ESO INS SLIT20 OBJ RA = 150.2657 / RA (J2000) of object                
HIERARCH ESO INS SLIT20 OBJ DEC = 2.320875 / DEC (J2000) of object              
HIERARCH ESO INS SLIT20 X = 28.405 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT20 Y = -61.437 / y co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT20 DIMX = 5.84 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT20 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT21 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT21 ID = 129940 / ID of SLIT i                              
HIERARCH ESO INS SLIT21 OBJ RA = 150.358563 / RA (J2000) of object              
HIERARCH ESO INS SLIT21 OBJ DEC = 2.317225 / DEC (J2000) of object              
HIERARCH ESO INS SLIT21 X = 35.264 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT21 Y = 132.797 / y co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT21 DIMX = 7.394 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT21 DIMY = 0.583 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT22 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT22 ID = 129028 / ID of SLIT i                              
HIERARCH ESO INS SLIT22 OBJ RA = 150.302108 / RA (J2000) of object              
HIERARCH ESO INS SLIT22 OBJ DEC = 2.314742 / DEC (J2000) of object              
HIERARCH ESO INS SLIT22 X = 41.623 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT22 Y = 14.682 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT22 DIMX = 4.886 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT22 DIMY = 0.581 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT23 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT23 ID = 128901 / ID of SLIT i                              
HIERARCH ESO INS SLIT23 OBJ RA = 150.318517 / RA (J2000) of object              
HIERARCH ESO INS SLIT23 OBJ DEC = 2.31175 / DEC (J2000) of object               
HIERARCH ESO INS SLIT23 X = 47.343 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT23 Y = 48.991 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT23 DIMX = 6.079 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT23 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT24 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT24 ID = 128682 / ID of SLIT i                              
HIERARCH ESO INS SLIT24 OBJ RA = 150.324967 / RA (J2000) of object              
HIERARCH ESO INS SLIT24 OBJ DEC = 2.309175 / DEC (J2000) of object              
HIERARCH ESO INS SLIT24 X = 53.185 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT24 Y = 62.484 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT24 DIMX = 5.126 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT24 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT25 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT25 ID = 128528 / ID of SLIT i                              
HIERARCH ESO INS SLIT25 OBJ RA = 150.324742 / RA (J2000) of object              
HIERARCH ESO INS SLIT25 OBJ DEC = 2.305881 / DEC (J2000) of object              
HIERARCH ESO INS SLIT25 X = 59.324 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT25 Y = 62.007 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT25 DIMX = 6.676 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT25 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT26 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT26 ID = 128053 / ID of SLIT i                              
HIERARCH ESO INS SLIT26 OBJ RA = 150.301421 / RA (J2000) of object              
HIERARCH ESO INS SLIT26 OBJ DEC = 2.302322 / DEC (J2000) of object              
HIERARCH ESO INS SLIT26 X = 66.473 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT26 Y = 13.246 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT26 DIMX = 7.152 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT26 DIMY = 0.581 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT27 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT27 ID = 128133 / ID of SLIT i                              
HIERARCH ESO INS SLIT27 OBJ RA = 150.313362 / RA (J2000) of object              
HIERARCH ESO INS SLIT27 OBJ DEC = 2.299283 / DEC (J2000) of object              
HIERARCH ESO INS SLIT27 X = 73.388 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT27 Y = 38.213 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT27 DIMX = 6.199 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT27 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT28 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT28 ID = 127508 / ID of SLIT i                              
HIERARCH ESO INS SLIT28 OBJ RA = 150.323738 / RA (J2000) of object              
HIERARCH ESO INS SLIT28 OBJ DEC = 2.294797 / DEC (J2000) of object              
HIERARCH ESO INS SLIT28 X = 81.321 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT28 Y = 59.917 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT28 DIMX = 9.182 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT28 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT29 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT29 ID = 127314 / ID of SLIT i                              
HIERARCH ESO INS SLIT29 OBJ RA = 150.296713 / RA (J2000) of object              
HIERARCH ESO INS SLIT29 OBJ DEC = 2.291669 / DEC (J2000) of object              
HIERARCH ESO INS SLIT29 X = 89.244 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT29 Y = 3.402 / y co-ordinate of SLIT i                     
HIERARCH ESO INS SLIT29 DIMX = 6.2 / Size in x direction of SLIT i              
HIERARCH ESO INS SLIT29 DIMY = 0.581 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT30 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT30 ID = 127180 / ID of SLIT i                              
HIERARCH ESO INS SLIT30 OBJ RA = 150.317087 / RA (J2000) of object              
HIERARCH ESO INS SLIT30 OBJ DEC = 2.288025 / DEC (J2000) of object              
HIERARCH ESO INS SLIT30 X = 96.283 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT30 Y = 46.016 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT30 DIMX = 7.394 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT30 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT31 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT31 ID = 127438 / ID of SLIT i                              
HIERARCH ESO INS SLIT31 OBJ RA = 150.321308 / RA (J2000) of object              
HIERARCH ESO INS SLIT31 OBJ DEC = 2.28495 / DEC (J2000) of object               
HIERARCH ESO INS SLIT31 X = 103.324 / x co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT31 Y = 54.842 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT31 DIMX = 6.203 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT31 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT32 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT32 ID = 126585 / ID of SLIT i                              
HIERARCH ESO INS SLIT32 OBJ RA = 150.324013 / RA (J2000) of object              
HIERARCH ESO INS SLIT32 OBJ DEC = 2.282417 / DEC (J2000) of object              
HIERARCH ESO INS SLIT32 X = 109.233 / x co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT32 Y = 60.508 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT32 DIMX = 5.13 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT32 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT33 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT33 ID = 126275 / ID of SLIT i                              
HIERARCH ESO INS SLIT33 OBJ RA = 150.301817 / RA (J2000) of object              
HIERARCH ESO INS SLIT33 OBJ DEC = 2.279189 / DEC (J2000) of object              
HIERARCH ESO INS SLIT33 X = 117.03 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT33 Y = 14.079 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT33 DIMX = 10.02 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT33 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS PIXSCALE = 0.205 / Pixel scale                                 
HIERARCH ESO INS MODE = 'MOS     ' / Instrument Observation Mode                
HIERARCH ESO INS FILT1 FIX = F / Filter mode degraded                           
HIERARCH ESO INS FILT1 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS FILT1 NAME = 'GG475   ' / Filter name                          
HIERARCH ESO INS FILT1 ID = 'vm-GG475-1.10' / Filter ESO ID                     
HIERARCH ESO INS MASK1 FIX = F / Mask mode degraded                             
HIERARCH ESO INS MASK1 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS MASK1 NO =  3 / Mask set number                                
HIERARCH ESO INS FOCU1 FIX = F / Focusing mode degraded                         
HIERARCH ESO INS FOCU1 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS FOCU1 VALUE = 6976.122982 / Camera lens position(micron)       
HIERARCH ESO INS FOCU3 FIX = F / Focusing mode degraded                         
HIERARCH ESO INS FOCU3 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS FOCU3 TEMP = 14.629684 / Temperature(C) for the beam           
HIERARCH ESO INS FOCU3 VALUE = 5181.099499 / Camera lens position(micron)       
HIERARCH ESO INS FOCU2 FIX = F / Focusing mode degraded                         
HIERARCH ESO INS FOCU2 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS FOCU2 TEMP = 14.841312 / Temperature(C) for the beam           
HIERARCH ESO INS FOCU2 VALUE = 6081.203659 / Camera lens position(micron)       
HIERARCH ESO INS MSHU1 FIX = F / Mask Shutter mode degraded                     
HIERARCH ESO INS MSHU1 SIM = T / Device Hardware Simulation                     
HIERARCH ESO INS MSHU1 MODE = 'OFF     ' / Mask Shutter mode                    
HIERARCH ESO INS MSHU1 POSH = 7.348085 / High Mask Shutters positions (mm)      
HIERARCH ESO INS MSHU1 POSL = 7.348085 / Low Mask Shutters position (mm)        
HIERARCH ESO INS GRIS1 FIX = F / Grism mode degraded                            
HIERARCH ESO INS GRIS1 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS GRIS1 NAME = 'HR_red  ' / Grism name                           
HIERARCH ESO INS GRIS1 ID = 'vm-HRredHolog-1.1' / Grism ESO ID                  
HIERARCH ESO INS AFC11 FIX = F / AFC mode degraded                              
HIERARCH ESO INS AFC11 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS AFC12 FIX = F / AFC mode degraded                              
HIERARCH ESO INS AFC12 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS FOCU4 FIX = F / Focusing mode degraded                         
HIERARCH ESO INS FOCU4 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS FOCU4 TEMP = 14.588951 / Temperature(C) for the beam           
HIERARCH ESO INS FOCU4 VALUE = 7108.120407 / Camera lens position(micron)       
HIERARCH ESO INS LAMP1 FIX = F / Lamp[i] mode degraded                          
HIERARCH ESO INS LAMP1 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS LAMP1 TIME = 0 / Lamp lightening time                          
HIERARCH ESO INS LAMP1 STATE = 'OFF     ' / State for this lamp in each box     
HIERARCH ESO INS LAMP1 NAME = 'QTH50   ' / Calibration lamp name                
HIERARCH ESO INS LAMP1 ID = 'vm-QTH50-1' / Calib. lamp ESO identification       
HIERARCH ESO INS LAMP2 FIX = F / Lamp[i] mode degraded                          
HIERARCH ESO INS LAMP2 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS LAMP2 TIME = 0 / Lamp lightening time                          
HIERARCH ESO INS LAMP2 STATE = 'OFF     ' / State for this lamp in each box     
HIERARCH ESO INS LAMP2 NAME = 'QTH10   ' / Calibration lamp name                
HIERARCH ESO INS LAMP2 ID = 'vm-QTH10-2' / Calib. lamp ESO identification       
HIERARCH ESO INS IFUM1 FIX = F / IFU Mask mode degraded                         
HIERARCH ESO INS IFUM1 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS IFUM1 MODE = 'OFF     ' / IFU Mask mode                        
HIERARCH ESO INS LAMP3 FIX = F / Lamp[i] mode degraded                          
HIERARCH ESO INS LAMP3 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS LAMP3 TIME = 0 / Lamp lightening time                          
HIERARCH ESO INS LAMP3 STATE = 'OFF     ' / State for this lamp in each box     
HIERARCH ESO INS LAMP3 NAME = 'He      ' / Calibration lamp name                
HIERARCH ESO INS LAMP3 ID = 'vm-He-3 ' / Calib. lamp ESO identification         
HIERARCH ESO INS IFUM2 FIX = F / IFU Mask mode degraded                         
HIERARCH ESO INS IFUM2 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS IFUM2 MODE = 'OFF     ' / IFU Mask mode                        
HIERARCH ESO INS LAMP4 FIX = F / Lamp[i] mode degraded                          
HIERARCH ESO INS LAMP4 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS LAMP4 TIME = 0 / Lamp lightening time                          
HIERARCH ESO INS LAMP4 STATE = 'OFF     ' / State for this lamp in each box     
HIERARCH ESO INS LAMP4 NAME = 'Ne      ' / Calibration lamp name                
HIERARCH ESO INS LAMP4 ID = 'vm-Ne-4 ' / Calib. lamp ESO identification         
HIERARCH ESO INS LAMP5 FIX = F / Lamp[i] mode degraded                          
HIERARCH ESO INS LAMP5 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS LAMP5 TIME = 0 / Lamp lightening time                          
HIERARCH ESO INS LAMP5 STATE = 'OFF     ' / State for this lamp in each box     
HIERARCH ESO INS LAMP5 NAME = 'Ar      ' / Calibration lamp name                
HIERARCH ESO INS LAMP5 ID = 'vm-Ar-5 ' / Calib. lamp ESO identification         
HIERARCH ESO INS LAMP6 FIX = F / Lamp[i] mode degraded                          
HIERARCH ESO INS LAMP6 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS LAMP6 TIME = 0 / Lamp lightening time                          
HIERARCH ESO INS LAMP6 STATE = 'OFF     ' / State for this lamp in each box     
HIERARCH ESO INS LAMP6 NAME = 'Afc     ' / Calibration lamp name                
HIERARCH ESO INS LAMP6 ID = 'vm-Afc-6' / Calib. lamp ESO identification         
HIERARCH ESO INS IFUE FIX =  F / IFU Elongator mode degraded                    
HIERARCH ESO INS IFUE SIM =  F / Device Hardware Simulation                     
HIERARCH ESO INS IFUE MAG = '1:1     ' / IFU Magnification 0.33(1:2) or 0.6     
HIERARCH ESO INS IFUS FIX =  F / IFU Shutter mode degraded                      
HIERARCH ESO INS IFUS SIM =  F / Device Hardware Simulation                     
HIERARCH ESO INS IFUS MODE = 'OFF     ' / IFU Shutter CLOSED(ON) or OPEN(OFF    
HIERARCH ESO DET CHIP1 INDEX = 1 / Chip index                                   
HIERARCH ESO DET CHIP1 ID = 'BRIAN   ' / Detector chip identification           
HIERARCH ESO DET CHIP1 NAME = 'EEV CCD-44-82' / Detector chip name              
HIERARCH ESO DET CHIP1 DATE = '2000-08-30' / Date of installation [YYYY-MM-DD]  
HIERARCH ESO DET CHIP1 X =   1 / X location in array                            
HIERARCH ESO DET CHIP1 Y =   1 / Y location in array                            
HIERARCH ESO DET CHIP1 NX = 2048 / # of pixels along X                          
HIERARCH ESO DET CHIP1 NY = 4096 / # of pixels along Y                          
HIERARCH ESO DET CHIP1 PSZX = 15. / Size of pixel in X                          
HIERARCH ESO DET CHIP1 PSZY = 15. / Size of pixel in Y                          
HIERARCH ESO DET CHIP1 XGAP = 0. / Gap between chips along x                    
HIERARCH ESO DET CHIP1 YGAP = 0. / Gap between chips along y                    
HIERARCH ESO DET OUT1 INDEX = 1 / Output index                                  
HIERARCH ESO DET OUT1 ID = 'L       ' / Output ID as from manufacturer          
HIERARCH ESO DET OUT1 NAME = 'L       ' / Description of output                 
HIERARCH ESO DET OUT1 CHIP = 1 / index of chip it belongs to                    
HIERARCH ESO DET OUT1 X =    1 / X location of output                           
HIERARCH ESO DET OUT1 Y =    1 / Y location of output                           
HIERARCH ESO DET OUT1 NX = 2048 / valid pixels along X                          
HIERARCH ESO DET OUT1 NY = 4096 / valid pixels along Y                          
HIERARCH ESO DET OUT1 CONAD = 0.55 / Conversion from ADUs to electrons          
HIERARCH ESO DET OUT1 RON = 3.21 / [e-] Readout noise per output (ele           
HIERARCH ESO DET OUT1 GAIN = 1.82 / Conversion from electrons to ADU            
HIERARCH ESO DET ID = 'CCD FIERA - Rev: 238334' / Detector system Id            
HIERARCH ESO DET NAME = 'Vimosb - vimosb' / Name of detector system             
HIERARCH ESO DET DATE = '2010-07-01' / Installation date                        
HIERARCH ESO DET DID = 'ESO-VLT-DIC.CCDDCS,ESO-VLT-DIC.FCDDCS' / Diction        
HIERARCH ESO DET BITS =     16 / Bits per pixel readout                         
HIERARCH ESO DET RA =       0. / Apparent 00:00:00.0 RA at start                
HIERARCH ESO DET DEC =      0. / Apparent 00:00:00.0 DEC at start               
HIERARCH ESO DET CHIPS =     1 / # of chips in detector array                   
HIERARCH ESO DET OUTPUTS =   1 / # of outputs                                   
HIERARCH ESO DET OUTREF =    0 / reference output                               
HIERARCH ESO DET WINDOWS =   1 / # of windows readout                           
HIERARCH ESO DET SOFW MODE = 'Normal  ' / CCD sw operational mode               
HIERARCH ESO DET EXP NO = 42197 / Unique exposure ID number                     
HIERARCH ESO DET EXP TYPE = 'Normal  ' / Exposure type                          
HIERARCH ESO DET EXP RDTTIME = 30.939 / image readout time                      
HIERARCH ESO DET EXP XFERTIM = 44.981 / image transfer time                     
HIERARCH ESO DET WIN1 ST =   T / If T, window enabled                           
HIERARCH ESO DET WIN1 STRX = 1 / Lower left pixel in X                          
HIERARCH ESO DET WIN1 STRY = 1 / Lower left pixel in Y                          
HIERARCH ESO DET WIN1 NX = 2148 / # of pixels along X                           
HIERARCH ESO DET WIN1 NY = 4096 / # of pixels along Y                           
HIERARCH ESO DET WIN1 BINX = 1 / Binning factor along X                         
HIERARCH ESO DET WIN1 BINY = 1 / Binning factor along Y                         
HIERARCH ESO DET WIN1 NDIT = 1 / # of subintegrations                           
HIERARCH ESO DET WIN1 UIT1 = 1200. / user defined subintegration time           
HIERARCH ESO DET WIN1 DIT1 = 1199.997383 / actual subintegration time           
HIERARCH ESO DET WIN1 DKTM = 1200.0992 / Dark current time                      
HIERARCH ESO DET READ MODE = 'normal  ' / Readout method                        
HIERARCH ESO DET READ SPEED = 'medium  ' / Readout speed                        
HIERARCH ESO DET READ CLOCK = '225kHz,1port,HG' / Readout clock pattern used    
HIERARCH ESO DET READ NFRAM = 1 / Number of readouts buffered in sin            
HIERARCH ESO DET FRAM ID =   1 / Image sequencial number                        
HIERARCH ESO DET FRAM TYPE = 'Normal  ' / Type of frame                         
HIERARCH ESO DET SHUT TYPE = 'seso    ' / type of shutter                       
HIERARCH ESO DET SHUT ID = 'VIMOS shutter' / Shutter unique identifier          
HIERARCH ESO DET SHUT TMOPEN = 0.074 / Time taken to open shutter               
HIERARCH ESO DET SHUT TMCLOS = 0.069 / Time taken to close shutter              
HIERARCH ESO DET TELE INT = 900. / Interval between two successive te           
HIERARCH ESO DET TELE NO =   7 / # of sources active                            
HIERARCH ESO DET TLM1 NAME = 'CCD Cold Plate' / Description of telemetry param  
HIERARCH ESO DET TLM1 ID = 'CCD_B1_T1' / ID of telemetry sensor                 
HIERARCH ESO DET TLM1 START = 140.1 / Telemetry value at read start             
HIERARCH ESO DET TLM1 END = 139.9 / Telemetry value at read completion          
HIERARCH ESO DET TLM2 NAME = 'CCD Cold Plate2' / Description of telemetry para  
HIERARCH ESO DET TLM2 ID = 'CCD_B1_T2' / ID of telemetry sensor                 
HIERARCH ESO DET TLM2 START = 151.6 / Telemetry value at read start             
HIERARCH ESO DET TLM2 END = 151.4 / Telemetry value at read completion          
HIERARCH ESO DET TLM3 NAME = 'CCD Cold Plate' / Description of telemetry param  
HIERARCH ESO DET TLM3 ID = 'CCD_B4_T1' / ID of telemetry sensor                 
HIERARCH ESO DET TLM3 START = 140. / Telemetry value at read start              
HIERARCH ESO DET TLM3 END = 140. / Telemetry value at read completion           
HIERARCH ESO DET TLM4 NAME = 'CCD Cold Plate2' / Description of telemetry para  
HIERARCH ESO DET TLM4 ID = 'CCD_B4_T2' / ID of telemetry sensor                 
HIERARCH ESO DET TLM4 START = 140.6 / Telemetry value at read start             
HIERARCH ESO DET TLM4 END = 140.6 / Telemetry value at read completion          
HIERARCH ESO DET TLM5 NAME = 'FBOX    ' / Description of telemetry param.       
HIERARCH ESO DET TLM5 ID = 'FBOX    ' / ID of telemetry sensor                  
HIERARCH ESO DET TLM5 START = 288.8 / Telemetry value at read start             
HIERARCH ESO DET TLM5 END = 288.7 / Telemetry value at read completion          
HIERARCH ESO DET TLM6 NAME = 'Vaccum_B1' / Description of telemetry param.      
HIERARCH ESO DET TLM6 ID = 'Vaccum_B1' / ID of telemetry sensor                 
HIERARCH ESO DET TLM6 START = 0. / Telemetry value at read start                
HIERARCH ESO DET TLM6 END = 0. / Telemetry value at read completion             
HIERARCH ESO DET TLM7 NAME = 'Vaccum_B4' / Description of telemetry param.      
HIERARCH ESO DET TLM7 ID = 'Vaccum_B4' / ID of telemetry sensor                 
HIERARCH ESO DET TLM7 START = 0. / Telemetry value at read start                
HIERARCH ESO DET TLM7 END = 0. / Telemetry value at read completion             
HIERARCH ESO PRO EXPTTOT = 4799.988                                             
HIERARCH ESO PRO ANCESTOR = 'VIMOS.2015-01-18T04:03:48.110.fits'                
HIERARCH ESO PRO DID = 'PRO-1.15' / Data dictionary for PRO                     
HIERARCH ESO PRO CATG = 'MOS_SCIENCE_SKY_EXTRACTED' / Category of pipeline produ
HIERARCH ESO PRO TYPE = 'REDUCED ' / Product type                               
HIERARCH ESO PRO TECH = 'MOS     ' / Observation technique                      
HIERARCH ESO PRO SCIENCE =   T / Scientific product if T                        
HIERARCH ESO PRO REC1 ID = 'vmmosscience' / Pipeline recipe (unique) identifier 
HIERARCH ESO PRO REC1 DRS ID = 'cpl-6.5.1' / Data Reduction System identifier   
HIERARCH ESO PRO REC1 PIPE ID = 'vimos-2.9.15' / Pipeline (unique) identifier   
HIERARCH ESO PRO REC1 RAW1 NAME = 'VIMOS.2015-01-18T04:03:48.110.fits' / File na
HIERARCH ESO PRO REC1 RAW1 CATG = 'MOS_SCIENCE' / Category of raw frame         
HIERARCH ESO PRO REC1 RAW2 NAME = 'VIMOS.2015-01-18T04:24:43.107.fits' / File na
HIERARCH ESO PRO REC1 RAW2 CATG = 'MOS_SCIENCE' / Category of raw frame         
HIERARCH ESO PRO REC1 RAW3 NAME = 'VIMOS.2015-01-18T04:45:38.102.fits' / File na
HIERARCH ESO PRO REC1 RAW3 CATG = 'MOS_SCIENCE' / Category of raw frame         
HIERARCH ESO PRO REC1 RAW4 NAME = 'VIMOS.2015-01-18T05:06:33.097.fits' / File na
HIERARCH ESO PRO REC1 RAW4 CATG = 'MOS_SCIENCE' / Category of raw frame         
HIERARCH ESO PRO DATANCOM =  4 / Number of combined frames                      
HIERARCH ESO PRO REC1 CAL1 NAME = 'M.VIMOS.2010-08-19T16:16:17.563.fits' / File 
HIERARCH ESO PRO REC1 CAL1 CATG = 'CONFIG_TABLE' / Category of calibration frame
HIERARCH ESO PRO REC1 CAL1 DATAMD5 = '827358dea0214d1587a3f99820ec5810' / MD5 si
HIERARCH ESO PRO REC1 CAL2 NAME = 'master_bias.fits' / File name of calibration 
HIERARCH ESO PRO REC1 CAL2 CATG = 'MASTER_BIAS' / Category of calibration frame 
HIERARCH ESO PRO REC1 CAL2 DATAMD5 = 'e61a66abf72ac06c6b41c8c3770404e3' / MD5 si
HIERARCH ESO PRO REC1 CAL3 NAME = 'mos_curv_coeff.fits' / File name of calibrati
HIERARCH ESO PRO REC1 CAL3 CATG = 'MOS_CURV_COEFF' / Category of calibration fra
HIERARCH ESO PRO REC1 CAL3 DATAMD5 = '8c4b8da9b95016a47162e7a44741b6e9' / MD5 si
HIERARCH ESO PRO REC1 CAL4 NAME = 'mos_disp_coeff.fits' / File name of calibrati
HIERARCH ESO PRO REC1 CAL4 CATG = 'MOS_DISP_COEFF' / Category of calibration fra
HIERARCH ESO PRO REC1 CAL4 DATAMD5 = '1c1a4eb0e8edb0d621c99d716b0c63ae' / MD5 si
HIERARCH ESO PRO REC1 CAL5 NAME = 'mos_master_screen_flat.fits' / File name of c
HIERARCH ESO PRO REC1 CAL5 CATG = 'MOS_MASTER_SCREEN_FLAT' / Category of calibra
HIERARCH ESO PRO REC1 CAL5 DATAMD5 = '0f2d03e3af5e30f7485a0a300f70b02f' / MD5 si
HIERARCH ESO PRO REC1 CAL6 NAME = 'mos_slit_location.fits' / File name of calibr
HIERARCH ESO PRO REC1 CAL6 CATG = 'MOS_SLIT_LOCATION' / Category of calibration 
HIERARCH ESO PRO REC1 CAL6 DATAMD5 = 'f53c43722c662a5763c5b94db9b2468e' / MD5 si
HIERARCH ESO PRO REC1 PARAM1 NAME = 'dispersion' / Resampling step (Angstrom/pix
HIERARCH ESO PRO REC1 PARAM1 VALUE = '0.6     ' / Default: 0                    
HIERARCH ESO PRO REC1 PARAM2 NAME = 'skyalign' / Polynomial order for sky lines 
HIERARCH ESO PRO REC1 PARAM2 VALUE = '0       ' / Default: 0                    
HIERARCH ESO PRO REC1 PARAM3 NAME = 'wcolumn ' / Name of sky line catalog table 
HIERARCH ESO PRO REC1 PARAM3 VALUE = 'WLEN    ' / Default: 'WLEN'               
HIERARCH ESO PRO REC1 PARAM4 NAME = 'startwavelength' / Start wavelength in spec
HIERARCH ESO PRO REC1 PARAM4 VALUE = '5800    ' / Default: 0                    
HIERARCH ESO PRO REC1 PARAM5 NAME = 'endwavelength' / End wavelength in spectral
HIERARCH ESO PRO REC1 PARAM5 VALUE = '9500    ' / Default: 0                    
HIERARCH ESO PRO REC1 PARAM6 NAME = 'reference' / Reference wavelength for calib
HIERARCH ESO PRO REC1 PARAM6 VALUE = '7948.18 ' / Default: 0                    
HIERARCH ESO PRO REC1 PARAM7 NAME = 'flux    ' / Apply flux conservation        
HIERARCH ESO PRO REC1 PARAM7 VALUE = 'true    ' / Default: true                 
HIERARCH ESO PRO REC1 PARAM8 NAME = 'flatfield' / Apply flat field              
HIERARCH ESO PRO REC1 PARAM8 VALUE = 'true    ' / Default: true                 
HIERARCH ESO PRO REC1 PARAM9 NAME = 'skyglobal' / Subtract global sky spectrum f
HIERARCH ESO PRO REC1 PARAM9 VALUE = 'false   ' / Default: false                
HIERARCH ESO PRO REC1 PARAM10 NAME = 'skymedian' / Sky subtraction from extracte
HIERARCH ESO PRO REC1 PARAM10 VALUE = 'false   ' / Default: false               
HIERARCH ESO PRO REC1 PARAM11 NAME = 'skylocal' / Sky subtraction from CCD slit 
HIERARCH ESO PRO REC1 PARAM11 VALUE = 'true    ' / Default: true                
HIERARCH ESO PRO REC1 PARAM12 NAME = 'cosmics ' / Eliminate cosmic rays hits (on
HIERARCH ESO PRO REC1 PARAM12 VALUE = 'true    ' / Default: true                
HIERARCH ESO PRO REC1 PARAM13 NAME = 'slit_margin' / Number of pixels to exclude
HIERARCH ESO PRO REC1 PARAM13 VALUE = '3       ' / Default: 3                   
HIERARCH ESO PRO REC1 PARAM14 NAME = 'ext_radius' / Maximum extraction radius fo
HIERARCH ESO PRO REC1 PARAM14 VALUE = '6       ' / Default: 6                   
HIERARCH ESO PRO REC1 PARAM15 NAME = 'cont_radius' / Minimum distance at which t
HIERARCH ESO PRO REC1 PARAM15 VALUE = '0       ' / Default: 0                   
HIERARCH ESO PRO REC1 PARAM16 NAME = 'ext_mode' / Object extraction method: 0 = 
HIERARCH ESO PRO REC1 PARAM16 VALUE = '1       ' / Default: 1                   
HIERARCH ESO PRO REC1 PARAM17 NAME = 'detection' / Object detection threshold (A
HIERARCH ESO PRO REC1 PARAM17 VALUE = '2       ' / Default: 2                   
HIERARCH ESO PRO REC1 PARAM18 NAME = 'time_normalise' / Normalise output spectra
HIERARCH ESO PRO REC1 PARAM18 VALUE = 'true    ' / Default: true                
HIERARCH ESO PRO REC1 PARAM19 NAME = 'anyframe' / Look for a standard star in an
HIERARCH ESO PRO REC1 PARAM19 VALUE = 'false   ' / Default: false               
HIERARCH ESO PRO REC1 PARAM20 NAME = 'response' / Order of polynomial modeling t
HIERARCH ESO PRO REC1 PARAM20 VALUE = '5       ' / Default: 5                   
HIERARCH ESO PRO REC1 PARAM21 NAME = 'alignment' / Type of alignment of dithered
HIERARCH ESO PRO REC1 PARAM21 VALUE = 'integer ' / Default: 'integer'           
HIERARCH ESO PRO REC1 PARAM22 NAME = 'stack_method' / Frames combination method 
HIERARCH ESO PRO REC1 PARAM22 VALUE = 'average ' / Default: 'average'           
HIERARCH ESO PRO REC1 PARAM23 NAME = 'minrejection' / Number of lowest values to
HIERARCH ESO PRO REC1 PARAM23 VALUE = '1       ' / Default: 1                   
HIERARCH ESO PRO REC1 PARAM24 NAME = 'maxrejection' / Number of highest values t
HIERARCH ESO PRO REC1 PARAM24 VALUE = '1       ' / Default: 1                   
HIERARCH ESO PRO REC1 PARAM25 NAME = 'klow    ' / Low threshold in ksigma method
HIERARCH ESO PRO REC1 PARAM25 VALUE = '3       ' / Default: 3                   
HIERARCH ESO PRO REC1 PARAM26 NAME = 'khigh   ' / High threshold in ksigma metho
HIERARCH ESO PRO REC1 PARAM26 VALUE = '3       ' / Default: 3                   
HIERARCH ESO PRO REC1 PARAM27 NAME = 'kiter   ' / Max number of iterations in ks
HIERARCH ESO PRO REC1 PARAM27 VALUE = '999     ' / Default: 999                 
HIERARCH ESO PRO REC1 PARAM28 NAME = 'dither  ' / Align dithered frames before s
HIERARCH ESO PRO REC1 PARAM28 VALUE = 'true    ' / Default: true                
HIERARCH ESO PRO REC1 PARAM29 NAME = 'compute ' / Compute offsets of dithered im
HIERARCH ESO PRO REC1 PARAM29 VALUE = 'false   ' / Default: false               
HIERARCH ESO PRO REC1 PARAM30 NAME = 'fringing' / Apply fringing correction (onl
HIERARCH ESO PRO REC1 PARAM30 VALUE = 'true    ' / Default: true                
HIERARCH ESO PRO REC1 PARAM31 NAME = 'offset  ' / Minimum required offset betwee
HIERARCH ESO PRO REC1 PARAM31 VALUE = '3       ' / Default: 3                   
HIERARCH ESO PRO REC1 PARAM32 NAME = 'qc      ' / Compute QC1 parameters        
HIERARCH ESO PRO REC1 PARAM32 VALUE = 'true    ' / Default: true                
HIERARCH ESO ADA ABSROT START = 21.56262 / Abs rot angle at exp start           
HIERARCH ESO ADA POSANG =  90. / Position angle at start                        
HIERARCH ESO ADA GUID STATUS = 'ON      ' / Status of autoguider                
HIERARCH ESO ADA GUID RA = 150.241537 / 10:00:57.9 Guide star RA J2000          
HIERARCH ESO ADA GUID DEC = 2.59418 / 02:35:39.0 Guide star DEC J2000           
HIERARCH ESO ADA ABSROT PPOS = 'NEG     ' / sign of probe position              
HIERARCH ESO ADA ABSROT END = 20.56639 / Abs rot angle at exp end               
HIERARCH ESO OCS CON QUAD =  1 / Quadrant Number                                
HIERARCH ESO OCS CMF FILE1 NAME = 'MOS_sky2ccd_1.cmf' / Calibration matrix file 
HIERARCH ESO OCS CMF FILE1 CRTEDATE = '2012-05-18T14:37:53' / Calibration matrix
HIERARCH ESO OCS CMF FILE2 NAME = 'MOS_wavecal_HR_red_1.cmf' / Calibration matri
HIERARCH ESO OCS CMF FILE2 CRTEDATE = '2005-10-20T17:18:57' / Calibration matrix
HIERARCH ESO OCS COMP ID = '226602 2012-03-18 17:49:45Z' / OS Software Versi    
HIERARCH ESO OCS DID = 'ESO-VLT-DIC.VIMOS_OS-0.2' / Data dictionary for         
HISTORY FTU-2_5_2/2015-01-18/REMOVE: XTENSION, "IMAGE"                          
HISTORY FTU-2_5_2/2015-01-18/REMOVE: EXTNAME, "BRIAN"                           
HISTORY FTU-2_5_2/2015-01-18/REMOVE: ORIGIN, "VIRMOS OS Software"               
HISTORY FTU-2_5_2/2015-01-18/REMOVE: CDELT1, 1.00000000                         
HISTORY FTU-2_5_2/2015-01-18/REMOVE: CDELT2, 1.00000000                         
HISTORY FTU-2_5_2/2015-01-18/CHANGE: DET.CHIPS, 1                               
HISTORY FTU-2_5_2/2015-01-18/CHANGE: DET.OUTPUTS, 1                             
HISTORY FTU-2_5_2/2015-01-18/CHANGE: DET.WIN1.NX, 2148                          
COMMENT FTU-2_5_2/2015-01-18T04:24:44/vimos.ftt                                 
CHECKSUM= 'NMTOQMTNNMTNNMTN'   / HDU checksum updated 2015-01-26T12:24:13       
DATASUM = '338047765'          / data unit checksum updated 2015-01-26T12:24:13 

Developer note#

I would appreciate the interactive tools here to zoom and pan through the 2D spectrum. Hoovering to know the precise wavelength of a feature would also be very useful. With that, the interactive tool could show automatically the calibration in wavelength reading it from the header.

plt.figure(figsize=[20,20])
plt.imshow(hdu2d[0].data)
plt.xlim(2000,3000) #spec is very big, plot just a bit
(2000.0, 3000.0)
../../_images/e3147c3f1c243e122359a3b70738b4c1256fd601d04ffc39a7dc74a53cc4cfcf.png

Now work with 1D spectrum#

Calibrate (in wavelength), inspect, and write in Spectrum1D object#

hdu1d = fits.open(file1d)
hdu1dwht = fits.open(file1dwht)
hdu1d.info()
Filename: ./mos_spectroscopy/observed/legac_M1_v3.7_spec1d_130902.fits
No.    Name      Ver    Type      Cards   Dimensions   Format
  0  PRIMARY       1 PrimaryHDU     715   (6166,)   float32   
hdu1d[0].header
SIMPLE  =                    T / file does conform to FITS standard             
BITPIX  =                  -32 / number of bits per data pixel                  
NAXIS   =                    1 / number of data axes                            
NAXIS1  =                 6166 / length of data axis 1                          
EXTEND  =                    T / FITS dataset may contain extensions            
COMMENT   FITS (Flexible Image Transport System) format is defined in 'Astronomy
COMMENT   and Astrophysics', volume 376, page 359; bibcode: 2001A&A...376..359H 
DATE    = '2015-01-26T12:24:07' / file creation date (YYYY-MM-DDThh:mm:ss UT)   
ORIGIN  = 'ESO     '           / European Southern Observatory                  
MJD-OBS =       57040.16930682 / Obs start 2015-01-18T04:03:48.109              
DATE-OBS= '2015-01-18T04:03:48.108' / Date of observation                       
EXPTIME =             1199.997 / Total integration time. 00:19:59.997           
TELESCOP= 'ESO-VLT-U3'         / ESO <TEL>                                      
RA      =           150.220555 / 10:00:52.9 RA (J2000) pointing                 
DEC     =              2.41104 / 02:24:39.7  DEC (J2000) pointing               
EQUINOX =                2000. / Standard FK5                                   
RADECSYS= 'FK5     '           / Coordinate reference frame                     
LST     =            25866.379 / 07:11:06.379 LST at start                      
UTC     =               14623. / 04:03:43.000 UT at start                       
OBSERVER= 'VIMOS   '           / Name of observer                               
INSTRUME= 'VIMOS   '           / Instrument used                                
PI-COI  = '555555555'          / Name of PI and COI                             
OBJECT  = 'COSMOS_M01'         / Original target                                
CRVAL1  =               5800.3 / 10:00:52.9, RA at ref pixel                    
CRPIX1  =                   1. / Reference pixel in X                           
CTYPE1  = 'LINEAR  '           / pixel coordinate system                        
CRVAL2  =                   1. / 02:24:39.7, DEC at ref pixel                   
CRPIX2  =                   1. / Reference pixel in Y                           
CTYPE2  = 'PIXEL   '           / pixel coordinate system                        
INHERIT =                    T / denotes the INHERIT keyword convention         
CD1_1   =                  0.6 / Translation matrix element                     
CD1_2   =                   0. / Translation matrix element                     
CD2_1   =                   0. / Translation matrix element                     
CD2_2   =                   1. / Translation matrix element                     
DATAMD5 = '77a1ef86719040e2fea8d1d96fef35e0' / MD5 checksum                     
PIPEFILE= 'mos_science_sky_extracted.fits' / Filename of data product           
ARCFILE = 'VIMOS.2015-01-18T04:03:48.110.fits' / Archive File Name              
HIERARCH ESO OBS DID = 'ESO-VLT-DIC.OBS-1.12' / OBS Dictionary                  
HIERARCH ESO OBS EXECTIME = 6835 / Expected execution time                      
HIERARCH ESO OBS GRP = '0       ' / linked blocks                               
HIERARCH ESO OBS ID =  1163772 / Observation block ID                           
HIERARCH ESO OBS NAME = 'legac_M01_2' / OB name                                 
HIERARCH ESO OBS OBSERVER = 'UNKNOWN ' / Observer Name                          
HIERARCH ESO OBS PI-COI ID = 73498 / ESO internal PI-COI ID                     
HIERARCH ESO OBS PI-COI NAME = 'UNKNOWN ' / PI-COI name                         
HIERARCH ESO OBS PROG ID = '194.A-2005(B)' / ESO program identification         
HIERARCH ESO OBS START = '2015-01-18T03:38:14' / OB start time                  
HIERARCH ESO OBS NTPL =      3 / Number of templates within OB                  
HIERARCH ESO OBS TPLNO =     2 / Template number within OB                      
HIERARCH ESO OBS AIRM =     2. / Req. max. airmass                              
HIERARCH ESO OBS AMBI FWHM = 1.3 / Req. max. seeing                             
HIERARCH ESO OBS AMBI TRANS = '2CLR    ' / Req. sky transparency                
HIERARCH ESO OBS ATM = 'NATM    ' / Req. Atmospheric Turbulence Model           
HIERARCH ESO OBS CONTAINER ID = -999 / Scheduling container ID                  
HIERARCH ESO OBS CONTAINER TYPE = 'U       ' / Scheduling container type        
HIERARCH ESO OBS MOON DIST = 120 / Req. min. angular dist. from moon            
HIERARCH ESO OBS MOON FLI = 0.3 / Req. max. fractional lunar illum.             
HIERARCH ESO OBS STREHLRATIO = 0. / Req. strehl ratio                           
HIERARCH ESO OBS TARG NAME = 'COSMOS_M01' / OB target name                      
HIERARCH ESO OBS TWILIGHT = -15 / Req. twilight                                 
HIERARCH ESO OBS WATERVAPOUR = 0. / Req. water vapour                           
HIERARCH ESO TPL FILE DIRNAME = '$INS_ROOT/$INS_USER/ADF' / Storage filename    
HIERARCH ESO TPL FILE KEEP = T / paramfile keep-flag                            
HIERARCH ESO TPL DID = 'ESO-VLT-DIC.TPL-1.9' / Data dictionary for TPL          
HIERARCH ESO TPL ID = 'VIMOS_mos_obs_Offset' / Template signature ID            
HIERARCH ESO TPL NAME = 'MOS Observation' / Template name                       
HIERARCH ESO TPL PRESEQ = 'VIMOS_mos_obs_Offset' / Sequencer script             
HIERARCH ESO TPL START = '2015-01-18T04:01:12' / TPL start time                 
HIERARCH ESO TPL VERSION = '@(#) $Revision: 77996 $' / Version of the templa    
HIERARCH ESO TPL NEXP =      4 / Number of exposures within templat             
HIERARCH ESO TPL EXPNO =     1 / Exposure number within template                
HIERARCH ESO TEL DID = 'ESO-VLT-DIC.TCS' / Data dictionary for TEL              
HIERARCH ESO TEL ID = 'v 261159' / TCS version number                           
HIERARCH ESO TEL DATE = '2000-01-01T00:00:00' / TCS installation date           
HIERARCH ESO TEL ALT =  40.643 / Alt angle at start                             
HIERARCH ESO TEL AZ =  243.107 / Az angle at start S=0,W=90                     
HIERARCH ESO TEL GEOELEV = 2648. / Elevation above sea level (m)                
HIERARCH ESO TEL GEOLAT = -24.6268 / Tel geo latitute (+=North)                 
HIERARCH ESO TEL GEOLON = -70.4045 / Tel geo longitude (+=East)                 
HIERARCH ESO TEL OPER = 'I, Condor' / Telescope Operator                        
HIERARCH ESO TEL FOCU ID = 'NB      ' / Telescope focus station ID              
HIERARCH ESO TEL FOCU LEN = 120. / Focal length                                 
HIERARCH ESO TEL FOCU SCALE = 1.718 / Focal scale                               
HIERARCH ESO TEL FOCU VALUE = -31.348 / M2 setting                              
HIERARCH ESO TEL PARANG START = -125.768 / Parallactic angle at start           
HIERARCH ESO TEL AIRM START = 1.533 / Airmass at start                          
HIERARCH ESO TEL AMBI FWHM START = 1.08 / Observatory Seeing queried from AS    
HIERARCH ESO TEL AMBI PRES START = 744.4 / Observatory ambient air pressure q   
HIERARCH ESO TEL AMBI WINDSP = 0.82 / Observatory ambient wind speed que        
HIERARCH ESO TEL AMBI WINDDIR = 325.5 / Observatory ambient wind direction      
HIERARCH ESO TEL AMBI RHUM = 3. / Observatory ambient relative humid            
HIERARCH ESO TEL AMBI TEMP = 15.77 / Observatory ambient temperature qu         
HIERARCH ESO TEL MOON RA = 265.331996 / 17:41:19.6 RA (J2000)                   
HIERARCH ESO TEL MOON DEC = -17.93387 / -17:56:01.9 DEC (J2000)                 
HIERARCH ESO TEL TH M1 TEMP = 13.63 / M1 superficial temperature                
HIERARCH ESO TEL TRAK STATUS = 'NORMAL  ' / Tracking status                     
HIERARCH ESO TEL DOME STATUS = 'FULLY-OPEN' / Dome status                       
HIERARCH ESO TEL CHOP ST =   F / True when chopping is active                   
HIERARCH ESO TEL TARG ALPHA = 100052.8 / Alpha coordinate for the target        
HIERARCH ESO TEL TARG DELTA = 22436. / Delta coordinate for the target          
HIERARCH ESO TEL TARG EPOCH = 2000. / Epoch                                     
HIERARCH ESO TEL TARG EPOCHSYSTEM = 'J       ' / Epoch system (default J=Julian)
HIERARCH ESO TEL TARG EQUINOX = 2000. / Equinox                                 
HIERARCH ESO TEL TARG PMA = 0. / Proper Motion Alpha                            
HIERARCH ESO TEL TARG PMD = 0. / Proper motion Delta                            
HIERARCH ESO TEL TARG RADVEL = 0. / Radial velocity                             
HIERARCH ESO TEL TARG PARALLAX = 0. / Parallax                                  
HIERARCH ESO TEL TARG COORDTYPE = 'M       ' / Coordinate type (M=mean A=apparen
HIERARCH ESO TEL PARANG END = -128.789 / Parallactic angle at end               
HIERARCH ESO TEL AIRM END = 1.421 / Airmass at end                              
HIERARCH ESO TEL AMBI FWHM END = 1.03 / Observatory Seeing queried from AS      
HIERARCH ESO TEL IA FWHM =  1. / Delivered seeing corrected by airm             
HIERARCH ESO TEL IA FWHMLIN = 1.01 / Delivered seeing on IA detector (l         
HIERARCH ESO TEL IA FWHMLINOBS = 1.25 / Delivered seeing on IA detector (l      
HIERARCH ESO TEL AMBI PRES END = 744.3 / Observatory ambient air pressure q     
HIERARCH ESO TEL AMBI TAU0 = 0.0039 / Average coherence time                    
HIERARCH ESO TEL GS1 ALPHA = 100057.969 / Guide Star alpha                      
HIERARCH ESO TEL GS1 DELTA = 23539.05 / Guide Star delta                        
HIERARCH ESO TEL GS1 PPOS = 'NEG     ' / Selected probe position (POS or NE     
HIERARCH ESO TEL TARG NAME = '        ' / Name of special target position       
HIERARCH ESO INS SWSIM = 'NORMAL  ' / Software simulation                       
HIERARCH ESO INS ID = 'VIMOS   ' / Instrument identification                    
HIERARCH ESO INS DID = 'ESO-VLT-DIC.VIMOS_ICS-242857' / Data dictionary         
HIERARCH ESO INS SOFT DATE = 'Mar 20 2014' / Instrument SW Release date         
HIERARCH ESO INS SOFT ID = '186976  ' / Version of ICS                          
HIERARCH ESO INS FOCU1 TEMP = 14.61792 / Temperature(C) for the beam            
HIERARCH ESO INS ADF TYPE = 'MOS     ' / Type of ADF                            
HIERARCH ESO INS ADF VERSION = 3.28 / Version of ADF                            
HIERARCH ESO INS ADF UNIT = 'MILLIMETER' / Unit of slit/object coordinates      
HIERARCH ESO INS ADF ID = 'legac_M1Q1_HR_Red_M1Q1.adp' / Name for the ADF       
HIERARCH ESO INS ADF ADMID = '2015-01-18T03:38:06-101-1.adm' / Name for the *.  
HIERARCH ESO INS ADF COMMENT = '        ' / Optional user comment for ADF       
HIERARCH ESO INS MASK1 ID = 103063 / Mask identification number                 
HIERARCH ESO INS ADF GRISM NAME = 'HR_Red  ' / Grism name                       
HIERARCH ESO INS ADF GRISM SPECTLEN = 4096. / Spectrum length in pixels         
HIERARCH ESO INS ADF SKYREG = 1.8 / Sky region in arcsec                        
HIERARCH ESO INS REF NO =    1 / Number of reference apertures                  
HIERARCH ESO INS ARC NO =    0 / Number of curved apertures                     
HIERARCH ESO INS SLIT NO =  33 / Number of all apertures but refere             
HIERARCH ESO INS SHU NO =    1 / Number of shutter positions                    
HIERARCH ESO INS SHU1 POSL = 0. / Lower shutter positions                       
HIERARCH ESO INS SHU1 POSH = 2440. / Upper shutter positions                    
HIERARCH ESO INS REF1 TYPE = 'SQUARE  ' / Type of REF i                         
HIERARCH ESO INS REF1 ID = 132308 / ID of REF i                                 
HIERARCH ESO INS REF1 OBJ RA = 150.345742 / RA (J2000) of object                
HIERARCH ESO INS REF1 OBJ DEC = 2.347825 / DEC (J2000) of object                
HIERARCH ESO INS REF1 X = -27.353 / x co-ordinate of REF i                      
HIERARCH ESO INS REF1 Y = 105.937 / y co-ordinate of REF i                      
HIERARCH ESO INS REF1 DIMX = 3.576 / Size in x direction of REF i               
HIERARCH ESO INS REF1 DIMY = 3.579 / Size in y direction of REF i               
HIERARCH ESO INS SLIT1 TYPE = 'RECTANGLE' / Type of SLIT i                      
HIERARCH ESO INS SLIT1 ID = 135357 / ID of SLIT i                               
HIERARCH ESO INS SLIT1 OBJ RA = 150.295375 / RA (J2000) of object               
HIERARCH ESO INS SLIT1 OBJ DEC = 2.387456 / DEC (J2000) of object               
HIERARCH ESO INS SLIT1 X = -114.091 / x co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT1 Y = 0.598 / y co-ordinate of SLIT i                      
HIERARCH ESO INS SLIT1 DIMX = 11.926 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT1 DIMY = 0.582 / Size in y direction of SLIT i             
HIERARCH ESO INS SLIT2 TYPE = 'RECTANGLE' / Type of SLIT i                      
HIERARCH ESO INS SLIT2 ID = 135418 / ID of SLIT i                               
HIERARCH ESO INS SLIT2 OBJ RA = 150.294267 / RA (J2000) of object               
HIERARCH ESO INS SLIT2 OBJ DEC = 2.384864 / DEC (J2000) of object               
HIERARCH ESO INS SLIT2 X = -105.266 / x co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT2 Y = -1.718 / y co-ordinate of SLIT i                     
HIERARCH ESO INS SLIT2 DIMX = 5.247 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT2 DIMY = 0.582 / Size in y direction of SLIT i             
HIERARCH ESO INS SLIT3 TYPE = 'RECTANGLE' / Type of SLIT i                      
HIERARCH ESO INS SLIT3 ID = 134929 / ID of SLIT i                               
HIERARCH ESO INS SLIT3 OBJ RA = 150.314113 / RA (J2000) of object               
HIERARCH ESO INS SLIT3 OBJ DEC = 2.381319 / DEC (J2000) of object               
HIERARCH ESO INS SLIT3 X = -98.833 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT3 Y = 39.782 / y co-ordinate of SLIT i                     
HIERARCH ESO INS SLIT3 DIMX = 7.155 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT3 DIMY = 0.582 / Size in y direction of SLIT i             
HIERARCH ESO INS SLIT4 TYPE = 'RECTANGLE' / Type of SLIT i                      
HIERARCH ESO INS SLIT4 ID = 134657 / ID of SLIT i                               
HIERARCH ESO INS SLIT4 OBJ RA = 150.333579 / RA (J2000) of object               
HIERARCH ESO INS SLIT4 OBJ DEC = 2.377983 / DEC (J2000) of object               
HIERARCH ESO INS SLIT4 X = -91.636 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT4 Y = 80.503 / y co-ordinate of SLIT i                     
HIERARCH ESO INS SLIT4 DIMX = 6.799 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT4 DIMY = 0.582 / Size in y direction of SLIT i             
HIERARCH ESO INS SLIT5 TYPE = 'RECTANGLE' / Type of SLIT i                      
HIERARCH ESO INS SLIT5 ID = 134371 / ID of SLIT i                               
HIERARCH ESO INS SLIT5 OBJ RA = 150.297429 / RA (J2000) of object               
HIERARCH ESO INS SLIT5 OBJ DEC = 2.374828 / DEC (J2000) of object               
HIERARCH ESO INS SLIT5 X = -84.818 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT5 Y = 4.892 / y co-ordinate of SLIT i                      
HIERARCH ESO INS SLIT5 DIMX = 6.318 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT5 DIMY = 0.581 / Size in y direction of SLIT i             
HIERARCH ESO INS SLIT6 TYPE = 'RECTANGLE' / Type of SLIT i                      
HIERARCH ESO INS SLIT6 ID = 134169 / ID of SLIT i                               
HIERARCH ESO INS SLIT6 OBJ RA = 150.29875 / RA (J2000) of object                
HIERARCH ESO INS SLIT6 OBJ DEC = 2.372181 / DEC (J2000) of object               
HIERARCH ESO INS SLIT6 X = -78.739 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT6 Y = 7.657 / y co-ordinate of SLIT i                      
HIERARCH ESO INS SLIT6 DIMX = 5.364 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT6 DIMY = 0.581 / Size in y direction of SLIT i             
HIERARCH ESO INS SLIT7 TYPE = 'RECTANGLE' / Type of SLIT i                      
HIERARCH ESO INS SLIT7 ID = 133783 / ID of SLIT i                               
HIERARCH ESO INS SLIT7 OBJ RA = 150.279029 / RA (J2000) of object               
HIERARCH ESO INS SLIT7 OBJ DEC = 2.368242 / DEC (J2000) of object               
HIERARCH ESO INS SLIT7 X = -71.828 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT7 Y = -33.576 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT7 DIMX = 7.987 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT7 DIMY = 0.582 / Size in y direction of SLIT i             
HIERARCH ESO INS SLIT8 TYPE = 'RECTANGLE' / Type of SLIT i                      
HIERARCH ESO INS SLIT8 ID = 133480 / ID of SLIT i                               
HIERARCH ESO INS SLIT8 OBJ RA = 150.338367 / RA (J2000) of object               
HIERARCH ESO INS SLIT8 OBJ DEC = 2.364314 / DEC (J2000) of object               
HIERARCH ESO INS SLIT8 X = -63.618 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT8 Y = 90.513 / y co-ordinate of SLIT i                     
HIERARCH ESO INS SLIT8 DIMX = 7.988 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT8 DIMY = 0.582 / Size in y direction of SLIT i             
HIERARCH ESO INS SLIT9 TYPE = 'RECTANGLE' / Type of SLIT i                      
HIERARCH ESO INS SLIT9 ID = 133163 / ID of SLIT i                               
HIERARCH ESO INS SLIT9 OBJ RA = 150.336546 / RA (J2000) of object               
HIERARCH ESO INS SLIT9 OBJ DEC = 2.36105 / DEC (J2000) of object                
HIERARCH ESO INS SLIT9 X = -56.106 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT9 Y = 86.692 / y co-ordinate of SLIT i                     
HIERARCH ESO INS SLIT9 DIMX = 6.557 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT9 DIMY = 0.582 / Size in y direction of SLIT i             
HIERARCH ESO INS SLIT10 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT10 ID = 133098 / ID of SLIT i                              
HIERARCH ESO INS SLIT10 OBJ RA = 150.292 / RA (J2000) of object                 
HIERARCH ESO INS SLIT10 OBJ DEC = 2.358644 / DEC (J2000) of object              
HIERARCH ESO INS SLIT10 X = -50.193 / x co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT10 Y = -6.454 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT10 DIMX = 4.767 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT10 DIMY = 0.581 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT11 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT11 ID = 132633 / ID of SLIT i                              
HIERARCH ESO INS SLIT11 OBJ RA = 150.276625 / RA (J2000) of object              
HIERARCH ESO INS SLIT11 OBJ DEC = 2.355419 / DEC (J2000) of object              
HIERARCH ESO INS SLIT11 X = -44.294 / x co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT11 Y = -38.594 / y co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT11 DIMX = 6.555 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT11 DIMY = 0.581 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT12 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT12 ID = 132419 / ID of SLIT i                              
HIERARCH ESO INS SLIT12 OBJ RA = 150.280325 / RA (J2000) of object              
HIERARCH ESO INS SLIT12 OBJ DEC = 2.351606 / DEC (J2000) of object              
HIERARCH ESO INS SLIT12 X = -35.117 / x co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT12 Y = -30.854 / y co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT12 DIMX = 11.322 / Size in x direction of SLIT i           
HIERARCH ESO INS SLIT12 DIMY = 0.581 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT13 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT13 ID = 132048 / ID of SLIT i                              
HIERARCH ESO INS SLIT13 OBJ RA = 150.254533 / RA (J2000) of object              
HIERARCH ESO INS SLIT13 OBJ DEC = 2.343506 / DEC (J2000) of object              
HIERARCH ESO INS SLIT13 X = -20.753 / x co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT13 Y = -84.791 / y co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT13 DIMX = 9.297 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT13 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT14 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT14 ID = 131393 / ID of SLIT i                              
HIERARCH ESO INS SLIT14 OBJ RA = 150.300879 / RA (J2000) of object              
HIERARCH ESO INS SLIT14 OBJ DEC = 2.340411 / DEC (J2000) of object              
HIERARCH ESO INS SLIT14 X = -12.715 / x co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT14 Y = 12.107 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT14 DIMX = 6.315 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT14 DIMY = 0.581 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT15 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT15 ID = 131127 / ID of SLIT i                              
HIERARCH ESO INS SLIT15 OBJ RA = 150.265312 / RA (J2000) of object              
HIERARCH ESO INS SLIT15 OBJ DEC = 2.337044 / DEC (J2000) of object              
HIERARCH ESO INS SLIT15 X = -5.918 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT15 Y = -62.244 / y co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT15 DIMX = 6.793 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT15 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT16 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT16 ID = 131013 / ID of SLIT i                              
HIERARCH ESO INS SLIT16 OBJ RA = 150.283083 / RA (J2000) of object              
HIERARCH ESO INS SLIT16 OBJ DEC = 2.334078 / DEC (J2000) of object              
HIERARCH ESO INS SLIT16 X = 0.693 / x co-ordinate of SLIT i                     
HIERARCH ESO INS SLIT16 Y = -25.086 / y co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT16 DIMX = 5.958 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT16 DIMY = 0.581 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT17 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT17 ID = 130497 / ID of SLIT i                              
HIERARCH ESO INS SLIT17 OBJ RA = 150.326658 / RA (J2000) of object              
HIERARCH ESO INS SLIT17 OBJ DEC = 2.330731 / DEC (J2000) of object              
HIERARCH ESO INS SLIT17 X = 7.242 / x co-ordinate of SLIT i                     
HIERARCH ESO INS SLIT17 Y = 66.006 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT17 DIMX = 6.673 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT17 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT18 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT18 ID = 130276 / ID of SLIT i                              
HIERARCH ESO INS SLIT18 OBJ RA = 150.348 / RA (J2000) of object                 
HIERARCH ESO INS SLIT18 OBJ DEC = 2.326914 / DEC (J2000) of object              
HIERARCH ESO INS SLIT18 X = 14.692 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT18 Y = 110.66 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT18 DIMX = 7.748 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT18 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT19 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT19 ID = 129957 / ID of SLIT i                              
HIERARCH ESO INS SLIT19 OBJ RA = 150.3202 / RA (J2000) of object                
HIERARCH ESO INS SLIT19 OBJ DEC = 2.323739 / DEC (J2000) of object              
HIERARCH ESO INS SLIT19 X = 22.018 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT19 Y = 52.509 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT19 DIMX = 6.435 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT19 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT20 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT20 ID = 129553 / ID of SLIT i                              
HIERARCH ESO INS SLIT20 OBJ RA = 150.2657 / RA (J2000) of object                
HIERARCH ESO INS SLIT20 OBJ DEC = 2.320875 / DEC (J2000) of object              
HIERARCH ESO INS SLIT20 X = 28.405 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT20 Y = -61.437 / y co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT20 DIMX = 5.84 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT20 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT21 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT21 ID = 129940 / ID of SLIT i                              
HIERARCH ESO INS SLIT21 OBJ RA = 150.358563 / RA (J2000) of object              
HIERARCH ESO INS SLIT21 OBJ DEC = 2.317225 / DEC (J2000) of object              
HIERARCH ESO INS SLIT21 X = 35.264 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT21 Y = 132.797 / y co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT21 DIMX = 7.394 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT21 DIMY = 0.583 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT22 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT22 ID = 129028 / ID of SLIT i                              
HIERARCH ESO INS SLIT22 OBJ RA = 150.302108 / RA (J2000) of object              
HIERARCH ESO INS SLIT22 OBJ DEC = 2.314742 / DEC (J2000) of object              
HIERARCH ESO INS SLIT22 X = 41.623 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT22 Y = 14.682 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT22 DIMX = 4.886 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT22 DIMY = 0.581 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT23 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT23 ID = 128901 / ID of SLIT i                              
HIERARCH ESO INS SLIT23 OBJ RA = 150.318517 / RA (J2000) of object              
HIERARCH ESO INS SLIT23 OBJ DEC = 2.31175 / DEC (J2000) of object               
HIERARCH ESO INS SLIT23 X = 47.343 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT23 Y = 48.991 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT23 DIMX = 6.079 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT23 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT24 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT24 ID = 128682 / ID of SLIT i                              
HIERARCH ESO INS SLIT24 OBJ RA = 150.324967 / RA (J2000) of object              
HIERARCH ESO INS SLIT24 OBJ DEC = 2.309175 / DEC (J2000) of object              
HIERARCH ESO INS SLIT24 X = 53.185 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT24 Y = 62.484 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT24 DIMX = 5.126 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT24 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT25 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT25 ID = 128528 / ID of SLIT i                              
HIERARCH ESO INS SLIT25 OBJ RA = 150.324742 / RA (J2000) of object              
HIERARCH ESO INS SLIT25 OBJ DEC = 2.305881 / DEC (J2000) of object              
HIERARCH ESO INS SLIT25 X = 59.324 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT25 Y = 62.007 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT25 DIMX = 6.676 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT25 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT26 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT26 ID = 128053 / ID of SLIT i                              
HIERARCH ESO INS SLIT26 OBJ RA = 150.301421 / RA (J2000) of object              
HIERARCH ESO INS SLIT26 OBJ DEC = 2.302322 / DEC (J2000) of object              
HIERARCH ESO INS SLIT26 X = 66.473 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT26 Y = 13.246 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT26 DIMX = 7.152 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT26 DIMY = 0.581 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT27 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT27 ID = 128133 / ID of SLIT i                              
HIERARCH ESO INS SLIT27 OBJ RA = 150.313362 / RA (J2000) of object              
HIERARCH ESO INS SLIT27 OBJ DEC = 2.299283 / DEC (J2000) of object              
HIERARCH ESO INS SLIT27 X = 73.388 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT27 Y = 38.213 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT27 DIMX = 6.199 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT27 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT28 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT28 ID = 127508 / ID of SLIT i                              
HIERARCH ESO INS SLIT28 OBJ RA = 150.323738 / RA (J2000) of object              
HIERARCH ESO INS SLIT28 OBJ DEC = 2.294797 / DEC (J2000) of object              
HIERARCH ESO INS SLIT28 X = 81.321 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT28 Y = 59.917 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT28 DIMX = 9.182 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT28 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT29 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT29 ID = 127314 / ID of SLIT i                              
HIERARCH ESO INS SLIT29 OBJ RA = 150.296713 / RA (J2000) of object              
HIERARCH ESO INS SLIT29 OBJ DEC = 2.291669 / DEC (J2000) of object              
HIERARCH ESO INS SLIT29 X = 89.244 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT29 Y = 3.402 / y co-ordinate of SLIT i                     
HIERARCH ESO INS SLIT29 DIMX = 6.2 / Size in x direction of SLIT i              
HIERARCH ESO INS SLIT29 DIMY = 0.581 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT30 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT30 ID = 127180 / ID of SLIT i                              
HIERARCH ESO INS SLIT30 OBJ RA = 150.317087 / RA (J2000) of object              
HIERARCH ESO INS SLIT30 OBJ DEC = 2.288025 / DEC (J2000) of object              
HIERARCH ESO INS SLIT30 X = 96.283 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT30 Y = 46.016 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT30 DIMX = 7.394 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT30 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT31 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT31 ID = 127438 / ID of SLIT i                              
HIERARCH ESO INS SLIT31 OBJ RA = 150.321308 / RA (J2000) of object              
HIERARCH ESO INS SLIT31 OBJ DEC = 2.28495 / DEC (J2000) of object               
HIERARCH ESO INS SLIT31 X = 103.324 / x co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT31 Y = 54.842 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT31 DIMX = 6.203 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT31 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT32 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT32 ID = 126585 / ID of SLIT i                              
HIERARCH ESO INS SLIT32 OBJ RA = 150.324013 / RA (J2000) of object              
HIERARCH ESO INS SLIT32 OBJ DEC = 2.282417 / DEC (J2000) of object              
HIERARCH ESO INS SLIT32 X = 109.233 / x co-ordinate of SLIT i                   
HIERARCH ESO INS SLIT32 Y = 60.508 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT32 DIMX = 5.13 / Size in x direction of SLIT i             
HIERARCH ESO INS SLIT32 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS SLIT33 TYPE = 'RECTANGLE' / Type of SLIT i                     
HIERARCH ESO INS SLIT33 ID = 126275 / ID of SLIT i                              
HIERARCH ESO INS SLIT33 OBJ RA = 150.301817 / RA (J2000) of object              
HIERARCH ESO INS SLIT33 OBJ DEC = 2.279189 / DEC (J2000) of object              
HIERARCH ESO INS SLIT33 X = 117.03 / x co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT33 Y = 14.079 / y co-ordinate of SLIT i                    
HIERARCH ESO INS SLIT33 DIMX = 10.02 / Size in x direction of SLIT i            
HIERARCH ESO INS SLIT33 DIMY = 0.582 / Size in y direction of SLIT i            
HIERARCH ESO INS PIXSCALE = 0.205 / Pixel scale                                 
HIERARCH ESO INS MODE = 'MOS     ' / Instrument Observation Mode                
HIERARCH ESO INS FILT1 FIX = F / Filter mode degraded                           
HIERARCH ESO INS FILT1 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS FILT1 NAME = 'GG475   ' / Filter name                          
HIERARCH ESO INS FILT1 ID = 'vm-GG475-1.10' / Filter ESO ID                     
HIERARCH ESO INS MASK1 FIX = F / Mask mode degraded                             
HIERARCH ESO INS MASK1 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS MASK1 NO =  3 / Mask set number                                
HIERARCH ESO INS FOCU1 FIX = F / Focusing mode degraded                         
HIERARCH ESO INS FOCU1 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS FOCU1 VALUE = 6976.122982 / Camera lens position(micron)       
HIERARCH ESO INS FOCU3 FIX = F / Focusing mode degraded                         
HIERARCH ESO INS FOCU3 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS FOCU3 TEMP = 14.629684 / Temperature(C) for the beam           
HIERARCH ESO INS FOCU3 VALUE = 5181.099499 / Camera lens position(micron)       
HIERARCH ESO INS FOCU2 FIX = F / Focusing mode degraded                         
HIERARCH ESO INS FOCU2 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS FOCU2 TEMP = 14.841312 / Temperature(C) for the beam           
HIERARCH ESO INS FOCU2 VALUE = 6081.203659 / Camera lens position(micron)       
HIERARCH ESO INS MSHU1 FIX = F / Mask Shutter mode degraded                     
HIERARCH ESO INS MSHU1 SIM = T / Device Hardware Simulation                     
HIERARCH ESO INS MSHU1 MODE = 'OFF     ' / Mask Shutter mode                    
HIERARCH ESO INS MSHU1 POSH = 7.348085 / High Mask Shutters positions (mm)      
HIERARCH ESO INS MSHU1 POSL = 7.348085 / Low Mask Shutters position (mm)        
HIERARCH ESO INS GRIS1 FIX = F / Grism mode degraded                            
HIERARCH ESO INS GRIS1 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS GRIS1 NAME = 'HR_red  ' / Grism name                           
HIERARCH ESO INS GRIS1 ID = 'vm-HRredHolog-1.1' / Grism ESO ID                  
HIERARCH ESO INS AFC11 FIX = F / AFC mode degraded                              
HIERARCH ESO INS AFC11 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS AFC12 FIX = F / AFC mode degraded                              
HIERARCH ESO INS AFC12 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS FOCU4 FIX = F / Focusing mode degraded                         
HIERARCH ESO INS FOCU4 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS FOCU4 TEMP = 14.588951 / Temperature(C) for the beam           
HIERARCH ESO INS FOCU4 VALUE = 7108.120407 / Camera lens position(micron)       
HIERARCH ESO INS LAMP1 FIX = F / Lamp[i] mode degraded                          
HIERARCH ESO INS LAMP1 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS LAMP1 TIME = 0 / Lamp lightening time                          
HIERARCH ESO INS LAMP1 STATE = 'OFF     ' / State for this lamp in each box     
HIERARCH ESO INS LAMP1 NAME = 'QTH50   ' / Calibration lamp name                
HIERARCH ESO INS LAMP1 ID = 'vm-QTH50-1' / Calib. lamp ESO identification       
HIERARCH ESO INS LAMP2 FIX = F / Lamp[i] mode degraded                          
HIERARCH ESO INS LAMP2 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS LAMP2 TIME = 0 / Lamp lightening time                          
HIERARCH ESO INS LAMP2 STATE = 'OFF     ' / State for this lamp in each box     
HIERARCH ESO INS LAMP2 NAME = 'QTH10   ' / Calibration lamp name                
HIERARCH ESO INS LAMP2 ID = 'vm-QTH10-2' / Calib. lamp ESO identification       
HIERARCH ESO INS IFUM1 FIX = F / IFU Mask mode degraded                         
HIERARCH ESO INS IFUM1 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS IFUM1 MODE = 'OFF     ' / IFU Mask mode                        
HIERARCH ESO INS LAMP3 FIX = F / Lamp[i] mode degraded                          
HIERARCH ESO INS LAMP3 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS LAMP3 TIME = 0 / Lamp lightening time                          
HIERARCH ESO INS LAMP3 STATE = 'OFF     ' / State for this lamp in each box     
HIERARCH ESO INS LAMP3 NAME = 'He      ' / Calibration lamp name                
HIERARCH ESO INS LAMP3 ID = 'vm-He-3 ' / Calib. lamp ESO identification         
HIERARCH ESO INS IFUM2 FIX = F / IFU Mask mode degraded                         
HIERARCH ESO INS IFUM2 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS IFUM2 MODE = 'OFF     ' / IFU Mask mode                        
HIERARCH ESO INS LAMP4 FIX = F / Lamp[i] mode degraded                          
HIERARCH ESO INS LAMP4 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS LAMP4 TIME = 0 / Lamp lightening time                          
HIERARCH ESO INS LAMP4 STATE = 'OFF     ' / State for this lamp in each box     
HIERARCH ESO INS LAMP4 NAME = 'Ne      ' / Calibration lamp name                
HIERARCH ESO INS LAMP4 ID = 'vm-Ne-4 ' / Calib. lamp ESO identification         
HIERARCH ESO INS LAMP5 FIX = F / Lamp[i] mode degraded                          
HIERARCH ESO INS LAMP5 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS LAMP5 TIME = 0 / Lamp lightening time                          
HIERARCH ESO INS LAMP5 STATE = 'OFF     ' / State for this lamp in each box     
HIERARCH ESO INS LAMP5 NAME = 'Ar      ' / Calibration lamp name                
HIERARCH ESO INS LAMP5 ID = 'vm-Ar-5 ' / Calib. lamp ESO identification         
HIERARCH ESO INS LAMP6 FIX = F / Lamp[i] mode degraded                          
HIERARCH ESO INS LAMP6 SIM = F / Device Hardware Simulation                     
HIERARCH ESO INS LAMP6 TIME = 0 / Lamp lightening time                          
HIERARCH ESO INS LAMP6 STATE = 'OFF     ' / State for this lamp in each box     
HIERARCH ESO INS LAMP6 NAME = 'Afc     ' / Calibration lamp name                
HIERARCH ESO INS LAMP6 ID = 'vm-Afc-6' / Calib. lamp ESO identification         
HIERARCH ESO INS IFUE FIX =  F / IFU Elongator mode degraded                    
HIERARCH ESO INS IFUE SIM =  F / Device Hardware Simulation                     
HIERARCH ESO INS IFUE MAG = '1:1     ' / IFU Magnification 0.33(1:2) or 0.6     
HIERARCH ESO INS IFUS FIX =  F / IFU Shutter mode degraded                      
HIERARCH ESO INS IFUS SIM =  F / Device Hardware Simulation                     
HIERARCH ESO INS IFUS MODE = 'OFF     ' / IFU Shutter CLOSED(ON) or OPEN(OFF    
HIERARCH ESO DET CHIP1 INDEX = 1 / Chip index                                   
HIERARCH ESO DET CHIP1 ID = 'BRIAN   ' / Detector chip identification           
HIERARCH ESO DET CHIP1 NAME = 'EEV CCD-44-82' / Detector chip name              
HIERARCH ESO DET CHIP1 DATE = '2000-08-30' / Date of installation [YYYY-MM-DD]  
HIERARCH ESO DET CHIP1 X =   1 / X location in array                            
HIERARCH ESO DET CHIP1 Y =   1 / Y location in array                            
HIERARCH ESO DET CHIP1 NX = 2048 / # of pixels along X                          
HIERARCH ESO DET CHIP1 NY = 4096 / # of pixels along Y                          
HIERARCH ESO DET CHIP1 PSZX = 15. / Size of pixel in X                          
HIERARCH ESO DET CHIP1 PSZY = 15. / Size of pixel in Y                          
HIERARCH ESO DET CHIP1 XGAP = 0. / Gap between chips along x                    
HIERARCH ESO DET CHIP1 YGAP = 0. / Gap between chips along y                    
HIERARCH ESO DET OUT1 INDEX = 1 / Output index                                  
HIERARCH ESO DET OUT1 ID = 'L       ' / Output ID as from manufacturer          
HIERARCH ESO DET OUT1 NAME = 'L       ' / Description of output                 
HIERARCH ESO DET OUT1 CHIP = 1 / index of chip it belongs to                    
HIERARCH ESO DET OUT1 X =    1 / X location of output                           
HIERARCH ESO DET OUT1 Y =    1 / Y location of output                           
HIERARCH ESO DET OUT1 NX = 2048 / valid pixels along X                          
HIERARCH ESO DET OUT1 NY = 4096 / valid pixels along Y                          
HIERARCH ESO DET OUT1 CONAD = 0.55 / Conversion from ADUs to electrons          
HIERARCH ESO DET OUT1 RON = 3.21 / [e-] Readout noise per output (ele           
HIERARCH ESO DET OUT1 GAIN = 1.82 / Conversion from electrons to ADU            
HIERARCH ESO DET ID = 'CCD FIERA - Rev: 238334' / Detector system Id            
HIERARCH ESO DET NAME = 'Vimosb - vimosb' / Name of detector system             
HIERARCH ESO DET DATE = '2010-07-01' / Installation date                        
HIERARCH ESO DET DID = 'ESO-VLT-DIC.CCDDCS,ESO-VLT-DIC.FCDDCS' / Diction        
HIERARCH ESO DET BITS =     16 / Bits per pixel readout                         
HIERARCH ESO DET RA =       0. / Apparent 00:00:00.0 RA at start                
HIERARCH ESO DET DEC =      0. / Apparent 00:00:00.0 DEC at start               
HIERARCH ESO DET CHIPS =     1 / # of chips in detector array                   
HIERARCH ESO DET OUTPUTS =   1 / # of outputs                                   
HIERARCH ESO DET OUTREF =    0 / reference output                               
HIERARCH ESO DET WINDOWS =   1 / # of windows readout                           
HIERARCH ESO DET SOFW MODE = 'Normal  ' / CCD sw operational mode               
HIERARCH ESO DET EXP NO = 42197 / Unique exposure ID number                     
HIERARCH ESO DET EXP TYPE = 'Normal  ' / Exposure type                          
HIERARCH ESO DET EXP RDTTIME = 30.939 / image readout time                      
HIERARCH ESO DET EXP XFERTIM = 44.981 / image transfer time                     
HIERARCH ESO DET WIN1 ST =   T / If T, window enabled                           
HIERARCH ESO DET WIN1 STRX = 1 / Lower left pixel in X                          
HIERARCH ESO DET WIN1 STRY = 1 / Lower left pixel in Y                          
HIERARCH ESO DET WIN1 NX = 2148 / # of pixels along X                           
HIERARCH ESO DET WIN1 NY = 4096 / # of pixels along Y                           
HIERARCH ESO DET WIN1 BINX = 1 / Binning factor along X                         
HIERARCH ESO DET WIN1 BINY = 1 / Binning factor along Y                         
HIERARCH ESO DET WIN1 NDIT = 1 / # of subintegrations                           
HIERARCH ESO DET WIN1 UIT1 = 1200. / user defined subintegration time           
HIERARCH ESO DET WIN1 DIT1 = 1199.997383 / actual subintegration time           
HIERARCH ESO DET WIN1 DKTM = 1200.0992 / Dark current time                      
HIERARCH ESO DET READ MODE = 'normal  ' / Readout method                        
HIERARCH ESO DET READ SPEED = 'medium  ' / Readout speed                        
HIERARCH ESO DET READ CLOCK = '225kHz,1port,HG' / Readout clock pattern used    
HIERARCH ESO DET READ NFRAM = 1 / Number of readouts buffered in sin            
HIERARCH ESO DET FRAM ID =   1 / Image sequencial number                        
HIERARCH ESO DET FRAM TYPE = 'Normal  ' / Type of frame                         
HIERARCH ESO DET SHUT TYPE = 'seso    ' / type of shutter                       
HIERARCH ESO DET SHUT ID = 'VIMOS shutter' / Shutter unique identifier          
HIERARCH ESO DET SHUT TMOPEN = 0.074 / Time taken to open shutter               
HIERARCH ESO DET SHUT TMCLOS = 0.069 / Time taken to close shutter              
HIERARCH ESO DET TELE INT = 900. / Interval between two successive te           
HIERARCH ESO DET TELE NO =   7 / # of sources active                            
HIERARCH ESO DET TLM1 NAME = 'CCD Cold Plate' / Description of telemetry param  
HIERARCH ESO DET TLM1 ID = 'CCD_B1_T1' / ID of telemetry sensor                 
HIERARCH ESO DET TLM1 START = 140.1 / Telemetry value at read start             
HIERARCH ESO DET TLM1 END = 139.9 / Telemetry value at read completion          
HIERARCH ESO DET TLM2 NAME = 'CCD Cold Plate2' / Description of telemetry para  
HIERARCH ESO DET TLM2 ID = 'CCD_B1_T2' / ID of telemetry sensor                 
HIERARCH ESO DET TLM2 START = 151.6 / Telemetry value at read start             
HIERARCH ESO DET TLM2 END = 151.4 / Telemetry value at read completion          
HIERARCH ESO DET TLM3 NAME = 'CCD Cold Plate' / Description of telemetry param  
HIERARCH ESO DET TLM3 ID = 'CCD_B4_T1' / ID of telemetry sensor                 
HIERARCH ESO DET TLM3 START = 140. / Telemetry value at read start              
HIERARCH ESO DET TLM3 END = 140. / Telemetry value at read completion           
HIERARCH ESO DET TLM4 NAME = 'CCD Cold Plate2' / Description of telemetry para  
HIERARCH ESO DET TLM4 ID = 'CCD_B4_T2' / ID of telemetry sensor                 
HIERARCH ESO DET TLM4 START = 140.6 / Telemetry value at read start             
HIERARCH ESO DET TLM4 END = 140.6 / Telemetry value at read completion          
HIERARCH ESO DET TLM5 NAME = 'FBOX    ' / Description of telemetry param.       
HIERARCH ESO DET TLM5 ID = 'FBOX    ' / ID of telemetry sensor                  
HIERARCH ESO DET TLM5 START = 288.8 / Telemetry value at read start             
HIERARCH ESO DET TLM5 END = 288.7 / Telemetry value at read completion          
HIERARCH ESO DET TLM6 NAME = 'Vaccum_B1' / Description of telemetry param.      
HIERARCH ESO DET TLM6 ID = 'Vaccum_B1' / ID of telemetry sensor                 
HIERARCH ESO DET TLM6 START = 0. / Telemetry value at read start                
HIERARCH ESO DET TLM6 END = 0. / Telemetry value at read completion             
HIERARCH ESO DET TLM7 NAME = 'Vaccum_B4' / Description of telemetry param.      
HIERARCH ESO DET TLM7 ID = 'Vaccum_B4' / ID of telemetry sensor                 
HIERARCH ESO DET TLM7 START = 0. / Telemetry value at read start                
HIERARCH ESO DET TLM7 END = 0. / Telemetry value at read completion             
HIERARCH ESO PRO EXPTTOT = 4799.988                                             
HIERARCH ESO PRO ANCESTOR = 'VIMOS.2015-01-18T04:03:48.110.fits'                
HIERARCH ESO PRO DID = 'PRO-1.15' / Data dictionary for PRO                     
HIERARCH ESO PRO CATG = 'MOS_SCIENCE_SKY_EXTRACTED' / Category of pipeline produ
HIERARCH ESO PRO TYPE = 'REDUCED ' / Product type                               
HIERARCH ESO PRO TECH = 'MOS     ' / Observation technique                      
HIERARCH ESO PRO SCIENCE =   T / Scientific product if T                        
HIERARCH ESO PRO REC1 ID = 'vmmosscience' / Pipeline recipe (unique) identifier 
HIERARCH ESO PRO REC1 DRS ID = 'cpl-6.5.1' / Data Reduction System identifier   
HIERARCH ESO PRO REC1 PIPE ID = 'vimos-2.9.15' / Pipeline (unique) identifier   
HIERARCH ESO PRO REC1 RAW1 NAME = 'VIMOS.2015-01-18T04:03:48.110.fits' / File na
HIERARCH ESO PRO REC1 RAW1 CATG = 'MOS_SCIENCE' / Category of raw frame         
HIERARCH ESO PRO REC1 RAW2 NAME = 'VIMOS.2015-01-18T04:24:43.107.fits' / File na
HIERARCH ESO PRO REC1 RAW2 CATG = 'MOS_SCIENCE' / Category of raw frame         
HIERARCH ESO PRO REC1 RAW3 NAME = 'VIMOS.2015-01-18T04:45:38.102.fits' / File na
HIERARCH ESO PRO REC1 RAW3 CATG = 'MOS_SCIENCE' / Category of raw frame         
HIERARCH ESO PRO REC1 RAW4 NAME = 'VIMOS.2015-01-18T05:06:33.097.fits' / File na
HIERARCH ESO PRO REC1 RAW4 CATG = 'MOS_SCIENCE' / Category of raw frame         
HIERARCH ESO PRO DATANCOM =  4 / Number of combined frames                      
HIERARCH ESO PRO REC1 CAL1 NAME = 'M.VIMOS.2010-08-19T16:16:17.563.fits' / File 
HIERARCH ESO PRO REC1 CAL1 CATG = 'CONFIG_TABLE' / Category of calibration frame
HIERARCH ESO PRO REC1 CAL1 DATAMD5 = '827358dea0214d1587a3f99820ec5810' / MD5 si
HIERARCH ESO PRO REC1 CAL2 NAME = 'master_bias.fits' / File name of calibration 
HIERARCH ESO PRO REC1 CAL2 CATG = 'MASTER_BIAS' / Category of calibration frame 
HIERARCH ESO PRO REC1 CAL2 DATAMD5 = 'e61a66abf72ac06c6b41c8c3770404e3' / MD5 si
HIERARCH ESO PRO REC1 CAL3 NAME = 'mos_curv_coeff.fits' / File name of calibrati
HIERARCH ESO PRO REC1 CAL3 CATG = 'MOS_CURV_COEFF' / Category of calibration fra
HIERARCH ESO PRO REC1 CAL3 DATAMD5 = '8c4b8da9b95016a47162e7a44741b6e9' / MD5 si
HIERARCH ESO PRO REC1 CAL4 NAME = 'mos_disp_coeff.fits' / File name of calibrati
HIERARCH ESO PRO REC1 CAL4 CATG = 'MOS_DISP_COEFF' / Category of calibration fra
HIERARCH ESO PRO REC1 CAL4 DATAMD5 = '1c1a4eb0e8edb0d621c99d716b0c63ae' / MD5 si
HIERARCH ESO PRO REC1 CAL5 NAME = 'mos_master_screen_flat.fits' / File name of c
HIERARCH ESO PRO REC1 CAL5 CATG = 'MOS_MASTER_SCREEN_FLAT' / Category of calibra
HIERARCH ESO PRO REC1 CAL5 DATAMD5 = '0f2d03e3af5e30f7485a0a300f70b02f' / MD5 si
HIERARCH ESO PRO REC1 CAL6 NAME = 'mos_slit_location.fits' / File name of calibr
HIERARCH ESO PRO REC1 CAL6 CATG = 'MOS_SLIT_LOCATION' / Category of calibration 
HIERARCH ESO PRO REC1 CAL6 DATAMD5 = 'f53c43722c662a5763c5b94db9b2468e' / MD5 si
HIERARCH ESO PRO REC1 PARAM1 NAME = 'dispersion' / Resampling step (Angstrom/pix
HIERARCH ESO PRO REC1 PARAM1 VALUE = '0.6     ' / Default: 0                    
HIERARCH ESO PRO REC1 PARAM2 NAME = 'skyalign' / Polynomial order for sky lines 
HIERARCH ESO PRO REC1 PARAM2 VALUE = '0       ' / Default: 0                    
HIERARCH ESO PRO REC1 PARAM3 NAME = 'wcolumn ' / Name of sky line catalog table 
HIERARCH ESO PRO REC1 PARAM3 VALUE = 'WLEN    ' / Default: 'WLEN'               
HIERARCH ESO PRO REC1 PARAM4 NAME = 'startwavelength' / Start wavelength in spec
HIERARCH ESO PRO REC1 PARAM4 VALUE = '5800    ' / Default: 0                    
HIERARCH ESO PRO REC1 PARAM5 NAME = 'endwavelength' / End wavelength in spectral
HIERARCH ESO PRO REC1 PARAM5 VALUE = '9500    ' / Default: 0                    
HIERARCH ESO PRO REC1 PARAM6 NAME = 'reference' / Reference wavelength for calib
HIERARCH ESO PRO REC1 PARAM6 VALUE = '7948.18 ' / Default: 0                    
HIERARCH ESO PRO REC1 PARAM7 NAME = 'flux    ' / Apply flux conservation        
HIERARCH ESO PRO REC1 PARAM7 VALUE = 'true    ' / Default: true                 
HIERARCH ESO PRO REC1 PARAM8 NAME = 'flatfield' / Apply flat field              
HIERARCH ESO PRO REC1 PARAM8 VALUE = 'true    ' / Default: true                 
HIERARCH ESO PRO REC1 PARAM9 NAME = 'skyglobal' / Subtract global sky spectrum f
HIERARCH ESO PRO REC1 PARAM9 VALUE = 'false   ' / Default: false                
HIERARCH ESO PRO REC1 PARAM10 NAME = 'skymedian' / Sky subtraction from extracte
HIERARCH ESO PRO REC1 PARAM10 VALUE = 'false   ' / Default: false               
HIERARCH ESO PRO REC1 PARAM11 NAME = 'skylocal' / Sky subtraction from CCD slit 
HIERARCH ESO PRO REC1 PARAM11 VALUE = 'true    ' / Default: true                
HIERARCH ESO PRO REC1 PARAM12 NAME = 'cosmics ' / Eliminate cosmic rays hits (on
HIERARCH ESO PRO REC1 PARAM12 VALUE = 'true    ' / Default: true                
HIERARCH ESO PRO REC1 PARAM13 NAME = 'slit_margin' / Number of pixels to exclude
HIERARCH ESO PRO REC1 PARAM13 VALUE = '3       ' / Default: 3                   
HIERARCH ESO PRO REC1 PARAM14 NAME = 'ext_radius' / Maximum extraction radius fo
HIERARCH ESO PRO REC1 PARAM14 VALUE = '6       ' / Default: 6                   
HIERARCH ESO PRO REC1 PARAM15 NAME = 'cont_radius' / Minimum distance at which t
HIERARCH ESO PRO REC1 PARAM15 VALUE = '0       ' / Default: 0                   
HIERARCH ESO PRO REC1 PARAM16 NAME = 'ext_mode' / Object extraction method: 0 = 
HIERARCH ESO PRO REC1 PARAM16 VALUE = '1       ' / Default: 1                   
HIERARCH ESO PRO REC1 PARAM17 NAME = 'detection' / Object detection threshold (A
HIERARCH ESO PRO REC1 PARAM17 VALUE = '2       ' / Default: 2                   
HIERARCH ESO PRO REC1 PARAM18 NAME = 'time_normalise' / Normalise output spectra
HIERARCH ESO PRO REC1 PARAM18 VALUE = 'true    ' / Default: true                
HIERARCH ESO PRO REC1 PARAM19 NAME = 'anyframe' / Look for a standard star in an
HIERARCH ESO PRO REC1 PARAM19 VALUE = 'false   ' / Default: false               
HIERARCH ESO PRO REC1 PARAM20 NAME = 'response' / Order of polynomial modeling t
HIERARCH ESO PRO REC1 PARAM20 VALUE = '5       ' / Default: 5                   
HIERARCH ESO PRO REC1 PARAM21 NAME = 'alignment' / Type of alignment of dithered
HIERARCH ESO PRO REC1 PARAM21 VALUE = 'integer ' / Default: 'integer'           
HIERARCH ESO PRO REC1 PARAM22 NAME = 'stack_method' / Frames combination method 
HIERARCH ESO PRO REC1 PARAM22 VALUE = 'average ' / Default: 'average'           
HIERARCH ESO PRO REC1 PARAM23 NAME = 'minrejection' / Number of lowest values to
HIERARCH ESO PRO REC1 PARAM23 VALUE = '1       ' / Default: 1                   
HIERARCH ESO PRO REC1 PARAM24 NAME = 'maxrejection' / Number of highest values t
HIERARCH ESO PRO REC1 PARAM24 VALUE = '1       ' / Default: 1                   
HIERARCH ESO PRO REC1 PARAM25 NAME = 'klow    ' / Low threshold in ksigma method
HIERARCH ESO PRO REC1 PARAM25 VALUE = '3       ' / Default: 3                   
HIERARCH ESO PRO REC1 PARAM26 NAME = 'khigh   ' / High threshold in ksigma metho
HIERARCH ESO PRO REC1 PARAM26 VALUE = '3       ' / Default: 3                   
HIERARCH ESO PRO REC1 PARAM27 NAME = 'kiter   ' / Max number of iterations in ks
HIERARCH ESO PRO REC1 PARAM27 VALUE = '999     ' / Default: 999                 
HIERARCH ESO PRO REC1 PARAM28 NAME = 'dither  ' / Align dithered frames before s
HIERARCH ESO PRO REC1 PARAM28 VALUE = 'true    ' / Default: true                
HIERARCH ESO PRO REC1 PARAM29 NAME = 'compute ' / Compute offsets of dithered im
HIERARCH ESO PRO REC1 PARAM29 VALUE = 'false   ' / Default: false               
HIERARCH ESO PRO REC1 PARAM30 NAME = 'fringing' / Apply fringing correction (onl
HIERARCH ESO PRO REC1 PARAM30 VALUE = 'true    ' / Default: true                
HIERARCH ESO PRO REC1 PARAM31 NAME = 'offset  ' / Minimum required offset betwee
HIERARCH ESO PRO REC1 PARAM31 VALUE = '3       ' / Default: 3                   
HIERARCH ESO PRO REC1 PARAM32 NAME = 'qc      ' / Compute QC1 parameters        
HIERARCH ESO PRO REC1 PARAM32 VALUE = 'true    ' / Default: true                
HIERARCH ESO ADA ABSROT START = 21.56262 / Abs rot angle at exp start           
HIERARCH ESO ADA POSANG =  90. / Position angle at start                        
HIERARCH ESO ADA GUID STATUS = 'ON      ' / Status of autoguider                
HIERARCH ESO ADA GUID RA = 150.241537 / 10:00:57.9 Guide star RA J2000          
HIERARCH ESO ADA GUID DEC = 2.59418 / 02:35:39.0 Guide star DEC J2000           
HIERARCH ESO ADA ABSROT PPOS = 'NEG     ' / sign of probe position              
HIERARCH ESO ADA ABSROT END = 20.56639 / Abs rot angle at exp end               
HIERARCH ESO OCS CON QUAD =  1 / Quadrant Number                                
HIERARCH ESO OCS CMF FILE1 NAME = 'MOS_sky2ccd_1.cmf' / Calibration matrix file 
HIERARCH ESO OCS CMF FILE1 CRTEDATE = '2012-05-18T14:37:53' / Calibration matrix
HIERARCH ESO OCS CMF FILE2 NAME = 'MOS_wavecal_HR_red_1.cmf' / Calibration matri
HIERARCH ESO OCS CMF FILE2 CRTEDATE = '2005-10-20T17:18:57' / Calibration matrix
HIERARCH ESO OCS COMP ID = '226602 2012-03-18 17:49:45Z' / OS Software Versi    
HIERARCH ESO OCS DID = 'ESO-VLT-DIC.VIMOS_OS-0.2' / Data dictionary for         
HISTORY FTU-2_5_2/2015-01-18/REMOVE: XTENSION, "IMAGE"                          
HISTORY FTU-2_5_2/2015-01-18/REMOVE: EXTNAME, "BRIAN"                           
HISTORY FTU-2_5_2/2015-01-18/REMOVE: ORIGIN, "VIRMOS OS Software"               
HISTORY FTU-2_5_2/2015-01-18/REMOVE: CDELT1, 1.00000000                         
HISTORY FTU-2_5_2/2015-01-18/REMOVE: CDELT2, 1.00000000                         
HISTORY FTU-2_5_2/2015-01-18/CHANGE: DET.CHIPS, 1                               
HISTORY FTU-2_5_2/2015-01-18/CHANGE: DET.OUTPUTS, 1                             
HISTORY FTU-2_5_2/2015-01-18/CHANGE: DET.WIN1.NX, 2148                          
COMMENT FTU-2_5_2/2015-01-18T04:24:44/vimos.ftt                                 
CHECKSUM= 'NMTOQMTNNMTNNMTN'   / HDU checksum updated 2015-01-26T12:24:13       
DATASUM = '338047765'          / data unit checksum updated 2015-01-26T12:24:13 

If i want to do it in Pandas.

flux = hdu1d[0].data.byteswap().newbyteorder()
wht = hdu1dwht[0].data.byteswap().newbyteorder()
unc = 1./ np.sqrt(wht)
wave = np.arange(flux.shape[0])*hdu1d[0].header['CD1_1'] + hdu1d[0].header['CRVAL1']

##for certain functions, I need to cut the spectrum where the weight is 0.
d = {'wavelength':wave, 'flux':flux, 'weight':wht, 'uncertainty':unc}
dataspec = pd.DataFrame(data=d)

##define subset where wht>0
dataspec_sub = dataspec[dataspec['weight'] > 0.].reset_index(drop=True)
dataspec_sub.head()
wavelength flux weight uncertainty
0 6229.3 14.757169 0.207031 2.197769
1 6229.9 26.817886 0.148122 2.598306
2 6230.5 23.955254 0.142191 2.651945
3 6231.1 22.239902 0.141403 2.659323
4 6231.7 22.505369 0.136912 2.702580

If I want to do it in astropy Tables.

flux = hdu1d[0].data
wht = hdu1dwht[0].data
unc = 1./ np.sqrt(wht)
wave = np.arange(flux.shape[0])*hdu1d[0].header['CD1_1'] + hdu1d[0].header['CRVAL1']

spec_unit = u.Unit('10^-19 erg s^-1 cm^-2 angstrom^-1')
dataspec = QTable([wave*u.angstrom, flux*spec_unit, wht, unc*spec_unit], 
                   names=('wavelength','flux','weight','uncertainty'))
dataspec_sub = dataspec[dataspec['weight']>0.]
dataspec_sub
QTable length=4059
wavelengthfluxweightuncertainty
Angstrom1e-19 erg / (Angstrom s cm2)1e-19 erg / (Angstrom s cm2)
float64float32float32float32
6229.314.7571687698364260.207031262.1977689266204834
6229.90000000000126.8178863525390620.148121982.5983057022094727
6230.523.955253601074220.142190642.651945114135742
6231.122.239902496337890.141402772.659322738647461
6231.722.5053691864013670.136912392.702580213546753
6232.321.2664356231689450.128615022.788393974304199
6232.90000000000121.2351989746093750.113255852.9714584350585938
6233.524.182773590087890.0968212043.2137696743011475
6234.132.750209808349610.0855739863.4184489250183105
............
8677.946.5547790527343750.087044923.3894426822662354
8678.545.34338378906250.085043873.429086923599243
8679.146.191535949707030.068361413.824674606323242
8679.749.6320571899414060.0585494824.132743835449219
8680.340.5991592407226560.049111784.51239538192749
8680.946.059841156005860.04548884.6886491775512695
8681.553.855808258056640.0498597554.478421211242676
8682.151.2714004516601560.0584816824.135138511657715
8682.751.619499206542970.0690493063.805575370788574
8683.351.206130981445310.089989483.3335282802581787
plt.figure(figsize=[10,6])
plt.fill_between(dataspec_sub['wavelength'],
                 dataspec_sub['flux']+dataspec_sub['uncertainty'],
                 dataspec_sub['flux']-dataspec_sub['uncertainty'],
                 color='gray',label="uncertainty")
plt.plot(dataspec_sub['wavelength'],dataspec_sub['flux'], label="data")
plt.xlabel("wavelength ({:latex})".format(dataspec_sub['wavelength'].unit))
plt.ylabel("flux ({:latex})".format(dataspec_sub['flux'].unit))
plt.title("Observed spectrum")
plt.legend()
plt.show()
../../_images/eb721e276b4bba24287aca60479706766416f09a64152061f6c7bfee1d47ac53.png

Go with specutils#

#write Spectrum1D object
spec1d = Spectrum1D(spectral_axis=dataspec_sub['wavelength'], 
                    flux=dataspec_sub['flux'], 
                    uncertainty=StdDevUncertainty(dataspec_sub['uncertainty']))

Developer note#

For supported datasets (like final JWST data products), this will be as simple as:

spec1d = Spectrum1d.read('datafile.fits')

Developer note#

Implemented but not yet released: snr_threshold, which will allow cutting the spectrum using that function.

Smooth to better inspect the features#

Developer note#

The uncertainty is not carried over, but would be useful. Same comment on the interactive tool as before: it would be very useful to be able to zoom, pan, hoover, etc. on the spectrum.

spec1d_gsmooth = gaussian_smooth(spec1d, stddev=5)
plt.figure(figsize=[10,6])
plt.plot(spec1d_gsmooth.spectral_axis,spec1d_gsmooth.flux)
plt.xlabel("wavelength ({:latex})".format(spec1d_gsmooth.spectral_axis.unit))
plt.ylabel("flux ({:latex})".format(spec1d_gsmooth.flux.unit))
plt.title("Smoothed observed spectrum")
plt.show()
../../_images/adf7db55e728ed7bad5402e5d7f633d025bb8a513c628d1575444a0337ec966a.png

Back to the non-smoothed spectrum to find lines#

Documentation says I need a continuum subtracted spectrum.

So fit continuum first#

cont_spec1d = fit_generic_continuum(spec1d)
cont_fit = cont_spec1d(spec1d.spectral_axis)
WARNING: Model is linear in parameters; consider using linear fitting methods. [astropy.modeling.fitting]
plt.figure(figsize=[10,6])
plt.plot(spec1d.spectral_axis, spec1d.flux, label="data")
plt.plot(spec1d.spectral_axis, cont_fit, label="modeled continuum")
plt.xlabel("wavelength ({:latex})".format(spec1d.spectral_axis.unit))
plt.ylabel("flux ({:latex})".format(spec1d.flux.unit))
plt.legend()
plt.title("Observed spectrum and fitted continuum")
plt.show()

plt.figure(figsize=[10,6])
plt.plot(spec1d.spectral_axis, spec1d.uncertainty.array, label="data")
plt.xlabel("wavelength ({:latex})".format(spec1d.spectral_axis.unit))
plt.ylabel("uncertainty ({:latex})".format(spec1d.uncertainty.unit))
plt.legend()
plt.title("Uncertianty of observed spectrum")
plt.show()
../../_images/bbbd09f15d1fac957675e16edf6e3ea4181c028046506f0b40c6508d33496f7d.png ../../_images/ff9cdf51228380c2a597da92926143a5ca1db068f6373b2a32559b02f2735d8a.png

Creating the continuum-subtracted spectrum#

Specutils will figure out what to do with the uncertainty!

spec1d_sub = spec1d - cont_fit
spec1d_sub
<Spectrum1D(flux=<Quantity [-9.11481202,  2.92418963,  0.03985127, ...,  2.46770635,
            2.81985914,  2.41054724] 1e-19 erg / (Angstrom s cm2)>, spectral_axis=<SpectralAxis 
   (observer to target:
      radial_velocity=0.0 km / s
      redshift=0.0)
  [6229.3, 6229.9, 6230.5, ..., 8682.1, 8682.7, 8683.3] Angstrom>, uncertainty=StdDevUncertainty([2.197769 , 2.5983057, 2.651945 , ..., 4.1351385,
                   3.8055754, 3.3335283]))>
plt.figure(figsize=[10,6])
plt.plot(spec1d_sub.spectral_axis, spec1d_sub.flux, label="data")
plt.xlabel("wavelength ({:latex})".format(spec1d_sub.spectral_axis.unit))
plt.ylabel("flux ({:latex})".format(spec1d_sub.flux.unit))
plt.legend()
plt.title("Continuum-subracted spectrum")
plt.show()

plt.figure(figsize=[10,6])
plt.plot(spec1d_sub.spectral_axis,spec1d_sub.uncertainty.array, label="data")
plt.xlabel("wavelength ({:latex})".format(spec1d_sub.spectral_axis.unit))
plt.ylabel("uncertainty ({:latex})".format(spec1d_sub.uncertainty.unit))
plt.legend()
plt.title("Uncertainty of continuum-subracted spectrum")
plt.show()
../../_images/012b67ca8619e9c1f2fe352cf5fda204897322d47a8ae6528cff3d9bf65973fc.png ../../_images/0dc93c82fbb22650a806865df95f11dc1fdcc90a3cf0ee5de4af640ed9b228a5.png

Now look for the lines#

lines = find_lines_derivative(spec1d_sub, flux_threshold=50)
lines
QTable length=4
line_centerline_typeline_center_index
Angstrom
float64str8int64
6549.1emission528
6553.900000000001emission536
6558.7emission544
8550.7emission3837
plt.figure(figsize=[10,6])
plt.plot(spec1d_sub.spectral_axis,spec1d_sub.flux,label="data")
plt.axvline(lines['line_center'][0].value, color="red", alpha=0.5, label='emission/absorption lines')
for line in lines:
    plt.axvline(line['line_center'].value, color='red',alpha=0.5)
plt.xlabel("wavelength ({:latex})".format(spec1d_sub.spectral_axis.unit))
plt.ylabel("flux ({:latex})".format(spec1d.flux.unit))
plt.legend()
plt.title("Continuum-subtracted spectrum and marked lines using find_lines_derivative")
plt.show()
../../_images/cfff7e1243022e591a2a70682c404239d8666cae10e0678fcc45374bb9d9092e.png

This works for cases where you understand the thresholds well, but doesn’t automate as well with noisy spectra.

Works better with find_lines_threshold#

lines = find_lines_threshold(spec1d_sub, noise_factor=6)
lines
QTable length=18
line_centerline_typeline_center_index
Angstrom
float64str10int64
6554.5emission537
7478.5emission2062
7488.1emission2072
8546.5emission3830
8551.3emission3838
6635.5absorption667
6671.5absorption727
6681.1absorption743
6739.9absorption841
6746.5absorption852
6913.3absorption1124
6916.9absorption1130
6975.7absorption1228
6979.9absorption1235
7205.5absorption1610
7217.5absorption1630
8079.700000000001absorption3055
8541.1absorption3821

Plot lines on the spectrum.

plt.figure(figsize=[10,6])
plt.plot(spec1d_sub.spectral_axis,spec1d_sub.flux,label="data")
plt.axvline(lines['line_center'][0].value, color="red", alpha=0.5, label='emission/absorption lines')
for line in lines:
    plt.axvline(line['line_center'].value, color='red', alpha=0.5)
plt.xlabel("wavelength ({:latex})".format(spec1d_sub.spectral_axis.unit))
plt.ylabel("flux ({:latex})".format(spec1d_sub.flux.unit))
plt.legend()
plt.title("Continuum-subtracted spectrum and marked lines using find_lines_threshold")
plt.show()
../../_images/38ed1be0dc9c2e0bd0efd344cab7426477c12538fa3dcad420c42baa76b7a27d.png

Developer note#

It would be useful to have a tool to cycle through the lines, show a zoom of the spectrum, and inspect how good the line identification is. For now I do it by hand on a single line.

plt.figure(figsize=[10,6])
plt.plot(spec1d_sub.spectral_axis,spec1d_sub.flux, label="data")
plt.scatter(spec1d_sub.spectral_axis,spec1d_sub.flux, label=None)
plt.axvline(lines['line_center'][0].value, color="red", alpha=0.5, label='[OII]')
for line in lines:
    plt.axvline(line['line_center'].value, alpha=0.5, color='red')
plt.xlim(6500,6600)
plt.xlabel("wavelength ({:latex})".format(spec1d_sub.spectral_axis.unit))
plt.ylabel("flux ({:latex})".format(spec1d_sub.flux.unit))
plt.legend()
plt.title("Continuum-subtracted spectrum zoomed on [OII]")
plt.show()
../../_images/fa954b0ec765493dd0981860b80b8d1e1c1b4e180647f3ecfe4c0e8b38be6e92.png

Measure line centroids and fluxes#

These too need spectra continuum subtracted.

#example with just one line
centroid(spec1d_sub, SpectralRegion(6540*u.AA, 6580*u.AA))
\[6554.415 \; \mathrm{\mathring{A}}\]
sline = centroid(spec1d_sub, SpectralRegion(6540*u.AA, 6580*u.AA))

plt.figure(figsize=[10,6])
plt.plot(spec1d_sub.spectral_axis,spec1d_sub.flux, label="data")
plt.scatter(spec1d_sub.spectral_axis,spec1d_sub.flux, label=None)
plt.axvline(sline.value, color='red', label="[OII]")
plt.axhline(0,color='black', label='flux = 0')
plt.xlim(6500,6600)
plt.xlabel("wavelength ({:latex})".format(spec1d_sub.spectral_axis.unit))
plt.ylabel("flux ({:latex})".format(spec1d_sub.flux.unit))
plt.legend()
plt.title("Continuum-subtracted spectrum zoomed on [OII]")
plt.show()
../../_images/11b044977023a626c0591230bb6893284df03f4d1d4894e87ccab5e3747527b2.png
line_flux(spec1d_sub, SpectralRegion(6540*u.AA, 6570*u.AA))  
\[965.51743 \; \mathrm{1 \times 10^{-19}\,\frac{erg}{s\,cm^{2}}}\]

Fit the line with a Gaussian#

Developer note#

Fitting lines with Gaussians is so common that it might make sense to have a one-line to iterate this over the list of lines, with sensible initialization of parameters and fit ranges for each line so that the fits generally work.

Also, the interactive tool here could allow one to select the continuum region and position of the line by clicking on the spectrum.

g_init = models.Gaussian1D(amplitude= 3 * 1e-19 * u.erg / u.s / u.cm**2 / u.AA, mean=6554*u.AA, stddev=2.*u.AA)
g_fit = fit_lines(spec1d_sub, g_init)
spec1d_fit = g_fit(spec1d_sub.wavelength)
g_fit
<Gaussian1D(amplitude=116.3316537 1e-19 erg / (Angstrom s cm2), mean=6554.93917566 Angstrom, stddev=3.71890257 Angstrom)>
#calculate the velocity dispertion from the stddev
vel = ((3.71890256/6554.415) * const.c.to('km/s').value)*u.km/u.s
print(vel)
170.0989240877931 km / s
plt.figure(figsize=[10,6])
plt.plot(spec1d_sub.wavelength,spec1d_sub.flux,label='data')
plt.plot(spec1d_sub.wavelength,spec1d_fit,color='darkorange',label='Gaussian fit')
plt.xlim(6500,6600)
plt.xlabel("wavelength ({:latex})".format(spec1d_sub.spectral_axis.unit))
plt.ylabel("flux ({:latex})".format(spec1d_sub.flux.unit))
plt.legend()
plt.title('Gaussian fit to the [OII] line')
plt.show()
../../_images/f5f16759986c2184ded03613872cc466933aae13a7430c94d08bdc0ca0942419.png

Measure the equivalent width of the lines#

This needs the spectrum continuum normalized.

spec1d_norm = spec1d / cont_fit
plt.figure(figsize=[10,6])
plt.plot(spec1d_norm.spectral_axis, spec1d_norm.flux, label='data')
plt.axhline(1,color='black',label='flux = 1')
plt.xlabel("wavelength ({:latex})".format(spec1d_norm.spectral_axis.unit))
plt.ylabel("flux (normalized)")
plt.xlim(6500,6600)
plt.legend()
plt.title("Continuum-normalized spectrum, zoomed on [OII]")
plt.show()

plt.figure(figsize=[10,6])
plt.plot(spec1d_norm.spectral_axis, spec1d_norm.uncertainty.array, label='data')
plt.xlabel("wavelength ({:latex})".format(spec1d_norm.spectral_axis.unit))
plt.ylabel("uncertainty (normalized)")
plt.xlim(6500,6600)
plt.legend()
plt.title("Uncertainty of continuum-normalized spectrum, zoomed on [OII]")
plt.show()
../../_images/23fb9e53c394f7453edfeadf62f666acce8275da1307b72c70adcce1aa379934.png ../../_images/77a1347aca5d21a74ab418681ab37585ec3898a24c45e42a72151e02ba5050ba.png
equivalent_width(spec1d_norm, regions=SpectralRegion(6540*u.AA, 6570*u.AA))
\[-28.243948 \; \mathrm{\mathring{A}}\]

Find the best-fitting template#

It needs a list of templates and the redshift of the observed galaxy. For the templates, I am using a set of model SEDs generated with Bruzual & Charlot stellar population models, emission lines, and dust attenuation as described in Pacifici et al. (2012).

Developer note#

Maybe there is a way to speed this up (maybe using astropy model_sets)? This fit is run with 100 models, but ideally, if we want to extract physical parameters from this, we would need at least 10,000 models.

A dictionary structure with meaningful keys (which can be, e.g., tuples of the relevant physical parameters) could be better than a list? It could make later analysis much clearer than having to map from the list indices back to the relevant parameters.

templatedir = './mos_spectroscopy/templates/'
zz = 0.758

templatelist = []
for i in range (1, 101):
    template_file = "{0}{1:05d}.dat".format(templatedir,i)
    template = ascii.read(template_file)
    temp1d = Spectrum1D(spectral_axis=template['col1']*u.AA,flux=template['col2']*u.erg/u.s/u.AA)
    templatelist.append(temp1d)
    
    
tm_results = template_comparison.template_match(observed_spectrum=spec1d, spectral_templates=templatelist, resample_method="flux_conserving", redshift=zz)
tm_results[0]
<Spectrum1D(flux=<Quantity [0.        , 0.        , 0.        , ..., 0.09480264, 0.09410727,
           0.09341189] 1e-19 erg / (Angstrom s cm2)>, spectral_axis=<SpectralAxis [   177.0306,    178.6128,    180.195 , ..., 170526.    , 171229.2   ,
   171932.4   ] Angstrom>)>
plt.figure(figsize=[10, 6])
plt.plot(spec1d.wavelength, spec1d.flux, label="data")
plt.plot(tm_results[0].wavelength, tm_results[0].flux,color='r',alpha=0.5,label='model')
plt.xlim(6000, 9000)
plt.xlabel("wavelength ({:latex})".format(spec1d_norm.spectral_axis.unit))
plt.ylabel("flux ({:latex})".format(spec1d_sub.flux.unit))
plt.legend()
plt.title("Observed spectrum and best-fitting model template")
plt.show()
../../_images/e4175dd71461945439d610e0e349d7f888b646f67a66a2b5ca4a3bc8ef9e5559.png

Potential next steps:#

- Automatic template *fitting* to get the redshift
- Measure if emission-line profile is consistent with PSF
- Measure line intensities in 2D
- Convert distances from pixels to kpc
- Run the line measurements on a set of lines
- Run the whole procedure on all galaxies detected on a mask