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

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 JTable无法向新创建的列添加值   java如何调整JEditorPane中编辑区域的大小?   Java通过反射确定未知数组中的数组大小   java Intellij Idea有时无法按其预期的方式构建应用程序   java Swing GUI带有IntelliJ错误“contentPane不能设置为null”从终端编译时   如何将这些通用方法调用从C#转换为Java   在null上找不到java属性或字段“index”   从Java HashMap获取整数值时是否需要调用intValue()方法?   java Android谷歌地图获取相机中的图像块   unix无法捕获JAVA中“who m”命令的输出   java,同时将邮件发送到“收件人”标题“我”中的多个收件人   在java中向链表添加未知数量的节点   无法为Heroku上的discord bot设置java端口   java使用Apache HttpClient进行选项请求   与元素类型“ApplicationName”关联的属性“Application Version”需要java Open quote   Android Studio Java中的两个变量求和