VESPA 0.6

Calculate astrophysical false positive probabilities for transiting exoplanet signals

VESPA
======
.. image:: https://travis-ci.com/timothydmorton/VESPA.svg?branch=master
:target: https://travis-ci.com/timothydmorton/VESPA

Validation of Exoplanet Signals using a Probabilistic Algorithm--- calculating false positive probabilities for transit signals

For usage and more info, `check out the documentation <http://vespa.rtfd.org>`_.

Installation
------------

To install, you can get the most recently released version from PyPI::

pip install vespa [--user]

Or you can clone the repository::

git clone https://github.com/timothydmorton/vespa.git
cd vespa
python setup.py install [--user]

The ``--user`` argument may be necessary if you don't have root privileges.

Depends on typical scientific packages (e.g. `numpy`, `scipy`, `pandas`),
as well as `isochrones <http://github.com/timothydmorton/isochrones>`_, and (in several corners of the code), another package of mine called `simpledist <http://github.com/timothydmorton/simpledist>`_. All dependencies *should* get resolved upon install, though this has only been tested under the anaconda Python distribution, which has all the scientific stuff already well-organized.

For best results, it is also recommended to have ``MultiNest`` and ``pymultinest`` installed. Without this, ``emcee`` will be used for stellar modeling, but the ``MulitNest`` results are a bit more trustworthy given the often multi-modal nature of stellar model fitting.

Basic Usage
-----------

The simplest way to run an FPP calculation straight out of the box is
as follows.

1. Make a text file containing the transit photometry in three columns: ``t_from_midtransit`` [days], ``flux`` [relative, where out-of-transit is normalized to unity], and ``flux_err``. The file should not have a header row (no titles); and can be either whitespace or comma-delimited (will be ingested by ``np.loadtxt``).

2. Make a ``star.ini`` file that contains the observed properties of the target star (photometric and/or spectroscopic, whatever is available)::

#provide spectroscopic properties if available
#Teff = 3503, 80 #value, uncertainty
#feh = 0.09, 0.09
#logg = 4.89, 0.1

#observed magnitudes of target star
# If uncertainty provided, will be used to fit StarModel
J = 9.763, 0.03
H = 9.135, 0.03
K = 8.899, 0.02
Kepler = 12.473

3. Make a ``fpp.ini`` file containing the following information::

name = k2oi #anything
ra = 11:30:14.510 #can be decimal form too
dec = +07:35:18.21

period = 32.988 #days
rprs = 0.0534 #Rp/Rstar best estimate
photfile = lc_k2oi.csv #contains transit photometry

[constraints]
maxrad = 12 # aperture radius [arcsec]
secthresh = 1e-4 # Maximum allowed depth of potential secondary eclipse

4. Run the following from the command line (from within the same folder that has ``star.ini`` and ``fpp.ini``)::

$ calcfpp

Or, if you put the files in a folder called ``mycandidate``, then you can run ``calcfpp mycandidate``::

This will run the calculation for you, creating result files, diagnostic plots, etc.
It should take 20-30 minutes. If you want to do a shorter
version to test, you can try ``calcfpp -n 1000`` (the default is 20000). The first
time you run it though, about half the time is doing the stellar modeling, so it will still
take a few minutes.


Attribution
-----------

If you use this code, please cite both the paper and the code.

Paper citation::

@ARTICLE{2012ApJ...761....6M,
author = {{Morton}, T.~D.},
title = "{An Efficient Automated Validation Procedure for Exoplanet Transit Candidates}",
journal = {\apj},
archivePrefix = "arXiv",
eprint = {1206.1568},
primaryClass = "astro-ph.EP",
keywords = {planetary systems, stars: statistics },
year = 2012,
month = dec,
volume = 761,
eid = {6},
pages = {6},
doi = {10.1088/0004-637X/761/1/6},
adsurl = {http://adsabs.harvard.edu/abs/2012ApJ...761....6M},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}

code::

@MISC{2015ascl.soft03011M,
author = {{Morton}, T.~D.},
title = "{VESPA: False positive probabilities calculator}",
howpublished = {Astrophysics Source Code Library},
year = 2015,
month = mar,
archivePrefix = "ascl",
eprint = {1503.011},
adsurl = {http://adsabs.harvard.edu/abs/2015ascl.soft03011M},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}