计算外行星信号传递的天体物理假阳性概率
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天体物理数据系统提供}
}
==
…图片::https://travis-ci.com/timothydmoton/vespa.svg?branch=master
:target:https://travis ci.com/timothydmoton/vespa
使用概率算法验证系外行星信号---计算传输信号的假阳性概率
有关用法和更多信息,请参阅文档<;http://vespa.rtfd.org>;`.
安装您可以从pypi获得最新发布的版本:
或者克隆存储库::
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天体物理数据系统提供}
}