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

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第三方库


热门话题
ByteArrayOutputStream的java解码属性   java S3 SDK在上载时更新单个对象,而不是创建新文件   java hibernate:无法从eclipse连接到DB   java如何在强制转换JComboBox之前检查其类型?   http从Java中的GETPOST请求方法捕获URI、资源名称,如开发人员工具中所示   java在Spring@Bean方法中返回接口的局限性   Java中的Web服务和客户端(使用Eclipse Apache Axis 2自底向上服务)某些代码会引发异常   java spring安全+rest不起作用   java将LinkedList添加到包含LinkedList的LinkedList并更改添加的LinkedList   java是否临时删除对象的属性?   java使用AnimatedGifEncoder类创建的gif图像的部分帧是不透明的   java如何高效地处理maven3时间戳快照?   java向集合对象添加另一项   java如何将动态参数传递给jquery函数   java使用libGdx桌面端口作为Android GLES20的仿真器