This code's aim is to performing parameter estimation for stellar flares in data from the Kepler satellite.
The Data class in this example will attempt to read in the TIME and
PDCSAP_FLUX light curve data from a Kepler
FITS file. It will set the epoch of the time stamps (read in as days) to be the first time stamp value and
remove any data containing NaNs or infinities. The light curve data will have the median value removed from it.
The aim is to be able to simultaneously fit a variable number of background/noise components (in this case
sinusoids, single-bin impulses and step changes in the underlying background level) and also a variable
number of flares. The maximum number of the various allowed model components (and some of their ranges)
can be set through a custom JSON configuration file, rather than having it hardwired in the code.
I have started to implement the flare burst model from the magnetron
code (by Huppenkothen, Brewer, Frean, Hogg & Murray) as described in Huppenkothen et al,
arXiv:1501.05251. This currently is unlikely to work though!
The DNest4 algorithm is described in Brewer et al, arXiv:0912.2380 and the documentation can be found within the code repository. The additional RJObject code is described in Brewer, arXiv:1411.3921.
- DNest4
- CFITSIO - on Debian/Ubuntu install with e.g.
sudo apt-get install libcfitsio3-dev - CCfits - on Debian/Ubuntu install with e.g.
sudo apt-get install libccfits-dev - Boost - in particulate the boost IOStreams and System libraries, on Debian/Ubuntu install with e.g.
sudo apt-get install libboost-dev - GSL - on Debian/Ubuntu (16.04) install with e.g.
sudo apt-get install libgsl-dev
Compilation of the code makes use of scons. To compile the code just run scons in the base directory
© 2013 Brendon J. Brewer, 2015 Matthew Pitkin
Licence: MIT License for software.