Commit 06176db680ce33e5f7f04b4dfaf4db6dbe7886f8
1 parent
63cb462c
Adjust isotropic spectrum according to results from Georg Weidenspointer's PhD thesis
Showing
1 changed file
with
10 additions
and
9 deletions
modules/comscripts/comobsmodel.py
| ... | ... | @@ -2,7 +2,7 @@ |
| 2 | 2 | # ========================================================================== |
| 3 | 3 | # Generate model for binned COMPTEL observations |
| 4 | 4 | # |
| 5 | -# Copyright (C) 2019-2021 Juergen Knoedlseder | |
| 5 | +# Copyright (C) 2019-2022 Juergen Knoedlseder | |
| 6 | 6 | # |
| 7 | 7 | # This program is free software: you can redistribute it and/or modify |
| 8 | 8 | # it under the terms of the GNU General Public License as published by |
| ... | ... | @@ -451,15 +451,16 @@ class comobsmodel(ctools.csobservation): |
| 451 | 451 | spatial = gammalib.GModelSpatialDiffuseConst(gammalib.fourpi) |
| 452 | 452 | |
| 453 | 453 | # Set spectral model using initial scaling factors that correspond |
| 454 | - # to Georg Weidenspointers extragalactic spectrum. Values integrated | |
| 455 | - # over the standard energy bands communicated by Werner Collmar and | |
| 456 | - # divided by width of standard energy bands. | |
| 457 | - # See his e-mail from 14 July 2021. | |
| 454 | + # to the spectrum | |
| 455 | + # | |
| 456 | + # I(E) = 1.12e-4 (E/5MeV)^-2.2 ph/cm2/s/sr/MeV | |
| 457 | + # | |
| 458 | + # determined by Georg Weidenspointer in this PhD thesis. | |
| 458 | 459 | spectral = gammalib.GModelSpectralNodes() |
| 459 | - spectral.append(gammalib.GEnergy(0.87, 'MeV'), 1.327e-3/0.25) | |
| 460 | - spectral.append(gammalib.GEnergy(1.73, 'MeV'), 2.358e-3/2.0) | |
| 461 | - spectral.append(gammalib.GEnergy(5.48, 'MeV'), 0.658e-3/7.0) | |
| 462 | - spectral.append(gammalib.GEnergy(17.32, 'MeV'), 0.148e-3/20.0) | |
| 460 | + spectral.append(gammalib.GEnergy(0.87, 'MeV'), 5.2481e-3) | |
| 461 | + spectral.append(gammalib.GEnergy(1.73, 'MeV'), 1.1567e-3) | |
| 462 | + spectral.append(gammalib.GEnergy(5.48, 'MeV'), 9.15e-5) | |
| 463 | + spectral.append(gammalib.GEnergy(17.32, 'MeV'), 7.2e-6) | |
| 463 | 464 | |
| 464 | 465 | # Set fix |
| 465 | 466 | fix = (self['iso'].string() == 'CONSTFIX') | ... | ... |