search

fermiPy / fermipy

Fermi-LAT Python Analysis Framework
http://fermipy.readthedocs.org/
BSD 3-Clause "New" or "Revised" License
51 stars 53 forks source link

Using FL8Y as input catalog #172

Closed astro-fermi closed 6 years ago

astro-fermi commented 6 years ago

I am trying to use 8 year list (FL8Y) as input source catalog. I add the following in the config file

catalogs :

However, I get the following error File "stacking_pre_processing_3.py", line 99, in main gta = GTAnalysis('%s.yaml' % srcname,logging={'verbosity' : 3}) File "/u/gl/vaidehi/.local/lib/python2.7/site-packages/fermipy-0.16.0-py2.7.egg/fermipy/gtanalysis.py", line 331, in init coordsys=configs[0]['binning']['coordsys']) File "/u/gl/vaidehi/.local/lib/python2.7/site-packages/fermipy-0.16.0-py2.7.egg/fermipy/roi_model.py", line 1669, in create config, kwargs) File "/u/gl/vaidehi/.local/lib/python2.7/site-packages/fermipy-0.16.0-py2.7.egg/fermipy/roi_model.py", line 1703, in create_from_source coordsys=coordsys, kwargs) File "/u/gl/vaidehi/.local/lib/python2.7/site-packages/fermipy-0.16.0-py2.7.egg/fermipy/roi_model.py", line 1689, in create_from_position roi = cls(config, skydir=skydir, coordsys=coordsys, kwargs) File "/u/gl/vaidehi/.local/lib/python2.7/site-packages/fermipy-0.16.0-py2.7.egg/fermipy/roi_model.py", line 1300, in init self.load(coordsys=coordsys, srcname=srcname) File "/u/gl/vaidehi/.local/lib/python2.7/site-packages/fermipy-0.16.0-py2.7.egg/fermipy/roi_model.py", line 1624, in load srcname=srcname) File "/u/gl/vaidehi/.local/lib/python2.7/site-packages/fermipy-0.16.0-py2.7.egg/fermipy/roi_model.py", line 1899, in load_fits_catalog self.load_existing_catalog(cat, kwargs) File "/u/gl/vaidehi/.local/lib/python2.7/site-packages/fermipy-0.16.0-py2.7.egg/fermipy/roi_model.py", line 1963, in load_existing_catalog catalog_dict) File "/u/gl/vaidehi/.local/lib/python2.7/site-packages/fermipy-0.16.0-py2.7.egg/fermipy/roi_model.py", line 241, in spectral_pars_from_catalog pars = get_function_defaults(cat['SpectrumType']) File "/u/gl/vaidehi/.local/lib/python2.7/site-packages/fermipy-0.16.0-py2.7.egg/fermipy/model_utils.py", line 37, in get_function_defaults fn_spec = get_function_spec(name) File "/u/gl/vaidehi/.local/lib/python2.7/site-packages/fermipy-0.16.0-py2.7.egg/fermipy/model_utils.py", line 63, in get_function_spec raise Exception('Invalid Function Name: %s' % name) Exception: Invalid Function Name: PowerLaw

Please help. Thanks.

woodmd commented 6 years ago

Could you point me to the location of the catalog file you're using?

astro-fermi commented 6 years ago

Here it is

/u/gl/vaidehi/gll_psc_8year_v3.fit

astro-fermi commented 6 years ago

I also tested using 'gll_pscP305uw8011_v2.xml' file available at the confluence page of FL8Y. However, I got another error ....... lots of lines ........ Exception: Invalid Function Name: PLSuperExpCutoff2

woodmd commented 6 years ago

@eacharles has done some recent work to add support for FL8Y that was just merged in #174. I haven't tested it myself but you might give it a try. Note that this is only available in master at the moment so you'll need to install from source to have this functionality.

woodmd commented 6 years ago

FL8Y should now be fully supported as of #179 . Note that FL8Y must be used in conjunction with the v18 extended archive. Here is a sample model configuration:

model:
    catalogs : ['FL8Y']
    extdir : '/path/to/ExtendedArchive_v18'

This functionality is only available in master at present. To install the development version of fermipy you can follow these instructions.