Goal Figure! - Githubissues

jradavenport commented 8 years ago

Here's the end-goal figure I want you to shoot for:

The "Flare Frequency Diagram" for your star with 2 sets of data plotted: the (newest) results from Appaloosa, and your by-hand FBEYE results.

The X-axis is the flare energy

The Y-axis is the cumulative flare rate (in number per day)

To make this figure in Python you'll need to do these steps (and probably more)

[x] Read the FBEYE results (.fbeye) file in to Python using np.loadtxt or similar
[x] Also read the original light curve in to Python the same way
[x] Also read the Appaloosa results (.flare) in to Python the same way
[x] Compute the total duration of the data in the light curve
[ ] Compute the cumulative flare energy distribution
[ ] Convert this cumulative distribution in to a rate, dividing by the duration

There is much more complication to this if we have long and short cadence data to combine...

cornetj2 commented 8 years ago

Progress has been made, my only problem is computing the cumulative flare energy distribution. Also, there is an issue with current rate calculation.

jradavenport commented 8 years ago

What's going on with the current rate calculation?

You get the (reverse) cumulative distribution by adding up the number of flares from smallest to biggest. Let's talk this week about this!

cornetj2 / FlareData