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

.. only:: html

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

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

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

.. _sphx_glr_auto_examples_plot_trojans.py:


Plotting Jupiter Trojans
========================

.. GENERATED FROM PYTHON SOURCE LINES 5-30

.. code-block:: Python


    import kete
    import numpy as np
    import matplotlib.pyplot as plt

    # Load orbit data and select just the jupiter trojans
    orbs = kete.mpc.fetch_known_orbit_data()
    subset = orbs[kete.population.jup_trojan(orbs.peri_dist, orbs.ecc)]

    # Construct the states and propagate them to a common epoch
    states = kete.mpc.table_to_states(subset)
    states = kete.propagate_n_body(states, states[0].jd)

    # Where is jupiter?
    jupiter = kete.spice.get_state("Jupiter", states[0].jd)
    # found it!

    # Compute the positions, and relative longitudinal distance from jupiter
    positions = np.transpose([[s.pos.x, s.pos.y] for s in states])

    longitudes = np.array([(s.pos.lon - jupiter.pos.lon + 180) % 360 - 180 for s in states])
    leading = longitudes > 0
    trailing = longitudes < 0

    # Plot the results







.. GENERATED FROM PYTHON SOURCE LINES 31-33

Histogram of Trojan counts
--------------------------

.. GENERATED FROM PYTHON SOURCE LINES 33-47

.. code-block:: Python

    plt.figure(figsize=(6, 4), dpi=150)
    bins = np.linspace(-180, 180, 120)
    plt.hist(longitudes[leading], bins=bins, histtype="step", color="C0")
    plt.hist(longitudes[trailing], bins=bins, histtype="step", color="C2")
    plt.title("Jupiter Trojans")
    plt.xlabel("Ecliptic Lon away from Jupiter (Deg)")
    plt.ylabel("Counts")

    print(
        "Fraction of trojans which are in the leading camp: "
        f"{sum(leading) / len(longitudes):0.2%}"
    )
    plt.show()




.. image-sg:: /auto_examples/images/sphx_glr_plot_trojans_001.png
   :alt: Jupiter Trojans
   :srcset: /auto_examples/images/sphx_glr_plot_trojans_001.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    Fraction of trojans which are in the leading camp: 63.44%




.. GENERATED FROM PYTHON SOURCE LINES 48-50

Top down view of the solar system
---------------------------------

.. GENERATED FROM PYTHON SOURCE LINES 50-63

.. code-block:: Python

    plt.figure(dpi=150)
    plt.scatter(*positions[:, leading], s=0.2, c="C0")
    plt.scatter(*positions[:, trailing], s=0.2, c="C2")
    plt.scatter(jupiter.pos.x, jupiter.pos.y, c="C4", label="Jupiter")
    plt.scatter(0, 0, c="Orange")
    plt.scatter([], [], c="C0", label="Leading")
    plt.scatter([], [], c="C2", label="Trailing")
    plt.gca().axis("equal")
    plt.legend()
    plt.xlabel("Distance (AU)")
    plt.ylabel("Distance (AU)")
    plt.title("Jupiter Trojans")
    plt.show()



.. image-sg:: /auto_examples/images/sphx_glr_plot_trojans_002.png
   :alt: Jupiter Trojans
   :srcset: /auto_examples/images/sphx_glr_plot_trojans_002.png
   :class: sphx-glr-single-img






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

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


.. _sphx_glr_download_auto_examples_plot_trojans.py:

.. only:: html

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

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

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

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

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

    .. container:: sphx-glr-download sphx-glr-download-zip

      :download:`Download zipped: plot_trojans.zip <plot_trojans.zip>`


.. only:: html

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

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