test_put.py 647 Bytes
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
from pykat import finesse
from pykat.detectors import *
from pykat.components import *
from pykat.commands import *
from pykat.structs import *

import numpy as np
import pylab as pl

code = """
l l1 1 0 0 n1
s s1 10 1 n1 n2
mod eom 10 0.1 1 am 0 n2 n3
"""

kat = finesse.kat()

kat.parseCommands(code)

20 21
kat.add(pd('pdp',1,'n3'))
kat.add(pd('pdm',1,'n3'))
22

23 24 25 26
kat.add(xaxis("lin", [0, 1000], kat.eom, "f", 100))

kat.pdp.f1.put(kat.xaxis.x)
kat.pdm.f1.put(kat.xaxis.mx)
27 28 29 30

out = kat.run(printout=0, printerr=0)

pl.figure()
31
pl.plot(out.x, out["pdp"], out.x, out["pdm"])
32 33 34 35
pl.xlabel(out.xlabel)
pl.ylabel("Intensity [W]")
pl.legend(out.ylabels)
pl.show()