compute_false_alarm_threshold

exovetter.modshift.compute_false_alarm_threshold(period_days, duration_hrs)

Compute the stat, significance needed to invalidate the null hypothesis

An event should be considered statistically significant if its peak in the convolved lightcurves is greater than the value computed by this function.

Note that this number is computed on a per-TCE basis. If you are looking at many TCEs you will need a stronger threshold. (For example, if you use this function to conclude that there is a less than 0.1% chance a given event is a false alarm due to Gaussian noise, you expect to see one such false alarm in 1,000 TCEs. See Coughlin et al. for the formula to ensure less than 1 false alarm over many TCEs.

Parameters:

period_daysfloat

Orbital period

duration_hrsfloat

Duration of transit in hours.

Returns:

fafloat

TODO What exactly is returned. Is this the 1 sigma false alarm threshold?