ekaterinailin / AltaiPony

Find flares in Kepler and TESS light curves. Notebooks for quickstart inside.
https://altaipony.readthedocs.io
MIT License
26 stars 10 forks source link

Write a function that checks each flare if it suffers from saturation. #35

Closed ekaterinailin closed 6 years ago

ekaterinailin commented 6 years ago

What needs to be created or improved?

Some function may go back to the TPF, get the aperture mask and see if the pixels inside the mask are saturated (>10094 or so). Then add a column to FlareLightCurve.flares with the average saturation in the mask.

Can you provide an example?

def measure_saturation(self):
    TPF = from_archive(self.targetid)
    relevanttpf = TPF.flux(in='some aperture')
    for i, row in self.flares.iterrows(): 
        [...]

What is the goal / expected behaviour?

We want to control for saturation to isolate possible departures from a single power law.

ekaterinailin commented 6 years ago

Now use FlareLightcurve.get_saturation().