test_plot.py 851 Bytes
Newer Older
Daniel Brown's avatar
Daniel Brown committed
1
2
3
4
5
from pykat import finesse
from pykat.detectors import *
from pykat.components import *
from pykat.commands import *
from pykat.structs import *
6

Daniel Brown's avatar
Daniel Brown committed
7
8
9
import numpy as np
import pylab as pl

10
11
12
13
14
15
16
17
18
19
code = """
l l1 1 0 0 n1
s s1 10 1 n1 n2
m m1 0.5 0.5 0 n2 n3
s s2 10 1 n3 n4
m m2 0.5 0.5 0 n4 n5
s s3 10 1 n5 n6
"""

kat = finesse.kat(kat_code=code)
Daniel Brown's avatar
Daniel Brown committed
20

21
kat.add(cavity('cav1','m1','n3','m2','n4'))
Daniel Brown's avatar
Daniel Brown committed
22

23
24
25
kat.add(photodiode('pd_cav','n4'))
kat.add(photodiode('pd_ref','n2'))
kat.add(photodiode('pd_trs','n5'))
Daniel Brown's avatar
Daniel Brown committed
26

27
kat.add(xaxis(Scale.linear, [0, 360], kat.m2, kat.m2.phi, 100))
Daniel Brown's avatar
Daniel Brown committed
28
29
30
31
32
33
34
35
36
37

kat.m1.Rcx = -1000.0
kat.m1.Rcy = -1000.0
kat.m2.Rcx =  1000.0
kat.m2.Rcy =  1000.0

kat.maxtem = 0

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

38
39
40
41
42
43
pl.figure()
pl.plot(run.x,run.y)
pl.xlabel(run.xlabel)
pl.ylabel("Intensity [W]")
pl.legend(run.ylabels)
pl.show()