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)

Joined stellar fits

Wrap-up function to call joint fits of stellar properties and profiles

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)

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)

Fit function: joined global stellar property

Calls corresponding model function for optimization

Parameters:

TBD

Returns:

TBD

mod_DIProp(param, args, inst, vis, n_coord)

Model function: global stellar property

Defines the model for global stellar property over a visit.

Parameters:

TBD

Returns:

TBD

joined_DIProp(param, args)

Model function: joined global stellar property

Defines the joined model for global stellar properties

Parameters:

TBD

Returns:

TBD

main_joined_IntrProp(rout_mode, fit_prop_dic, gen_dic, system_param, theo_dic, plot_dic, coord_dic, data_dic)

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)

Fit function: joined local stellar property

Calls corresponding model function for optimization

Parameters:

TBD

Returns:

TBD

joined_IntrProp(param, args)

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)

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

FIT_joined_IntrProf(param, x_tab, args=None)

Fit function: joined intrinsic stellar profiles

Calls corresponding model function for optimization

Parameters:

TBD

Returns:

TBD

joined_IntrProf(param, fixed_args)

Model function: joined intrinsic stellar profiles

Defines the joined model for intrinsic stellar profiles.

Parameters:

TBD

Returns:

TBD

main_joined_DiffProf(rout_mode, data_dic, gen_dic, system_param, fit_prop_dic, theo_dic, plot_dic, coord_dic)

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

FIT_joined_DiffProf(param, x_tab, args=None)

Fit function: joined differential stellar profiles

Calls corresponding model function for optimization

Parameters:

TBD

Returns:

TBD

joined_DiffProf(param, fixed_args)

Model function: joined differential profiles

Defines the joined model for differential profiles. This is done in three steps

1. 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.

2. 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}\]

3. 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

calc_chi_Prof(mod_dic, args)

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 intrinsic 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

outputs_Prof(inst, vis, coeff_line_dic, mod_prop_dic, args, param)

Fit function: outputs

Defines outputs for profile fit functions

Parameters:

TBD

Returns:

TBD