adding kat files for the 4 'builddown' configs and single_test.py to check these quickly

playground/anna/PI_paper/kat_files/single_test.py

+from pykat import finesse
+from pykat.commands import *
+import pykat.gw_detectors.ifo as ifo
+import matplotlib.pyplot as plt
+import tune_tools as tt
+from tune import testing
+
+LIGO1 = ifo.aLIGO(debug=True)
+kat = finesse.kat()
+
+kat.loadKatFile("tuned/PI-Cav.kat")
+
+LIGO1.APOP_f2 = ifo.port("APOP_f2","nAPOP", LIGO1.f2, phase=13)
+LIGO1.aMICH = ifo.DOF("aMICH", LIGO1.APOP_f2, "Q", ["ITMX", "ETMX", "ITMY", "ETMY"], [1,1,-1,-1], 100.0)
+
+dofs = [LIGO1.CARM]#, LIGO1.DARM, LIGO1.aMICH, LIGO1.PRCL, LIGO1.SRCL]
+locks = [kat.ETMX_lock]#, kat.ETMY_lock, kat.MICH_lock, kat.ITMY_lock, kat.PRCL_lock, kat.SRCL_lock]
+
+testing(kat,dofs,locks,prc=False,opt='ITMX')
+##testing(kat,dofs,locks,prc=True,xlim=1e-2,opt="BS")
+
+END = input('press enter to finish and close all plots')#to override auto-closing.
+
+
+
+#TODO: tunings on ITM ETM still offset from zero.
\ No newline at end of file

playground/anna/PI_paper/kat_files/tune.py

@@ -22,14 +22,6 @@ def main():
 maxtem {0}
 phase {1}
     """.format(kat.maxtem, kat.phase)
-    #######################
-    #creating auxMICH port that works for non-recycled MICH monitoring
-    LIGO1.APOP_f2 = ifo.port("APOP_f2","nAPOP", LIGO1.f2, phase=13)
-    LIGO1.aMICH = ifo.DOF("aMICH", LIGO1.APOP_f2, "Q", ["ITMX", "ETMX", "ITMY", "ETMY"], [1,1,-1,-1], 100.0)
-    #######################
-    #for testing
-    dofs = [LIGO1.CARM, LIGO1.DARM, LIGO1.aMICH, LIGO1.PRCL, LIGO1.SRCL]
-    locks = [kat.ETMX_lock, kat.ETMY_lock, kat.MICH_lock, kat.ITMY_lock, kat.PRCL_lock, kat.SRCL_lock]
 
     # #Making lock blocks
     # ######################## #high masses to ensure TF calcs work in DRMi_tuning
@@ -69,8 +61,20 @@ phase {1}
     # o.close()
     # ## now *manually* edit the file to use aMICH as MICH_lock rather than MICH.
 
+    ##for testing
+    #######################
+    #creating auxMICH port that works for non-recycled MICH monitoring
+    LIGO1.APOP_f2 = ifo.port("APOP_f2","nAPOP", LIGO1.f2, phase=13)
+    LIGO1.aMICH = ifo.DOF("aMICH", LIGO1.APOP_f2, "Q", ["ITMX", "ETMX", "ITMY", "ETMY"], [1,1,-1,-1], 100.0)
+    #######################
+    #for testing
+    dofs = [LIGO1.CARM, LIGO1.DARM, LIGO1.aMICH, LIGO1.PRCL, LIGO1.SRCL]
+    locks = [kat.ETMX_lock, kat.ETMY_lock, kat.MICH_lock, kat.ITMY_lock, kat.PRCL_lock, kat.SRCL_lock]
     
+    # kat = finesse.kat()
+    # kat.loadKatFile("tuned/PI-DRMi.kat")
     testing(kat,dofs,locks)
+    
     ###############  cut down  #####################
     # ######################## # DR --> PRMi
     #removing kat lines
@@ -92,7 +96,7 @@ phase {1}
     remove SRCL lock commands and puts
     kat.L0.P = {0}
     kat.DARM_lock.gain = {1}
-    kat.POP_f2_Q.phase1 = {2}
+    kat.APOP_f2.phase1 = {2}
     ---------------------------------------------------
     """.format(kat.L0.P,kat.DARM_lock.gain,kat.APOP_f2.phase1))
     testing(kat,dofs[:-1],locks[:-1])

playground/anna/PI_paper/kat_files/tune_tools.py

@@ -81,13 +81,13 @@ def test_locks(_kat,locks,PRC=False,optic="BS"):
     if optic =="BS":
         # test lock with BS tuning
         code = """
-        xaxis BS phi lin 0 2 100 
+        xaxis* BS phi lin 0 2 100 
         yaxis abs
         #showiterate -1
         """
     else:
         code = """
-        xaxis {0} phi lin 0 2 100 
+        xaxis* {0} phi lin 0 2 100 
         yaxis abs
         #showiterate -1
         """.format(optic)

playground/anna/PI_paper/kat_files/tuned/PI-Cav.kat

+
+%%% FTblock NO_BLOCK
+# reintroduced: errsigs, locks, powers (partial)
+%%% FTend NO_BLOCK
+
+%%% FTblock header
+#--------------------------------------------------------------------------
+# an aLIGO file to introduce a PI to one arm
+# based on
+# aLIGO_IFO_AWC_tuning_DCoff_maxtem4.kat (details in the original file header)
+#
+# updated to re-include locks
+# (errsig, powers and locks blocks, plus EOMs 1 and 2 (9 and 45kHz))
+#
+# - R. Abbott. Advanced LIGO length sensing and control final design. LIGO Document T1000298, 2010.
+# - Advanced LIGO, CQG, 2015
+#   http://iopscience.iop.org/article/10.1088/0264-9381/32/7/074001/meta#cqg507871s4-8
+# - G. Billingsley. LIGO core optics reference page. https://galaxy.ligo. caltech.edu/optics/
+#
+# Anna Green, Andreas Freise, 28th November 2016
+#--------------------------------------------------------------------------
+#
+# References used for aLIGO parameters
+# [1] R. Abbott. Advanced LIGO length sensing and control final design. LIGO Document T1000298, 2010.
+# [2] Advanced LIGO, CQG, 2015
+#     http://iopscience.iop.org/article/10.1088/0264-9381/32/7/074001/meta#cqg507871s4-8
+# [3] G. Billingsley. LIGO core optics reference page. https://galaxy.ligo.caltech.edu/optics/
+# [4] D. Martynov, `Lock Acquisition and Sensitivity Analysis of Advanced
+#     LIGO Interferometers', PhD thesis, Caltech, 2015
+#     http://thesis.library.caltech.edu/8899/1/DenisMartynovThesis.pdf
+# [4] A. Staley, `Locking the Advanced LIGO Gravitational Wave Detector:
+#     with a focus on the Arm Length Stabilization Technique', PhD thesis,
+#     Columbia University, 2015
+#     https://academiccommons.columbia.edu/catalog/ac:189457
+#
+# Changes since 12/2016:
+# - 06.12.2016, initial corrections compared to previous files, Anna and Andreas
+#
+#   Bug fixes
+#   - AR coatings of ETMs and ITMs were using mirror components with non-zero
+#     reflectivity. This is wrng as it causes residual resoances with the HR
+#     surfaces and ignores the wedge in the substrate. Instead the mirror component
+#     should have losses but no reflectivity. Alternatively beamsplitter components
+#     with non-zero reflectivity can be used to detect the AR reflections
+#   - Same as above for PRM and SRM AR coatings
+#   - change order of `cav' commands, arm cavities need to come first
+#
+#   Change of parameter values
+#   - BS substrate length were assymetric and did not match substrate thickness,
+#     changed to symmetric 0.0687m, using pykat.GWdetectors.aLIGO.BSpath(0.06)
+#   - changed SRM transmission from 35% to 20%, matching the aLIGO 2015 paper
+#   - changed arm length to L=3994.4692, see [4], Table 10.2
+#
+#   Changes in file structure
+#   - removed length calculation fucntions for PRC and SRC sub-length
+#     and added length values into space components
+#   - removed maxtem 4 command from HOMs block
+#
+#   Naming changes
+#   - changed naming for thermal lens components
+#   - changed naming of PRMHR to PRM and node names in PRM substrate
+#   - changed naming of SRMHR to SRM and node names in SRM substrate
+#   - changed name of ITMHR to ITM and ETMHR to ETM
+#   - changed names of arm spaces from LXarm, LYarm to LX, LY
+#		- changed name of HOMs block to cavities
+#		- changed PR2 node4 name from dump11 to nAPOP (auxilliary POP)
+#
+#
+#--------------------------------------------------------------------------
+%%% FTend header
+
+%%% FTblock laser
+###########################################################################
+# Laser and input optics
+l L0 2001.999 0 n0
+s lmod1 1 n0 n1
+# modulators for core interferometer sensing - Advanced LIGO, CQG, 2015
+# http://iopscience.iop.org/article/10.1088/0264-9381/32/7/074001/meta#cqg507871s4-8
+# 9MHz (CARM, PRC, SRC loops)
+const f1 9099471
+const f2 45497355
+mod mod1 $f1 0.18 1 pm n1 n2
+s lmod2 1 n2 n3
+# 45MHz (MICH, SRC loops)
+mod mod2 $f2 0.18 1 pm n3 n4
+s lmod3 1 n4 nREFL # TODO
+###########################################################################
+%%% FTend laser
+
+%%% FTblock PRC
+###########################################################################
+# PRM
+# AR surface
+m2 PRMAR 0 0 0 nREFL nPRMs1
+#bs2 PRMAR 35u 4.5u 0 0 nREFL nPRMAR1 nPRMs1 nPRMAR2
+# Substrate
+s sPRMsub1 0.0737 $nsilica nPRMs1 nPRMs2
+# HR surface
+m1 PRM 1 0 $phi_PRM nPRMs2 nPRM2
+attr PRM Rc 11.009
+# Distance between PRM and PR2
+s lp1 16.6107 nPRM2 nPR2a
+# PR2
+bs1 PR2 250u $Mloss 0 -0.79 nPR2a nPR2b nPOP nAPOP
+attr PR2 Rc -4.545
+# Distance from PR2 to PR3
+s lp2 16.1647 nPR2b nPR3a
+# PR3
+bs1 PR3 0 $Mloss 0 0.615 nPR3a nPR3b dump12 dump13
+attr PR3 Rc 36.027
+# Distance from PR3
+s lp3 19.5381 nPR3b nPRBS
+
+###########################################################################
+%%% FTend PRC
+
+%%% FTblock BS
+###########################################################################
+# BS beamsplitter
+##------------------------------------------------------------
+## BS
+##                             ^
+##                  to IMY     |
+##                             |      ,'-.
+##                             |     +    `.
+##                        nYBS |   ,'       :'
+##      nPR3b                  |  +i1      +
+##         ---------------->    ,:._  i2 ,'
+##    from the PRC       nPRBS + \  `-. + nXBS
+##                           ,' i3\   ,' --------------->
+##                          +      \ +     to IMX
+##                        ,'     i4.'
+##                       `._      ..
+##                          `._ ,' |nSRBS
+##                             -   |
+##                                 |to the SRC
+##                                 |
+##                                 v
+##------------------------------------------------------------
+bs1 BS 1 0 $phi_BS 45 nPRBS nYBS nBSi1 nBSi3
+s BSsub1 0.0687 $nsilica nBSi1 nBSi2
+s BSsub2 0.0687 $nsilica nBSi3 nBSi4
+bs2 BSAR1 50u 0 0 -29.195 nBSi2 dump14 nXBS nPOX
+bs2 BSAR2 50u 0 0 29.195 nBSi4 dump15 nSRBS dump16
+
+###########################################################################
+%%% FTend BS
+
+%%% FTblock Yarm
+###########################################################################
+# Distance from beam splitter to Y arm input mirror
+s ly1 5.0126 nYBS nITMY1a
+
+# Thermal lens correction
+lens ITMY_therm $th_f nITMY1a nITM_th1
+s ITMY_th1 0 nITM_th1 nITM_th2
+lens ITMY_lens inf nITM_th2 nITM_th3
+s ITMY_th2 0 nITM_th3 nITMY1
+
+# Y arm input mirror
+m2 ITMYAR 0 0 0 nITMY1 nITMYs1
+#bs2 ITMYAR 20u 0 0 0 nITMY1 nITMYAR1 nITMYs1 nITMYAR2
+s ITMYsub 0.2 $nsilica nITMYs1 nITMYs2
+m1 ITMY 1 0 $phi_ITMY nITMYs2 nITMY2
+attr ITMY Rc -1934
+
+# Y arm length
+s LY $Larm nITMY2 nETMY1
+
+# Y arm end mirror
+m1 ETMY 1 0 $phi_ETMY nETMY1 nETMYs1
+s ETMYsub 0.2 $nsilica nETMYs1 nETMYs2
+m2 ETMYAR 0 0 0 nETMYs2 nPTY
+#bs2 ETMYAR 500u 0 0 0 nETMYs2 nETMYAR1 nPTY nETMYAR1
+attr ETMY Rc 2245
+attr ETMY mass 40
+attr ITMY mass 40
+
+###########################################################################
+%%% FTend Yarm
+
+%%% FTblock Xarm
+###########################################################################
+
+# Distance from beam splitter to X arm input mirror
+s lx1 4.993 nXBS nITMX1a
+
+# Thermal lens correction
+lens ITMX_therm $th_f nITMX1a nITMX_th1
+s ITMX_th1 0 nITMX_th1 nITMX_th2
+lens ITMX_lens inf nITMX_th2 nITMX_th3
+s ITMX_th2 0 nITMX_th3 nITMX1
+
+# X arm input mirror
+m2 ITMXAR 0 20u 0 nITMX1 nITMXs1
+#bs2 ITMXAR 20u 0 0 0 nITMX1 nITMXAR1 nITMXs1 nITMXAR2
+s ITMXsub 0.2 $nsilica nITMXs1 nITMXs2
+m1 ITMX 0.014 $Mloss $phi_ITMX nITMXs2 nITMX2
+attr ITMX Rc -1934
+
+# X arm length
+s LX $Larm nITMX2 nETMX1
+
+# X arm end mirror
+m1 ETMX 5u $Mloss $phi_ETMX nETMX1 nETMXs1
+s ETMXsub 0.2 $nsilica nETMXs1 nETMXs2
+#bs2 ETMXAR 500u 0 0 0 nETMXs2 nETMXAR1 nPTX nETMXAR1
+m2 ETMXAR 0 500u 0 nETMXs2 nPTX
+attr ETMX Rc 2245
+attr ETMX mass 40
+attr ITMX mass 40
+
+###########################################################################
+%%% FTend Xarm
+
+%%% FTblock SRC
+###########################################################################
+# Distance to SR3
+s ls3 19.3661 nSRBS nSR3b
+# SR3
+bs1 SR3 0 $Mloss 0 0.785 nSR3b nSR3a dump17 dump18
+attr SR3 Rc 35.972841
+# Distance from SR3 to SR2
+s ls2 15.4435 nSR3a nSR2b
+# SR2
+bs1 SR2 0 $Mloss 0 -0.87 nSR2b nSR2a dump19 dump20
+attr SR2 Rc -6.406
+# Distance from SR2 to SRM
+s ls1 15.7586 nSR2a nSRM1
+# Signal recycling mirror SRM-08
+m1 SRM 1 0 $phi_SRM nSRM1 nSRMs1
+s SRMsub 0.0749 $nsilica nSRMs1 nSRMs2
+m2 SRMAR 0 0 0 nSRMs2 nSRM2
+#bs2 SRMAR 50n 0 0 0 nSRMs2 nSRMAR1 nSRM2 nSRMAR2
+attr SRM Rc -5.6938
+
+###########################################################################
+%%% FTend SRC
+
+%%% FTblock cavities
+###########################################################################
+#cav cavIMC MC2 nMC2in MC2 nMC2refl
+cav cavXARM ITMX nITMX2 ETMX nETMX1
+cav cavYARM ITMY nITMY2 ETMY nETMY1
+cav cavPRX PRM nPRM2 ITMX nITMXs2
+cav cavPRY PRM nPRM2 ITMY nITMYs2
+cav cavSRX SRM nSRM1 ITMX nITMXs2
+cav cavSRY SRM nSRM1 ITMY nITMYs2
+# removed OMC
+
+###########################################################################
+%%% FTend cavities
+
+%%% FTblock constants
+###########################################################################
+const nsilica 1.44963098985906
+const nTGG 1.954
+const nCalcite 1.65846
+const Mloss 37.5u
+const Larm 3994.4692
+const fM 36.397884M
+const fP 54.596826M
+const th_f 34.5k
+###########################################################################
+%%% FTend constants
+
+%%% FTblock powers
+###########################################################################
+pd P_DC_AS  nSRM2
+# pd P_DC_OMC nOMC_HROC_trans
+pd PIMCtrans nREFL*
+pd Px nITMX2
+pd Py nITMY2
+pd Pprc nPRM2
+pd Psrc nSRM1*
+ad prc0 0 nPRM2
+ad prcf1 $f1 nPRM2
+ad prcf2 $f2 nPRM2
+
+#ad aoc0 0 nOMC_HROC_trans
+ad asc0 0 nSRM2
+ad asf1 $f1 nSRM2
+ad asf2 $f2 nSRM2
+ad src0 0 nSRM1*
+ad srcf1 $f1 nSRM1*
+ad srcf2 $f2 nSRM1*
+###########################################################################
+%%% FTend powers
+
+%%% FTblock tunings
+###########################################################################
+const phi_ITMX -0.0012
+const phi_ITMY 0.0012
+const phi_ETMX -0.0046012985
+const phi_ETMY 0.0046015197
+const phi_BS   0.0
+const phi_PRM  0.003
+const phi_SRM  -90.0012
+###########################################################################
+%%% FTend tunings
+
+%%% FTblock errsigs
+###########################################################################
+##updated to use nAPOP port for MICH locks rather than nPOP (default)
+pd AS_DC nSRM2
+pd1 REFL_f1_I 9099471.0 101 nREFL
+pd1 POP_f1_I 9099471.0 101 nPOP
+pd1 APOP_f2 45497355.0 93.0 nAPOP
+#pd1 POP_f2_Q 45497355.0 103.0 nPOP
+pd1 REFL_f2_I 45497355.0 14 nREFL
+noplot AS_DC
+noplot REFL_f1_I
+noplot POP_f1_I
+#noplot POP_f2_Q
+noplot APOP_f2
+noplot REFL_f2_I
+###########################################################################
+%%% FTend errsigs
+
+%%% FTblock locks
+###########################################################################
+set CARM_err REFL_f1_I re
+lock ETMX_lock $CARM_err  0.00039  0.00043
+put* ETMX phi $ETMX_lock
+noplot ETMX_lock
+###########################################################################
+
+%%% FTend locks
+
+maxtem 0
+phase 2
+    
\ No newline at end of file

playground/anna/PI_paper/kat_files/tuned/PI-DRMi.kat

@@ -301,11 +301,12 @@ const phi_SRM  -90.0012
 
 %%% FTblock errsigs
 ###########################################################################
+##updated to use nAPOP port for MICH locks rather than nPOP (default)
 pd AS_DC nSRM2
 pd1 REFL_f1_I 9099471.0 101 nREFL
 pd1 POP_f1_I 9099471.0 101 nPOP
-#pd1 POP_f2_Q 45497355.0 103.0 nPOP
 pd1 APOP_f2 45497355.0 103.0 nAPOP
+#pd1 POP_f2_Q 45497355.0 103.0 nPOP
 pd1 REFL_f2_I 45497355.0 14 nREFL
 noplot AS_DC
 noplot REFL_f1_I
@@ -318,6 +319,7 @@ noplot REFL_f2_I
 
 %%% FTblock locks
 ###########################################################################
+##updated to lock using auxMICH lock rather than MICH_lock
 set AS_f2_I_re AS_DC re
 set CARM_err REFL_f1_I re
 set PRCL_err POP_f1_I re
@@ -329,8 +331,8 @@ func DARM_err = $AS_f2_I_re - 1.56140670384883
 lock DARM_lock $DARM_err   0.0035   0.0048
 lock CARM_lock $CARM_err  0.00039  0.00043 
 lock PRCL_lock $PRCL_err     -2.8  0.00012
-#lock MICH_lock $MICH_err       43  3.9e-05 
 lock MICH_lock $aMICH_err     43  3.9e-05
+#lock MICH_lock $MICH_err       43  3.9e-05 
 lock SRCL_lock $SRCL_err      2.6   0.0016
 noplot ITMY_lock
 func ITMY_lock = (-1.0) * $MICH_lock 

playground/anna/PI_paper/kat_files/tuned/PI-FPMi.kat

+
+%%% FTblock NO_BLOCK
+# reintroduced: errsigs, locks, powers (partial)
+%%% FTend NO_BLOCK
+
+%%% FTblock header
+#--------------------------------------------------------------------------
+# an aLIGO file to introduce a PI to one arm
+# based on
+# aLIGO_IFO_AWC_tuning_DCoff_maxtem4.kat (details in the original file header)
+#
+# updated to re-include locks
+# (errsig, powers and locks blocks, plus EOMs 1 and 2 (9 and 45kHz))
+#
+# - R. Abbott. Advanced LIGO length sensing and control final design. LIGO Document T1000298, 2010.
+# - Advanced LIGO, CQG, 2015
+#   http://iopscience.iop.org/article/10.1088/0264-9381/32/7/074001/meta#cqg507871s4-8
+# - G. Billingsley. LIGO core optics reference page. https://galaxy.ligo. caltech.edu/optics/
+#
+# Anna Green, Andreas Freise, 28th November 2016
+#--------------------------------------------------------------------------
+#
+# References used for aLIGO parameters
+# [1] R. Abbott. Advanced LIGO length sensing and control final design. LIGO Document T1000298, 2010.
+# [2] Advanced LIGO, CQG, 2015
+#     http://iopscience.iop.org/article/10.1088/0264-9381/32/7/074001/meta#cqg507871s4-8
+# [3] G. Billingsley. LIGO core optics reference page. https://galaxy.ligo.caltech.edu/optics/
+# [4] D. Martynov, `Lock Acquisition and Sensitivity Analysis of Advanced
+#     LIGO Interferometers', PhD thesis, Caltech, 2015
+#     http://thesis.library.caltech.edu/8899/1/DenisMartynovThesis.pdf
+# [4] A. Staley, `Locking the Advanced LIGO Gravitational Wave Detector:
+#     with a focus on the Arm Length Stabilization Technique', PhD thesis,
+#     Columbia University, 2015
+#     https://academiccommons.columbia.edu/catalog/ac:189457
+#
+# Changes since 12/2016:
+# - 06.12.2016, initial corrections compared to previous files, Anna and Andreas
+#
+#   Bug fixes
+#   - AR coatings of ETMs and ITMs were using mirror components with non-zero
+#     reflectivity. This is wrng as it causes residual resoances with the HR
+#     surfaces and ignores the wedge in the substrate. Instead the mirror component
+#     should have losses but no reflectivity. Alternatively beamsplitter components
+#     with non-zero reflectivity can be used to detect the AR reflections
+#   - Same as above for PRM and SRM AR coatings
+#   - change order of `cav' commands, arm cavities need to come first
+#
+#   Change of parameter values
+#   - BS substrate length were assymetric and did not match substrate thickness,
+#     changed to symmetric 0.0687m, using pykat.GWdetectors.aLIGO.BSpath(0.06)
+#   - changed SRM transmission from 35% to 20%, matching the aLIGO 2015 paper
+#   - changed arm length to L=3994.4692, see [4], Table 10.2
+#
+#   Changes in file structure
+#   - removed length calculation fucntions for PRC and SRC sub-length
+#     and added length values into space components
+#   - removed maxtem 4 command from HOMs block
+#
+#   Naming changes
+#   - changed naming for thermal lens components
+#   - changed naming of PRMHR to PRM and node names in PRM substrate
+#   - changed naming of SRMHR to SRM and node names in SRM substrate
+#   - changed name of ITMHR to ITM and ETMHR to ETM
+#   - changed names of arm spaces from LXarm, LYarm to LX, LY
+#		- changed name of HOMs block to cavities
+#		- changed PR2 node4 name from dump11 to nAPOP (auxilliary POP)
+#
+#
+#--------------------------------------------------------------------------
+%%% FTend header
+
+%%% FTblock laser
+###########################################################################
+# Laser and input optics
+l L0 4003.998 0 n0
+s lmod1 1 n0 n1
+# modulators for core interferometer sensing - Advanced LIGO, CQG, 2015
+# http://iopscience.iop.org/article/10.1088/0264-9381/32/7/074001/meta#cqg507871s4-8
+# 9MHz (CARM, PRC, SRC loops)
+const f1 9099471
+const f2 45497355
+mod mod1 $f1 0.18 1 pm n1 n2
+s lmod2 1 n2 n3
+# 45MHz (MICH, SRC loops)
+mod mod2 $f2 0.18 1 pm n3 n4
+s lmod3 1 n4 nREFL # TODO
+###########################################################################
+%%% FTend laser
+
+%%% FTblock PRC
+###########################################################################
+# PRM
+# AR surface
+m2 PRMAR 0 0 0 nREFL nPRMs1
+#bs2 PRMAR 35u 4.5u 0 0 nREFL nPRMAR1 nPRMs1 nPRMAR2
+# Substrate
+s sPRMsub1 0.0737 $nsilica nPRMs1 nPRMs2
+# HR surface
+m1 PRM 1 0 $phi_PRM nPRMs2 nPRM2
+attr PRM Rc 11.009
+# Distance between PRM and PR2
+s lp1 16.6107 nPRM2 nPR2a
+# PR2
+bs1 PR2 250u $Mloss 0 -0.79 nPR2a nPR2b nPOP nAPOP
+attr PR2 Rc -4.545
+# Distance from PR2 to PR3
+s lp2 16.1647 nPR2b nPR3a
+# PR3
+bs1 PR3 0 $Mloss 0 0.615 nPR3a nPR3b dump12 dump13
+attr PR3 Rc 36.027
+# Distance from PR3
+s lp3 19.5381 nPR3b nPRBS
+
+###########################################################################
+%%% FTend PRC
+
+%%% FTblock BS
+###########################################################################
+# BS beamsplitter
+##------------------------------------------------------------
+## BS
+##                             ^
+##                  to IMY     |
+##                             |      ,'-.
+##                             |     +    `.
+##                        nYBS |   ,'       :'
+##      nPR3b                  |  +i1      +
+##         ---------------->    ,:._  i2 ,'
+##    from the PRC       nPRBS + \  `-. + nXBS
+##                           ,' i3\   ,' --------------->
+##                          +      \ +     to IMX
+##                        ,'     i4.'
+##                       `._      ..
+##                          `._ ,' |nSRBS
+##                             -   |
+##                                 |to the SRC
+##                                 |
+##                                 v
+##------------------------------------------------------------
+bs1 BS 0.5 $Mloss $phi_BS 45 nPRBS nYBS nBSi1 nBSi3
+s BSsub1 0.0687 $nsilica nBSi1 nBSi2
+s BSsub2 0.0687 $nsilica nBSi3 nBSi4
+bs2 BSAR1 50u 0 0 -29.195 nBSi2 dump14 nXBS nPOX
+bs2 BSAR2 50u 0 0 29.195 nBSi4 dump15 nSRBS dump16
+
+###########################################################################
+%%% FTend BS
+
+%%% FTblock Yarm
+###########################################################################
+# Distance from beam splitter to Y arm input mirror
+s ly1 5.0126 nYBS nITMY1a
+
+# Thermal lens correction
+lens ITMY_therm $th_f nITMY1a nITM_th1
+s ITMY_th1 0 nITM_th1 nITM_th2
+lens ITMY_lens inf nITM_th2 nITM_th3
+s ITMY_th2 0 nITM_th3 nITMY1
+
+# Y arm input mirror
+m2 ITMYAR 0 20u 0 nITMY1 nITMYs1
+#bs2 ITMYAR 20u 0 0 0 nITMY1 nITMYAR1 nITMYs1 nITMYAR2
+s ITMYsub 0.2 $nsilica nITMYs1 nITMYs2
+m1 ITMY 0.014 $Mloss $phi_ITMY nITMYs2 nITMY2
+attr ITMY Rc -1934
+
+# Y arm length
+s LY $Larm nITMY2 nETMY1
+
+# Y arm end mirror
+m1 ETMY 5u $Mloss $phi_ETMY nETMY1 nETMYs1
+s ETMYsub 0.2 $nsilica nETMYs1 nETMYs2
+m2 ETMYAR 0 500u 0 nETMYs2 nPTY
+#bs2 ETMYAR 500u 0 0 0 nETMYs2 nETMYAR1 nPTY nETMYAR1
+attr ETMY Rc 2245
+attr ETMY mass 40
+attr ITMY mass 40
+
+###########################################################################
+%%% FTend Yarm
+
+%%% FTblock Xarm
+###########################################################################
+
+# Distance from beam splitter to X arm input mirror
+s lx1 4.993 nXBS nITMX1a
+
+# Thermal lens correction
+lens ITMX_therm $th_f nITMX1a nITMX_th1
+s ITMX_th1 0 nITMX_th1 nITMX_th2
+lens ITMX_lens inf nITMX_th2 nITMX_th3
+s ITMX_th2 0 nITMX_th3 nITMX1
+
+# X arm input mirror
+m2 ITMXAR 0 20u 0 nITMX1 nITMXs1
+#bs2 ITMXAR 20u 0 0 0 nITMX1 nITMXAR1 nITMXs1 nITMXAR2
+s ITMXsub 0.2 $nsilica nITMXs1 nITMXs2
+m1 ITMX 0.014 $Mloss $phi_ITMX nITMXs2 nITMX2
+attr ITMX Rc -1934
+
+# X arm length
+s LX $Larm nITMX2 nETMX1
+
+# X arm end mirror
+m1 ETMX 5u $Mloss $phi_ETMX nETMX1 nETMXs1
+s ETMXsub 0.2 $nsilica nETMXs1 nETMXs2
+#bs2 ETMXAR 500u 0 0 0 nETMXs2 nETMXAR1 nPTX nETMXAR1
+m2 ETMXAR 0 500u 0 nETMXs2 nPTX
+attr ETMX Rc 2245
+attr ETMX mass 40
+attr ITMX mass 40
+
+###########################################################################
+%%% FTend Xarm
+
+%%% FTblock SRC
+###########################################################################
+# Distance to SR3
+s ls3 19.3661 nSRBS nSR3b
+# SR3
+bs1 SR3 0 $Mloss 0 0.785 nSR3b nSR3a dump17 dump18
+attr SR3 Rc 35.972841
+# Distance from SR3 to SR2
+s ls2 15.4435 nSR3a nSR2b
+# SR2
+bs1 SR2 0 $Mloss 0 -0.87 nSR2b nSR2a dump19 dump20
+attr SR2 Rc -6.406
+# Distance from SR2 to SRM
+s ls1 15.7586 nSR2a nSRM1
+# Signal recycling mirror SRM-08
+m1 SRM 1 0 $phi_SRM nSRM1 nSRMs1
+s SRMsub 0.0749 $nsilica nSRMs1 nSRMs2
+m2 SRMAR 0 0 0 nSRMs2 nSRM2
+#bs2 SRMAR 50n 0 0 0 nSRMs2 nSRMAR1 nSRM2 nSRMAR2
+attr SRM Rc -5.6938
+
+###########################################################################
+%%% FTend SRC
+
+%%% FTblock cavities
+###########################################################################
+#cav cavIMC MC2 nMC2in MC2 nMC2refl
+cav cavXARM ITMX nITMX2 ETMX nETMX1
+cav cavYARM ITMY nITMY2 ETMY nETMY1
+cav cavPRX PRM nPRM2 ITMX nITMXs2
+cav cavPRY PRM nPRM2 ITMY nITMYs2
+cav cavSRX SRM nSRM1 ITMX nITMXs2
+cav cavSRY SRM nSRM1 ITMY nITMYs2
+# removed OMC
+
+###########################################################################
+%%% FTend cavities
+
+%%% FTblock constants
+###########################################################################
+const nsilica 1.44963098985906
+const nTGG 1.954
+const nCalcite 1.65846
+const Mloss 37.5u
+const Larm 3994.4692
+const fM 36.397884M
+const fP 54.596826M
+const th_f 34.5k
+###########################################################################
+%%% FTend constants
+
+%%% FTblock powers
+###########################################################################
+pd P_DC_AS  nSRM2
+# pd P_DC_OMC nOMC_HROC_trans
+pd PIMCtrans nREFL*
+pd Px nITMX2
+pd Py nITMY2
+pd Pprc nPRM2
+pd Psrc nSRM1*
+ad prc0 0 nPRM2
+ad prcf1 $f1 nPRM2
+ad prcf2 $f2 nPRM2
+
+#ad aoc0 0 nOMC_HROC_trans
+ad asc0 0 nSRM2
+ad asf1 $f1 nSRM2
+ad asf2 $f2 nSRM2
+ad src0 0 nSRM1*
+ad srcf1 $f1 nSRM1*
+ad srcf2 $f2 nSRM1*
+###########################################################################
+%%% FTend powers
+
+%%% FTblock tunings
+###########################################################################
+const phi_ITMX -0.0012
+const phi_ITMY 0.0012
+const phi_ETMX -0.0046012985
+const phi_ETMY 0.0046015197
+const phi_BS   0.0
+const phi_PRM  0.003
+const phi_SRM  -90.0012
+###########################################################################
+%%% FTend tunings
+
+%%% FTblock errsigs
+###########################################################################
+##updated to use nAPOP port for MICH locks rather than nPOP (default)
+pd AS_DC nSRM2
+pd1 REFL_f1_I 9099471.0 101 nREFL
+pd1 POP_f1_I 9099471.0 101 nPOP
+pd1 APOP_f2 45497355.0 93.0 nAPOP
+#pd1 POP_f2_Q 45497355.0 103.0 nPOP
+pd1 REFL_f2_I 45497355.0 14 nREFL
+noplot AS_DC
+noplot REFL_f1_I
+noplot POP_f1_I
+#noplot POP_f2_Q
+noplot APOP_f2
+noplot REFL_f2_I
+###########################################################################
+%%% FTend errsigs
+
+%%% FTblock locks
+###########################################################################
+##updated to lock using auxMICH lock rather than MICH_lock
+set AS_f2_I_re AS_DC re
+set CARM_err REFL_f1_I re
+#set MICH_err POP_f2_Q re
+set aMICH_err APOP_f2 re
+func DARM_err = $AS_f2_I_re - 1.56140670384883
+
+lock DARM_lock $DARM_err  0.00035   0.0048
+lock CARM_lock $CARM_err  0.00039  0.00043 
+lock MICH_lock $aMICH_err     430  3.9e-07
+#lock MICH_lock $MICH_err       43  3.9e-05 
+noplot ITMY_lock
+func ITMY_lock = (-1.0) * $MICH_lock 
+func ETMX_lock = $CARM_lock + $MICH_lock + $DARM_lock
+func ETMY_lock = $CARM_lock - $MICH_lock - $DARM_lock
+
+put* ITMX    phi     $MICH_lock
+put* ITMY    phi     $ITMY_lock
+put* ETMX    phi     $ETMX_lock
+put* ETMY    phi     $ETMY_lock
+
+noplot MICH_lock
+noplot DARM_lock
+noplot CARM_lock
+noplot ETMX_lock
+noplot ETMY_lock
+###########################################################################
+
+%%% FTend locks
+
+maxtem 0
+phase 2
+    
\ No newline at end of file

playground/anna/PI_paper/kat_files/tuned/PI-PRMi.kat

@@ -17,7 +17,7 @@
 #   http://iopscience.iop.org/article/10.1088/0264-9381/32/7/074001/meta#cqg507871s4-8
 # - G. Billingsley. LIGO core optics reference page. https://galaxy.ligo. caltech.edu/optics/
 #
-# Anna Green, 28th November 2016
+# Anna Green, Andreas Freise, 28th November 2016
 #--------------------------------------------------------------------------
 #
 # References used for aLIGO parameters
@@ -54,7 +54,7 @@
 #   Changes in file structure
 #   - removed length calculation fucntions for PRC and SRC sub-length
 #     and added length values into space components
-#   -
+#   - removed maxtem 4 command from HOMs block
 #
 #   Naming changes
 #   - changed naming for thermal lens components
@@ -62,6 +62,8 @@
 #   - changed naming of SRMHR to SRM and node names in SRM substrate
 #   - changed name of ITMHR to ITM and ETMHR to ETM
 #   - changed names of arm spaces from LXarm, LYarm to LX, LY
+#		- changed name of HOMs block to cavities
+#		- changed PR2 node4 name from dump11 to nAPOP (auxilliary POP)
 #
 #
 #--------------------------------------------------------------------------
@@ -99,7 +101,7 @@ attr PRM Rc 11.009
 # Distance between PRM and PR2
 s lp1 16.6107 nPRM2 nPR2a
 # PR2
-bs1 PR2 250u $Mloss 0 -0.79 nPR2a nPR2b nPOP dump11
+bs1 PR2 250u $Mloss 0 -0.79 nPR2a nPR2b nPOP nAPOP
 attr PR2 Rc -4.545
 # Distance from PR2 to PR3
 s lp2 16.1647 nPR2b nPR3a
@@ -299,32 +301,37 @@ const phi_SRM  -90.0012
 
 %%% FTblock errsigs
 ###########################################################################
+##updated to use nAPOP port for MICH locks rather than nPOP (default)
 pd AS_DC nSRM2
 pd1 REFL_f1_I 9099471.0 101 nREFL
 pd1 POP_f1_I 9099471.0 101 nPOP
-pd1 POP_f2_Q 45497355.0 93.0 nPOP
+pd1 APOP_f2 45497355.0 93.0 nAPOP
+#pd1 POP_f2_Q 45497355.0 103.0 nPOP
 pd1 REFL_f2_I 45497355.0 14 nREFL
 noplot AS_DC
 noplot REFL_f1_I
 noplot POP_f1_I
-noplot POP_f2_Q
+#noplot POP_f2_Q
+noplot APOP_f2
 noplot REFL_f2_I
 ###########################################################################
 %%% FTend errsigs
 
 %%% FTblock locks
 ###########################################################################
+##updated to lock using auxMICH lock rather than MICH_lock
 set AS_f2_I_re AS_DC re
 set CARM_err REFL_f1_I re
 set PRCL_err POP_f1_I re
-set MICH_err POP_f2_Q re
-set SRCL_err REFL_f2_I re
+#set MICH_err POP_f2_Q re
+set aMICH_err APOP_f2 re
 func DARM_err = $AS_f2_I_re - 1.56140670384883
 
 lock DARM_lock $DARM_err  0.00035   0.0048
 lock CARM_lock $CARM_err  0.00039  0.00043 
 lock PRCL_lock $PRCL_err     -2.8  0.00012
-lock MICH_lock $MICH_err       43  3.9e-05 
+lock MICH_lock $aMICH_err     43  3.9e-05
+#lock MICH_lock $MICH_err       43  3.9e-05 
 noplot ITMY_lock
 func ITMY_lock = (-1.0) * $MICH_lock 
 func ETMX_lock = $CARM_lock + $MICH_lock + $DARM_lock