#!/usr/bin/env python
# coding: utf-8

# In[1]:


import finesse
finesse.configure(plotting=True)

kat = finesse.Model()
kat.parse(
    """
    # Laser with power 1W
    laser L0 P=1
    s s0 L0.p1 eom1.p1

    # Phase modulator at 40kHz, with a modulation index of 0.05
    mod eom1 f=40k midx=.05 order=5 mod_type=pm

    # Three amplitude detectors, looking at the first three positive modulation
    # sidebands
    ad bessel1 eom1.p2.o f=40k
    ad bessel2 eom1.p2.o f=80k
    ad bessel3 eom1.p2.o f=120k

    # Sweep the modulation index of the modulator from 0 to 10
    xaxis(eom1.midx, lin, 0, 10, 300)
    """
)


# In[2]:


out = kat.run()
out.plot();