FERMI LAT data analysis using fermipy

[NithyaSathyan]

Data Analysis Question

How do I...

Not sure if a certain command exists in fermipy or how to use it? Looking for examples or best practices? Ask here!

To reproduce
Please let us know which data file, instrument response function etc. you are using, and include your model(s) and analysis script.

[NithyaSathyan]

I am doing FERMI LAT data analysis using fermipy. Recently i had to reinstall fermipy. However, after the installation, I encounter the following error message after running the analysis script.

File "/home/dell/Documents/1H1914-194/xx.py", line 3, in
from gammapy.irf.io import IRF_DL3_AXES_SPECIFICATION
File "/home/dell/miniconda3/envs/fermipy/lib/python3.9/site-packages/gammapy/irf/init.py", line 6, in
from .background import Background2D, Background3D
File "/home/dell/miniconda3/envs/fermipy/lib/python3.9/site-packages/gammapy/irf/background.py", line 10, in
from .core import IRF
File "/home/dell/miniconda3/envs/fermipy/lib/python3.9/site-packages/gammapy/irf/core.py", line 17, in
from .io import IRF_DL3_HDU_SPECIFICATION, IRF_MAP_HDU_SPECIFICATION, gadf_is_pointlike
File "/home/dell/miniconda3/envs/fermipy/lib/python3.9/site-packages/gammapy/irf/io.py", line 4, in
from gammapy.data.hdu_index_table import HDUIndexTable
File "/home/dell/miniconda3/envs/fermipy/lib/python3.9/site-packages/gammapy/data/init.py", line 3, in
from .data_store import DataStore
File "/home/dell/miniconda3/envs/fermipy/lib/python3.9/site-packages/gammapy/data/data_store.py", line 15, in
from .observations import Observation, ObservationChecker, Observations
File "/home/dell/miniconda3/envs/fermipy/lib/python3.9/site-packages/gammapy/data/observations.py", line 22, in
from .event_list import EventList, EventListChecker
File "/home/dell/miniconda3/envs/fermipy/lib/python3.9/site-packages/gammapy/data/event_list.py", line 8, in
from astropy.coordinates.angle_utilities import angular_separation
ModuleNotFoundError: No module named 'astropy.coordinates.angle_utilities'

The analysis script is given below.

from fermipy.gtanalysis import GTAnalysis
gta = GTAnalysis('config.yaml',logging={'verbosity': 3})
gta.setup()
gta.free_sources(distance=3.0,pars='norm')
gta.free_source('galdiff')
gta.free_source('isodiff')
gta.free_sources(minmax_ts=[10,None],pars='norm')
gta.free_source('4FGL J1917.7-1921')
gta.fit()
gta.write_roi('fit_model')
import numpy as np
sed=gta.sed('4FGL J1917.7-1921', make_plots=True)

Please help me solve this issue. Thanks in advance

[nyctophile-aayu]

@NithyaSathyan

Have you tried installing the astropy using conda after activating the conda environment for fermipy? Also, try installing astropy with pip. Hopefully that should work.

-Ayush

[NithyaSathyan]

I reinstalled fermipy using the code below: mamba create --name fermipy -c conda-forge -c fermi python=3.9 "fermitools>=2.2.0" healpy numpy=1.20 astropy=5.1 gammapy=0.20 mamba activate fermipy pip install fermipy Now, everything works fine. Thanks for the reply.

…

On Fri, 12 Jan, 2024, 12:04 pm nyctophile-aayu, ***@***.***> wrote: @NithyaSathyan <https://github.com/NithyaSathyan> Have you tried installing the astropy using conda after activating the conda environment for fermipy? Also, try installing astropy with pip. Hopefully that should work. -Ayush — Reply to this email directly, view it on GitHub <#561 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/BE6FP4IGK5ROSXZOUJGOYKDYODKPVAVCNFSM6AAAAABBKD7PESVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTQOBYGUYDSNBSGQ> . You are receiving this because you were mentioned.Message ID: ***@***.***>

[NithyaSathyan]

Hii.. I can see data points with TS (Test Statistic) values less than 10 after running the analysis using an analysis script with the line gta.free_sources(minmax_ts=[10,None],pars='norm'). Why is it so? On Fri, Jan 12, 2024 at 12:16 PM NITHYA SATHYAN ***@***.***> wrote:

…

I reinstalled fermipy using the code below: mamba create --name fermipy -c conda-forge -c fermi python=3.9 "fermitools>=2.2.0" healpy numpy=1.20 astropy=5.1 gammapy=0.20 mamba activate fermipy pip install fermipy Now, everything works fine. Thanks for the reply. On Fri, 12 Jan, 2024, 12:04 pm nyctophile-aayu, ***@***.***> wrote: > @NithyaSathyan <https://github.com/NithyaSathyan> > > Have you tried installing the astropy using conda after activating the > conda environment for fermipy? Also, try installing astropy with pip. > Hopefully that should work. > > -Ayush > > — > Reply to this email directly, view it on GitHub > <#561 (comment)>, > or unsubscribe > <https://github.com/notifications/unsubscribe-auth/BE6FP4IGK5ROSXZOUJGOYKDYODKPVAVCNFSM6AAAAABBKD7PESVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTQOBYGUYDSNBSGQ> > . > You are receiving this because you were mentioned.Message ID: > ***@***.***> >

[omodei]

The minmax_ts option in gta.free_sources is used to free the sources in the ROI with TS>10, not to set a minimum threshold for datapoint/upper limits.