.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "user-guide/model-gallery/spectral/plot_logparabola2.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code. or to run this example in your browser via Binder .. rst-class:: sphx-glr-example-title .. _sphx_glr_user-guide_model-gallery_spectral_plot_logparabola2.py: .. _logparabola2-spectral-model: Log parabola 2 spectral model =========================== This model parametrises a log parabola spectrum, where the energy scale of the exponente, :math:`E_s`, and the reference energy, :math:`E_0`, can be different. It is defined by the following equation: .. math:: \phi(E) = \phi_0 \left( \frac{E}{E_0} \right) ^ { - \alpha - \beta \log{ \left( \frac{E}{E_s} \right) } } Note that :math:`log` refers to the natural logarithm. If you have parametrization based on :math:`log_{10}` you can use the :func:`~gammapy.modeling.models.LogParabola2SpectralModel.from_log10` method. .. GENERATED FROM PYTHON SOURCE LINES 23-26 Example plot ------------ Here is an example plot of the model: .. GENERATED FROM PYTHON SOURCE LINES 26-42 .. code-block:: Python from astropy import units as u import matplotlib.pyplot as plt from gammapy.modeling.models import LogParabola2SpectralModel, Models, SkyModel energy_bounds = [0.1, 100] * u.TeV model = LogParabola2SpectralModel( alpha=2.3, amplitude="1e-12 cm-2 s-1 TeV-1", reference=10 * u.TeV, beta=0.5, escale=1 * u.TeV, ) model.plot(energy_bounds) plt.grid(which="both") .. image-sg:: /user-guide/model-gallery/spectral/images/sphx_glr_plot_logparabola2_001.png :alt: plot logparabola2 :srcset: /user-guide/model-gallery/spectral/images/sphx_glr_plot_logparabola2_001.png :class: sphx-glr-single-img .. GENERATED FROM PYTHON SOURCE LINES 43-46 YAML representation ------------------- Here is an example YAML file using the model: .. GENERATED FROM PYTHON SOURCE LINES 46-51 .. code-block:: Python model = SkyModel(spectral_model=model, name="log-parabola2-model") models = Models([model]) print(models.to_yaml()) .. rst-class:: sphx-glr-script-out .. code-block:: none components: - name: log-parabola2-model type: SkyModel spectral: type: LogParabola2SpectralModel parameters: - name: amplitude value: 1.0e-12 unit: TeV-1 s-1 cm-2 - name: reference value: 10.0 unit: TeV - name: alpha value: 2.3 - name: beta value: 0.5 - name: escale value: 1.0 unit: TeV metadata: creator: Gammapy 2.0.dev2512+g0733504a8 date: '2025-11-30T05:22:18.018928' origin: null .. _sphx_glr_download_user-guide_model-gallery_spectral_plot_logparabola2.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: binder-badge .. image:: images/binder_badge_logo.svg :target: https://mybinder.org/v2/gh/gammapy/gammapy-webpage/main?urlpath=lab/tree/notebooks/dev/user-guide/model-gallery/spectral/plot_logparabola2.ipynb :alt: Launch binder :width: 150 px .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: plot_logparabola2.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: plot_logparabola2.py ` .. container:: sphx-glr-download sphx-glr-download-zip :download:`Download zipped: plot_logparabola2.zip ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_