run_ls

tessilator.lc_analysis.run_ls(cln, n_sca=10, p_min_thresh=0.05, p_max_thresh=100.0, samples_per_peak=10)

Run Lomb-Scargle periodogram and return a dictionary of results.

Parameters:

clndict

A dictionary containing the lightcurve data. The keys must include | “time0” -> The time coordinate relative to the first data point | “nflux” -> The detrended, cleaned, normalised flux values | “enflux” -> The uncertainty for each value of nflux | “lc_part” -> An running index describing the various contiguous sections

p_min_threshfloat, optional, default=0.05

The minimum period (in days) to be calculated.

p_max_threshfloat, optional, default=100.

The maximum period (in days) to be calculated.

samples_per_peakint, optional, default=10

The number of samples to measure in each periodogram peak.

Returns:

LS_dictdict

A dictionary of parameters calculated from the periodogram analysis. These are: | “median_MAD_nLC” : The median and median absolute deviation of the normalised lightcurve. | “jump_flag” : A flag determining if the lightcurve has sharp jumps in flux. | “period” : A list of period values from the periodogram analysis. | “power” : A list of power values from the periodogram analysis. | “period_best” : The period corrseponding to the highest power output. | “power_best” : The highest power output. | “time” : The time coordinate corresponding to the normalised lightcurve. | “y_fit_LS” : The best fit sinusoidal function. | “AIC_sine” : The Aikaike Information Criterion value of the best-fit sinusoid | “AIC_line” : The Aikaike Information Criterion value of the best-fit linear function. | “FAPs” : The power output for the false alarm probability values of 0.1, 1 and 10% | “Gauss_fit_peak_parameters” : Parameters for the Gaussian fit to the highest power peak | “Gauss_fit_peak_y_values” : The corresponding y-values for the Gaussian fit | “period_around_peak” : The period values covered by the Gaussian fit | “power_around_peak” : The power values across the period range covered by the Gaussian fit | “period_not_peak” : The period values not covered by the Gaussian fit | “power_not_peak” : The power values across the period range not covered by the Gaussian fit | “period_second” : The period of the second highest peak. | “power_second” : The power of the second highest peak. | “phase_fit_x” : The time co-ordinates from the best-fit sinusoid to the phase-folded lightcurve. | “phase_fit_y” : The normalised flux co-ordinates from the best-fit sinusoid to the phase-folded lightcurve. | “phase_x” : The time co-ordinates from the phase-folded lightcurve. | “phase_y” : The normalised flux co-ordinates from the phase-folded lightcurve. | “phase_chisq” : The chi-square fit between the phase-folded lightcurve and the sinusoidal fit. | “phase_col” : The cycle number for each data point. | “pops_vals” : The best-fit parameters from the sinusoidal fit to the phase-folded lightcurve. | “pops_cov” : The corresponding co-variance matrix from the “pops_val” parameters. | “phase_scatter” : The typical scatter in flux around the best-fit. | “frac_phase_outliers” : The fraction of data points that are more than 3 median absolute deviation values from the best-fit. | “Ndata” : The number of data points used in the periodogram analysis.