antaress.ANTARESS_analysis.ANTARESS_joined_star module#
- joined_Star_ana(glob_fit_dic, system_param, theo_dic, data_dic, gen_dic, plot_dic, coord_dic, data_prop)[source]#
Joined stellar fits
Wrap-up function to call joint fits of stellar properties and profiles
- Parameters:
TBD
- Returns:
TBD
- MAIN_single_anaprof(vis_mode, data_type, data_dic, gen_dic, inst, vis, coord_dic, theo_dic, plot_dic, system_param, data_prop)[source]#
Main analysis routine for single profiles
Initializes and calls analysis of line profiles
profiles in their input rest frame, original exposures, for all formats
profiles in their input or star (if aligned) rest frame, original exposures, converted from 2D->1D
profiles in their input or star (if aligned) rest frame, binned exposures, all formats
profiles in the star rest frame, binned exposures, all formats
Formats are CCF profiles, 1D spectra, or a specific 2D spectral order. Each profile is analyzed independently from the others. The routine can also be used to fit the unocculted disk-integrated stellar profile using best estimates for the intrinsic stellar profiles (measured, theoretical, imported)
- Parameters:
TBD
- Returns:
TBD
- single_prof(param, fixed_args)[source]#
Generic single profile model function
Computes a single spectral line profiles (by opposition to joint routines computing multiple profiles). The profile is built using a line profile model and continuum model. If requested, it is then convolved at the instrumental resolution, converted from RV to spectral space, and resampled. The function is defined in this file to avoid circular imports.
The line profile model functions must return several arguments (or a list of one) with the model as first argument. The profile is calculated over a continuous table to allow for convolution and use of covariance matrix (fitted pixels are accounted for directly in the minimization routine).
- Parameters:
TBD
- Returns:
TBD
- main_joined_DIProp(rout_mode, fit_prop_dic, gen_dic, system_param, theo_dic, plot_dic, coord_dic, data_dic, data_prop)[source]#
Joined disk-integrated stellar property fits
Main routine to fit time-series of properties derived from disk-integrated profiles as a function of various parameters, to search for systematic trends
Individual visits can still be fitted, but the use of a joint model over instruments and visits allows better characterizing stellar variations.
Results of the analysis are saved, to be used in the detrending module.
- Parameters:
TBD
- Returns:
TBD
- FIT_joined_DIProp(param, x_tab, args=None)[source]#
Fit function: joined global stellar property
Calls corresponding model function for optimization
- Parameters:
TBD
- Returns:
TBD
- mod_DIProp(param, system_param, fixed_args, inst, vis, n_coord)[source]#
Model function: global stellar property
Defines the model for global stellar property over a visit.
- Parameters:
TBD
- Returns:
TBD
- joined_DIProp(param, fixed_args)[source]#
Model function: joined global stellar property
Defines the joined model for global stellar properties
- Parameters:
TBD
- Returns:
TBD
- joined_DIProf(param, fixed_args)[source]#
Model function: joined disk-integrated profiles
Defines the joined model for disk-integrated profiles.
The stellar disk-integrated profile is computed in the star rest frame.
- Parameters:
TBD
- Returns:
TBD
- joined_DIProf_exp(isub_exp, iexp, fib_exp, inst, vis, args, param_val, mod_dic, key_fluxderiv, key_compderiv, output_dic, sp_polcont_coeffs, intr_prof_dic, fixed_args)[source]#
Model function: joined disk-integrated profile per exposure
Defines the joined model for disk-integrated profile in a given exposure.
- Parameters:
TBD
- Returns:
TBD
- main_joined_DiffProf(rout_mode, data_dic, gen_dic, system_param, fit_prop_dic, theo_dic, plot_dic, coord_dic)[source]#
Joined differential profiles fits
Main routine to fit a given stellar surface property from planet-occulted regions with a joined model over instruments and visits.
- Parameters:
TBD
- Returns:
TBD
- joined_DiffProf(param, fixed_args)[source]#
Model function: joined differential profiles
Defines the joined model for differential profiles. This is done in three steps
We calculate all DI profiles of the star (fitted exposures + exposures that contributed to the master-out), and we scale them at the same value as after the Broadband flux Scaling module. A given out-of-transit profile corresponds to
\[\begin{split}F_\mathrm{DI}(t_\mathrm{out}) &= \sum_{k} dF_\mathrm{quiet}(k) + \sum_{i(t)} dF_\mathrm{quiet}(i(t)) + \sum_{j(t)} dF_\mathrm{active region}(j(t)) \\\end{split}\]Where stellar cells are either quiet or within active regions, and k indicates cells that are never within active regions. Other cells nature depend on time because the active region and planet are moving. A given in-transit profile corresponds to
\[\begin{split}F_\mathrm{DI}(t_\mathrm{in}) &= \sum_{k} dF_\mathrm{quiet}(k) + \sum_{i(t),nopl} dF_\mathrm{quiet}(i(t)) + \sum_{i(t),pl} dF_\mathrm{quiet}(i(t)) + \sum_{j(t),nopl} dF_\mathrm{active region}(j(t)) + \sum_{j(t),pl} dF_\mathrm{active region}(j(t)) \\ &= \sum_{k} dF_\mathrm{quiet}(k) + \sum_{i(t),nopl} dF_\mathrm{quiet}(i(t)) + \sum_{j(t),nopl} dF_\mathrm{active region}(j(t)) \\\end{split}\]Where the pl cells are occulted by the planet at time t and thus null.
We compute the master out, with same weights as those used in the corresponding module. The master out writes (neglecting weights) as
\[\begin{split}F_\mathrm{DI}(out) = \sum_{k} dF_\mathrm{quiet}(k) + <_{t_\mathrm{out}} \sum_{i(t)} dF_\mathrm{quiet}(i(t)) + \sum_{j(t)} dF_\mathrm{active region}(j(t)) > \\\end{split}\]We extract differential profiles as \(F_\mathrm{diff} = F_\mathrm{out} - F_\mathrm{sc}\), corresponding to
\[\begin{split}F_\mathrm{diff}(t_\mathrm{out}) &= <_{t_\mathrm{out}} \sum_{i(t)} dF_\mathrm{quiet}(i(t)) + \sum_{j(t)} dF_\mathrm{active region}(j(t)) > - \sum_{i(t)} dF_\mathrm{quiet}(i(t)) + \sum_{j(t)} dF_\mathrm{active region}(j(t)) \\ F_\mathrm{diff}(t_\mathrm{in}) &= <_{t_\mathrm{out}} \sum_{i(t)} dF_\mathrm{quiet}(i(t)) + \sum_{j(t)} dF_\mathrm{active region}(j(t)) - \sum_{i(t),nopl} dF_\mathrm{quiet}(i(t)) - \sum_{j(t),nopl} dF_\mathrm{active region}(j(t))>\end{split}\]If the active region is fixed, then
\[\begin{split}F_\mathrm{diff}(t_\mathrm{out}) &= 0 \\ F_\mathrm{diff}(t_\mathrm{in}) &= \sum_{i(t)} dF_\mathrm{quiet}(i(t)) + \sum_{j} dF_\mathrm{active region}(j) - \sum_{i(t),nopl} dF_\mathrm{quiet}(i(t)) - \sum_{j,nopl} dF_\mathrm{active region}(j) &= \sum_{i(t),pl} dF_\mathrm{quiet}(i(t)) + \sum_{j,pl} dF_\mathrm{active region}(j)\end{split}\]In that case in-transit differential profiles can be fitted directly with the joined_IntrProf routine, calculating the quiet or active regions profiles of planet-occulted cells.
- Parameters:
TBD
- Returns:
TBD
- joined_DiffProf_exp(inst, vis, iexp, isub_exp, isub_ord, fixed_args, master_out_dic, master_out_tab_inst, st_prof_grids_exp, st_prof_grids, sp_line_model, base_DI_prof, cond_def_mast_vis, cond_undef_weights_mast, glob_weight_all, contrib_profs2mast)[source]#
Model function: joined differential profile per exposure
Processing exposures and preparing master
- Parameters:
TBD
- Returns:
TBD
- joined_DiffProf_corr(inst, vis, master_out_dic, fixed_args, mod_dic, param)[source]#
Model function: joined differential profile creation
Computed differential profile model
- Parameters:
TBD
- Returns:
TBD
- main_joined_IntrProp(rout_mode, fit_prop_dic, gen_dic, system_param, theo_dic, plot_dic, coord_dic, data_dic)[source]#
Joined intrinsic stellar property fits
Main routine to fit a given stellar surface property from planet-occulted regions with a joined model over instruments and visits.
- Parameters:
TBD
- Returns:
TBD
- FIT_joined_IntrProp(param, x_tab, args=None)[source]#
Fit function: joined local stellar property
Calls corresponding model function for optimization
- Parameters:
TBD
- Returns:
TBD
- joined_IntrProp(param, fixed_args)[source]#
Model function: joined stellar property
Defines the joined model for stellar properties
- Parameters:
TBD
- Returns:
TBD
- main_joined_IntrProf(rout_mode, data_dic, gen_dic, system_param, fit_prop_dic, theo_dic, plot_dic, coord_dic)[source]#
Joined stellar profile fits
Main routine to fit intrinsic stellar profiles from planet-occulted regions with a joined model over instruments and visits.
Profile description
We use analytical models, measured profiles, or theoretical models to describe the intrinsic profiles.
Positions of the profiles along the transit chord are linked by the stellar surface RV model.
Shapes of the profiles, when analytical, are linked across the transit chord by polynomial laws as a function of a chosen dimension. Polynomial coefficients can depend on the visit and their associated instrument, to account for possible variations in the line shape between visits
Stellar line profiles are defined before instrumental convolution, so that data from all instruments and visits can be fitted together
Beware that the intrinsic and disk-integrated profiles have the same continuum, but that it is not necessarily unity Thus the continuum of analytical and theoretical model profiles must be let free to vary
- Parameters:
TBD
- Returns:
TBD