Question

我正在尝试使用FermiTools和FermiPy python软件包对一些Fermi LAT数据进行一些分析。最终目标是将我的数据仅限制在选定的目标上，并创建一个用于进一步分析的光曲线，但我还没有开始。

我已经根据以下链接设置了环境并安装了FermiPy：

https://github.com/fermi-lat/Fermitools-conda/wiki/Quickstart-Guide https://docs.gammapy.org/dev/install/index.html#quick-install-guide

我有三个详细描述LAT事件数据的PH Fits文件和一个SC spacecraft Fits文件。我已经这样创建了一个YAML配置文件，并将其放置在工作目录中：

data :
  evfile : L181204201833D8AF139E12_PH00.fits, 
L181204201833D8AF139E12_PH01.fits, L181204201833D8AF139E12_PH02.fits
  scfile : L181204201833D8AF139E12_SC00.fits

gtlike:
  edisp : True
  irfs : 'P8R2_SOURCE_V6'
  edisp_disable : ['isodiff','galdiff']

fileio:
  outdir : null

selection :
  emin : 100
  emax : 316227.76
  zmax    : 180
  evclass : 128
  evtype  : 3
  filter  : null
  ra: '21.9811292171612'
  dec: '-29.774410928364'

lightcurve :
  nbins : 10

model :
  src_roiwidth : 3.0
  galdiff  : 'gll_iem_v06.fits'
  isodiff  : 'iso_P8R2_SOURCE_V6_v06.txt'

我正在运行以下内容：

from fermipy.gtanalysis import GTAnalysis

gta = GTAnalysis('config.yaml')
gta.setup()

该脚本在安装命令过程中的某个地方出现了IndexError：列表超出范围问题。寻找回溯，似乎应该生成一些实时多维数据集适合文件，但事实并非如此。然后，它尝试在文件列表中调用该文件，但失败。

Traceback (most recent call last):

  File "<ipython-input-2-46cd43b49eae>", line 1, in <module>
    runfile('/home/jo/Desktop/Astrophysics_Proj/Proj_Files/FermiPy_Lightcurve.py', wdir='/home/jo/Desktop/Astrophysics_Proj/Proj_Files')

  File "/home/jo/anaconda2/envs/fermi/lib/python2.7/site-packages/spyder_kernels/customize/spydercustomize.py", line 704, in runfile
    execfile(filename, namespace)

  File "/home/jo/anaconda2/envs/fermi/lib/python2.7/site-packages/spyder_kernels/customize/spydercustomize.py", line 100, in execfile
    builtins.execfile(filename, *where)

  File "/home/jo/Desktop/Astrophysics_Proj/Proj_Files/FermiPy_Lightcurve.py", line 15, in <module>
    gta.setup(use_local_ltcube=True)

  File "/home/jo/anaconda2/envs/fermi/lib/python2.7/site-packages/fermipy/gtanalysis.py", line 1039, in setup
    c.setup(overwrite=overwrite)

  File "/home/jo/anaconda2/envs/fermi/lib/python2.7/site-packages/fermipy/gtanalysis.py", line 5004, in setup
    self._ltc = LTCube.create(self.files['ltcube'])

  File "/home/jo/anaconda2/envs/fermi/lib/python2.7/site-packages/fermipy/ltcube.py", line 192, in create
    ltc = cls.create_from_fits(files[0])

IndexError: list index out of range

如果有人对从这里出发有什么建议，那就太好了！如果您想复制我的问题，也许尝试下载一些更小，更简单的fits文件来代替我的文件（它们很大）。一位朋友能够使用以下FITS文件复制该问题：

https://fermi.gsfc.nasa.gov/cgi-bin/ssc/LAT/QueryResults.cgi?id=L181207223916D8AF139E12

您只需要更改config.yaml文件中的文件名，然后对这些行进行旋转即可。

编辑：可以在配置文件中添加一个“ ltcube：”部分，然后自己提供该文件，但是我专门尝试生成它，因为我没有。

FermiPy和FermiTools实时多维数据集生成错误

0 个答案: