Mexican Hat Power Spectra (MHPS)

Based in work of Patricia Arevalo Power Spectra.

Installing MHPS

Run:

pip install -e .

Using MHPS

import pandas as pd
import mhps

df = pd.read_csv(".../ZTF18aaajeru.csv", index_col=0)
mag = df.magpsf_corr
magerr = df.sigmapsf_corr
time = df.index
t1 = 100
t2 = 10

#Other parameters
# dt = 3.0
# mag0 = 19.0
# epsilon = 1.0
ratio,low,high,non_zero, PN_flag = mhps.statistics(mag.values,
                                                   magerr.values,
                                                   time.values,
                                                   t1,t2)

Known issues

Some times the algorithm gets an ZeroDivisionException this will return all values as numpy NaN (np.nan).