.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "examples_gallery/plot_indirect_detection_GCE.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_examples_gallery_plot_indirect_detection_GCE.py>`
        to download the full example code

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_examples_gallery_plot_indirect_detection_GCE.py:


Indirect detection calculation for GCE
======================================

This script calculates the gamma-ray flux for the Galactic center excess.

.. GENERATED FROM PYTHON SOURCE LINES 7-13

.. code-block:: default


    import matplotlib.pyplot as plt
    import numpy as np
    from singletscalar_dm import *
    from scipy.interpolate import interp1d








.. GENERATED FROM PYTHON SOURCE LINES 14-16

In order to calculate the relic density it is possible to use the function `interpolate_Omega`, which uses the code DRAKE near the Higgs resonance and MicrOMEGAs elsewhere.
For examplte, let's compute the relic density, expressed as :math:`\Omega h^2` for the QCDA model.

.. GENERATED FROM PYTHON SOURCE LINES 16-37

.. code-block:: default


    table_best  = np.loadtxt(import_data_file('GCE_BestFitModel_flux_Inner40x40_masked_disk.dat'))
    table_up  = np.loadtxt(import_data_file('GCE_band_up.dat'))
    table_down  = np.loadtxt(import_data_file('GCE_band_down.dat'))
    func_min = interp1d(table_down[:,0],table_down[:,1])
    func_max = interp1d(table_up[:,0],table_up[:,1])
    energy_C_vec = np.logspace(np.log10(0.315),np.log10(34.0),100)
    energy_Cerr_vec = table_best[:,0]
    flux_C_min = func_min(energy_C_vec)
    flux_C_max = func_max(energy_C_vec)
    flux_Cerr_min = func_min(energy_Cerr_vec)
    flux_Cerr_max = func_max(energy_Cerr_vec)
    flux_av = table_best[:,1]

    DMmass=62.485
    lambda_hs=2e-4
    energy_vec = np.logspace(-1.,2.,100)
    flux_DM = np.zeros(len(energy_vec))
    for i in range(len(energy_vec)):
        flux_DM[i] = flux_DM_prompt(energy_vec[i],DMmass,lambda_hs,False)








.. GENERATED FROM PYTHON SOURCE LINES 38-39

Now we create the plot.

.. GENERATED FROM PYTHON SOURCE LINES 39-57

.. code-block:: default


    fig = plt.figure(figsize=(8,6))
    plt.plot(energy_vec,flux_DM*np.power(energy_vec,2.), lw=5.0, ls='--', color='blue', label=r'SHP model')
    plt.fill_between(energy_C_vec,flux_C_min,flux_C_max,alpha=0.5,color="grey")
    plt.errorbar(table_best[:,0], table_best[:,1], yerr=[table_best[:,1]-table_best[:,2],table_best[:,3]-table_best[:,1]], fmt='*', color='black',label=r'Cholis+2022')
    plt.ylabel(r'$E^2 \frac{dN}{dE}$ [GeV/cm$^2$/s/sr]', fontsize=18)
    plt.xlabel(r'$E$ [GeV]', fontsize=18)
    plt.axis([0.1,100,8.e-8,3e-6])
    plt.xticks(fontsize=18)
    plt.yticks(fontsize=18)
    plt.tick_params('both', length=7, width=2, which='major')
    plt.tick_params('both', length=5, width=2, which='minor')
    plt.grid(True)
    plt.yscale('log')
    plt.xscale('log')
    plt.legend(loc=2,prop={'size':16},numpoints=1, scatterpoints=1, ncol=1)
    fig.tight_layout(pad=0.5)
    plt.show()



.. image-sg:: /examples_gallery/images/sphx_glr_plot_indirect_detection_GCE_001.png
   :alt: plot indirect detection GCE
   :srcset: /examples_gallery/images/sphx_glr_plot_indirect_detection_GCE_001.png
   :class: sphx-glr-single-img






.. rst-class:: sphx-glr-timing

   **Total running time of the script:** ( 0 minutes  33.187 seconds)


.. _sphx_glr_download_examples_gallery_plot_indirect_detection_GCE.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example




    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: plot_indirect_detection_GCE.py <plot_indirect_detection_GCE.py>`

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: plot_indirect_detection_GCE.ipynb <plot_indirect_detection_GCE.ipynb>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_