Process a science observation.
For detailed information please refer to the GIRAFFE pipeline user manual.
It is available at http://www.eso.org/pipelines.
Create an object for the recipe giscience.
import cpl
giscience = cpl.Recipe("giscience")
Enable bias removal (bool; default: True) [default=True].
Bias removal method (str; default: ‘PROFILE’) [default=”PROFILE”].
Bias areas to use (Xl0:Xr0:Yl0:Yr0, ... ,Xln:Xrn:Yln:Yrn) (str; default: ‘5:40:0:4095’) [default=”5:40:0:4095”].
Sigma Clipping: sigma threshold factor (float; default: 2.5) [default=2.5].
Sigma Clipping: maximum number of iterations (long; default: 5) [default=5].
Sigma Clipping: minimum fraction of points accepted/total [0.0..1.0] (float; default: 0.8) [default=0.8].
Order of X polynomial fit (CURVE method only) (long; default: 1) [default=1].
Order of Y polynomial fit (CURVE method only) (long; default: 1) [default=1].
Sampling step along X (CURVE method only) (long; default: 1) [default=1].
Sampling step along Y (CURVE method only) (long; default: 1) [default=1].
Extraction method: ‘SUM’, ‘HORNE’ or ‘OPTIMAL’ (str; default: ‘SUM’) [default=”SUM”].
New bias sigma (RON) value for bias and dark corrected image (float; default: -1.0) [default=-1.0].
PSF profile model: `psfexp’, `psfexp2’ (str; default: ‘psfexp2’) [default=”psfexp2”].
Sigma clippging threshold used for rejecting data points during PSF fitting (Horne’s sigma). It is used to reject bad pixels and cosmics. (float; default: 7.0) [default=7.0].
Maximum number of iterations used for fitting the PSF profile. (long; default: 2) [default=2].
Horne extraction method: Number of extra pixels added to the fiber half-width. (long; default: 2) [default=2].
Horne extraction method: Minimum number of points used for the profile fit. It sets the lower limit of data points for the pixel rejection. (long; default: 3) [default=3].
Optimal extraction method: Minimum fraction of the data points used for fitting the fiber profiles. It sets the lower limit for the pixel rejection. (float; default: 0.9) [default=0.9].
Optimal extraction method: Factor by which the fiber PSF half width is multiplied. Adjacent spectra within this area are assumed to affect the spectrum being extracted. (float; default: 3.0) [default=3.0].
Optimal extraction method: Order of the polynomial background model, which is fitted for each wavelength bin along the spatial direction. (long; default: 2) [default=2].
Controls the flat field correction. (bool; default: True) [default=True].
Controls the fiber to fiber transmission correction. (bool; default: False) [default=False].
Method to use : `linear’ or `spline’ (str; default: ‘linear’) [default=”linear”].
Use x residuals during rebinning? `true’/`false’ (bool; default: True) [default=True].
Lambda step size, only used if scaling method is ‘linear’ (float; default: 0.005) [default=0.005].
Scaling method: `log’ or `linear’ (str; default: ‘linear’) [default=”linear”].
Size of output rebinned spectra, 0 means calculate size based on wavelength range and lambda stepsize (long; default: 0) [default=0].
Rebinning range: `setup’ or `common’ (str; default: ‘setup’) [default=”setup”].
Enable simultaneous wavelength calibration correction. (bool; default: True) [default=True].
Slit geometry calibration maximum number of iterations. (long; default: 1) [default=1].
Maximum allowed pixel value. To be effective it must be larger than 0. (float; default: 10000.0) [default=10000.0].
Cross-correlation step. (float; default: -0.005) [default=-0.005].
Restricts the cross-correlation to the given domain. (str; default: ‘0.,0.’) [default=”0.,0.”].
Delta RV limits of the cross-correlation window in km/s. (str; default: ‘-200.,200.’) [default=”-200.,200.”].
Maximum number of iterations used for the RV determination. (long; default: 3) [default=3].
Data window width factor. The FWHM times this value determines the window width. (float; default: 1.5) [default=1.5].
Peak model fit maximum number of iterations. (long; default: 50) [default=50].
Cross-correlation peak fit maximum number of tests (long; default: 7) [default=7].
Cross-correlation peak fit minimum chi-square difference. (float; default: 0.0001) [default=0.0001].
Minimum wavelength for image reconstruction (float; default: 0.0) [default=0.0].
Maximum wavelength for image reconstruction (float; default: 0.0) [default=0.0].
Turns data cube creation on and off (bool; default: True) [default=True].
Selects the file format for cubes (str; default: ‘single’) [default=”single”].
The following code snippet shows the default settings for the available parameters.
import cpl
giscience = cpl.Recipe("giscience")
giscience.param.remove_bias = True
giscience.param.bsremove_method = "PROFILE"
giscience.param.bsremove_areas = "5:40:0:4095"
giscience.param.bsremove_sigma = 2.5
giscience.param.bsremove_niter = 5
giscience.param.bsremove_mfrac = 0.8
giscience.param.bsremove_xorder = 1
giscience.param.bsremove_yorder = 1
giscience.param.bsremove_xstep = 1
giscience.param.bsremove_ystep = 1
giscience.param.extr_method = "SUM"
giscience.param.extr_ron = -1.0
giscience.param.extr_psfmodel = "psfexp2"
giscience.param.extr_psfsigma = 7.0
giscience.param.extr_psfniter = 2
giscience.param.extr_hewidth = 2
giscience.param.extr_hmingood = 3
giscience.param.extr_omfrac = 0.9
giscience.param.extr_owfactor = 3.0
giscience.param.extr_obkgorder = 2
giscience.param.flat_apply = True
giscience.param.transmission_apply = False
giscience.param.rbin_method = "linear"
giscience.param.rbin_xresid = True
giscience.param.rbin_lstep = 0.005
giscience.param.rbin_scmethod = "linear"
giscience.param.rbin_size = 0
giscience.param.rbin_range = "setup"
giscience.param.siwc_apply = True
giscience.param.scal_cniter = 1
giscience.param.scal_zmax = 10000.0
giscience.param.scal_cstep = -0.005
giscience.param.scal_cdomain = "0.,0."
giscience.param.scal_rvlimits = "-200.,200."
giscience.param.scal_rvniter = 3
giscience.param.scal_rvwfactor = 1.5
giscience.param.scal_pfniter = 50
giscience.param.scal_pfntest = 7
giscience.param.scal_pfdchisq = 0.0001
giscience.param.recon_min = 0.0
giscience.param.recon_max = 0.0
giscience.param.recon_cube = True
giscience.param.recon_format = "single"
You may also set or overwrite some or all parameters by the recipe parameter param, as shown in the following example:
import cpl
giscience = cpl.Recipe("giscience")
[...]
res = giscience( ..., param = {"remove_bias":True, "bsremove_method":"PROFILE"})
See also
cpl.Recipe for more information about the recipe object.
Please report any problems to Giraffe Pipeline. Alternatively, you may send a report to the ESO User Support Department.
This file is part of the GIRAFFE Instrument Pipeline Copyright (C) 2002-2014 European Southern Observatory
This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version.
This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
Code author: Giraffe Pipeline <usd-help@eso.org>