Select Git revision
createEleutheriaDagSubmitFile.py
-
Grant Meadors authored
Made changes so as to show graphs of MICH, PRC and their product with transfer functions in the frequency domain so that we can better anticipate what pre-processing measures may prove necessary.
Grant Meadors authoredMade changes so as to show graphs of MICH, PRC and their product with transfer functions in the frequency domain so that we can better anticipate what pre-processing measures may prove necessary.
test_aperture.py 709 B
from pykat.utilities.optics.gaussian_beams import gauss_param
from pykat import finesse
from pykat.commands import xaxis
import pylab as pl
import numpy as np
import math
code = """
l l1 1 0 0 n1
s s1 10 1 n1 n2
m m1 1 0 0 n2 n3
pd refl n2
xaxis m1 r_ap lin 0.1e-3 2e-3 10
"""
kat = finesse.kat()
kat.parseCommands(code)
maxtem = np.arange(0, 4)
kat.m1.n2.q = gauss_param(w0=1e-3, z=0)
kat.verbose = False
for tem in maxtem:
print "Calculating maxtem ", tem, "..."
kat.maxtem = tem
r = kat.run()
pl.plot(r.x/1e-3, r.y, label="maxtem={0}".format(tem))
pl.ylabel("Reflected Power [W]")
pl.xlabel("Mirror aperture [mm]")
pl.legend()
pl.show()