Skip to content
Snippets Groups Projects
Select Git revision
  • a76f88c16370e3db396376d1b684d24aa7f00c89
  • master default
2 results

test_aperture.py

Blame
  • Forked from finesse / pykat
    492 commits behind the upstream repository.
    test_aperture.py 711 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
    gauss g1 m1 n2 1e-3 0
    xaxis m1 r_ap lin 0.1e-3 2e-3 10
    """
    
    kat = finesse.kat()
    kat.parseCommands(code)
    
    maxtem = np.linspace(0, 1, 2)
    
    #kat.m1.n2.q = gauss_param(w0=1e-3, z=0)
    
    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()