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

.. only:: html

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

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

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

.. _sphx_glr_gallery_validation_internal_phase_function_detector.py:


Goniometric Coupling vs S1 S2 Comparison
========================================

.. GENERATED FROM PYTHON SOURCE LINES 6-96



.. image-sg:: /gallery/validation/internal/images/sphx_glr_phase_function_detector_001.png
   :alt: Polarization 90 degree
   :srcset: /gallery/validation/internal/images/sphx_glr_phase_function_detector_001.png
   :class: sphx-glr-single-img





.. code-block:: Python


    # Standard library imports
    import numpy as np
    import matplotlib.pyplot as plt

    # PyMieSim imports
    from PyMieSim.experiment.detector import Photodiode
    from PyMieSim.experiment.scatterer import Sphere as ExperimentSphere
    from PyMieSim.experiment.source import Gaussian as ExperimentGaussian
    from PyMieSim.experiment import Setup
    from PyMieSim.units import degree, watt, AU, RIU, micrometer

    from PyMieSim.single.scatterer import Sphere as SingleSphere
    from PyMieSim.single.source import Gaussian as SingleGaussian
    from MPSPlots.styles import mps

    # Setup parameters
    scatterer_diameter = 0.3 * micrometer  # Diameter of the scatterer in meters
    scatterer_index = 1.4 * RIU  # Refractive index of the scatterer
    source_wavelength = 1.2 * micrometer  # Wavelength of the source in meters

    # Experiment source and scatterer setup
    source = ExperimentGaussian(
        wavelength=1.2 * micrometer,
        polarization=[0, 90] * degree,
        optical_power=1 * watt,
        NA=0.2 * AU
    )

    scatterer = ExperimentSphere(
        diameter=scatterer_diameter,
        property=scatterer_index,
        medium_property=1.0 * RIU,
        source=source
    )

    # Detector setup
    detector = Photodiode(
        NA=[0.1] * AU,
        phi_offset=np.linspace(-180, 180, 100) * degree,
        gamma_offset=0.0 * degree,
        sampling=1000 * AU,
        polarization_filter=None
    )

    # Configure experiment
    experiment = Setup(scatterer=scatterer, source=source, detector=detector)

    # Gather data
    dataframe = experiment.get('coupling', drop_unique_level=True)
    # dataframe.index = /= 180 / np.pi
    dataframe['coupling'] /= dataframe['coupling'].max()  # Normalize data

    # Single scatterer simulation for S1 and S2
    single_source = SingleGaussian(
        wavelength=source_wavelength,
        polarization=90 * degree,
        optical_power=1 * watt,
        NA=0.2 * AU
    )

    single_scatterer = SingleSphere(
        diameter=scatterer_diameter,
        source=single_source,
        property=scatterer_index,
        medium_property=1.0 * RIU
    )

    s1s2 = single_scatterer.get_s1s2()
    phi, s1, s2 = s1s2.phi, np.abs(s1s2.S1)**2, np.abs(s1s2.S2)**2
    s1 /= s1.max()  # Normalize S1 data
    s2 /= s2.max()  # Normalize S2 data


    with plt.style.context(mps):
        figure, ax0 = plt.subplots(1, 1, subplot_kw=dict(projection='polar'))


    df = dataframe.unstack('source:polarization').pint.dequantify().reset_index().pint.quantify()
    df['detector:phi_offset'] /= (180 / np.pi)

    df.plot(x='detector:phi_offset', y='coupling', ax=ax0, linewidth=3, title='Polarization 90 degree')

    ax0.plot(np.deg2rad(phi), s1, color='white', linestyle='--', linewidth=1, label='Computed S1')

    ax0.plot(np.deg2rad(phi), s2, color='cyan', linestyle='--', linewidth=1, label='Computed S2')

    ax0.grid()

    plt.show()


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

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


.. _sphx_glr_download_gallery_validation_internal_phase_function_detector.py:

.. only:: html

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

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

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

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

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

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

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


.. only:: html

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

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