xsh_compute_response.c File Reference

#include <gsl/gsl_spline.h>
#include <gsl/gsl_errno.h>
#include <gsl/gsl_bspline.h>
#include <gsl/gsl_multifit.h>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>
#include <gsl/gsl_statistics.h>
#include <xsh_dfs.h>
#include <xsh_error.h>
#include <xsh_msg.h>
#include <cpl.h>
#include <string.h>
#include <time.h>
#include <xsh_utils_table.h>
#include <xsh_data_star_flux.h>
#include <xsh_data_atmos_ext.h>
#include <xsh_data_spectrum.h>
#include <xsh_pfits.h>
#include <xsh_utils.h>
#include <xsh_utils_wrappers.h>
#include <xsh_drl.h>
#include <xsh_utils_efficiency.h>
#include <xsh_efficiency_response.h>
#include <xsh_utils_response.h>
#include <xsh_utils_vector.h>

Go to the source code of this file.

Macros
#define	INTERPOL_WSTEP_NM   2
#define	FILTER_MEDIAN_HSIZE   3

Functions
static int	find_lambda_idx (double lambda, double *wave, int from, int to, double step)
static void	find_lambda_idx_limit (double min, double max, double *spectrum, int from, int to, double step, int *if0, int *if1)
static cpl_error_code	xsh_interpolate_atm_ext (xsh_star_flux_list *star_list, cpl_table *atmos_ext_tab, xsh_instrument *instrument, double **atmos_lambda, double **atmos_K)
static cpl_error_code	xsh_interpolate_high_abs_regions (xsh_star_flux_list *star_list, xsh_star_flux_list *resp_list, HIGH_ABS_REGION *phigh)
static cpl_error_code	xsh_response_crea_ascii (xsh_star_flux_list *resp_list, xsh_star_flux_list *star_list, double *lambda_spectrum, double *flux_spectrum, double *flux_added)
static cpl_error_code	xsh_flux_integrate_and_corr_for_badpix (int npix_in_interval, double *flux_spectrum, int *qual_spectrum, int if0, int if1, int i, double **flux_added, int *npixels, int *nbad)
static cpl_error_code	xsh_spectrum_correct (xsh_star_flux_list *star_list, xsh_spectrum *spectrum, double *atmos_K, XSH_ARM the_arm, double airmass, double exptime, double gain, xsh_star_flux_list **obj_list)
static xsh_star_flux_list *	xsh_response_calculate (xsh_star_flux_list *star_list, xsh_star_flux_list **obj_list, xsh_spectrum *spectrum, cpl_table *atmos_ext_tab, XSH_ARM the_arm, double *atmos_K, double airmass, double exptime, double gain)
static xsh_star_flux_list *	do_compute (xsh_star_flux_list *star_list, xsh_star_flux_list **obj_list, xsh_spectrum *spectrum, cpl_table *atmos_ext_tab, HIGH_ABS_REGION *phigh, double airmass, double exptime, double gain, xsh_instrument *instrument)
cpl_error_code	xsh_response_merge_obj_std_info (cpl_frame *result, xsh_star_flux_list *star_list, xsh_star_flux_list *obj_list)
static xsh_star_flux_list *	xsh_obs_std_correct (cpl_frame *obs_std_star, xsh_star_flux_list *ref_std_star_list, cpl_frame *atmos_ext, HIGH_ABS_REGION *phigh, xsh_instrument *instrument)
static cpl_error_code	xsh_std_star_spectra_to_vector_range (xsh_star_flux_list *obs_std_star_list, const double wmin, const double wmax, cpl_vector **vec_wave, cpl_vector **vec_flux)
cpl_error_code	xsh_sort_double_pairs (double *u1, double *u2, cpl_size n)
	Sort an array u1 of doubles, and permute an array u2 in the same way as u1 is permuted.
static cpl_error_code	xsh_select_line_core (cpl_vector *wave, cpl_vector *flux, const double wguess, const double wrange, cpl_vector **xfit, cpl_vector **yfit)
static double	xsh_std_star_spectra_correlate (xsh_star_flux_list *obs_std_star_list, xsh_star_flux_list *ref_std_star_list, xsh_rv_ref_wave_param *w)
static xsh_star_flux_list *	xsh_bspline_fit_interpol (xsh_star_flux_list *ref_std_star_list, cpl_table *resp_fit_points, HIGH_ABS_REGION *phigh, xsh_instrument *inst)
static xsh_star_flux_list *	xsh_bspline_fit_smooth (xsh_star_flux_list *ref_std_star_list, HIGH_ABS_REGION *phigh, xsh_instrument *inst)
double *	xsh_bspline_interpolate_data_at_pos (double *w_data, double *f_data, const int n_data, double *w_pos, const int n_pos)
double *	xsh_bspline_fit_smooth_data (double *wave, double *flux, const int size, HIGH_ABS_REGION *phigh, xsh_instrument *inst, const int fit_region)
double *	xsh_bspline_fit_smooth_data2 (double *wave, double *flux, const int size, HIGH_ABS_REGION *phigh, xsh_instrument *inst, const int fit_region)
cpl_frame *	xsh_compute_response2 (cpl_frame *obs_std_star, cpl_frame *flux_std_star_cat, cpl_frame *atmos_ext, cpl_frame *high_abs, cpl_frame *resp_fit_points, cpl_frame *tell_mod_cat, xsh_instrument *instrument, double exptime, const int tell_corr)
cpl_frame *	xsh_compute_response (cpl_frame *spectrum_frame, cpl_frame *flux_std_star_cat_frame, cpl_frame *atmos_ext_frame, cpl_frame *high_abs_frame, xsh_instrument *instrument, double exptime)
cpl_frame *	xsh_compute_response_ord (cpl_frame *spectrum_frame, cpl_frame *flux_std_star_cat_frame, cpl_frame *atmos_ext_frame, cpl_frame *high_abs_win_frame, xsh_instrument *instrument, double exptime)

Macro Definition Documentation

FILTER_MEDIAN_HSIZE

#define FILTER_MEDIAN_HSIZE   3

Definition at line 81 of file xsh_compute_response.c.

INTERPOL_WSTEP_NM

#define INTERPOL_WSTEP_NM   2

Definition at line 80 of file xsh_compute_response.c.

Function Documentation

do_compute()

static xsh_star_flux_list * do_compute ( xsh_star_flux_list *  star_list, xsh_star_flux_list **  obj_list, xsh_spectrum *  spectrum, cpl_table *  atmos_ext_tab, HIGH_ABS_REGION *  phigh, double  airmass, double  exptime, double  gain, xsh_instrument *  instrument  )

static

Integrate the flux on the same sampling that the star_list, taking into account the bad pixels, apply atmos correction, and divide by the star_list table.

Parameters

star_list	Table of ref star flux
star_list	Table of obs obj flux
spectrum	1./CPL/DOC/cpl_um_2.0.1.psD merged spectrum of standard star
atmos_ext_table	Atmos Extinction table
airmass	observed object airmass mean
exptime	Exposure time
gain	detector's gain
instrument	instrument arm and lamp setting

Returns

The computed response frame (table)

Definition at line 1280 of file xsh_compute_response.c.

References check, exptime, instrument, xsh_star_flux_list::lambda, HIGH_ABS_REGION::lambda_max, HIGH_ABS_REGION::lambda_min, xsh_star_flux_list::size, XSH_ASSURE_NOT_NULL, XSH_FREE, xsh_instrument_get_arm(), xsh_interpolate_atm_ext(), xsh_interpolate_high_abs_regions(), xsh_msg_dbg_low, xsh_msg_dbg_medium, and xsh_response_calculate().

Referenced by xsh_compute_response(), and xsh_compute_response_ord().

find_lambda_idx()

static int find_lambda_idx ( double  lambda, double *  wave, int  from, int  to, double  step  )

static

Definition at line 96 of file xsh_compute_response.c.

References step.

Referenced by xsh_interpolate_high_abs_regions().

find_lambda_idx_limit()

static void find_lambda_idx_limit ( double  min, double  max, double *  spectrum, int  from, int  to, double  step, int *  if0, int *  if1  )

static

Definition at line 197 of file xsh_compute_response.c.

References max, step, xsh_msg_dbg_high, and xsh_msg_dbg_low.

Referenced by xsh_response_calculate().

xsh_bspline_fit_interpol()

static xsh_star_flux_list * xsh_bspline_fit_interpol ( xsh_star_flux_list *  ref_std_star_list, cpl_table *  resp_fit_points, HIGH_ABS_REGION *  phigh, xsh_instrument *  inst  )

static

Definition at line 2244 of file xsh_compute_response.c.

References check, xsh_star_flux_list::flux, xsh_star_flux_list::lambda, HIGH_ABS_REGION::lambda_max, HIGH_ABS_REGION::lambda_min, xsh_star_flux_list::size, XSH_ARM_NIR, XSH_ARM_UVB, XSH_ARM_VIS, xsh_bspline_interpolate_data_at_pos(), XSH_CALLOC, XSH_FREE, xsh_free_table(), xsh_instrument_get_arm(), xsh_msg, xsh_star_flux_list_duplicate(), xsh_star_flux_list_get_flux(), and xsh_star_flux_list_get_lambda().

Referenced by xsh_compute_response2().

xsh_bspline_fit_smooth()

static xsh_star_flux_list * xsh_bspline_fit_smooth ( xsh_star_flux_list *  ref_std_star_list, HIGH_ABS_REGION *  phigh, xsh_instrument *  inst  )

static

Definition at line 2488 of file xsh_compute_response.c.

References xsh_star_flux_list::flux, xsh_star_flux_list::lambda, HIGH_ABS_REGION::lambda_max, HIGH_ABS_REGION::lambda_min, n, order, sigma, xsh_star_flux_list::size, x, XSH_ARM_NIR, XSH_ARM_UVB, XSH_ARM_VIS, xsh_free_table(), xsh_instrument_get_arm(), xsh_msg, xsh_star_flux_list_create(), and y.

Referenced by xsh_compute_response2().

xsh_bspline_fit_smooth_data()

double * xsh_bspline_fit_smooth_data	(	double *	wave,
		double *	flux,
		const int	size,
		HIGH_ABS_REGION *	phigh,
		xsh_instrument *	inst,
		const int	fit_region
	)

Definition at line 2715 of file xsh_compute_response.c.

References HIGH_ABS_REGION::lambda_max, HIGH_ABS_REGION::lambda_min, n, order, sigma, size, x, XSH_ARM_NIR, XSH_ARM_UVB, XSH_ARM_VIS, xsh_free_table(), xsh_instrument_get_arm(), xsh_msg, and y.

xsh_bspline_fit_smooth_data2()

double * xsh_bspline_fit_smooth_data2	(	double *	wave,
		double *	flux,
		const int	size,
		HIGH_ABS_REGION *	phigh,
		xsh_instrument *	inst,
		const int	fit_region
	)

Definition at line 2911 of file xsh_compute_response.c.

References n, order, sigma, size, x, XSH_ARM_NIR, XSH_ARM_UVB, XSH_ARM_VIS, xsh_free_table(), xsh_instrument_get_arm(), xsh_msg, and y.

xsh_bspline_interpolate_data_at_pos()

double * xsh_bspline_interpolate_data_at_pos	(	double *	w_data,
		double *	f_data,
		const int	n_data,
		double *	w_pos,
		const int	n_pos
	)

Definition at line 2664 of file xsh_compute_response.c.

References xsh_msg.

Referenced by xsh_bspline_fit_interpol(), xsh_fit_spline1(), xsh_fit_spline2(), and xsh_telluric_model_eval().

xsh_compute_response()

cpl_frame * xsh_compute_response	(	cpl_frame *	spectrum_frame,
		cpl_frame *	flux_std_star_cat_frame,
		cpl_frame *	atmos_ext_frame,
		cpl_frame *	high_abs_frame,
		xsh_instrument *	instrument,
		double	exptime
	)

Compute the instrument response using the rectified frame and the Standard star flux frame. The input should include also the Master Response Function and the Athmospheric Extinction Curve. Also should be produced the Efficiency frame.

Parameters

spectrum_frame	Rectified frame (1D spectrum)
flux_std_star_cat_frame	Standard star flux catalog frame
atmos_ext_frame	Atmospheric Extinction Frame
instrument	Instrument strcture
exptime	Exposure time

Returns

Frame of the computed response

< computed response frame

< Std star flux

< Std star flux

< 1d spectrum

< 1d spectrum

Definition at line 3523 of file xsh_compute_response.c.

References check, check_msg, do_compute(), exptime, instrument, INTERPOL_WSTEP_NM, STAR_MATCH_DEPSILON, xsh_add_temporary_file(), XSH_ARM_NIR, XSH_ASSURE_NOT_NULL, XSH_ATMOS_EXT, xsh_atmos_ext_list_free(), xsh_atmos_ext_list_load(), xsh_fill_high_abs_regions(), xsh_frame_product(), xsh_frame_sci_get_ra_dec_airmass(), xsh_free_frame(), xsh_free_propertylist(), xsh_free_table(), xsh_instrument_arm_tostring(), xsh_instrument_get_arm(), xsh_instrument_mode_tostring(), xsh_msg, xsh_msg_warning, xsh_parse_catalog_std_stars(), xsh_pfits_get_cdelt1(), xsh_pfits_get_crval1(), xsh_pfits_get_gain(), xsh_pfits_get_naxis1(), xsh_response_merge_obj_std_info(), xsh_spectrum_free(), xsh_spectrum_interpolate(), xsh_spectrum_load(), xsh_spectrum_resample(), xsh_spectrum_save(), xsh_star_flux_list_free(), xsh_star_flux_list_load(), and xsh_star_flux_list_save().

Referenced by main(), and xsh_util_integrate().

xsh_compute_response2()

cpl_frame * xsh_compute_response2	(	cpl_frame *	obs_std_star,
		cpl_frame *	flux_std_star_cat,
		cpl_frame *	atmos_ext,
		cpl_frame *	high_abs,
		cpl_frame *	resp_fit_points,
		cpl_frame *	tell_mod_cat,
		xsh_instrument *	instrument,
		double	exptime,
		const int	tell_corr
	)

Compute the instrument response using the rectified frame and the Standard star flux frame. The input should include also the Master Response Function and the Athmospheric Extinction Curve. Also should be produced the Efficiency frame.

Parameters

obbs_std_star	Rectified frame (1D spectrum)
flux_std_star_cat	Standard star flux catalog frame
atmos_ext	Atmospheric Extinction Frame
instrument	Instrument structure
exptime	Exposure time
tell_corr	apply telluric correction

Returns

Frame of the computed response

Definition at line 3145 of file xsh_compute_response.c.

References check, xsh_star_flux_list::flux, instrument, xsh_star_flux_list::lambda, xsh_star_flux_list::size, STAR_MATCH_DEPSILON, xsh_rv_ref_wave_param::wguess, xsh_add_temporary_file(), XSH_ARM_UVB, XSH_ASSURE_NOT_NULL, xsh_bspline_fit_interpol(), xsh_bspline_fit_smooth(), xsh_fill_high_abs_regions(), xsh_frame_product(), xsh_frame_sci_get_ra_dec_airmass(), xsh_frame_spectrum_save(), xsh_free_frame(), xsh_free_propertylist(), xsh_free_table(), xsh_get_basename(), xsh_instrument_arm_tostring(), xsh_instrument_get_arm(), xsh_instrument_mode_tostring(), xsh_msg, xsh_msg_warning, xsh_obs_std_correct(), xsh_parse_catalog_std_stars(), xsh_pfits_get_cdelt1(), xsh_pfits_get_crval1(), xsh_pfits_get_naxis1(), xsh_response_merge_obj_std_info(), xsh_rv_ref_wave_init(), xsh_rv_ref_wave_param_create(), xsh_rv_ref_wave_param_destroy(), xsh_spectrum_free(), xsh_spectrum_interpolate_linear(), xsh_spectrum_load(), xsh_star_flux_list_divide(), xsh_star_flux_list_duplicate(), xsh_star_flux_list_extrapolate_wave_end(), xsh_star_flux_list_filter_median(), xsh_star_flux_list_free(), xsh_star_flux_list_load(), xsh_star_flux_list_load_spectrum(), xsh_star_flux_list_save(), xsh_star_flux_list_to_frame(), xsh_std_star_spectra_correlate(), xsh_table_shift_rv(), and xsh_telluric_model_eval().

Referenced by xsh_respon_slit_nod(), xsh_respon_slit_offset(), xsh_respon_slit_stare(), and xsh_util_compute_response().

xsh_compute_response_ord()

cpl_frame * xsh_compute_response_ord	(	cpl_frame *	spectrum_frame,
		cpl_frame *	flux_std_star_cat_frame,
		cpl_frame *	atmos_ext_frame,
		cpl_frame *	high_abs_win_frame,
		xsh_instrument *	instrument,
		double	exptime
	)

Compute the instrument response using the rectified frame and the Standard star flux frame. The input should include also the Master Response Function and the Athmospheric Extinction Curve. Also should be produced the Efficiency frame.

Parameters

spectrum_frame	Rectified frame (1D spectrum)
flux_std_star_cat_frame	Standard star flux catalog frame
atmos_ext_frame	Atmospheric Extinction Frame
instrument	Instrument strcture
exptime	Exposure time

Returns

Frame of the computed response

< computed response frame

< Std star flux

< obj star flux

< 1d spectrum

< 1d spectrum

Definition at line 3709 of file xsh_compute_response.c.

References check, check_msg, do_compute(), exptime, instrument, INTERPOL_WSTEP_NM, STAR_MATCH_DEPSILON, xsh_add_temporary_file(), XSH_ARM_NIR, XSH_ASSURE_NOT_NULL, xsh_atmos_ext_list_free(), xsh_atmos_ext_list_load(), xsh_fill_high_abs_regions(), xsh_frame_product(), xsh_frame_sci_get_ra_dec_airmass(), xsh_free_frame(), xsh_free_propertylist(), xsh_free_table(), xsh_instrument_arm_tostring(), xsh_instrument_get_arm(), xsh_instrument_mode_tostring(), xsh_msg, xsh_msg_warning, xsh_parse_catalog_std_stars(), xsh_pfits_get_cdelt1(), xsh_pfits_get_crval1(), xsh_pfits_get_gain(), xsh_pfits_get_naxis1(), xsh_spectrum_free(), xsh_spectrum_interpolate(), xsh_spectrum_load_order(), xsh_spectrum_save_order(), xsh_star_flux_list_free(), xsh_star_flux_list_load(), and xsh_star_flux_list_save_order().

Referenced by xsh_respon_slit_nod(), xsh_respon_slit_offset(), and xsh_respon_slit_stare().

xsh_flux_integrate_and_corr_for_badpix()

static cpl_error_code xsh_flux_integrate_and_corr_for_badpix ( int  npix_in_interval, double *  flux_spectrum, int *  qual_spectrum, int  if0, int  if1, int  i, double **  flux_added, int *  npixels, int *  nbad  )

static

Definition at line 595 of file xsh_compute_response.c.

References QFLAG_GOOD_PIXEL, and XSH_ASSURE_NOT_NULL_MSG.

Referenced by xsh_response_calculate().

xsh_interpolate_atm_ext()

static cpl_error_code xsh_interpolate_atm_ext ( xsh_star_flux_list *  star_list, cpl_table *  atmos_ext_tab, xsh_instrument *  instrument, double **  atmos_lambda, double **  atmos_K  )

static

Interpolate atms extinction at same sampling step as ref std star

Parameters

	star_list	Table of ref star flux
[in]	atmos_ext_tab	atmospheric extinction table
[in]	instrument	instrument arm and lamp setting
[out]	atmos_lambda	interpolated wave array
[out]	atmos_K	interpolated atmospheric extinction K

Returns

cpl error code

Definition at line 237 of file xsh_compute_response.c.

References instrument, xsh_star_flux_list::lambda, xsh_star_flux_list::size, XSH_ARM_NIR, XSH_ASSURE_NOT_NULL_MSG, XSH_ATMOS_EXT_LIST_COLNAME_K, XSH_ATMOS_EXT_LIST_COLNAME_OLD, XSH_CALLOC, xsh_data_interpolate(), xsh_instrument_get_arm(), and xsh_msg_warning.

Referenced by do_compute(), and xsh_obs_std_correct().

xsh_interpolate_high_abs_regions()

static cpl_error_code xsh_interpolate_high_abs_regions ( xsh_star_flux_list *  star_list, xsh_star_flux_list *  resp_list, HIGH_ABS_REGION *  phigh  )

static

Interpolate high absorbtion regions

Parameters

[in]	star_list	Table of ref star flux
[out]	resp_list	Table of response values
[in]	HIGH_ABS_REGION	pointer to hight abs regions

Returns

cpl error code

Definition at line 308 of file xsh_compute_response.c.

References find_lambda_idx(), xsh_star_flux_list::flux, xsh_star_flux_list::lambda, HIGH_ABS_REGION::lambda_max, HIGH_ABS_REGION::lambda_min, xsh_star_flux_list::size, XSH_ASSURE_NOT_NULL_MSG, xsh_msg, xsh_msg_dbg_low, xsh_msg_dbg_medium, and xsh_msg_debug.

Referenced by do_compute().

xsh_obs_std_correct()

static xsh_star_flux_list * xsh_obs_std_correct ( cpl_frame *  obs_std_star, xsh_star_flux_list *  ref_std_star_list, cpl_frame *  atmos_ext, HIGH_ABS_REGION *  phigh, xsh_instrument *  instrument  )

static

Corrects the observed spectrum by gain, exptime, spectral bin size, atmospheric extinction, airmass. at the same wavelengths as the ones at which it is sampled the ref std star

Parameters

obs_std_star	Rectified frame (1D spectrum)
ref_std_star_list	Ref std star structure
atmos_ext	Atmospheric Extinction Frame
phigh	pointer to high abs region structure
instrument	Instrument structure

Returns

structure with corrected observed std starspectrum

< 1d spectrum

Definition at line 1423 of file xsh_compute_response.c.

References check, check_msg, exptime, instrument, XSH_ARM_NIR, XSH_ATMOS_EXT, xsh_atmos_ext_list_free(), xsh_atmos_ext_list_load(), xsh_frame_sci_get_ra_dec_airmass(), XSH_FREE, xsh_free_propertylist(), xsh_free_table(), xsh_instrument_arm_tostring(), xsh_instrument_get_arm(), xsh_interpolate_atm_ext(), xsh_msg_warning, xsh_pfits_get_airm_mean(), xsh_pfits_get_exptime(), xsh_pfits_get_gain(), xsh_spectrum_correct(), xsh_spectrum_free(), and xsh_spectrum_load().

Referenced by xsh_compute_response2().

xsh_response_calculate()

static xsh_star_flux_list * xsh_response_calculate ( xsh_star_flux_list *  star_list, xsh_star_flux_list **  obj_list, xsh_spectrum *  spectrum, cpl_table *  atmos_ext_tab, XSH_ARM  the_arm, double *  atmos_K, double  airmass, double  exptime, double  gain  )

static

Definition at line 843 of file xsh_compute_response.c.

References binx, biny, check, exptime, FILTER_MEDIAN_HSIZE, find_lambda_idx_limit(), xsh_star_flux_list::flux, xsh_star_flux_list::lambda, max, xsh_star_flux_list::size, xsh_add_temporary_file(), XSH_ARM_NIR, XSH_CALLOC, xsh_debug_level_get(), XSH_DEBUG_LEVEL_LOW, xsh_flux_integrate_and_corr_for_badpix(), XSH_FREE, xsh_free_frame(), xsh_free_table(), xsh_msg_dbg_low, xsh_msg_dbg_medium, xsh_response_crea_ascii(), xsh_spectrum_get_flux(), xsh_spectrum_get_lambda_max(), xsh_spectrum_get_lambda_min(), xsh_spectrum_get_lambda_step(), xsh_spectrum_get_qual(), xsh_spectrum_get_size(), xsh_star_flux_list_create(), xsh_star_flux_list_filter_median(), and xsh_star_flux_list_save().

Referenced by do_compute().

xsh_response_crea_ascii()

static cpl_error_code xsh_response_crea_ascii ( xsh_star_flux_list *  resp_list, xsh_star_flux_list *  star_list, double *  lambda_spectrum, double *  flux_spectrum, double *  flux_added  )

static

Integrate the flux on the same sampling that the star_list, taking into account the bad pixels, apply atmos correction, and divide by the star_list table.

Parameters

resp_list	Table of response values
star_list	Table of ref star flux
lambda_spectrum	array with lambda values
flux_spectrum	array with flux values
flux_added	array with integrated flux values

Returns

cpl error code

Definition at line 543 of file xsh_compute_response.c.

References xsh_star_flux_list::flux, xsh_star_flux_list::lambda, xsh_star_flux_list::size, and XSH_ASSURE_NOT_NULL_MSG.

Referenced by xsh_response_calculate().

xsh_response_merge_obj_std_info()

cpl_error_code xsh_response_merge_obj_std_info	(	cpl_frame *	result,
		xsh_star_flux_list *	star_list,
		xsh_star_flux_list *	obj_list
	)

Definition at line 1362 of file xsh_compute_response.c.

References check, xsh_star_flux_list::flux, xsh_free_propertylist(), and xsh_free_table().

Referenced by xsh_compute_response(), and xsh_compute_response2().

xsh_select_line_core()

static cpl_error_code xsh_select_line_core ( cpl_vector *  wave, cpl_vector *  flux, const double  wguess, const double  wrange, cpl_vector **  xfit, cpl_vector **  yfit  )

static

Definition at line 1674 of file xsh_compute_response.c.

References check, and size.

Referenced by xsh_std_star_spectra_correlate().

xsh_sort_double_pairs()

cpl_error_code xsh_sort_double_pairs	(	double *	u1,
		double *	u2,
		cpl_size	n
	)

Sort an array u1 of doubles, and permute an array u2 in the same way as u1 is permuted.

Parameters

[in,out]	u1	Pointer to the first array.
[in,out]	u2	Pointer to the second array.
	n	The common length of both arrays.

Returns

CPL_ERROR_NONE or the appropriate error code.

Definition at line 1644 of file xsh_compute_response.c.

References n.

Referenced by xsh_std_star_spectra_correlate().

xsh_spectrum_correct()

static cpl_error_code xsh_spectrum_correct ( xsh_star_flux_list *  star_list, xsh_spectrum *  spectrum, double *  atmos_K, XSH_ARM  the_arm, double  airmass, double  exptime, double  gain, xsh_star_flux_list **  obj_list  )

static

This function corrects a spectrum for contributes due to atmospheric extinction, airmass,exptime,gain

Parameters

star_list	Table of ref star flux
spectrum	Table of observed star flux
atmos_K	atmospheric extinction curve
the_arm	observing arm
airmass	airmass value
exptime	exptime (dit*ndit) value
gain	gain value
obj_list	Table of observed star flux corrected by atmospheric extinction, airmass,exptime,gain

Returns

cpl error code

Definition at line 659 of file xsh_compute_response.c.

References binx, biny, check, exptime, xsh_star_flux_list::flux, xsh_spectrum::flux_header, xsh_star_flux_list::lambda, xsh_star_flux_list::size, size, xsh_add_temporary_file(), XSH_ARM_NIR, XSH_CALLOC, XSH_FREE, xsh_free_frame(), xsh_free_table(), xsh_msg, xsh_msg_dbg_low, xsh_pfits_get_binx(), xsh_pfits_get_biny(), xsh_spectrum_get_flux(), xsh_spectrum_get_lambda_max(), xsh_spectrum_get_lambda_min(), xsh_spectrum_get_lambda_step(), xsh_spectrum_get_size(), xsh_star_flux_list_create(), and xsh_star_flux_list_save().

Referenced by xsh_obs_std_correct().

xsh_std_star_spectra_correlate()

static double xsh_std_star_spectra_correlate ( xsh_star_flux_list *  obs_std_star_list, xsh_star_flux_list *  ref_std_star_list, xsh_rv_ref_wave_param *  w  )

static

Definition at line 1743 of file xsh_compute_response.c.

References check, xsh_star_flux_list::header, xsh_rv_ref_wave_param::ipol_hbox, xsh_rv_ref_wave_param::ipol_wmax_min, xsh_rv_ref_wave_param::ipol_wmin_max, xsh_rv_ref_wave_param::range_wmax, xsh_rv_ref_wave_param::range_wmin, xsh_rv_ref_wave_param::wguess, xsh_free_polynomial(), xsh_free_vector(), xsh_msg, xsh_pfits_get_cdelt1(), xsh_pfits_get_crval1(), xsh_pfits_get_naxis1(), xsh_polynomial_fit_1d_create(), xsh_select_line_core(), xsh_sort_double_pairs(), xsh_std_star_spectra_to_vector_range(), and xsh_vector_fit_slope().

Referenced by xsh_compute_response2().

xsh_std_star_spectra_to_vector_range()

static cpl_error_code xsh_std_star_spectra_to_vector_range ( xsh_star_flux_list *  obs_std_star_list, const double  wmin, const double  wmax, cpl_vector **  vec_wave, cpl_vector **  vec_flux  )

static

Definition at line 1594 of file xsh_compute_response.c.

References xsh_star_flux_list::flux, xsh_star_flux_list::lambda, xsh_star_flux_list::size, and size.

Referenced by xsh_std_star_spectra_correlate().