antaress.ANTARESS_analysis.ANTARESS_ana_comm module#
- par_formatting_inst_vis(p_start, fixed_args, elem_prof_inst)[source]#
Parameter formatting: automatic definitions
Defines complementary information for parameters based on their name
- Parameters:
TBD
- Returns:
TBD
- prior_check(par, priors_par, params, args)[source]#
Prior check
Checks that priors are physically valid.
- Parameters:
TBD
- Returns:
None
- prior_def(par, params, args)[source]#
Prior default
Defines default priors
- Parameters:
TBD
- Returns:
None
- MCMC_walkers_check(par, fit_dic, params, args)[source]#
MCMC walkers check
Checks that walkers initialization range is valid.
- Parameters:
par (str) – Parameter to check.
fit_dic (dict) – Dictionary containing fitting information.
params (dict) – Parameter values.
args (dict) – Additional arguments.
- Returns:
None
- model_par_names(par)[source]#
Naming function
Returns name and unit of input variable for plot display.
- Parameters:
None
- Returns:
name_par (str) – parameter name
Example
>>> model_par_names(x) x_name
- model_par_names_txt(par)[source]#
Naming function
Returns plain text name of input variable for plot and log display.
- Parameters:
None
- Returns:
name_par (str) – parameter name
Example
>>> model_par_names(x) x_name
- model_par_units(par)[source]#
Unit function
Returns plain text unit of input variable.
- Parameters:
None
- Returns:
unit_par (str) – parameter unit
Example
>>> model_par_units(x) x_name
- init_model_par(mod_prop, prop_loc, rout_mode, fixed_args, p_start, fit_dic, coord_dic)[source]#
Fits: model parameter initialization.
Initializes parameters that are required by the model function but were left undefined by the user.
- Parameters:
TBD
- Returns:
TBD
- init_fit_routines(rout_mode, gen_dic, system_param, theo_dic, data_dic, fit_prop_dic, prop_loc)[source]#
Joined fits: general initialization.
Initializes properties for the joined fits to stellar and planetary lines and properties.
- Parameters:
TBD
- Returns:
TBD
- init_joined_routines_inst(inst, fit_dic, fixed_args, verbose=True)[source]#
Joined fits: instrument initialization.
Initializes properties for the joined fits to stellar and planetary lines.
- Parameters:
inst (str) – Instrument considered.
fit_prop_dic (dict) – Dictionary containing the parameters of the fitting process.
fixed_args (dict) – Dictionary containing the arguments that will be passed to the fitting function.
- Returns:
None
- init_joined_routines_vis(inst, vis, fit_dic, fixed_args)[source]#
Joined fits: visit initialization.
Initializes properties for the joined fits to stellar and planetary lines.
- Parameters:
inst (str) – Instrument considered.
vis (str) – Visit considered.
fit_dic (dict) – Dictionary containing the parameters of the fitting process.
fixed_args (dict) – Dictionary containing the arguments that will be passed to the fitting function.
- Returns:
None
- init_joined_routines_vis_fit(rout_mode, inst, vis, fit_dic, fixed_args, data_vis, gen_dic, data_dic, coord_dic, theo_dic, data_prop, plot_dic)[source]#
Joined fits: fitted visit initialization.
Initializes properties for the joined fits to stellar and planetary lines.
- Parameters:
TBD
- Returns:
TBD
- PCA_init(fit_dic, inst, vis, fixed_args, data_dic)[source]#
Joined fits: PC initialization
Initializes PC analysis for profile series fits.
- Parameters:
TBD
- Returns:
TBD
- PCA_init_prof(fixed_args, inst, vis, isub, iexp_eff, fit_dic)[source]#
**Joined fits: PC profile initialization **
Initializes PC analysis for profile series fits, in a given profile
- Parameters:
TBD
- Returns:
TBD
- trim_prof(fit_range, fixed_args, edge_bins, inst, vis, isub_ord, fiber, nexp, conv_model=True)[source]#
Fit function: trimming
Defines minimum trimming range based on fitted ranges
- Parameters:
TBD
- Returns:
TBD
- calc_chi_Prof(mod_dic, args)[source]#
Profiles fit function: merit grid
Calculates merit table for optimization.
Because exposures are specific to each visit, defined on different bins, and stored as objects we define the output table as :
chi = concatenate( exp, (obs(exp)-mod(exp))/err(exp)) )
Or the equivalent with the covariance matrix, so that the merit function will compare chi to a table of same size filled with 0 and with errors of 1 in the residual() function (where the condition to use covariance has been set to False for this purpose)
Observed profiles may have gaps, but due to the convolution the model must be calculated over the continuous table and then limited to fitted bins
- Parameters:
TBD
- Returns:
TBD
- sub_calc_chi_Prof(chi, meas_prof, mod_prof, cond_fit, meas_cov, use_cov)[source]#
Profiles fit function: profile merit grid with covariance
Sub-function of calc_chi_Prof()
- Parameters:
TBD
- Returns:
TBD
- com_joint_fits(rout_mode, prop_loc, fit_dic, fixed_args, gen_dic, data_dic, theo_dic, mod_prop, fit_save, coord_dic)[source]#
Wrap-up for time-series fits.
Performs joint fits to time-series.
- Parameters:
TBD
- Returns:
TBD
- FIT_joined_Prof(param, x_tab, args=None)[source]#
Joint profiles fit function: wrap-up
Calls specific model function for optimization of joint profile fitting
- Parameters:
TBD
- Returns:
TBD
- outputs_Prof(inst, vis, iord, coeff_line_dic, mod_prop_dic, fixed_args, param, coeff_line)[source]#
Joint profiles fit function: outputs
Defines outputs for profile fit functions
- Parameters:
fixed_args (dict) – fixed arguments that are not modified within this routine.
- Returns:
TBD
- com_single_fits(fixed_args_all, isub_exp, iexp, inst, data_dic, vis, fit_prop_dic, gen_dic, verbose, cond_def_fit, data_type_gen, idx_force_det, theo_dic, priors, model_prop, fit_dic, flux_cont)[source]#
Single profile analysis routine
Performs fits and measurements of line profiles.
- Parameters:
TBD
- Returns:
TBD
- FIT_single_Prof(param, x_tab, args=None)[source]#
Single profile fit function: wrap-up
Calls specific model function for optimization of single profile fitting
- Parameters:
TBD
- Returns:
TBD
- prior_sini_geom(p_step_loc, fixed_args, prior_func_prop)[source]#
Prior: sin(istar)
Calculates \(\log(p)\) for stellar inclination \(i_\star\).
Assumes isotropic distribution of stellar orientations (equivalent to draw a uniform distribution on \(\cos(i_\star)\)) that favor stellar inclination closer to 90:math:^circ
\[\begin{split}&p(x) = \sin(i_\star(x))/2 \\ &\ln(p(x)) = \ln(0.5 \sin(i_\star(x)))\end{split}\]Where p is normalized by the integral of \(\sin(i_\star)\) over the definition space of \(i_\star\) (0:180:math:^circ), which is 2.
- Parameters:
TBD
- Returns:
TBD
- prior_cosi(p_step_loc, fixed_args, prior_func_prop)[source]#
Prior: cos(ip)
Calculates \(\log(p)\) for orbital inclination \(i_p\).
Prior imposes that \(0<i_p<90`^\circ\)
- Parameters:
TBD
- Returns:
TBD
- prior_sini(p_step_loc, fixed_args, prior_func_prop)[source]#
Prior: sin(ip)
Calculates \(\log(p)\) for orbital inclination \(i_p\).
Imposes Gaussian prior on \(\sin(i_p)\).
- Parameters:
TBD
- Returns:
TBD
- prior_b(p_step_loc, fixed_args, prior_func_prop)[source]#
Prior: b
Calculates \(\log(p)\) for impact parameter.
Prior imposes that \(b<1\).
Used when orbital inclination \(i_p\) and/or semi-major axis \(a/R_\star\) are varied. If only \(i_p\) or \(a/R_\star\) is free to vary, the other property should still be defined as a constant parameter for this prior to be used.
- Parameters:
TBD
- Returns:
TBD
- prior_DR(p_step_loc, fixed_args, prior_func_prop)[source]#
Prior: DR
Calculates \(\log(p)\) for differential rotation coefficients \(\alpha_\mathrm{rot}\) and \(\beta_\mathrm{rot}\).
If both \(\alpha_\mathrm{rot}\) and \(\beta_\mathrm{rot}\) are used, their sum cannot be larger than 1 if the star rotates in the same direction at all latitudes.
- Parameters:
TBD
- Returns:
TBD
- prior_vsini(p_step_loc, fixed_args, prior_func_prop)[source]#
Prior: veq sin(istar)
Calculates \(\log(p)\) for sky-projected stellar rotational velocity \(v_\mathrm{eq} \sin(i_\star)\).
Imposes Gaussian prior on \(v_\mathrm{eq} \sin(i_\star)\). Relevant when \(v_\mathrm{eq}\) and \(i_\star\) are fitted independently, otherwise set a prior on \(v_\mathrm{eq}\) alone.
- Parameters:
TBD
- Returns:
TBD
- prior_vsini_deriv(p_step_loc, fixed_args, prior_func_prop)[source]#
Prior: veq sin(istar)
Calculates \(\log(p)\) for sky-projected stellar rotational velocity \(v_\mathrm{eq} \sin(i_\star)\).
Imposes Gaussian prior on \(v_\mathrm{eq} \sin(i_\star)\). Relevant when \(P_\mathrm{eq}\) and \(i_\star\) are fitted independently, otherwise set a prior on \(v_\mathrm{eq}\) alone.
- Parameters:
TBD
- Returns:
TBD
- prior_contrast(p_step_loc, fixed_args, prior_func_prop)[source]#
Prior: intrinsic contrast
Calculates \(\log(p)\) for the contrast \(C_\mathrm{Intr}\) of intrinsic stellar lines.
Prior imposes that \(0<C_\mathrm{Intr}<1\). The contrast of intrinsic (and a fortiori measured) stellar lines should be physical. See details on the calculation of line contrast in joined_intr_prof() (for now assumes a single planet per visit transits).
- Parameters:
TBD
- Returns:
TBD
- prior_FWHM_vsini(p_step_loc, args, prior_func_prop)[source]#
Prior: intrinsic FWHM
Calculates \(\log(p)\) for the mathrm{FWHM} \(\mathrm{FWHM}_\mathrm{Intr}\) of intrinsic stellar lines.
Prior imposes that \(\mathrm{FWHM}_\mathrm{Intr}^2 + \mathrm{FWHM}_\mathrm{Inst}^2 + v_\mathrm{eq} \sin(i_\star)^2<\mathrm{FWHM}_\mathrm{DI}\). We assume that the intrinsic line cannot be larger than the disk-integrated one, after rotational broadening. This sets an upper limit on \(\mathrm{FWHM}_\mathrm{Intr}\), considering there can be other broadening contributions See details on the calculation of line contrast in joined_intr_prof() (for now assumes intrinsic line must have constant width).
- Parameters:
TBD
- Returns:
TBD
- com_joint_postproc(p_final, fixed_args, fit_dic, merged_chain, gen_dic)[source]#
Wrap-up for post-processing of time-series fits.
Performs joint fits to time-series.
- Parameters:
TBD
- Returns:
TBD
- conv_cosistar(deriv_prop, fixed_args_in, fit_dic_in, p_final_in, merged_chain_in)[source]#
Parameter conversion: cos(i_star)
Converts results from \(\chi^2\), mcmc, or ns fitting: from \(\cos(i_\star)\) to \(i_\star\).
- Parameters:
TBD
- Returns:
TBD
- conv_CF_intr_meas(deriv_prop, inst_list, inst_vis_list, fixed_args, merged_chain, gen_dic, p_final, fit_dic)[source]#
Parameter conversion: line contrast and FWHM
Converts results from \(\chi^2\), mcmc, or ns fitting: from intrinsic width and contrast to measured-like values.
- Parameters:
TBD
- Returns:
TBD