Created
May 24, 2019 12:59
-
-
Save debaice/681f1d357da26af86ec8cf5ca3eb4ad2 to your computer and use it in GitHub Desktop.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
(gammapy-dev) Debanjans-MacBook-Pro:CTA debanjanbose$ ipython | |
Python 3.7.0 | packaged by conda-forge | (default, Nov 12 2018, 12:34:36) | |
Type 'copyright', 'credits' or 'license' for more information | |
IPython 7.5.0 -- An enhanced Interactive Python. Type '?' for help. | |
In [1]: from gammapy.data import DataStore | |
In [2]: DataStore.from_events_files ("M82_events_18.fits") | |
--------------------------------------------------------------------------- | |
TypeError Traceback (most recent call last) | |
<ipython-input-2-262488ac4bbb> in <module> | |
----> 1 DataStore.from_events_files ("M82_events_18.fits") | |
~/gammapy-dev/gammapy/gammapy/data/data_store.py in from_events_files(cls, paths) | |
188 data_store.obs_table.write("obs-index.fits.gz") | |
189 """ | |
--> 190 return DataStoreMaker(paths).run() | |
191 | |
192 @staticmethod | |
~/gammapy-dev/gammapy/gammapy/data/data_store.py in __init__(self, paths) | |
401 def __init__(self, paths): | |
402 if isinstance(paths, (str, Path)): | |
--> 403 raise TypeError("Need list of paths, not a single string or Path object.") | |
404 | |
405 self.paths = [make_path(path) for path in paths] | |
TypeError: Need list of paths, not a single string or Path object. | |
In [3]: DataStore.from_events_files (["M82_events_18.fits"]) | |
Out[3]: <gammapy.data.data_store.DataStore at 0x11307c208> | |
In [4]: d=DataStore.from_events_files (["M82_events_18.fits"]) | |
In [5]: d.obs_table | |
Out[5]: | |
<ObservationTable length=1> | |
OBS_ID RA_PNT DEC_PNT GLON_PNT GLAT_PNT ZEN_PNT ... OBJECT TELESCOP CALDB IRF EVENTS_FILENAME EVENT_COUNT | |
deg deg deg deg deg ... | |
int64 float64 float64 float64 float64 float64 ... str1 str3 str13 str13 str18 int64 | |
------ ------- ------- ------------------ ------------------ ------- ... ------ -------- ------------- ------------- ------------------ ----------- | |
1 148.97 69.679 141.40980733845007 40.567806671343114 0.0 ... CTA NOT AVAILABLE NOT AVAILABLE M82_events_18.fits 13664312 | |
In [6]: d.hdu_table | |
Out[6]: | |
<HDUIndexTable length=6> | |
HDU_TYPE HDU_CLASS HDU_NAME OBS_ID FILE_DIR FILE_NAME | |
str6 str10 str21 int64 str47 str18 | |
-------- ---------- --------------------- ------ ----------------------------------------------- ------------------ | |
events events EVENTS 1 . M82_events_18.fits | |
gti gti GTI 1 . M82_events_18.fits | |
aeff aeff_2d EFFECTIVE AREA 1 $CALDB/data/CTA/NOT AVAILABLE/bcf/NOT AVAILABLE irf_file.fits | |
edisp edisp_2d ENERGY DISPERSION 1 $CALDB/data/CTA/NOT AVAILABLE/bcf/NOT AVAILABLE irf_file.fits | |
psf psf_3gauss POINT SPREAD FUNCTION 1 $CALDB/data/CTA/NOT AVAILABLE/bcf/NOT AVAILABLE irf_file.fits | |
bkg bkg_3d BACKGROUND 1 $CALDB/data/CTA/NOT AVAILABLE/bcf/NOT AVAILABLE irf_file.fits | |
In [7]: d.obs | |
Out[7]: <bound method DataStore.obs of <gammapy.data.data_store.DataStore object at 0x11307c7b8>> | |
In [8]: d.obs(1) | |
Out[8]: <gammapy.data.observations.DataStoreObservation at 0x11346c8d0> | |
In [9]: obs=d.obs(1) | |
In [10]: obs.events | |
Out[10]: <gammapy.data.event_list.EventList at 0x1130db630> | |
In [11]: obs.psf | |
--------------------------------------------------------------------------- | |
FileNotFoundError Traceback (most recent call last) | |
<ipython-input-11-02d00e5a608c> in <module> | |
----> 1 obs.psf | |
~/gammapy-dev/gammapy/gammapy/data/observations.py in psf(self) | |
128 def psf(self): | |
129 """Load point spread function object.""" | |
--> 130 return self.load(hdu_type="psf") | |
131 | |
132 @property | |
~/gammapy-dev/gammapy/gammapy/data/observations.py in load(self, hdu_type, hdu_class) | |
93 """ | |
94 location = self.location(hdu_type=hdu_type, hdu_class=hdu_class) | |
---> 95 return location.load() | |
96 | |
97 @property | |
~/gammapy-dev/gammapy/gammapy/data/hdu_index_table.py in load(self) | |
98 from ..irf import EnergyDependentMultiGaussPSF | |
99 | |
--> 100 return EnergyDependentMultiGaussPSF.read(filename, hdu=hdu) | |
101 elif hdu_class == "psf_king": | |
102 from ..irf import PSFKing | |
~/gammapy-dev/gammapy/gammapy/irf/psf_gauss.py in read(cls, filename, hdu) | |
122 """ | |
123 filename = make_path(filename) | |
--> 124 with fits.open(str(filename), memmap=False) as hdulist: | |
125 psf = cls.from_fits(hdulist[hdu]) | |
126 | |
~/anaconda3/envs/gammapy-dev/lib/python3.7/site-packages/astropy/io/fits/hdu/hdulist.py in fitsopen(name, mode, memmap, save_backup, cache, lazy_load_hdus, **kwargs) | |
149 | |
150 return HDUList.fromfile(name, mode, memmap, save_backup, cache, | |
--> 151 lazy_load_hdus, **kwargs) | |
152 | |
153 | |
~/anaconda3/envs/gammapy-dev/lib/python3.7/site-packages/astropy/io/fits/hdu/hdulist.py in fromfile(cls, fileobj, mode, memmap, save_backup, cache, lazy_load_hdus, **kwargs) | |
388 return cls._readfrom(fileobj=fileobj, mode=mode, memmap=memmap, | |
389 save_backup=save_backup, cache=cache, | |
--> 390 lazy_load_hdus=lazy_load_hdus, **kwargs) | |
391 | |
392 @classmethod | |
~/anaconda3/envs/gammapy-dev/lib/python3.7/site-packages/astropy/io/fits/hdu/hdulist.py in _readfrom(cls, fileobj, data, mode, memmap, save_backup, cache, lazy_load_hdus, **kwargs) | |
1037 if not isinstance(fileobj, _File): | |
1038 # instantiate a FITS file object (ffo) | |
-> 1039 fileobj = _File(fileobj, mode=mode, memmap=memmap, cache=cache) | |
1040 # The Astropy mode is determined by the _File initializer if the | |
1041 # supplied mode was None | |
~/anaconda3/envs/gammapy-dev/lib/python3.7/site-packages/astropy/utils/decorators.py in wrapper(*args, **kwargs) | |
519 kwargs[new_name[i]] = value | |
520 | |
--> 521 return function(*args, **kwargs) | |
522 | |
523 return wrapper | |
~/anaconda3/envs/gammapy-dev/lib/python3.7/site-packages/astropy/io/fits/file.py in __init__(self, fileobj, mode, memmap, overwrite, cache) | |
176 self._open_fileobj(fileobj, mode, overwrite) | |
177 elif isinstance(fileobj, str): | |
--> 178 self._open_filename(fileobj, mode, overwrite) | |
179 else: | |
180 self._open_filelike(fileobj, mode, overwrite) | |
~/anaconda3/envs/gammapy-dev/lib/python3.7/site-packages/astropy/io/fits/file.py in _open_filename(self, filename, mode, overwrite) | |
553 | |
554 if not self._try_read_compressed(self.name, magic, mode, ext=ext): | |
--> 555 self._file = fileobj_open(self.name, IO_FITS_MODES[mode]) | |
556 self.close_on_error = True | |
557 | |
~/anaconda3/envs/gammapy-dev/lib/python3.7/site-packages/astropy/io/fits/util.py in fileobj_open(filename, mode) | |
386 """ | |
387 | |
--> 388 return open(filename, mode, buffering=0) | |
389 | |
390 | |
FileNotFoundError: [Errno 2] No such file or directory: '/Users/debanjanbose/fermi/ScienceTools-v11r5p3-fssc-20180124-x86_64-apple-darwin16.7.0/x86_64-apple-darwin16.7.0/refdata/fermi/caldb/CALDB/data/CTA/NOT AVAILABLE/bcf/NOT AVAILABLE/irf_file.fits' | |
In [12]: |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment