Skip to content

Commit

Permalink
Merge pull request gammapy#5460 from aaguasca/info_fermi_catalog
Browse files Browse the repository at this point in the history 
Fix get index of a PLSuperExpCutoff spectral model (fermi catalog)
  • Loading branch information
@registerrier
registerrier authored Oct 17, 2024
2 parents 0b2c129 + 1545062 commit 2510115
Show file tree
Hide file tree
Showing 3 changed files with 101 additions and 9 deletions.
11 changes: 8 additions & 3 deletions gammapy/catalog/fermi.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -407,9 +407,14 @@ def _info_spectral_fit(self):
)

fmt = "{:<45s} : {:.3f} +- {:.3f}\n"
ss += fmt.format(
"Spectral index", d[tag + "_Index"], d["Unc_" + tag + "_Index"]
)
if f"{tag}_ExpfactorS" in d:
ss += fmt.format(
"Spectral index", d[tag + "_IndexS"], d["Unc_" + tag + "_IndexS"]
)
else:
ss += fmt.format(
"Spectral index", d[tag + "_Index"], d["Unc_" + tag + "_Index"]
)

fmt = "{:<45s} : {:.3} +- {:.3} {}\n"
ss += fmt.format(
Expand Down
70 changes: 70 additions & 0 deletions gammapy/catalog/tests/data/4fgl-dr4_J0534.5+2200.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,70 @@

*** Basic info ***

Catalog row index (zero-based) : 1391
Source name : 4FGL J0534.5+2200
Extended name :
Associations : PSR J0534+2200, Crab IC field, Crab pulsar, 3FGL J0534.5+2201, 3FHL J0534.5+2201, 2AGL J0534+2205, EGR J0534+2159
ASSOC_PROB_BAY : 1.000
ASSOC_PROB_LR : 1.000
Class1 : PSR
Class2 :
TeVCat flag : P

*** Other info ***

Significance (100 MeV - 1 TeV) : 25.186
Npred : 98594.1

Other flags : 784

*** Position info ***

RA : 83.637 deg
DEC : 22.015 deg
GLON : 184.559 deg
GLAT : -5.781 deg

Semimajor (68%) : 0.0045 deg
Semiminor (68%) : 0.0044 deg
Position angle (68%) : -71.52 deg
Semimajor (95%) : 0.0073 deg
Semiminor (95%) : 0.0072 deg
Position angle (95%) : -71.52 deg
ROI number : 270

*** Spectral info ***

Spectrum type : PLSuperExpCutoff
Detection significance (100 MeV - 1 TeV) : 25.186
Exponential factor : 0.2337 +- 0.0188
Super-exponential cutoff index : 0.5219 +- 0.0946
Significance curvature : 27.6
Pivot energy : 1343 MeV
Spectral index : 2.274 +- 0.008
Flux Density at pivot energy : 1.1e-10 +- 6.99e-13 cm-2 MeV-1 s-1
Integral flux (1 - 100 GeV) : 1.54e-07 +- 8.44e-10 cm-2 s-1
Energy flux (100 MeV - 100 GeV) : 1.46e-09 +- 2.16e-11 erg cm-2 s-1

*** Spectral points ***

e_min e_max flux flux_errn flux_errp e2dnde e2dnde_errn e2dnde_errp is_ul flux_ul e2dnde_ul sqrt_ts
MeV MeV 1 / (s cm2) 1 / (s cm2) 1 / (s cm2) erg / (s cm2) erg / (s cm2) erg / (s cm2) 1 / (s cm2) erg / (s cm2)
---------- ----------- ----------- ----------- ----------- ------------- ------------- ------------- ----- ----------- ------------- -------
50.000 100.000 2.365e-06 2.901e-07 2.916e-07 3.800e-10 4.662e-11 4.686e-11 False nan nan 8.309
100.000 300.000 1.595e-06 3.088e-08 3.088e-08 3.840e-10 7.436e-12 7.436e-12 False nan nan 60.667
300.000 1000.000 5.546e-07 5.472e-09 5.472e-09 3.760e-10 3.710e-12 3.710e-12 False nan nan 138.346
1000.000 3000.000 1.255e-07 8.804e-10 8.804e-10 2.918e-10 2.047e-12 2.047e-12 False nan nan 155.573
3000.000 10000.000 2.464e-08 8.789e-10 8.789e-10 1.552e-10 5.536e-12 5.536e-12 False nan nan 9.059
10000.000 30000.000 2.149e-09 4.120e-10 2.330e-10 4.498e-11 8.622e-12 4.877e-12 False nan nan 3.116
30000.000 100000.000 1.542e-16 nan 4.113e-12 8.085e-18 nan 2.157e-13 True 8.227e-12 4.314e-13 0.000
100000.000 1000000.000 6.125e-18 nan 3.875e-12 5.675e-19 nan 3.590e-13 True 7.750e-12 7.180e-13 0.000
*** Lightcurve info ***

Lightcurve measured in the energy band: 100 MeV - 100 GeV

Variability index : 156.729
Significance peak (100 MeV - 100 GeV) : 154.788
Integral flux peak (100 MeV - 100 GeV) : 2.54e-06 +- 2.55e-08 cm^-2 s^-1
Time peak : 3.18e+08 s (Mission elapsed time)
Peak interval : 3.65e+02 day
29 changes: 23 additions & 6 deletions gammapy/catalog/tests/test_fermi.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -113,6 +113,17 @@
),
]

SOURCES_4FGL_DR4 = [
dict(
idx=1391,
name="4FGL J0534.5+2200",
str_ref_file="data/4fgl-dr4_J0534.5+2200.txt",
spec_type=SuperExpCutoffPowerLaw4FGLDR3SpectralModel,
dnde=u.Quantity(1.1048e-07, "cm-2 s-1 GeV-1"),
dnde_err=u.Quantity(6.9934e-10, "cm-2 s-1 GeV-1"),
)
]

SOURCES_3FGL = [
dict(
idx=0,
Expand Down Expand Up @@ -186,9 +197,10 @@


@requires_data()
def test_4FGL_DR4():
@pytest.mark.parametrize("ref", SOURCES_4FGL_DR4, ids=lambda _: _["name"])
def test_4FGL_DR4(ref):
cat = SourceCatalog4FGL("$GAMMAPY_DATA/catalogs/fermi/gll_psc_v32.fit.gz")
source = cat["4FGL J0534.5+2200"]
source = cat[ref["name"]]
model = source.spectral_model()
fp = source.flux_points
not_ul = ~fp.is_ul.data.squeeze()
Expand All @@ -199,6 +211,11 @@ def test_4FGL_DR4():
models = cat.to_models()
assert len(models) == len(cat.table)

actual = str(cat[ref["idx"]])
with open(get_pkg_data_filename(ref["str_ref_file"])) as fh:
expected = fh.read()
assert actual == modify_unit_order_astropy_5_3(expected)


@requires_data()
class TestFermi4FGLObject:
Expand Down Expand Up @@ -363,9 +380,9 @@ def test_lightcurve_dr1(self):
assert_allclose(table["flux_errn"][0], 4.437058e-8, rtol=1e-3)

def test_lightcurve_dr4(self):
dr2 = SourceCatalog4FGL("$GAMMAPY_DATA/catalogs/fermi/gll_psc_v32.fit.gz")
source_dr2 = dr2[self.source_name]
table = source_dr2.lightcurve(interval="1-year").to_table(
dr4 = SourceCatalog4FGL("$GAMMAPY_DATA/catalogs/fermi/gll_psc_v32.fit.gz")
source_dr4 = dr4[self.source_name]
table = source_dr4.lightcurve(interval="1-year").to_table(
format="lightcurve", sed_type="flux"
)

Expand All @@ -379,7 +396,7 @@ def test_lightcurve_dr4(self):
assert_allclose(table["flux_errn"][0], 2.298336e-08, rtol=1e-3)

with pytest.raises(ValueError):
source_dr2.lightcurve(interval="2-month")
source_dr4.lightcurve(interval="2-month")


@requires_data()
Expand Down

0 comments on commit 2510115

Please sign in to comment.