计算外行星信号传递的天体物理假阳性概率

VESPA的Python项目详细描述


vespa
==
…图片::https://travis-ci.com/timothydmoton/vespa.svg?branch=master
:target:https://travis ci.com/timothydmoton/vespa


使用概率算法验证系外行星信号---计算传输信号的假阳性概率

有关用法和更多信息,请参阅文档<;http://vespa.rtfd.org>;`.

安装您可以从pypi获得最新发布的版本:



或者克隆存储库::

cd vespa
python setup.py install[--user]

特权。

取决于典型的科学软件包(例如“numpy”、“scipy”、“pandas”)、
以及“等时线”<;http://github.com/timothydmoton/isochrones>;`,以及(在代码的几个角落中)我的另一个软件包“simpledit<;http://github.com/timothydmoton/simpledit>;`。所有依赖项*都应该*在安装时得到解决,尽管这只是在anaconda python发行版下测试过的,该发行版已经将所有的科学内容组织得很好了。

为了获得最佳结果,还建议安装“multinest”和“pymultinest”。没有这个,``emcee``将用于恒星建模,但鉴于恒星模型拟合的多模态特性,``mulitnest`结果更可信一些。

制作一个文本文件,将透射光度分为三列:``t_from嫒midtranit```[days]、``flux```[relative,where out-transit is normalized to unity]和``flux嫒err``。文件不应有标题行(没有标题);可以是空格或逗号分隔的(将由“np.loadtxt”接收)。


2。制作一个“star.ini”文件,其中包含目标恒星的观测属性(光度和/或光谱,任何可用的):

目标星的大小
如果提供不确定性,将用于拟合starmodel
j=9.763,0.03
h=9.135,0.03
k=8.899,0.02
kepler=12.473


3。制作一个包含以下信息的“fpp.ini”文件::


name=k2oi;任何内容
ra=11:30:14.510也可以是十进制形式
dec=+07:35:18.21


period=32.988天
rprs=0.0534\rp/rstar best estimate
photfile=lc\k2oi.csv渡越光度计
BR/>〔约束〕Br/> Max rad=12×孔径半径[ARCSEC]
SeCdjHy= 1E-4最大潜在潜蚀深度
BR/> 4。从命令行运行以下命令(在具有“star.ini”和“fpp.ini”的同一文件夹中):

$calcfpp

等等,应该需要20-30分钟。如果您想测试一个较短的
版本,可以尝试使用“calcfpp-n 1000”(默认值为20000)。第一次运行时,大约有一半的时间是在进行恒星建模,因此仍然需要几分钟。





attribution
--


如果使用此代码,请同时引用论文和代码。


论文引用::

@article{2012apj…761….6m,
author={{{morton},t.~d.},
title=“{一个有效的外行星过境候选者自动验证程序”,
journal={\apj},
archiveprefix=“arxiv”,
eprint={1206.1568},
primaryclass=“astro-ph.ep”,
keywords={行星系统,恒星:统计},
year=2012,
month=dec,
volume=761,
eid={6},
pages={6},
doi={10.1088/0004-637X/761/1/6/6},
adurl={http://adabs.harvarvard.edu/abs/2012APJ/761…..6m},
adsnote={由SAO/NASA天体物理学数据系统提供



code:





@misc{2015MISC{20152015MISC{201520152015MISC{2015201520152015ascl.soft03011m,
作者={{morton},t.~d.},
title={vespa:false positive probability calculator},
howpublished={astrophysics源代码库},
year=2015,
month=mar,
archiveprefix=“ascl”,
eprint={1503.011},
adsurl={http://adsabs.harvard.edu/abs/2015ascl.soft03011m},
adsnote={由sao/nasa天体物理数据系统提供}
}

欢迎加入QQ群-->: 979659372 Python中文网_新手群

推荐PyPI第三方库


热门话题
java置换相关算法   在Java中读取/打开文本文件   java为什么这段代码不适用于CodeEval?   java如何将日历转换为JSON   从txt文件中读取字符串并将其存储到java中的字符数组中   字符编码Java ResourceBundles umlauts搞砸了   java为什么hashmap会根据总大小而不是填充的存储桶调整大小   java如何将Excel单元格中的数字字符串读取为字符串(而不是数字)?   java Guava的LocalCache无法使用,为什么?   java有没有办法强制JVM在单个处理器或内核上运行   java Eclipse不安装软件   将字节转换为java字符串(可能是汉字)   Java正则表达式:提取函数名   JavaTestNG:如何从多个类中指定测试方法顺序?