.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "auto_examples/ESR/esr_example.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_ESR_esr_example.py>`
        to download the full example code

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

.. _sphx_glr_auto_examples_ESR_esr_example.py:


Bruker cw ESR Data
==================

Load cw Bruker ESR data, both 1D and 2D.

Check out the
:ref:`Simple UV <sphx_glr_auto_examples_UV_Cary_simple.py>`
example to understand how
pySpecData locates the file here.

.. GENERATED FROM PYTHON SOURCE LINES 12-16

.. code-block:: Python


    import matplotlib.pyplot as plt
    import pyspecdata as psd








.. GENERATED FROM PYTHON SOURCE LINES 17-20

Load some 1D ESR data with harmonic + phase info.
The data is initial organized into two dimensions -- `harmonic` and $B_0$.


.. GENERATED FROM PYTHON SOURCE LINES 20-33

.. code-block:: Python


    d = psd.find_file(
        "S175R1a.*DHPC.*today.*200304.*",
        exp_type="francklab_esr/Sam",
        zenodo="21287247",
    )
    print(d.shape)
    print(
        "here, we see the harmonic axis contains both harmonic and phase info",
        repr(d.getaxis("harmonic")),
    )
    d.chunk_auto("harmonic", "phase")





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

 .. code-block:: none

    [(2048, '$B_0$'), (10, 'harmonic')]
    here, we see the harmonic axis contains both harmonic and phase info array([(1,  0), (1, 90), (2,  0), (2, 90), (3,  0), (3, 90), (4,  0),
           (4, 90), (5,  0), (5, 90)],
          dtype=[('harmonic', '<i8'), ('phase', '<i8')])

    array([[[-2.94546828e-01, -3.68288604e-02],
            [-3.08371764e-03,  6.13238812e-03],
            [-4.49149462e-03,  7.18220516e-03],
            [ 6.73636716e-04, -2.11555613e-03],
            [ 6.88192341e-03,  5.08947960e-03]],

           [[-2.94215753e-01, -3.57182955e-02],
            [-3.51026174e-03,  4.29212270e-03],
            [-4.31634015e-03,  8.39053228e-03],
            [ 2.66558776e-04, -2.56135657e-03],
            [ 7.00552768e-03,  5.10833628e-03]],

           [[-2.94953691e-01, -3.34736934e-02],
            [-2.56422486e-03,  1.12643511e-03],
            [-2.53831179e-03,  8.33304986e-03],
            [ 6.55038547e-04, -2.93743236e-03],
            [ 5.11506815e-03,  6.28893360e-03]],

           ...,

           [[-2.80658800e-01, -3.28446011e-02],
            [ 4.01519075e-03, -2.80185439e-03],
            [ 1.54833082e-03,  4.65094659e-03],
            [ 5.32599721e-03, -8.04821331e-03],
            [-5.21515998e-03, -1.34590515e-03]],

           [[-2.80400928e-01, -3.25151160e-02],
            [ 2.18324464e-03, -9.15192415e-04],
            [-1.20783223e-04,  6.36937529e-03],
            [ 5.44008378e-03, -9.96095938e-03],
            [-2.31595131e-03, -1.64052881e-03]],

           [[-2.80445226e-01, -2.97726370e-02],
            [ 1.72233837e-03, -9.90763420e-05],
            [-2.16840206e-05,  6.24455227e-03],
            [ 5.35087800e-03, -9.71951497e-03],
            [ 5.82405591e-04, -3.06429791e-03]]],
          shape=(2048, 5, 2), dtype='>f8')
    	dimlabels=['$B_0$', 'harmonic', 'phase']
    	axes={`$B_0$':array([3436.        , 3436.07324219, 3436.14648438, ..., 3585.78027362,
           3585.85351581, 3585.926758  ], shape=(2048,))
    			+/-None,
    		`harmonic':array([1, 2, 3, 4, 5])
    			+/-None,
    		`phase':array([ 0, 90])
    			+/-None}




.. GENERATED FROM PYTHON SOURCE LINES 34-36

`chunk_auto` breaks the `harmonic` dimensions since it was labeled with an
axis that had 2 fields.

.. GENERATED FROM PYTHON SOURCE LINES 36-44

.. code-block:: Python


    print(d.shape)

    plt.figure(1)
    psd.plot(d["phase", 0], alpha=0.5)
    psd.plot(d["phase", 1], ":", alpha=0.5)
    plt.title("1D Data with Multiple Harmonics")




.. image-sg:: /auto_examples/ESR/images/sphx_glr_esr_example_001.png
   :alt: 1D Data with Multiple Harmonics
   :srcset: /auto_examples/ESR/images/sphx_glr_esr_example_001.png, /auto_examples/ESR/images/sphx_glr_esr_example_001_2_00x.png 2.00x
   :class: sphx-glr-single-img


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

 .. code-block:: none

    [(2048, '$B_0$'), (5, 'harmonic'), (2, 'phase')]

    Text(0.5, 1.0, '1D Data with Multiple Harmonics')



.. GENERATED FROM PYTHON SOURCE LINES 45-46

Next, let's load some power-saturation data

.. GENERATED FROM PYTHON SOURCE LINES 46-62

.. code-block:: Python


    d = psd.find_file(
        ".*Power.*Sat_200303.*",
        exp_type="francklab_esr/Sam",
        zenodo="21287247",
    )
    d.chunk_auto("harmonic", "phase")
    plt.figure(2)
    psd.image(
        d["harmonic", 0]["phase", 0]
        .C.set_axis("Microwave Power", "#")
        .set_units("Microwave Power", "scan #")
    )
    plt.title("2D Power Saturation")
    plt.gca().set_aspect("auto")
    plt.show()



.. image-sg:: /auto_examples/ESR/images/sphx_glr_esr_example_002.png
   :alt: 2D Power Saturation
   :srcset: /auto_examples/ESR/images/sphx_glr_esr_example_002.png, /auto_examples/ESR/images/sphx_glr_esr_example_002_2_00x.png 2.00x
   :class: sphx-glr-single-img






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

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


.. _sphx_glr_download_auto_examples_ESR_esr_example.py:

.. only:: html

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

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

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

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

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


.. only:: html

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

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