Skip to content
Snippets Groups Projects
Select Git revision
  • 7a4f6ba8cb6654ceae89e56b094d0cd8ba3552f0
  • master default protected
2 results

RDGW150914_ptemcee_diff1.py

Blame
  • Forked from Xisco Jimenez Forteza / RDStackingProject
    Source project has a limited visibility.
    test_plot.py 851 B
    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
    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)
    
    kat.add(cavity('cav1','m1','n3','m2','n4'))
    
    kat.add(photodiode('pd_cav','n4'))
    kat.add(photodiode('pd_ref','n2'))
    kat.add(photodiode('pd_trs','n5'))
    
    kat.add(xaxis(Scale.linear, [0, 360], kat.m2, kat.m2.phi, 100))
    
    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)
    
    pl.figure()
    pl.plot(run.x,run.y)
    pl.xlabel(run.xlabel)
    pl.ylabel("Intensity [W]")
    pl.legend(run.ylabels)
    pl.show()