[Anchal, Yuta]

We are still in the progress of re-commissioning Yarm ASS.
Today, we tried to adjust output matrix by measuring the sensing matrix at DC.
Turning on yaw loops kind of works, but pitch does not. It seems like there is too much coupling in pitch to yaw.
We might need to adjust the coil output matrix of ITMY and ETMY to go further, and/or try measuring the sensing matrix including pitch - yaw coupling.

What we did:
 - Confirmed that turning on TT1 and TT2 loops (max-transmission loops) work fine. When we intentionally misalign TT1/2, the ASS loops correct it. So, we moved on to measure the sensing matrix of A2L paths, instead of using theoretical matrix caluclated from cavity geometry we used last week (40m/16909).
 - Instead of +/-1's, we put +/-2's in the ITMY coil output matrix to balance the actuation between ETMY and ITMY to take into account that ITMY is now using only two coils for actuating pitch and yaw (40m/16899).
 - Measured the change in C1:ASS-YARM_(E|I)TM_(PIT|YAW)_L_DEMOD_I_OUT16 error signals when offset was added to C1:SUS-(E|I)TMY_ASC(PIT|YAW)_OFFSET. We assumed pitch-yaw coupling is small enough here. Below was the result.

                            ETM PIT error  ITM PIT error
ETM PIT OFFSET of +100cnts: -3.0cnts       -2.99cnts
ITM PIT OFFSET of +100cnts: -11.94cnts      -5.38cnts

                            ETM PIT error  ITM PIT error
ETM YAW OFFSET of +100cnts: -3.42cnts      -16.93cnts
ITM YAW OFFSET of +10 cnts: +1.41cnts      +0.543cnts

 - Inverted the matrix to get A2L part of C1:ASS-YARM_OUT_MTRX. Attachment #1 is the current configuration so far.
 - With this, we could close all yaw loops when pitch loops were not on. But vise versa didn't work.
 - Anyway, we aligned the IFO by centering the beams on test masses by our eyes and centered all the oplevs (Attachment #2).

Next:
 - Do coil balancing to reduce pitch-yaw coupling
 - Measure sensing matrix also for pitch-yaw coupling
 - Xarm ASS is also not working now. We need to do similar steps also for Xarm

[JC, Anchal, Yuta]

We are working on resonant frequency idendification from the free swing test done last weekend.
Table below is the resonant frequencies identified, and attached are the plots of peak identification for some of our new suspensions.
To identify the resonant frequencies, the kicks were done in each degrees of freedom so that we can assume, for example, SUSPOS will be mostly excited when kicked in POS and the heighest peak is at the POS resonant frequency.
For PR3, AS1 and ETMY, the resonant frequency idendification needs to be done in the order of POS, PIT, YAW, SIDE and identified frequencies need to be removed when finding a peak.
Other than that, the identification was done without any prior assumptions on the suspensions.
For ITMY, ETMY, PR2, PR3, AS1, AS4, yaw has lower resonant frequencies than pitch, as opposed to other suspensions.
For LO1, POS and PIT frequencies might be swapped because LLCOIL is not working (40m/16898) and POS/PIT kicks both might be excited SUSPOS/PIT.
LO1 coil output matrix was temporarily modified so that we use only two coils for POS/PIT/YAW excitation (Attachment #7), as we did for ITMY (40m/16899).

The scripts for the free swinging test and analysis live in /Git/40m/scripts/SUS/InMatCalc

     POS    PIT    YAW    SIDE
BS   0.990  0.748  0.794  0.959 
ITMY 0.987  0.739  0.634  0.948 fPIT > fYAW
ETMY 0.979  0.816  0.649  0.954 fPIT > fYAW
ITMX 0.978  0.586  0.758  0.959 
ETMX 0.962  0.725  0.847  1.000 
PRM  0.939  0.541  0.742  0.990 
PR3  1.019  0.885  0.751  0.989 fPIT > fYAW
PR2  0.996  0.816  0.724  0.999 fPIT > fYAW
SRM  0.969  0.533  0.815  0.985 
SR2  0.978  0.720  0.776  0.997 
LO1  0.926  1.011  0.669  0.993 POS AND PIT MIGHT BE SWAPPED
LO2  0.964  0.998  0.995  0.990 WRONG DUE TO STUCK (40m/16913)
AS1  1.028  0.832  0.668  0.988 fPIT > fYAW
AS4  1.015  0.800  0.659  0.991 fPIT > fYAW
MC1  0.967  0.678  0.797  0.995 
MC2  0.968  0.748  0.815  0.990 
MC3  0.978  0.770  0.841  0.969 

We tracked the issue of LO1 LLCOIL not actuating LO1, and found that the DB9 cable from the coil driver to the sat amp was loose.
I tightened the screws and now it is working.
Never ever connect cables without screwing the connectors tightly!

What I did:
 - Measured the resistance and the inductance of each coil with BK PRECISION LCR meter, as I did for ITMY (Attachment #1, 40m/16896). The result is the following and it shows that LLCOIL is there.

Feedthru connector: LO1 1
Pin 3-15 / R = 16.0Ω / L = 3.27 mH (UL)
Pin 7-19 / R = 15.8Ω / L = 3.27 mH (UR)
Pin11-23 / R = 15.7Ω / L = 3.27 mH (LL)

Feedthru connector: LO1 2
Pin 3-15 / N/A
Pin 7-19 / R = 15.6Ω / L = 3.22 mH (SD)
Pin11-23 / R = 15.9Ω / L = 3.30 mH (LR)

 - Swapped the DB25 cable which goes to the feedthru LO1 1 and feedthru LO1 2. LLCOIL could be drived from LR coil driver and LRCOIL could not be drived from LL coil driver. SD and UR worked fine with the swap. This means that there is something wrong with the LL driving.
 - Went to see the rack which have coil drivers and sat amp for LO1, and immediately found that the DB9 cable was loose (Attachment #2). Tightened them and the issue was fixed.
 - C1:SUS-LO1_TO_COIL matrix gains were reverted to default values (Attachment #3).

[Anchal, Yuta]

We fixed the issue of ITMY ULCOIL not driving ITMY by replacing one of the 64pin ribbon cable in the satellite amplifier box.
We thought the coil driver and the sat amp box are OK by checking the voltage change at the output of the sat amp box by giving an offset to UL coil driver, but it was not giving a current change, probably due to too much contact resistance in the cables.
It was sneaky because it was not completely disconnected.

All the coils for our suspensions are now working!

What we did:
 - Using breakout boards, the output current of sat amp box was measured using FLUKE multimeter. It turned out that UL is not giving measurable current. We also confirmed that UR coil driver can drive UL by re-directing the current from UR coil driver to UL. This means that the UL magnet was not de-magnetized!
 - Measured the coil resistance from at the coil driver output and found that UL coil seen from there has too high resistance which cannot be measured with the multimeter, whereas UR coil was measured to be ~30 Ohms.
 - Went back to the feedthru and measured the resistance of UL coil. Upto the output of the Satellite Amp Terimator, the resistance was measured to be ~16 Ohms, but not at the input of the Satellite Amp Terimator (Attachment #1,2).
 - It turned out that #16 pin of 64pin ribbon cable in between the Satellite Amp Terimator (LIGO-D990021) and the Satellite Amp board (LIGO-D961289) at the Satellite Amp Terimator side was not good (Attachment #3).
 - Replaced the cable and confirmed that ULCOIL can kick ITMY (Attachment #4).
 - C1:SUS-ITMY_TO_COIL matrix was reverted to default values.

Next:
 - We might have to re-commission Yarm ASS again since pitch-yaw coupling have changed. -> EDIT: Checked that it works (except for ITM PIT L), including offloading offsets (writeASS_offsets.py), 18:30 local.
 - Now that LO1 LLCOIL issue is solved and LO2 stuck is solved, we should do the free swing test again to identify the resonant frequencies.
 - OSEM sensor diagonalization (input matrix), coil balancing (and F2A)

[Paco, Tomislav, Yuta]

Somehow, when we were trying to measure WFS open loop transfer functions, PMC unlocked many times for the past two hours and PMC transmission got low.
PMC iput beam was aligned again, and IMC WFS DC offsets and RF offsets were adjusted.
PMC transmission is now C1:PSL-PMC_PMCTRANSPD~0.75, and IMC transmission is C1:IOO-MC_TRANS_SUM~1.4e4.
Actually, IMC transmission once reached 1.5e4 at 06-16-2022 20:01 UTC with PMC transmission of 0.75 (see Attached). There might be a better alignment.

Following what we have done in 2013 (40m/8242), actuator calibration was done using MICH.

AS55_Q in MICH : 9.74e8 counts/m
BS   : 26.08e-9 /f^2 m/counts
ITMX : 5.29e-9 /f^2 m/counts
ITMY : 4.74e-9 /f^2 m/counts

Optical gain is 25% lower than the measurement in June 6 (40m/16892), probably because our estimate was too rough then and also we now have ~15% lower IMC transmission.
Actuator gains are 2-30% higher than the measurement in 2013.

MICH error signal calibration:
 C1:LSC-AS55_Q_ERR was calibrated by taking data with C1:LSC-ASDC_OUT, when Michelson was aligned and free swinging (Attachment #1).
 AS55_Q and ASDC were X-Y plotted and fitted with ellipse to get an amplitude of AS55_Q to be 82.51 counts (Attachment #2).
 4*pi*A/lambda gives you 9.74e8 counts/m, where meters are in terms of difference between BS to ITMX length and BS to ITMY length.
 Jupyter notebook: https://git.ligo.org/40m/scripts/-/blob/main/CAL/MICH/MICHOpticalGainCalibration.ipynb

Openloop transfer function for actuator calibration:
 C1:LSC-MICH_GAIN was lowered to -1 (instead of -6), and some of filters are turned off to make the MICH UGF to be ~10.
 Also, ellip("LowPass",4,1,40,50) was added to C1:LSC-MICH_A filter bank to cut the feedback above 50 Hz, so that the loop does not suppress the measurement.
 The configuration is in Attachment #3.

Actuator calibration of BS, ITMX, ITMY:
 With this MICH OLG, transfer functions from C1:LSC-BS,ITMX,ITMY_EXC to C1:LSC-AS55_Q_ERR were measured.
 AS55_Q was calibrated to meters using the calibration factor above, and fitted the transfer function with 1/f^2 in 70-150 Hz range to get the actuator efficiency mentioned above (Attachement #4).
 Thus, meters in this calibration is in terms of ITM POS motion (not in BS POS motion).
 Jupyter notebook: https://git.ligo.org/40m/scripts/-/blob/main/CAL/MICH/MICHActuatorCalibration.ipynb

MICH displacement noise:
 Measured values were added to c1cal model as follows.
  C1:CAL-MICH_CINV FM2: 1/9.74e8 = 1.03e-9
  C1:CAL-MICH_A FM2: 2.608e-8 (it was 2.07e-8 from 2013!)
  C1:CAL-MICH_A_GAIN = 0.5 to take into account of C1:LSC-OUTPUT_MTRX_8_2=0.5 in the LSC output matrix for BS
 Spectrum of C1:CAL-MICH_W_OUT (now calibrated in nm) with configuration in Attachment #5 was taken.
 Attachement #6 is the result. I also took the spectrum with PSL shutter off to measure the sensing noise. The sensing noise limits our sensitivity above ~40 Hz at 5e-11 m/rtHz.

After discussing with Anchal, we decided to route BHD related PD signals directly to ADC of c1sus2, which handles our new suspensions including LO1, LO2, AS1, AS4, so that we can control them directly.
BHD related PD signals will be sent to c1lsc for DARM control.

Re-cabling was done, and now they are online at C1:X07-MADC1_EPICS_CH16 (DC PD A) and CH17 (DC PD B) with 15ft DB9 cable.
Here, DC PD A is the transmission of BHD BS for AS beam, and DC PD B is the reflection of BHD BS for AS beam (see attached photo).

[Anchal, Yuta]

RTS models for BHD homodyne phase control (c1hpc) and angular control (c1bac) are created and added to c1sus2.
c1su2 and c1lsc models were modified accordingly.
We still have issues with IPC PCIE connection sending DCPD A and B signals to c1lsc and DC error 0x2000 in c1su2 model.

c1hpc (host: c1sus2) Attachment #1
 This model is for homodyne phase control.
 It can dither LO1, LO2, AS1, AS4 in POS and demodulate mixture of DCPD A/B signals for the phase control to feedback to those optics.
 It also sends DCPD A/B signals to c1lsc via cdsIPCx_PCIE.
 Dither and controls signals are sent to the optics via cdsIPCx_SHMEM.

c1bac (host: c1sus2)
 This model is for BHD angular control.
 It is basically the same as c1hpc, but it is for PIT and YAW dithering of LO1, LO2, AS1, AS4.

c1su2 (host: c1su2) Attachment #2
 LSC and ASCPIT/YAW feedback signals from c1hpc and c1bac via shared memory were added to send them to corresponding optics.
 Somehow Mux/Demux didn't work to send SHMEM signals inside the subsystem in the Simulink model (this works for ADC, but probably not for IPC stuff?), and we had hard time make-install-ing this model.

c1lsc (host: c1lsc) Attachment #3
 DCPD A/B signals from c1hpc via PCIE were added for our new error signals for LSC.

Starting and restarting the models
 After having some troubles make-install-ing modified models (be careful of goto and from tags!), we stopped all the models in c1sus, c1ioo, c1lsc, c1sus2 and started all of them, including new c1hpc and c1bac models.
 This somehow created RFM errors in c1scx and c1scy.
 So, we proceeded to do the same step we did in 40m/16887 and 40m/15646, now including c1sus2 for the restart.
 Initial attempt made c1lsc, c1sus, c1ioo mostly red, so scripts/cds/rebootC1LSC.sh was run again on pianosa.
 RFM issues for c1scx and c1scy were solved.
 Shared memory within c1sus2 seems to be working, but sending DCPD A/B signals from c1hpc to c1lsc is not working (see Attachement #4).

Next:
 - Fix C1:HPC-LSC_DCPC_A/B issue
 - Make/modify MEDM screens

[Koji, Yuta]

I found that Yarm cannot be locked today. Both POY11 and POYDC were not there when Yarm was aligned, and ITMY needed to be highly misaligned to get POYDC.
POY beam also could not be found at ITMY table.
Koji suggested to use AS55 instead to lock Yarm. We did it (AS55_I_ERR, C1:LSC-YARM_GAIN=-0.002) and manually ASS-ed to get Yarm aligned (ASS with AS55 somehow didn't work).
After that, we checked ITMY table and found that POY beam was clipped at an iris which was closed!
I opened it and now Yarm locks with POY11 again. ASS works.
PMC was also aligned.

C1:PSL-PMC_PMCTRANSPD ~0.74
C1:IOO-MC_TRANS_SUM ~14000
C1:LSC-TRY_OUT ~0.7
C1:LSC-TRX_OUT ~0.8

[Anchal, Yuta]

We found that MC1 local damping loop signs were revereted to the state before our standardization on June 7th (40m/16898), but the WFS output matrix was not reverted.
This caused the sign flip in the feedback to MC1, which caused the IMC WFS issue.
This probably happened when we were restarting the models after RTS modeling (40m/16935). We might have used wrong snap files for burt-restoring.

We went back to the snapshot taken at 09:19 June 21, 2022 and now the IMC WFS is working,

[Yehonathan, Yuta]
EDITED by YM on 22:11 June 27, 2022 to correct for a factor of two in the modulation index

Since we have measured optical gain in MICH to be an order of magnitude less compared with Yehonathan's FINESSE model (40m/16923), we measured the power at AS55 RF PD, and measured the modulation depths using Yarm cavity scan.
We found that 50/50 beam splitter which splits AS55 path into RF PD and RF QPD was not included in the FINESSE model. Measured modulation index were as follows:

TEM00 peak height: 0.6226 +/- 0.0237
RF11 peak height: 0.0067 +/- 0.0007
RF55 peak height: 0.0081 +/- 0.0014
RF11 modulation index: 0.208 +/- 0.012
RF55 modulation index: 0.229 +/- 0.020
RF11 modulation index: 0.104 +/- 0.006
RF55 modulation index: 0.114 +/- 0.010

Here, modulation depth m is defined in E=E_0*exp(i*(w*t+m*sin(w_m*t))), and m m/2 equals to square of the intensity ratio between sidebands and TEM00.

Power measurement at AS55 RF PD:
 - ITMY and ITMX single bounce reflection was measured to be 50-60 uW at the front of AS55 RFPD.
 - In the FINESSE model, it was expected to be ~110 uW with 0.8 W input to PRM (0.8 W * 5%(PRM) * 50%(BS) * 50%(BS) * 10%(SRM) * 10%(AS2) gives 100 uW)
 - In AP table, AS55 beam was split into two paths with 50/50 beam splitter, one for AS55 RF PD and one for AS WFS and AS110. This will be included in the FINESSE model.

Modulation depth measurement using Yarm cavity scan:
 - Aligned Yarm using ASS, and unlocked Yarm to get the 2sec scan data of C1:LSC-TRY_OUT_DQ, C1:LSC-POY11_I_ERR_DQ, C1:LSC-AS55_I_ERR_DQ.
 - TRY data was used to get TEM00 peak heights
 - POY11/AS55 data was used to find RF11/RF55 sideband peaks, and height was measured at TRY (see attached).
 - If we define m to be E=E_0*exp(i*(w*t+m*sin(w_m*t))), the amplitude of TEM00 I_00 is proportional to J_0(m) and the amplitude of upper/lower sideband I_f1 is proportional to J_1(m), where J_n(m) is the bessel function of the first kind.
 - m can be calculated using 2*sqrt(I_f1 / I_00).
 - Results were shown above. Error is calculated from the standard deviation of multiple measurements with multiple peaks,
 - The code for doing this lives in https://git.ligo.org/40m/measurements/-/blob/main/LSC/YARM/modulationIndex.ipynb

Discussion:
 - Power at AS55 account for the factor of 2, In the FINESSE model, modulation index of 0.3 was used (could be m=0.3/2 or m=0.3; needs check). These combined can explain a factor of 3 at least (or 6).
 - Gautam's measurement in Jan 2021 (40m/15769) gives almost double modulation index, but I'm not sure what is the definition Gautam used. It agrees with Gautam's measurement in Jan 2021.

[JC, Yuta]

To circumvent IPC error sending BHD DC PD signals from c1sus2 to c1lsc, DB9 cable from BHD DC PD box sent to c1sus2 is now split and sent also to c1lsc.
They are now available in both

c1sus2 ADC1
C1:X07-MADC1_EPICS_CH16 (DC PD A) and CH17 (DC PD B)

c1lsc ADC1
C1:X04-MADC1_EPICS_CH4 (DC PD A) and CH5 (DC PD B)

Next:
 - Add battery powered SR560 to decouple c1sus2 and c1lsc to avoid the ground loop

[Anchal, Paco, Yuta]

We tried to lock FPMI with REFL55 and AS55 this week, but no success yet.
FPMI locks with POX11, POY11 and ASDC for MICH stably, but handing over to 55's couldn't be done yet.

What we did:
 - REFL55: Increased the whitening gain to 24dB. Demodulation phase tuned to minimize MICH signal in I when both arms are locked with POX and POY. REFL55 is noisier than AS55. Demodulation phase and amplitude of the signal seem to drift a lot also. Might need investigation.
 - AS55: Demodulation phase tuned to minimize MICH signal in I when both arms are locked with POX and POY. Whitening gain is 24dB.
 - Script for demodulation phase tuning lives in https://git.ligo.org/40m/scripts/-/blob/main/RFPD/getPhaseAngle.py
 - Locking MICH with REFL55 Q: Kicks BS much and not so stable probably because of noisy REFL55. Offtet also needs to be adjusted to lock MICH to dark fringe.
 - BS coil balancing: When MICH is "locked" with REFL55 Q, TRX drops rapidly and AS fringe gets worse, indicating BS coil balancing is not good. We balanced the coils by dithering POS with different coil output matrix gains to minimize oplev PIT and YAW output manually using LOCKINs.
 - Locking MICH with ASDC: Works nicely. Offset is set to -0.1 in MICH filter and reduced to -0.03 after lock acquisition.
 - ETMX/ETMY actuation balancing: We found that feedback signal to ETMX and ETMY at LSC output is unbalanced when locking with POX and POY. We dithered MC2 at 71 Hz, and checked feedback signals when Xarm/Yarm are locked to find out actuation efficiency imbalance. A gain of 2.9874 is put into C1:LSC-ETMX filter to balance ETMX/ETMY. I think we need to check this factor carefully again.
 - TRX and TRY: We normalized TRX and TRY to give 1 when arms are aligned. Before doing this, we also checked the alignment of TRX and TRY DC PDs (also reduced green scattering for TRY). Together with ETMX/ETMY balancing, this helped making filter gains the same for POX and POY lock to be 0.02 (See, also 40m/16888).
 - Single arm with REFL55/AS55: We checked that single arm locking with both REFL55_I and AS55_Q works. Single arm locking feeding back to MC2 also worked.
 - Handing over to REFL55/AS55: After locking Xarm and Yarm using POX to ETMX and POY to ETMY, MICH is locked with ASDC to BS. Handing over to REFL55_I for CARM using ETMX+ETMY and AS55_Q for DARM using -ETMX+ETMY was not successful. Changing an actuator for CARM to MC2 also didn't work. There might be an unstable point when turning off XARM/YARM filter modules and switching on DARM/CARM filter modules with a ramp time. We also need to re-investigate correct gains and signs for DARM and CARM. (Right now, gains are 0.02 for POX and POY, -0.02 for DARM with AS55_Q (-ETMX+ETMY), -0.02 for CARM with REFL55_I with MC2 are the best we found so far)
 
Next:
 - Measure ETMX and ETMY actuation efficiencies with Xarm/Yarm to balance the output matrix for DARM.
 - Measure optical gains of POX11, POY11, AS55 and REFL55 when FPMI is locked with POX/POY/ASDC to find out correct filter gains for them.
 - Make sure to measure OLTFs when doing above to correct for loop gains.
 - Lock CARM with POY11 to MC2, DARM with POX11 to ETMX. Use input matrix to hand over instead of changing filter modules from XARM/YARM to DARM/CARM.
 - Try using ALS to lock FPMI.

(This is a complete restore of elog 40m/16970 from July 5, 2022 at 14:34)

ETMX and ETMY actuators were calibrated using single arm lock by taking the actuation efficiency ratio between ITMs. Below is the result.

ETMX :  2.65e-9 /f^2 m/counts (0.5007 times ITMX)
ETMY : 10.91e-9 /f^2 m/counts (2.3017 times ITMY)

Motivation:
- ETMX and ETMY actuators seemed to be unbalanced when locking DARM (see 40m/16968)

What we did:
- Reverted to C1:LSC-ETMX_GAIN = 1
- XARM was locked using POX11_I_ERR (42dB whitening gain, 132.95 deg for demod phase) with ETMX and C1:LSC-XARM_GAIN=0.06
- YARM was locked using POY11_I_ERR (18dB whitening gain, -66.00 deg for demod phase) with ETMX and C1:LSC-YARM_GAIN=0.02
- OLTFs for each was measured to be Attachment #1; UGF was ~180 Hz for XARM, ~200 Hz for YARM.
- Measured TF from C1:LSC-(E|I)TM(X|Y)_EXC to C1:LSC-(X|Y)ARM_IN1 (see Attachment #2)
- Took the ratio between ITM actuation and ETM actuation to calculate ETM actuation. For ITM actuation, we used the value measured using MICH (see 40m/16929). The average of the ratio in the frequency range 70-150 Hz was used.

Files:
- Measurement files live in https://git.ligo.org/40m/measurements/-/tree/main/LSC/XARM and YARM
- Script for calculation lives at https://git.ligo.org/40m/scripts/-/blob/main/CAL/ARM/ETMActuatorCalibration.ipynb

Discussion:
- ETMX actuation is 4.12 times less compared with ETMY. This is more or less consistent with what we measured in 40m/16968, but we didn't do loop-correction at that time.
- We should check if this imbalance is as expected or not.

Summary of actuation calibration so far:
BS   : 26.08e-9 /f^2 m/counts (see 40m/16929)
ITMX :  5.29e-9 /f^2 m/counts (see 40m/16929)
ITMY :  4.74e-9 /f^2 m/counts (see 40m/16929)
ETMX :  2.65e-9 /f^2 m/counts (0.5007 times ITMX)
ETMY : 10.91e-9 /f^2 m/counts (2.3017 times ITMY)
 

(This is also a restore of elog 40m/16971 from Jul 5, 2022 at 17:36)

MC2 actuator calibration was also done using Yarm in the same way as we did in 40m/16970 (now 40m/16977).
The result is the following;
MC2 : -14.17e-9 /f^2 m/counts in arm length (-2.9905 times ITMY)
MC2 :   5.06e-9 /f^2 m/counts in IMC length
MC2 :  1.06e+05 /f^2 Hz/counts in IR laser frequency

What we did:
- Measured TF from C1:LSC-MC2_EXC to C1:LSC-YARM_IN1 during YARM lock using ETMY (see Attachment #1). Note that the sign of MC2 actuation and ITMY actuation is flipped.
- Took the ratio between ITM actuation and MC2 actuation to calculate MC2 actuation. For ITM actuation, we used the value measured using MICH (see 40m/16929). The average of the ratio in the frequency range 70-150 Hz was used (see Attachment #2).
- The actuation efficiency in meters in arm length was converted into meters in IMC length by multiplying it by IMCLength/ArmLength, where IMCLength=13.5 m is half of IMC round-trip length, ArmLength=37.79 m is the arm length.
- The actuation efficiency in meters in arm length was converted into Hz in IR laser frequency by multiplying it by LaserFreq/ArmLength, where LaserFreq=1064 nm / c is the laser frequency.

Files:
- Measurement files live in https://git.ligo.org/40m/measurements/-/tree/main/LSC/YARM
- Script for calculation lives at https://git.ligo.org/40m/scripts/-/blob/main/CAL/ARM/ETMActuatorCalibration.ipynb

Summary of actuation calibration so far:
BS   : 26.08e-9 /f^2 m/counts (see 40m/16929)
ITMX :  5.29e-9 /f^2 m/counts (see 40m/16929)
ITMY :  4.74e-9 /f^2 m/counts (see 40m/16929)
ETMX :  2.65e-9 /f^2 m/counts (0.5007 times ITMX)
ETMY : 10.91e-9 /f^2 m/counts (2.3017 times ITMY)
MC2 : -14.17e-9 /f^2 m/counts in arm length (-2.9905 times ITMY)
MC2 :   5.06e-9 /f^2 m/counts in IMC length
 

NOTE ADDED by YM on July 7, 2022

To account for the gain imbalance in ETMX, ETMY, MC2, LSC violin filter gains were set to:
C1:LSC-ETMX_GAIN = 4.12
C1:LSC-MC2_GAIN = -0.77
This is a temporary solution to make ETMX and MC2 actuation efficiencies from LSC in terms of arm length to be the same as ETMY 10.91e-9 /f^2 m/counts.

I think it is better to make C1:LSC-ETMX_GAIN = 1, and put 4.12 in C1:SUS-ETMX_TO_COIL gains. We need to adjust local damping gains and XARM ASS afterwards.
As for MC2, it is better to put -0.77 in LSC output matrix, since this balancing depends on LSC topology.

[Paco, Koji, Yuta]

We managed to lock MICH using REFL55_Q by setting the demodulation phases and offsets right.
The following is the current FPMI locking configuration we achieved so far.

DARM: POX11_I / gain 0.007 / 0.5*ETMX-0.5*ETMY (or 1*ETMX) / UGF of ~100 Hz
CARM: POY11_I / gain 0.018 / 1*MC2 / UGF of ~200 Hz
MICH: REFL55_Q / gain -10 / 0.5*BS / UGF of ~30 Hz

Transitioning DARM error signal from POX11_I to 0.5*POX11_I+0.5*POY11_I was possible with FM4 filter off in DARM filter bank, but not to AS55_Q yet.

REFL55 and AS55 demodulation phase tuning:
 - We found that both AS55 and REFL55 are contaminated by large non-MICH signal, by making a ASDC vs RF plot (see 40m/16929).
 - After both arms are locked with POX and POY, MICH was locked with AS55_Q. ASDC was minimized by putting an offset to MICH filter.
 - With this, REFL55 offsets were zeroed and demodulation phase was tuned to minimize REFL55_Q.
 - Locked MICH with REFL55_Q, and did the same thing for AS55_Q.
 - Resulting ASDC vs RF plots were attached. REFL55_Q now looks great, but REFL55_I and AS55 are noisy (due to signals from the arms?).

Jupyter notebook: https://git.ligo.org/40m/scripts/-/blob/main/CAL/MICH/MICHOpticalGainCalibration.ipynb

Sensing matrix:
 - With FPMI locked using POX/POY, DARM and CARM lines were injected at around 300 Hz to measure the sensing gains. For line injection, C1:CAL-SENSMAT was used, but for the demodulation we used a script. The following is the result.

 Sensors              DARM (ETMX)         CARM (MC2)        
C1:LSC-AS55_I_ERR    3.10e+00 (-34.1143 deg)    1.09e+01 (-14.907 deg)    
C1:LSC-AS55_Q_ERR    9.96e-01 (-33.9848 deg)    3.30e+00 (-27.9468 deg)    
C1:LSC-REFL55_I_ERR    6.75e+00 (-33.7723 deg)    2.92e+01 (-34.0958 deg)    
C1:LSC-REFL55_Q_ERR    7.07e-01 (-33.4296 deg)    3.08e+00 (-33.4437 deg)    
C1:LSC-POX11_I_ERR    3.97e+00 (-33.9164 deg)    1.51e+01 (-30.7586 deg)    
C1:LSC-POY11_I_ERR    6.25e-02 (-20.3946 deg)    3.59e+00 (38.4207 deg)

Jupyter notebook: https://git.ligo.org/40m/scripts/-/blob/main/CAL/SensingMatrix/MeasureSensMat.ipynb

 - By taking the ratios of POX11_I and AS55_Q for DARM, POY11_I and REFL55_I for CARM, we tried to find the correct gains for REFL55 and AS55 for DARM and CARM. x3.96 more gain for AS55_Q than POX11_I and x0.123 less gain for REFL55_I than POY11_I.

Next:
 - Try locking the arms with no triggering, and then try locking FPMI with REFL/AS without triggering. No FM4 for this, since FM4 kills gain margin.
 - Lock single arm with AS55_Q and make a noise budget. Make sure to misalign ITMX(Y) completely when locking Y(X)arm.
 - Lock single arm with REFL55_I and make a noise budget.
 - Repeat Xarm noise budget with Yarm locked with POY11_I and MC2 (40m/16975).
 - Check IMC to reduce frequency noise (40m/17001)

To balance the actuation on ETMX and ETMY, x4.12 was aded to C1:SUS-ETMX_(UL|UR|LR|LL|SD)COIL FM1. OSEM damping filter gains, oplev loop gains, and alignment offsets were divided by this factor.
C1:LSC-ETMX_GAIN is now 1.

To do:
 - Balance ETM and ITM. It should make ASS more sensible.
 - Re-commission Xarm ASS and Yarm ASS.

[Paco, Anchal-remote-support, Yuta]

We added fast channels to BHD DC PDs.
C1:LSC-DCPD_(A|B)_IN1 are now available, but C1:HPC-DCPD_(A|B)_IN1 still gives us zero.

c1hpc situation -> not good
 - We can see the slow signal at C1:X07-MADC1_EPICS_CH16 (DC PD A) and CH17 (DC PD B)
 - C1:HPC-DCPD_(A|B)_IN1 is there, but zero.
 - We have modified c1hpc model to add DCPD_(A|B) filters in front of the input matrix (see Attachment #1).
 - After modifying the model, we run
ssh c1sus2
rtcds make c1hpc
rtcds install c1hpc
ssh fb1
sudo systemctl restart daqd_*

 - After this, we got 0x2000 error. So, we ran the following. This removed 0x2000 error, but DCPD signals are still zero. They are also not available in C1HPC-MONITOR_ADC1.adl screen (see Attachment #3).
ssh c1sus2
rtcds restart c1hpc

c1lsc situation -> good
 - We could see the slow signal at C1:X04-MADC1_EPICS_CH4 (DC PD A) and CH5 (DC PD B), and also C1:LSC-DCPD_(A|B)_NORM after making C1:LSC-DCPD_(A|B)_POW_NORM=1. The ADC channel and DCPD channel are exactly the same.
 - After confirming the above, we modified the c1lsc model to add DCPD_(A|B) filters in front of the input matrix (see Attachment #2).
 - After modifying the model, we run
ssh c1lsc
rtcds make c1lsc
rtcds install c1lsc
ssh fb1
sudo systemctl restart daqd_*

 - After this, we also got 0x2000 error. We also noticed that, for example, C1:X04-MADC0_EPICS_CH31 and C1:LSC-ASDC_INMON are different, which used to be the same (ASDC_INMON was largely attenuated).
 - In the end, we run the following to remove 0x2000 error, but it crashed c1lsc, as well as c1sus, c1ioo.
ssh c1lsc
rtcds restart c1lsc

 - So, we did rebootC1LSC.sh. This made c1lsc, c1ioo and c1sus as green as before, except for RFM issue in TRX/TRY, like we saw in June. We followed the steps in 40m/16887 to hard reboot c1iscex/c1iscey and ran rebootC1LSC.sh again. This made C1CDS_FE_STATUS.adl screen as green as before (see Attachment #3).

 - Fast channels C1:LSC-DCPD_(A|B)_IN1 are now available. They are also available in C1LSC-MONITOR_ADC1.adl screen (see Attachment #3).

[Paco, Yuta]

We measured contrast of Michelson fringe in both arms locked and mis-aligned. It was around 90%.
We also measured the contrast of ITM single bounce vs LO beam using BHD DC PDs. It was around 43%.
The measurement depends on the alignment and how to measure the maximum and minimum of the fringe. ITM-LO fringe was also not stable because motions of AS/LO mirrors are large. More tuning necessary.

Background
 - As measured in elog 40m/17012, we see a lot of CARM in AS, which indicates large contrast defect.
 - We want to check mode-matching of LO beam to AS beam.

BHD DC PD conditioning
 - We added DCPD_A and DCPD_B to /opt/rtcds/caltech/c1/scripts/LSC/LSCoffsets3 script, which zeros the offsets when shutters are closed.
 - We also set C1:LSC-DCPD_(A|B)_GAIN = -1 since they are inverted.

Contrast measurement
 - Contrast was measured using channels ['C1:LSC-ASDC_OUT','C1:LSC-POPDC_OUT','C1:LSC-REFLDC_OUT','C1:LSC-DCPD_A_OUT','C1:LSC-DCPD_B_OUT']. For LO, only DCPD_(A|B) are used.
 - We took 15%-percentile (40% for ITM-LO fringe) from the maximum and minimum of the data, and took the median to estimate the maximum value and the minimum value (see Attachment).
 - Contrast = (Imax - Imin) / (Imax + Imin)
 - We measured three times in each configuration to estimate the standard error.
 - Jupyter notebook: https://git.ligo.org/40m/scripts/-/blob/main/CAL/BHD/measureContrast.ipynb

Results
Both arms locked, MICH fringe (15% percentile)
Contrast measured by C1:LSC-ASDC_OUT is 89.75 +/- 0.17 %
Contrast measured by C1:LSC-POPDC_OUT is 79.41 +/- 0.86 %
Contrast measured by C1:LSC-REFLDC_OUT is 97.34 +/- 0.34 %
Contrast measured by C1:LSC-DCPD_A_OUT is 95.41 +/- 1.55 %
Contrast measured by C1:LSC-DCPD_B_OUT is 89.76 +/- 1.49 %
Contrast measured by all is 90.34 +/- 1.68 %

Both arms mis-aligned, MICH fringe (15% percentile)
Contrast measured by C1:LSC-ASDC_OUT is 89.32 +/- 0.57 %
Contrast measured by C1:LSC-POPDC_OUT is 94.55 +/- 0.62 %
Contrast measured by C1:LSC-REFLDC_OUT is 97.95 +/- 1.37 %
Contrast measured by C1:LSC-DCPD_A_OUT is 96.40 +/- 1.04 %
Contrast measured by C1:LSC-DCPD_B_OUT is 90.98 +/- 1.07 %
Contrast measured by all is 93.84 +/- 0.94 %

ITMY-LO fringe (40% percentile)
Contrast measured by C1:LSC-DCPD_A_OUT is 45.51 +/- 0.45 %
Contrast measured by C1:LSC-DCPD_B_OUT is 38.69 +/- 0.43 %
Contrast measured by all is 42.10 +/- 1.03 %

ITMX-LO fringe (40% percentile)
Contrast measured by C1:LSC-DCPD_A_OUT is 46.65 +/- 0.65 %
Contrast measured by C1:LSC-DCPD_B_OUT is 39.82 +/- 0.51 %
Contrast measured by all is 43.24 +/- 1.45 %

Discussion
 - As you can see from the attachment, REFLDC is noisy and over estimating the contrast. ASDC is reliable. We need to tune the threshold to measure the maximum value and minimum value. We should also use the mode instead of median.
 - Contrast depends very much on the alignment. We didn't tweak too much today.
 - ITM-LO fringe was not stable, probably due to too much motion in AS1, AS4, LO1, LO2. Their damping needs to be re-tuned.

Next:
 - Model FPMI sensing matrix with measured contrast defect
 - Estimate AS-LO mode-mismatch using the measured contrast
 - Lock ITM-LO fringe using DCPD_(A|B) as error signal, and ITM or LO1/2 as actuator
 - Lock MICH with DCPD_(A|B), and with LO beam
 - Get better contrast data with better alignment and better AS1, AS4, LO1, LO2 damping

[Tega, Yuta]

We have added C1:HPC-DCPD_A_ERR and C1:HPC-DCPD_B_ERR testpoints, which can be used as A+B, A-B etc.
Restarting c1hpc crashed c1sus2, and also made c1lsc/ioo/sus models red.
We run /opt/rtcds/caltech/c1/Git/40m/scripts/cds/restartAllModels.sh to restart all the machines. It worked perfectly without manually pressing power buttons! Wow!

We have also edited /opt/rtcds/caltech/c1/medm/c1su2/C1SU2_WATCHDOGS.adl so that it will use new /opt/rtcds/caltech/c1/Git/40m/scripts/SUS/medm/resetFromWatchdogTrip.sh instead of old /opt/rtcds/caltech/c1/scripts/SUS/damprestore.py.

[Paco, Yehonathan, Yuta]

We measured power and counts at BHD DCPDs with LO beam only and ITM single bounce.
We found that LO beam power is ~7 times less than the expected.
We also confirmed that AS beam is clipped somewhere inside vacuum and have 20-50% less power compared with the expected.
LO/AS beams going to DCPD A and B also have power imbalance by 30-40%.

What we did:
 - Run LSCoffsets.py to zero the offsets. I modified the old script so that it can handle new BHD PDs. Also, a bug was fixed (it didn't take into account the gains in filer modules, so INMON is now used instead of OUT16 for calculating offsets).
https://git.ligo.org/40m/scripts/-/blob/main/RFPD/LSCoffsets.py
 - Measured powers and counts in BHD DCPDs at ITMY table with LO beam only and ITMX/ITMY single bounce.
 - During the measurement, we found that power into DCPD A and DCPD B are quite different. One of the reason was a lens and an iris right after the viewport for A path. We removed both of them. Also, only A path have a pickoff which picks off ~20% of light to BHD camera (called SRMF; 40m/16880).
 - We also found that LO beam shape is ugly. ITM single bounce beam from X and Y have similar clipping (see Attached photos). We tried to reduce clipping with various suspensions (LO1, LO2, AS1, AS4, SR2, SRM, BS, ITMX, PR2), but was not possible by moving only single suspension.

Result:
 - Result of counts and power measurements are as follows. Power was measured right in front of DCPD, and also right after the viewport to estimate the loss in the in-air paths. Note that LSC channels have gain of 1, but HPC channels have gain of -1.9 for DCPD_A and -1 for DCPD_B.

                       Blocked  LO       ITMX      ITMY    
C1:LSC-DCPD_A_OUT16    -0.01    -17.89    -91.62    -86.21    
C1:LSC-DCPD_B_OUT16    +0.06    -17.72   -131.83   -131.98    
C1:HPC-DCPD_A_OUT16    +0.07    +34.12   +174.63   +164.24    
C1:HPC-DCPD_B_OUT16    +0.13    +17.60   +131.31   +131.49    
Power at DCPD_A        19 uW    63 uW    278 uW    290 uW    
Power at DCPD_B        19 uW    65 uW    393 uW    404 uW    
Power at viewport A    -- uW    82 uW    350 uW    337 uW    
Power at viewport B    -- uW    64 uW    436 uW    431 uW

DCPD calibration:
 - From the measurements above, counts/W in IN1 can be calculated as follows. Offset of 19 uW is substracted from the measured power to take into account for background light.

C1:LSC-DCPD_A_IN1     -3.59e+05 counts/W
C1:LSC-DCPD_B_IN1     -3.61e+05 counts/W
C1:HPC-DCPD_A_IN1     -3.60e+05 counts/W
C1:HPC-DCPD_B_IN1     -3.57e+05 counts/W

Discussion:
 - DCPD calibration shows that DCPD to ADC itself is quite balanced within 1%. A factor of 1.8-1.9 seen was from unbalanced light between A path and B path (40m/17037).
 - Power expected for ITM single bounce to one of DCPDs is ~520 uW, but was 350-430 uW as measured right after the viewport. Power at A is significantly less than that for B. Note that power at AS55 was as expected (40m/16952). Also, clipping cannot be reduced by moving suspensions. These could mean that clipping is happening after AS2.

950 mW * 0.9 (IMC transmission?) * 5.637%(PRM) * 97.8%(PR2) * 50%(BS) * 98.6%(ITM) * 50%(BS) * 10%(SRM) * 90%(AS2) * 50%(BHDBS) = 520 uW

 - Power expected for LO beam to one of DCPDs is ~530 uW, but was 60-80 uW as measured right after the viewport. Power at A is significantly more than that for B, which is opposite for ITM single bounce. This could mean that something is happening at BHDBS? We are not sure why the power is so low. Are we seeing some ghost beam? For PR2 transmission, 22000 ppm was used for calculation, from 40m/16541.

950 mW * 0.9 (IMC transmission?) * 5.637%(PRM) * 2.2%(PR2) * 50%(BHDBS) = 530 uW

 - As far as we remember, beam shapes were not as bad when we closed out the chambers...

Next:
 - Check if measured power explains the visitivity of LO-ITM single bounce (40m/17020)
 - If not, what is the mode mismatch? Is it possible to explain the mode mismatch with deviations from designed mode-matching telescope?
 - Measure POP power to see if PR2 actually have T=2.2%
 - Play with LO1 and LO2 to invesitate LO beam shape and power
 - Check coherence between LO/AS power fluctuations with suspension motions
 - What is the expected counts/W for these DCPDs?
 - Balance the optical paths in ITMX table for A and B (same lenses, same mirrors)
 - Install better lens in front of camera

Last week, we could find an alignment which realizes LO beam and AS beam both unclipped, but it was not consistent with an alignment which realize BHD fringe (40m/17046).
Today, we tweaked the alignment of SR2, AS1, AS4 to have BHD fringe with reduced LO and AS beam clipping.
AS beams on AP table and BHD both still look clipped, but much better now.
Ideally, SR2 and AS1 will unclip AS beam, and LO1, LO2, AS4 would make BHD fringe, but it is hard right now since LO beam seem to have little room and LO2 have little actuation range.
BHD optics on ITMY table, including camera, and AS55/ASDC were realigned after the aglinment work (Note that DCPD_A path have a pick-off for camera path, and this pick-off mirror have quite significant incident angle dependence of R/T ratio).

Current alignment scheme:
Current alignment scheme I figured out is the following.
 - Check Y green. If it is transmitted at good spot on GTRY camera, Yarm is OK. If not, tweak ITMY/ETMY. alignment.
 - Mis-align AS4, align TT1, TT2, LO1 to have DCPD_A_OUT of ~130 and DCPD_B_OUT of ~125.
 - Align PR3, PR2 to maximize TRY_OUT to ~1.05.
 - Tweak ITMY/ETMY if the beam spot on them are not good.
 - Align BS, ITMX to have good MICH fringe and TRX_OUT to ~1.1.
 - Tweak ITMX/ETMX if the beam spot on them are not good.
 - Misalign ETMY, ETMX, ITMY to have LO-ITMX fringe in BHD DCPDs, and align AS beam with SR2 and AS4 differentially, with ratio of AS4/SR2=3.6.

DC PD values in various configurations:
Both arms locked with POX/POY, MICH free, PRM/SRM misaligned
                          Mean     Max      Min
C1:IOO-MC_TRANS_SUM :     14088.57 13947.52 14167.04
C1:LSC-ASDC_OUT16 :           0.16    -0.02     0.34
C1:LSC-POPDC_OUT16 :        369.34   -74.88   854.34
C1:LSC-REFLDC_OUT16 :         0.03    -0.00     0.06
C1:LSC-TRY_OUT16 :            1.00     0.95     1.04
C1:LSC-TRX_OUT16 :            1.07     1.04     1.08

Only LO beam to BHD DCPDs
                          Mean     Max      Min
C1:IOO-MC_TRANS_SUM :     14121.32 14057.71 14159.38
C1:HPC-DCPD_A_OUT16 :       129.80   128.37   130.68 (Consistent with, 40m/17046. Power as expected within 20%. Squashed shape)
C1:HPC-DCPD_B_OUT16 :       123.42   121.92   124.48

ITMX single bounce (ITMY, ETMX, ETMY, PRM, SRM, LO misalgined)
                          Mean     Max      Min
C1:IOO-MC_TRANS_SUM :     14105.13 14000.89 14171.91
C1:HPC-DCPD_A_OUT16 :        92.54    91.45    93.30 (Consistent with 40m/17040, Power as expected within 40%. Clipped to the left in camera)
C1:HPC-DCPD_B_OUT16 :       137.70   136.55   138.53 (Note that DCPD_A/B ratio is different from LO, due to BHD BS R/T unbalance; 40m/17044)
C1:LSC-ASDC_OUT16 :           0.10     0.09     0.10 (Power as expected 40m/16952. Clipped to the right in camera)
C1:LSC-POPDC_OUT16 :        309.19   288.93   327.10 (Power as expected within 30% 40m/17042.)
C1:LSC-REFLDC_OUT16 :         0.02     0.01     0.02

ITMY single bounce (ITMX, ETMX, ETMY, PRM, SRM, LO misalgined)
                          Mean     Max      Min
C1:IOO-MC_TRANS_SUM :     14112.09 14025.37 14154.51
C1:HPC-DCPD_A_OUT16 :        92.58    92.01    93.26
C1:HPC-DCPD_B_OUT16 :       137.68   136.81   138.27
C1:LSC-ASDC_OUT16 :           0.10     0.09     0.10
C1:LSC-POPDC_OUT16 :        308.48   290.49   319.73
C1:LSC-REFLDC_OUT16 :         0.02     0.01     0.02

MICH fringe only (ETMX, ETMY, PRM, SRM, LO misalgined)
                          Mean     Max      Min
C1:IOO-MC_TRANS_SUM :     14090.34 13979.15 14143.86
C1:HPC-DCPD_A_OUT16 :       325.60    91.92   714.57
C1:HPC-DCPD_B_OUT16 :       400.27    18.37   762.57
C1:LSC-ASDC_OUT16 :           0.19    -0.05     0.41
C1:LSC-POPDC_OUT16 :        595.66  -119.21  1334.11
C1:LSC-REFLDC_OUT16 :         0.03    -0.01     0.07

LO-ITMX fringe only (ITMY, ETMX, ETMY, PRM, SRM misalgined)
                          Mean     Max      Min
C1:IOO-MC_TRANS_SUM :     14062.58 13968.05 14113.67
C1:HPC-DCPD_A_OUT16 :       224.31    89.57   371.66
C1:HPC-DCPD_B_OUT16 :       259.74    85.37   421.86

Next:
 - Measure contrast (40m/17020) and estimate mode-matching of LO-AS again (40m/17041)
 - Now that we have better LO-AS fringe, lock LO phase in MICH (40m/17037)
 - Now that Dolphin issue was fixed, try double-demodulation to lock LO phase

[Anchal, Yehonathan, Yuta]

We did the constrast measurement with the method same as 40m/17020.
Contrast between ITM single bounce and LO beam increased to 74% (we had 43% before unclipping LO beam in 40m/17056).
From equations in 40m/17041 and measured ITM sigle bounce power (93 or 138 counts @ BHD DCPD) and LO power (130 or 124 counts @ BHD DCPD) from 40m/17056,  expected visibility for perfectly mode-matched case is 99%.
Measured constrast of 74% indicate mode-matching of 56%.

Both arms locked, MICH fringe (20% percentile)
Contrast measured by C1:LSC-ASDC_OUT is 80.66 +/- 0.20 %
Contrast measured by C1:LSC-POPDC_OUT is 92.27 +/- 0.66 %
Contrast measured by C1:LSC-REFLDC_OUT is 89.59 +/- 0.84 %
Contrast measured by all is 87.51 +/- 1.69 %

Both arms misaligned, MICH fringe (20% percentile)
Contrast measured by C1:LSC-ASDC_OUT is 82.50 +/- 0.61 %
Contrast measured by C1:LSC-POPDC_OUT is 94.18 +/- 0.26 %
Contrast measured by C1:LSC-REFLDC_OUT is 92.78 +/- 0.19 %
Contrast measured by all is 89.82 +/- 1.75 %

ITMX-LO fringe (40% percentile)
Contrast measured by C1:HPC-DCPD_A_OUT is 73.93 +/- 1.52 %
Contrast measured by C1:HPC-DCPD_B_OUT is 73.56 +/- 1.22 %
Contrast measured by all is 73.74 +/- 0.98 %

ITMY-LO fringe (40% percentile)
Contrast measured by C1:HPC-DCPD_A_OUT is 73.45 +/- 0.61 %
Contrast measured by C1:HPC-DCPD_B_OUT is 75.27 +/- 0.50 %
Contrast measured by all is 74.36 +/- 0.54 %

[Tega, Anchal, Yuta]

We resored FPMI locking settings. Below is the summary of locking configurations tonight.
To ease the lock acquisition, the step to feedback POX11_I to ETMX and POY11_I to MC2 before POX and POY mixing was necessary tonight.

CARM (YARM):
 - 0.5 * POX11_I + 0.5 * POY11_I handed to 0.5 * REFL55_I
 - YARM filter module, FM4,5 for acquisition, FM1,2,3,6,8 triggered, C1:LSC-YARM_GAIN = 0.012
 - Actuation on -0.77 * MC2
 - UGF ~ 250 Hz

DARM (XARM):
 - 0.5 * POX11_I - 0.5 * POY11_I handed to 4.6 * AS55_Q (it was 2.5 in 40m/17012)
 - XARM filter module, FM5 for acquisition (no FM4), FM1,2,3,6,8 triggered, C1:LSC-XARM_GAIN = 0.015
 - Actuation on 0.5 * ETMX - 0.5 * ETMY
 - UGF ~ 120 Hz

MICH:
 - 1 * REFL55_Q (turned on after XARM and YARM acquisition)
 - MICH filter module, FM4,5,8 for acquisition, FM2,3 triggered, C1:LSC-MICH_GAIN = +40
 - Actuation on 0.5 * BS
 - UGF ~ 100 Hz

Measured sensing matrix:
Sensing Matrix with the following demodulation phases
{'AS55': 200.41785156862835, 'REFL55': 93.7514468401475, 'POX11': 105.08325063571438, 'POY11': -11.343909976281823}
Sensors              DARM                    CARM                   MICH
C1:LSC-AS55_I_ERR_DQ 5.27e-02 (-154.105 deg) 2.83e-01 (132.395 deg) 1.17e-04 (-40.1051 deg)
C1:LSC-AS55_Q_ERR_DQ 3.99e-02 (-151.048 deg) 1.42e-02 (125.504 deg) 1.41e-04 (-2.42846 deg)
C1:LSC-REFL55_I_ERR_DQ 5.59e-02 (77.6871 deg) 1.15e+00 (-44.589 deg) 3.55e-04 (69.2585 deg)
C1:LSC-REFL55_Q_ERR_DQ 1.84e-03 (16.3186 deg) 3.35e-03 (125.67 deg) 4.59e-05 (4.18718 deg)
C1:LSC-POX11_I_ERR_DQ 1.54e-01 (-157.852 deg) 6.07e-01 (-42.1078 deg) 5.55e-05 (73.3963 deg)
C1:LSC-POX11_Q_ERR_DQ 6.83e-05 (-148.591 deg) 6.37e-04 (121.983 deg) 1.35e-06 (43.7201 deg)
C1:LSC-POY11_I_ERR_DQ 1.85e-01 (36.1624 deg) 5.73e-01 (-43.1776 deg) 2.12e-04 (82.16 deg)
C1:LSC-POY11_Q_ERR_DQ 2.16e-05 (130.937 deg) 6.38e-05 (-173.194 deg) 1.40e-06 (47.5416 deg)

FPMI locked periods:
  - 1344129143 - 1344129520
  - 1344131106 - 1344131305
  - 1344133503 - 1344134020

Next:
- Restore CM servo for CARM

FPMI related
- Better suspension damping HIGH
 - Investigate ITMX input matrix diagonalization (40m/16931)
 - Output matrix diagonalization
 * FPMI lock is not stable, only lasts a few minutes for so. MICH fringe is too fast; 5-10 fringes/sec in the evening.
- Noise budget HIGH
 - Calibrate error signals (actually already done with sensing matrix measurement 40m/17069)
 - Make a sensitivity curve using error and feedback signals (actuator calibration 40m/16978)
 * See if optical gain and actuation efficiency makes sense. REFL55 error signal amplitude is sensitive to cable connections.
- FPMI locking
 - Use CARM/DARM filters, not XARM/YARM filters
 - Remove FM4 belly
 - Automate lock acquisition procedure
- Initial alignment scheme
 - Investigate which suspension drifts much
 - Scheme compatible with BHD alignment
 * These days, we have to align almost from scratch every morning. Empirically, TT2 seems to recover LO alignment and PR2/3 seems to recover Yarm alignment (40m/17056). Xarm seems to be stable.
- ALS
 - Install alignment PZTs for Yarm
 - Restore ALS CARM and DARM
 * Green seems to be useful also for initial alignment of IR to see if arms drifted or not (40m/17056).
- ASS
 - Suspension output matrix diagonalization to minimize pitch-yaw coupling (current output matrix is pitch-yaw coupled 40m/16915)
 - Balance ITM and ETM actuation first so that ASS loops will be understandable (40m/17014)
- Suspension calibrations
 - Calibrate oplevs
 - Calibrate SUSPOS/PIT/YAW/SIDE signals (40m/16898)
 * We need better understanding of suspension motions. Also good for A2L noise budgeting.
- CARM servo with Common Mode Board
 - Do it with single arm first

BHD related
- Better suspension damping HIGH
 - Invesitage LO2 input matrix diagonalization (40m/16931)
 - Output matrix diagonalization (almost all new suspensions 40m/17073)
 * BHD fringe speed is too fast (~100 fringes/sec?), LO phase locking saturates (40m/17037).
- LO phase locking
 - With better suspensions
 - Measure open loop transfer function
 - Try dither lock with dithering LO or AS with MICH offset (single modulation)
 - Modify c1hpc/c1lsc so that it can modulate BS and do double demodulation, and try double demodulation
- Noise Budget HIGH
 - Calibrate MICH error signal and AS-LO fringe
 - Calibrate LO1, LO2, AS1, AS4 actuation using ITM single bounce - LO fringe
 - Check BHD DCPD signal chain (DCPD making negative output when fringes are too fast; 40m/17067)
 - Make a sensitivity curve using error and feedback signals
- AS-LO mode-matching 
 - Model what could be causing funny LO shape
 - Model if having low mode-matching is bad or not
 * Measured mode-matching of 56% sounds too low to explain with errors in mode-matching telescope (40m/16859, 40m/17067).

IMC related
- WFS loops too fast (40m/17061)
- Noise Budget
- Investigate MC3 damping (40m/17073)
- MC2 length control path

I measured the sideband frequencies for PMC and IMC lock (to use it for Mariner PMC and IMC design).

PMC: 35.498912(2) MHz
IMC: 29.485038(2) MHz

Details:
 - Mini-Circuits UFC-6000 was used. The spec sheet says the frequency accuracy in 1-40 MHz is +/- 2 Hz.
 - "29.5 MHz OUT" port of 40m Frequency Generation Unit (LIGO-T1000461) was connected to UFC-6000 to measure IMC sideband frequency.
 - "LO TO SERVO" port of Crystal Frequency Ref (LIGO-D980353) was connected to UFC-6000 to measure PMC sideband frequency.
 - It seems like IMC sideband frequency changed from 29.485 MHz to 29.491 MHz back in 2011 (40m/4621). We are back to 29.485 MHz. I'm not sure what happened after this.

[Paco, Yuta]

We are still struggling with locking LO phase in MICH or ITM single bounce with BH55 with audio dither.
Without audio dither, BH55 can be used to lock.

What works:
 - LO phase locking with ITMX single bounce, using BH55_Q
  - BH55_Q configuration: 45 dB whitening gain, with whitening filter on.
  - C1:LSC-BH55_PHASE_R=147.621 deg gives most signal in BH55_Q.
  - LO phase can be locked using BH55_Q, C1:HPC-LO_PHASE_GAIN=-0.5 (bright fringe for A, dark for B), feeding back to LO1 gives UGF of ~80Hz (funny structure in ~20 Hz region; see Attachment #1)

 - LO phase locking with ITMX single bounce, using BHDC_DIFF
  - BHDC B/A = 1.57 (gain balanced with C1:HPC-IN_MTRX)
  - LO phase can be locked using BHDC_DIFF, C1:HPC-LO_PHASE_GAIN=-0.4 (mid-fringe lock), feeding back to LO1 gives UGF of ~50 Hz (see Attachment #2).

 - LO phase locking with MICH locked with AS55_Q, using BH55_Q
  - AS55_Q configuration: 24 dB whitening gain, with whitening filter off
  - C1:LSC-AS55_PHASE_R=-150 deg gives most signal in AS55_Q
  - MICH can be locked using AS55_Q, C1:LSC-MICH_GAIN=-10, C1:LSC-MICH_OFFSET=30 (slightly off from AS dark fringe), feeding back to 0.5*BS gives UGF of ~100Hz (see Attachment #3)
  - LO phase can be locked using BH55_Q, C1:HPC-LO_PHASE_GAIN=-0.8 (bright fringe for A, dark for B), feeding back to LO1 gives UGF of ~45Hz (see Attachment #4)

 - LO phase locking with MICH locked with AS55_Q, using BHDC_DIFF
  - LO phase can be locked using BHDC_DIFF, C1:HPC-LO_PHASE_GAIN=1 (mid-fringe lock), feeding back to LO1. Not a very stable lock.

What does not work:
 - LO phase locking using BH55_Q demodulated at LO1 (or AS1) dither frequency, neither in ITMX sigle bounce or MICH locked with/without offset using AS55_Q
  - C1:HPC-AS1_POS_OSC_FREQ=142.7 Hz, C1:HPC-AS1_POS_OSC_CLKGAIN=3000, C1:HPC-BH55_Q_AS1_DEMOD_PHASE=-15 deg, BLP30 is used.
  - Attachment #5 shows error signals when LO phase is locked with BH55_Q. BHDC_DIFF and BH55_Q_AS1_DEMOD_I having some coherence is a good indication, but we cannot lock LO phase with BH55_Q_AS1_DEMOD_I yet.
  - Also, injection at 13.14 Hz with an amplitude of 300 for AS1 can be seen in both BH55_Q and BH55_Q_AS1_DEMOD_I (26 Hz peak for BHDC_DIFF, as it is quadratic, as expected), which means that BH55_Q_AS1_DEMOD_I is seeing something.

Next:
 - Check actuation TFs for LO1, LO2, AS1 too see if there are any funny structures at ~ 20 Hz.
 - LO phase locking might require at least ~50 Hz of UGF. Use higher audio dither frequency so that we can increase the control bandwidth.
 - Check analog filtering situation for BHDC_A and BHDC_B signals (they go minus when fringes are moving fast)

[Paco, Yuta]

We removed splitter to route BHD DC PD signals to c1hpc and c1lsc. This was necessary to circumvent IPC error, but this is no longer necessary. Now BHD DC PD signals are ADC-ed with c1hpc, and sent to c1lsc via IPC.
We also found that BHD DC PD signals have whitening filters as described in LIGO-T2000500 (Readout board is LIGO-D1400384).
We added unwhitening filter zpk([151.9;3388],[13.81],1,"n") to C1:HPC_BHDC_A and B, based on measured whitening stage gain (see Sec 3.1 of characterization reoprt in LIGO-T2000500).
This solved the signals leaking to minus (40m/17068).

Next:
 - Modify c1hpc model to send BHD DCPD signals to c1lsc after unwhitening. (Note added on Nov 15: The same unwhitening filter is also added to C1:LSC-DCPD_A and B for now. See attached.)
 - Redo visibility measurements,

[Paco, Yuta]

MICH was locked with BHD DCPD A-B signal with LO phase controlled.
Locking procedure and configuration was as follows (see Attachment #1).

1. Lock MICH with AS55_Q, with C1:LSC-MICH_GAIN=-3, FM4, FM5, FM8, FM10 (boost filters are turned off to have more phase margin).

2. Lock LO PHASE with BH55_Q, with C1:HPC-LO_PHASE_GAIN=6, FM5, FM8, feeding back to AS1.
  - C1:LSC-BH55_PHASE_R=136.136 deg was tuned to minimize I when AS-LO is fringing with MICH locked with an offset of 50 (we first thought 136.136 deg - 90 deg is better from 40m/17216, but today, 136.136 deg seems to work better; Reason needs to be investigated).
  - We are supposed to use C1:HPC-BH55_Q_AS1_DEMOD_I_OUT to control the LO phase to give maximum MICH signal on BHD_DIFF (40m/17170), but somehow BH55_Q without audio dither was OK to get MICH signal. Line injection at 211.1 Hz on BS was seen in BHDC_DIFF (and AS55_Q), even if we use BH55_Q to lock LO PHASE (see Attachment #2; MICH_B is BHDC_DIFF and MICH_A is AS55_Q) or BH55_Q_AS1_DEMOD_I to lock LO PHASE (with both signs). Reason needs to be investigated.
  - Audio dither was done using AS1 with excitation of 15000 counts at 281.79 Hz. C1:HPC-BH55_Q_AS1_DEMOD_PHASE=60 deg was tuned to minimize Q with injection of line at 13 Hz using LO1.

3. Handed over MICH lock from AS55_Q to 0.66 * C1:LSC-DCPD_A - 1 * C1:LSC-DCPD_B. This was done by using C1:LSC-MICH_A and MICH_B gains. C1:LSC-MICH_A_GAIN=1 was handed over to C1:LSC-MICH_B_GAIN=-1.
  - 0.66 * A - B was tuned so that BHDC_DIFF will be zero (as it supposed to be with MICH offset of zero).
  - AS55_Q and BHDC_DIFF had roughly the same optical gain at 211.1 Hz (actually, BHDC_DIFF had higher optical gain; see Attachment #2), so we used MICH_A_GAIN=1 and C1:LSC-MICH_B_GAIN=-1
  - After handing over of BHDC_DIFF, OLTF was measured. UGF was ~70 Hz (Attachment #3).

Next:
  - Investigate how to get optimal LO phase. With BH55_Q or BH55_Q + audio dithering? How to optimize demod phases?
  - How do we balance DCPD A and B? What is the effect of BHD BS being 44:56 not 50:50?
  - Measure amount of MICH signal in BHDC_DIFF with different LO phases.
  - Improve SNR in BH55.
  - It will be much simpler if we send BHDC_SUM and BHDC_DIFF to c1lsc from c1hpc, instead of sending un-unwhitened BHDC_A and B.

BHD fringe visibility was measured again with unwhitening filters on on BHDC_A and B, which removed signal leakage to zero (40m/17265).
The result didn't change much from previous measurement (40m/17067) thanks to using the 'mode' of signal to calculate visibility.
Measured constrast of 74% indicate mode-matching AS beam to LO beam of 56%.

ITMX-LO fringe (10% percentile)
Contrast measured by C1:HPC-BHDC_A_OUT is 74.46 +/- 0.07 %
Contrast measured by C1:HPC-BHDC_B_OUT is 74.25 +/- 0.07 %
Contrast measured by all is 74.35 +/- 0.07 %

ITMY-LO fringe (10% percentile)
Contrast measured by C1:HPC-BHDC_A_OUT is 74.01 +/- 0.10 %
Contrast measured by C1:HPC-BHDC_B_OUT is 73.85 +/- 0.09 %
Contrast measured by all is 73.93 +/- 0.08 %

Errors are from standard deviation of 3 measurements.
The notebook lives in /opt/rtcds/caltech/c1/Git/40m/scripts/CAL/BHD/measureContrast.ipynb

Optical gains of AS55 and BH55 are calibrated for BHD MICH.

LO-ITM single bounce:
 With LO-ITM signle bounce fringe, optical gain of BH55_Q is measured using a method similar to MICH calibration in AS55 (40m/16929).
 Demodulation phase for BH55 is tuned to minimize I when LO-ITM is freeswinging (using getPhaseAngle.py).
 (Notebook: /opt/rtcds/caltech/c1/Git/40m/scripts/CAL/BHD/BHDOpticalGainCalibration.ipynb)
 Results are the following:

 LO-ITMY fringe: 7.84e9 counts/m (demod phase 147.1 +/- 0.3 deg) See Attachment #1
 LO-ITMX fringe: 8.44e9 counts/m (demod phase 149.6 +/- 0.4 deg) See Attachment #1

 Difference in the optimal demodulation phase 2.5 +/- 0.5 deg agrees with half of Schnupp asymmetry, as expected (40m/17007).
 Difference in the optical gain for LO-ITMY and LO-ITMX is probably from statistical fluctuation.

BHD MICH:
 Sensing matrix was measured by injecting a line at BS (300 counts @ 211.1 Hz), LO1 (5000 counts @ 287.1 Hz) and AS1 (5000 counts @ 281.79 Hz), when MICH is locked with AS55_Q and LO PHASE is locked with BH55_Q (both with no offset).
 Using the sensing matrix, demodulation phase was tuned to minimize I phase for MICH signal in AS55 and LO1 signal in BH55.
 After the demodulation phase tuning. sensing matrix was measured to be the following.
 See, also Attachment #3 for injected peaks. I phase signal is successfully suppressed by at least an order of magnitude.
 (Notebook: /opt/rtcds/caltech/c1/Git/40m/scripts/CAL/SensingMatrix/MeasureSensMatBHD.ipynb)

Sensing Matrix with the following demodulation phases (counts/counts)
{'AS55': -160.15695076011946, 'BH55': 154.13916838400047}
      Sensors           MICH @211.1 Hz          LO1 @287.1 Hz            AS1 @281.79 Hz           
C1:LSC-AS55_I_ERR_DQ    1.22e-05 (120.53 deg)   7.24e-07 (85.64 deg)     1.26e-06 (40.42 deg)    
C1:LSC-AS55_Q_ERR_DQ    2.95e-03 (-101.62 deg)  1.24e-06 (-80.43 deg)    1.69e-06 (152.31 deg)    
C1:LSC-BH55_I_ERR_DQ    1.28e-03 (80.95 deg)    3.44e-06 (109.31 deg)    2.22e-06 (154.40 deg)    
C1:LSC-BH55_Q_ERR_DQ    7.44e-03 (77.38 deg)    2.56e-04 (-59.85 deg)    2.42e-04 (6.40 deg)    
C1:HPC-BHDC_DIFF_OUT    2.21e-03 (82.45 deg)    4.37e-05 (121.87 deg)    3.61e-05 (-169.09 deg)

 Using BS actuation efficiency of 26.08e-9 /f^2 m/counts (40m/16929), optical gain for AS55_Q and BHDC_DIFF for MICH is

2.95e-03 / (26.08e-9/(211.1**2)) = 5.04e9 counts/m (AS55_Q for MICH)
2.21e-03 / (26.08e-9/(211.1**2)) = 3.78e9 counts/m (BHDC_DIFF for MICH)

 For AS55_Q, this is a factor of 4~5 higher than the previous measurement from free swing (40m/16929). Why?
 Free swing measurement was done again, and this gave 1.24e9 counts/m, which is consistent with the previous measurement (see Attachment #3).

 Using LO1 and AS1 actuation efficiencies of 3.14e-8 /f^2 m/counts (40m/17206), optical gains for BH55_Q for LO1 and AS1 are

2.56e-04 / (3.14e-8/(287.1**2)) = 6.72e8 counts/m (BH55_Q for LO1)
2.42e-04 / (3.14e-8/(281.79**2)) = 6.12e8 counts/m (BH55_Q for AS1)

Next:
 - Compare them with expected values
 - Measure them with different locking points (different LO phases, MICH offsets)
 - Investigate why MICH optical gain in AS55 is 4~5 times higher than free swing measurement (use different modulation frequency?)

Summary of actuation calibration so far (counts from C1:LSC-xx_EXC or C1:SUS-xx_LSC_EXC):
BS   : 26.08e-9 /f^2 m/counts (see 40m/16929)
ITMX :  5.29e-9 /f^2 m/counts (see 40m/16929)
ITMY :  4.74e-9 /f^2 m/counts (see 40m/16929)
ETMX : 10.91e-9 /f^2 m/counts (see 40m/16977 and 40m/17014)
ETMY : 10.91e-9 /f^2 m/counts (see 40m/16977)
MC2 : -14.17e-9 /f^2 m/counts in arm length (see 40m/16978)
MC2 :   5.06e-9 /f^2 m/counts in IMC length (see 40m/16978)
LO1 : 3.14e-8 / f^2 m/counts (see 40m/17206)
LO2 : 2.52e-8 / f^2 m/counts (see 40m/17206)
AS1 : 3.14e-8 / f^2 m/counts (see 40m/17206)
AS4 : 2.38e-8 / f^2 m/counts (see 40m/17206)

[Paco, Yuta]

We found that MICH UGF was unexpectedly high, ~200 Hz, in the measurement yesterday, which makes the closed loop gain to be more than one at MICH line injection at 211.1 Hz.
We did optical gain calibrations for AS55, BH55 and BHDC_DIFF in BHD MICH again with UGF at around 10 Hz.
This solved the inconsistent result with free swing calibration.

What we did:
 Did the same measurement for BHD MICH as written in 40m/17274, but with MICH UGF of ~10 Hz and LO PHASE UGF of ~15 Hz (see OLTFs in Attachment #1, and filter configurations in Attachment #2).
 Updated sensing matrix is as follows

Sensing Matrix with the following demodulation phases (counts/counts)
{'AS55': -163.52789698340882, 'BH55': 152.7860744565449}
      Sensors        	MICH @211.1 Hz       	LO1 @287.1 Hz       	AS1 @281.79 Hz       	
C1:LSC-AS55_I_ERR_DQ	1.85e-05 (-118.82 deg)	3.31e-07 (-32.19 deg)	7.86e-07 (112.27 deg)	
C1:LSC-AS55_Q_ERR_DQ	7.32e-04 (59.57 deg)	1.19e-06 (158.17 deg)	9.07e-07 (-92.25 deg)	
C1:LSC-BH55_I_ERR_DQ	5.02e-04 (-123.21 deg)	1.79e-05 (-26.73 deg)	1.76e-05 (-120.23 deg)	
C1:LSC-BH55_Q_ERR_DQ	1.75e-03 (59.57 deg)	2.71e-04 (-22.64 deg)	2.56e-04 (-114.37 deg)	
C1:HPC-BHDC_DIFF_OUT	1.00e-03 (-115.93 deg)	3.09e-05 (-14.99 deg)	2.84e-05 (-110.23 deg)	

 Using BS actuation efficiency of 26.08e-9 /f^2 m/counts (40m/16929), optical gain for AS55_Q and BHDC_DIFF for MICH is

7.32e-03 / (26.08e-9/(211.1**2)) = 1.25e9 counts/m (AS55_Q for MICH) This is consistent with freeswing measurement (1.24e9 m/counts) 40m/17274
1.00e-03 / (26.08e-9/(211.1**2)) = 1.71e9 counts/m (BHDC_DIFF for MICH)

 Using LO1 and AS1 actuation efficiencies of 3.14e-8 /f^2 m/counts (40m/17206), optical gains for BH55_Q for LO1 and AS1 are

2.71e-04 / (3.14e-8/(287.1**2)) = 7.12e8 counts/m (BH55_Q for LO1)
2.56e-04 / (3.14e-8/(281.79**2)) = 6.47e8 counts/m (BH55_Q for AS1)

  (Notebook: /opt/rtcds/caltech/c1/Git/40m/scripts/CAL/SensingMatrix/MeasureSensMatBHD.ipynb)

Next:
 - Compare them with expected values
 - Measure them with different locking points (different LO phases, MICH offsets; LO phase can be calibrated using optical gain calibration of BH55_Q)

MICH optical gain was measured with different LO phases over ~90 degrees.
Zero crossing of BH55_Q_ERR seems to be roughly 55 degrees away from optimal LO phase.

What we did:
 - Locked MICH with AS55_Q with no offset, with UGF at ~10 Hz (same as configuration in 40m/17279).
 - Injected BS calibration line at amptilude of 300 counts at 211.1 Hz.
 - Locked LO Phase with BH55_Q with different offsets added at C1:HPC-LO_PHASE_OFFSET.
 - Measured sensing matrix at that frequency. Counts are calibrated into meters using actuator efficiencies as described in 40m/17279.
 - LO phase was obtained using a DC value of BH55_Q. This was calibrated into degrees from the following:

Amplitude of LO-AS fringe in BH55_Q was calculated to be

A = BH55optgain*lamb/(4*pi) = 60 counts

where BH55optgain is 7.12e8 counts/m, which is optical gain of BH55_Q for LO1 measured in 40m/17279.
(Actually, BH55_Q goes upto ~ +/-200 counts in time series data, but maybe 60 is the nominal fringe amplitude, considering alignment fluctuations and fluctuation in AS darkness? Note that, no offset in BH55_Q is assumed in this calculation, but AM etc can create an offset.) 
LO phase can be obtained by

LOphase = arcsin(BH55_Q/A)

where BH55_Q a DC value (10 sec average) of BH55_Q.

Result:
 Attachment #1 is uncalibrated plot C1:HPC-LO_PHASE_OFFSET of around +/- 50 was the maximum we could add, and more offset gave unstable lock.
 Attachment #2 is calibrated plot. AS55_Q does not depend on LO phase, as expected. BH55_Q and BHDC_DIFF depend on LO phase as expected. BH55_I and AS55_I stay at low level, as expected (this means that our RF demodulation phase is OK).
 Dotted gray line is an eyeball fit of expected curve (40m/17170) to fool your eyes.
 This tells you that we are roughly 55 deg away from LO phase which gives maximum MICH signal for BHDC_DIFF.
 Error bar in x-axis is from standard deviation of BH55_Q fluctuations. Error in y axis is probably ~20% at maximum.
 Notebook: /opt/rtcds/caltech/c1/Git/40m/scripts/CAL/SensingMatrix/MeasureSensMatBHD.ipynb

Next:
 - Repeat the measurement with
   - MICH locked with higher UGF, with notch at 211.1 Hz, for more robust AS dark fringe
   - DCPD A and B balanced at 211.1 Hz (null MICH signal for BHDC_SUM to balance?)
   - Measure optical gain also for BHDC_SUM and BH55_Q demodulated at audio dither
   - Lock LO phase at different sign so that we can sweep LO phase over ~180 deg
   - Sign-sensitive optical gain measurement (demodulation with BS motion necessary)
 - Compare with expected values from simulations
 - Why do we have 55 degrees offset? Expected offset is 90 degrees...
   - Check if there is any RAM in 55 MHz in the input beam by measuring AM with ITM single bounce

[Paco, Yuta]

We realized that bandstop filters ("violin" filters) we implemented in 40m/17206 had pass band gain of -1dB.
gain(1,"dB") was added to all the filters (see Attachment #1 for gain adjusted violin filters for AS1).
We also realized that audio dither frequency we chose to generate BH55+audio dither error signal and to measure sensing matrix at ~280 Hz was too close to violin filters.
These will affect calibrations by upto ~60%.
For example, actuation gains should be actually

• LO1 = 3.14e-8 / f^2 m / cts * 3dB = 4.44e-8 / f^2 m / cts (3 violin filters)
• LO2 = 2.52e-8 / f^2 m / cts * 0dB = 2.52e-8 / f^2 m / cts (no violin filters)
• AS1 = 3.14e-8 / f^2 m / cts * 3dB = 4.44e-8 / f^2 m / cts (3 violin filters)
• AS4 = 2.38e-8 / f^2 m / cts * 4dB = 3.36e-8 / f^2 m / cts (3 violin filters+ bandstop at 96.7 Hz)

Next:
 - Redo actuator calibrations for LO1, LO2, AS1, AS4
 - Redo sensing matrix measurements with different audio dither frequencies for LO1 and AS1

As there is some confusion in actuator calibration, we have done the measurement again from scratch.
Results are the following.
New values for LO1, LO2, AS1, AS4 are obtained from free swinging ITMY-LO, so it should be more robust.

BS   : 26.54e-9 /f^2 m/counts
ITMX :  4.93e-9 /f^2 m/counts
ITMY :  4.90e-9 /f^2 m/counts
LO1  : 26.34e-9 /f^2 m/counts
LO2  :  9.81e-9 /f^2 m/counts
AS1  : 23.35e-9 /f^2 m/counts
AS4  : 24.07e-9 /f^2 m/counts

BS, ITMX, and ITMY actuator calibration:
 Followed the procedure in 40m/16929.
 Calibrated AS55_Q using X-Y plot to be 9.72e8 counts/m (Attachment #1), locked MICH with UGF of 10 Hz, and measured the transfer function from C1:LSC-BS,ITMX,ITMY_EXC to C1:LSC-AS55_Q_ERR.
 The result is Attachment #2. They are consistent with 40m/16929.

LO1, LO2, AS1, and AS4 actuator calibration:
 Followed similar steps with ITMY-LO fringe.
 Calibrated BH55_Q using X-Y plot to be 7.40e9 counts/m (Attachment #3), locked ITMY-LO with UGF of ~15 Hz (Attachment #4), and measured the transfer function from C1:SUS-LO1,LO2,AS1,AS4_LSC_EXC to C1:LSC-BH55_Q_ERR.
 The result is Attachment #5. They are inconsistent with 40m/17284, but this one should be more robust (see discussions below).

LO1, LO2, AS1, and AS4 actuator calibration by taking the ratio between ITMY:
 We have also followed the steps in 40m/17206 to calibrate BHD actuators.
 This method does not depend on BH55_Q optical gain calibration, but depends on ITMY calibration.
 Measured OLTFs for ITMY-LO fringe locking is Attachment #6, and actuator ratio with respect to ITMY is Attachment #7. In this measurement, Bandstop filter at 96.7 Hz for AS4 was turned off, and gain was lowered by a factor of 2 to avoid AS4 oscillating.
 This gives

LO1  : 116.81e-9 /f^2 m/counts
LO2  :  51.69e-9 /f^2 m/counts
AS1  : 101.48e-9 /f^2 m/counts
AS4  : 117.84e-9 /f^2 m/counts

 These are not consistent with 40m/17284, and larger by a factor of ~2-3.
 These are also not consistent with the values from free swinging measurement, and are larger by a factor of ~4-5.
 I guess there are some gains missing when comparing ITMY loop in c1lsc and other loops in c1hpc.

MICH optical gain with a sign was measured with different LO phases over ~180 degrees, with updated calibration and higher MICH UGF.
Zero crossing of BH55_Q_ERR seems to be 68 degrees away from optimal LO phase.

Calibrated sensing matrix:
 - Locked MICH with AS55_Q at dark fringe, with UGF of ~200 Hz. Notch at 311.1 Hz was turned on.
 - Locked LO PHASE with BH55_Q, with UGF of ~10 Hz (C1:HPC-LO_PHASE_GAIN=-2, using LO1).
 - Measured the sensing matrix as written in 40m/17279, but with different dither frequencies to avoid violin mode frequencies and to match with already-installed notch filters.
 - Sensing matrix was calibrated into meters using actuator gains measured in 40m/17285
 - Sign was added by comparing the phase with C1:SUS-xx_LSC_OUT. If they are 90-270 deg apart, minus sign was added to the sensing matrix.
 - Resuts are as follows. At least important green ones are consistent with previous measurements (40m/17279).

Calibrated sensing matrix with the following demodulation phases (counts/m)
{'AS55': -164.1726747789845, 'BH55': 169.57651332419115}
      Sensors            MICH @311.1 Hz           LO1 @147.1 Hz           AS1 @141.79 Hz           
C1:LSC-AS55_I_ERR_DQ    2.72e+06 (84.89 deg)    6.41e+05 (14.60 deg)    -1.98e+05 (206.79 deg)    
C1:LSC-AS55_Q_ERR_DQ    -1.20e+09 (-228.85 deg)    -1.43e+06 (-106.51 deg)    1.41e+06 (29.21 deg)    
C1:LSC-BH55_I_ERR_DQ    -2.45e+09 (-230.64 deg)    -6.57e+07 (167.84 deg)    7.28e+07 (-16.56 deg)    
C1:LSC-BH55_Q_ERR_DQ    7.81e+09 (-48.64 deg)    -7.34e+08 (159.70 deg)    8.06e+08 (-10.08 deg)    
C1:HPC-BHDC_DIFF_OUT    -9.91e+08 (-224.55 deg)    -1.13e+08 (164.14 deg)    1.26e+08 (-3.95 deg)    
C1:HPC-BHDC_SUM_OUT    -6.84e+06 (-104.69 deg)    1.50e+07 (-8.71 deg)    -1.79e+07 (173.80 deg)    
LO phase from C1:LSC-BH55_Q_ERR_avg 4.98e-03 +/- 1.65e+01 deg

Estimating LO phase:
 - Using 7.34e+08 counts/m, which is an optical gain of BH55_Q for LO1, LO phase can be estimated as follows.

A = BH55optgain*lamb/(4*pi) = 62 counts
LOphase = arcsin(BH55_Q/A)

 - When C1:HPC-LO_PHASE_GAIN is plus, LOphase was calculated with the following to take into account of the sign flip in the controls.

LOphase = 180 - arcsin(BH55_Q/A)

Balancing A-B:
 - BHDC_A and BHDC_B were balanced to give null MICH signal in BHDC_SUM at 311.1 Hz. This gave BHDC_DIFF = 0.919*A - B.
 - It seems like this balancing gain changes over time by ~30%.

Result:
 - Attachment #1 is uncalibrated MICH optical gain in different LO phases, and Attachment #2 is the calbirated one. Basically the same with 40m/17282, but with updated calibration and sign considerations.
 - In addition to the previous measurements, we can see that BHD_SUM is not dependent on LO phase (small dependence probably from not perfect A and B balancing).
 - 0 deg of LO phase means that it is a zero crossing of BH55_Q with a slope that LO PHASE loop can be closed with a minus C1:HPC-LO_PHASE_GAIN, feeding back to LO1.
 - Dotted and dashed gray lines are from scipy.optimize.curve_fit using the following fitting function (not an eyeball fit this time!).

def fitfunc(x, a,b,c):
    return a*np.sin(np.deg2rad(x-b))+c

 - Fitting results show that we are -22 deg away from our intuition that BH55_Q crosses zero when BHDC_DIFF give no MICH signal (68 degrees away from optimal LO phase).
 - Fitting results also show that BH55_Q sensitivity to MICH crosses zero when BHD_DIFF sensitivity to MICH maximizes. This suggests that BH55+MICH dither can be used to lock LO phase to optimal LO phase.

Notebook: /opt/rtcds/caltech/c1/Git/40m/scripts/CAL/SensingMatrix/MeasureSensMatBHD.ipynb

Next:
 - Compare with expected values from simulations
 - Why do we have -22 deg?
   - Check if there is any RAM in 55 MHz in the input beam by measuring AM with ITM single bounce (quick measurement shows it is small)
   - Unbalanced BHD BS?
   - Contribution from 55 MHz sidebands from LO beating with 55 MHz sidebands from AS?
       - Lock LO phase using audio dither only (demodulate BHDC_DIFF?).

[Paco, Yuta]

MICH displacement sensitivity was compared under AS55_Q locking and BHD_DIFF locking.
Sensitivity with BHD was better by more than an order of magnitude due to smaller sensing noise.
During the measurement, LO phase fluctuation was ~13 deg RMS.

Locking configurations:
 - MICH was first locked with AS55_Q, no offset, and then handed over to BHD_DIFF after LO phase locked. FM2, FM3, FM4, FM5, FM6, FM8, FM10 on, C1:LSC-MICH_GAIN=-3 gave UGF of around 80 Hz.
 - LO PHASE was locked with BH55_Q, no offset. FM5, FM8 on, C1:HPC-LO_PHASE_GAIN=-2 feeding back to LO1 gave UGF of around 40 Hz.
 - Attachment #1 shows the OLTFs.

Sensitivity estimate:
 - Sensitivity was estimated using measured actuator gains and optical gains. Following numbers are used.

C1:LSC-AS55_Q_ERR to MICH 1.08e-9 counts/m (measured at 311.1 Hz today)
C1:HPC-BHDC_DIFF to MICH 1.91e-9 counts/m (measured at 311.1 Hz today)
BS   : 26.54e-9 /f^2 m/counts (40m/17285)
LO1  : 26.34e-9 /f^2 m/counts (40m/17285)

 These numbers were also reflected to C1:CAL-MICH_CINV and C1:CAL-MICH_A.
 C1:CAL-MICH_A_GAIN = 0.5 was used to take into account of LSC output matrix of MICH to BS being C1:LSC-OUTPUT_MTRX_8_2=0.5.

 - Attachment #2 shows the displacement spectrum of MICH (top) and LO PHASE (bottom). Brown MICH curve is when locked with AS55_Q and black MICH curve is when locked with BHD_DIFF. RMS of original and in-loop LO PHASE was estimated to be

 Original LO phase noise: 393 nm RMS (266 deg RMS)
 In-loop LO phase noise: 19.4 nm RMS (13 deg RMS)

Next:
 - Improve LO phase loops to reduce LO phase noise
 - Estimate LO phase noise contribution to MICH sensitivity

[Anchal, Yuta]

To send BHD signals from c1hpc after unwhitening and taking sum/diff, c1hpc and c1lsc models are modified.
PDDC_DOF_MTRX medm screen was modified to reflect this change.
We don't need to unwhiten and take sum/diff again in c1lsc model anymore
 

[Anchal, Paco, Yuta]

Attachment #1 is the same plot as 40m/17303 but with MICH sensitivity for ASDC and AS55 also included (in this measurement, BH55 demodulation phase was set to 140.07 deg to minimize I fringe).
Y-axis is now calibrated in to counts/m using BS actuation efficiency 26.54e-9 /f^2 m/counts (40m/17285) at 311.1 Hz.
2nd X-axis is calibrated into MICH offset using the measured AS55_Q value and it's MICH sensitivity, 8.81e8 counts/m (this is somehow ~10% less than our usual value 40m/17294).
ASDC have similar dependence with BHDC_SUM on MICH offset, as expected.
AS55_Q have little dependence with MICH offset on MICH offset, as expected.

This plot tells you that even a small MICH offset at nm level can create MICH sensitivity for BHDC_DIFF, even if we control LO phase to have BH55_Q to be zero, as MICH offset shifts zero crossing of BH55_Q for LO phase.

Notebook: /opt/rtcds/caltech/c1/Git/40m/scripts/CAL/BHD/BH_DIFFSens_pydemod.ipynb

Attachment #2 is the same plot, but BH55 demodulation phase was tuned to 227.569 deg to have no MICH signal in BH55_Q (a.k.a measurement (c)).
In this case, LO phase will be always controlled at 0 deg (90 deg away from optimal), even if we change the MICH offset, as BH55_Q will not be sensitive to MICH.
In this plot, BHD_DIFF have little sensitivity to MICH, irrelevant of MICH offset, as expected.
MICH sensitivity for BH55_I is also constant, which indicate that LO phase is constant over this measurement, as expected.

Notebook: /opt/rtcds/caltech/c1/Git/40m/scripts/CAL/BHD/BH_DIFFSens_pydemod.ipynb

Attachment #3 is the same plot, but BH55 demodulation phase was tuned to 70 deg.
This demodulation phase was tuned within 5 deg to maximize MICH signal in BHD_DIFF with large MICH offset (20).
In this case, LO phase will be always controlled at 90 deg (optimal), even if we change the MICH offset, as BH55_Q will not be sensitve to LO carrier x AS sideband component of the LO phase signal.
In this plot, BHD_DIFF have high sensitivity to MICH, irrelevant of MICH offset (at around zero MICH offset it is hard to see because LO_PHASE lock cannot hold lock, as there will be little LO phase signal in BH55_Q, and measurement error is high for BHD_DIFF and BH55 signals).
MICH sensitivity for BH55_I and BH55_Q is roughly constant, which indicate that LO phase is constant over this measurement, as expected.

These plots indicate that BH55 demodulated at MICH dither frequency can be used to control LO phase robustly at 90 deg, under unknown or zero MICH offset.


Notebook: /opt/rtcds/caltech/c1/Git/40m/scripts/CAL/BHD/BH_DIFFSens_pydemod.ipynb

LO phase delay:
 From these measurements of demodulation phases, I guess we can say that phase delay for 55 MHz in LO path with respect to MICH path (length difference in PR2->LO->BHDBS and PR2->ITMs->AS->BHDBS) is

2*(227.569-70(5)-90)-90 = 45(10) deg

 This means that the length difference is (omegam=5*2*pi*11.066195 MHz)

c * np.deg2rad(45(10)+360) / omegam = 6.1(2) m   (360 deg is added to make it close to the design)

  Is this consistent with our design? (According to Yehonathan, it is 12.02 m - 5.23 m = 6.79 m)

  Attachment #4 illustrates signals in BH55.

Next:
 - Lock LO PHASE with BH55 demodulated at MICH dither frequency (RF+audio double demodulation), and repeat the same measurement
 - Finer measurement at small MICH offsets (~1nm) to see how much MICH offset we have
 - Repeat the same measurement with BH55_Q demodulation phase tuned everytime we change the MICH offset to maximize LO phase sensitivity in BH55_Q (a.k.a measurement (b)).
 - What is the best way to tune BH55 demodulation phase?

[Anchal, Yuta]

We have modified c1cal model to support sensing matrix measurements for BHD PDs on Friday last week.
c1cal model now can inject dither to LO1, LO2, AS1, and AS4, and demodulate BH55_I, BH55_Q, BHDC_SUM and BHDC_DIFF signals.
Related models, c1lsc, c1hpc, and c1sus2 are also modifed accordingly.
MEDM screens are also edited accordingly.
Attachments highlight the modifications.

[Yehonathan, Yuta]

We found that moving AS1 in yaw improves power on ASDC and AS55.
We compensated this move with AS4 and SR2 to keep the BHD fringe (ITM single bounce and LO beam fringes ~600 counts in amplitude at BH55).
We have also aligned BHD CCD camera to avoid clipping on a lens just before the camera (all the other optics on ITMY table remain untouched).
After the alignment, MICH BHD sensing matrix were measured with new C1CAL model (40m/17339) under the following conditions.
 - Locked MICH with AS55_Q at dark fringe. Notch at 311.1 Hz was turned on.
 - Locked LO PHASE with BH55_Q with C1:HPC-LO_PHASE_GAIN=-2, using LO1.

Sensing matrix with the following demodulation phases (counts/m)
{'AS55': -161.16488964312092, 'BH55': 162.57275834049358}
Sensors      BS @311.1 Hz           LO1 @147.1 Hz           AS1 @141.79 Hz           
AS55_I       (-0.19+/-1.45)e+07    (-0.26+/-2.43)e+06    (+0.35+/-2.39)e+06    
AS55_Q       (-1.74+/-0.02)e+09    (+1.61+/-8.31)e+06    (+1.08+/-8.59)e+06    
BH55_I       (+3.01+/-0.17)e+09    (+3.20+/-9.59)e+07    (-3.67+/-9.46)e+07    
BH55_Q       (-6.77+/-0.45)e+09    (+1.09+/-0.17)e+09    (-1.22+/-0.18)e+09    
BHDC_DIFF    (-8.41+/-4.81)e+08    (-1.26+/-0.94)e+08    (+1.38+/-1.03)e+08    
BHDC_SUM     (-2.75+/-9.14)e+07    (+1.18+/-1.13)e+07    (-0.97+/-1.02)e+07    

AS55_Q optical gain to MICH and BH55_Q optical gain to LO phase was improved by ~45%, compared with previous measurements (see 40m/17287).
The value for AS55_Q is consistent with the free swing measurement as attached.
SENSMAT part of c1cal seems to be working fine.

[Yehonathan, Yuta]

Sensing matrix measurements at different LO phases were performed under LO phase locked to both BH55_Q and BH55_Q+MICH dither.
We confirmed that BH55_Q+MICHdither can lock LO phase to around maximum MICH sensitivity for BHD_DIFF.

Locking configuratons
 - MICH was lockied using AS55_Q feeding back to BS, at dark fringe. Notch at 311.1 Hz was turned on. C1:LSC-MICH_GAIN=-6 (lowered to reduce BS DAC saturation).
 - LO PHASE was locked using BH55_Q, feeding back to LO1. FM2, FM5, FM8 on. C1:HPC-LO_PHASE_GAIN=+/-2.
 - LO PHASE was also locked using BH55_Q+MICHdither. BS was dithered with C1:HPC-BS_POS_OSC_CLKGAIN=4000 at 281.768 Hz (2nd notch of ELP80 used for demodulation). Feeding back to LO1. FM5, FM8 on (no LF boost). C1:HPC-LO_PHASE_GAIN=+/-20.
  -- Note that we could not increase the dither amplitude more as BS DAC starts to saturate (we are using BS for MICH loop, sensing matrix measurement, and audio dither; see 40m/17343).

Sensing martix measurements
 - Lines are injected to BS @ 311.1 Hz with amplitude of 1000, LO1 @ 147.1 Hz and AS1 @ 141.79 Hz with amplitude of 5000.

Estimating LO phase
 - Estimation of LO phase was done in the same way described in 40m/17287. We used measured sensitivity of BH55_Q for LO1 at BH55_Q zero crossing (-1.42e9 counts/m) to estimate LO phase offset from BH55_Q zero crossing.
 - In BH55_Q+MICHdither case, LO phase was flipped using the following equation when C1:HPC-LO_PHASE_GAIN is minus (to have consistend LO phase dependence with BH55_Q locking. NEEDS CHECK).

LOphase = 180 - arcsin(BH55_Q/A)

Result
 - Attachment #1 shows the sensitivity of AS55, BH55, BHDC_DIFF/SUM to BS (upper panel), LO1 (middle) and AS1 (lower), under LO phase locked to BH55_Q. The upper plot is the same plot as 40m/17287. As we can see, "0 deg" in the x-axis is not the optimal phase for BHDC_DIFF to have maximum MICH sensitivity. "0 deg" is the optimal point in terms of BH55_Q sensitivity to LO1/AS1, as we tuned the demodulation phase to maximize it.
 - Attachment #2 shows the same plot, under LO phase locked to BH55_Q+MICH dither. Sensitivity of BH55_Q to MICH crosses zero at round these measurements, as we are zero-ing it with this locking scheme. Around these LO phases, sensitivity of BHDC_DIFF to MICH is maximized as expected. Also, sensitivity of BHDC_DIFF to LO1/AS1 is minimized, as expected (assuming residual MICH offset and contrast defect are small).
 - Attachment #3 is the combined data from #1 and #2. Data points from BH55_Q locking are marked with "o" and those from BH55_Q+MICH dither locking are marked with "x" (they have larger uncertainties in LO phase). Both measurements are somewhat inconsistent in some channels (BS to BHDC_DIFF and LO1/AS1 to BH55_Q). Needs further investigation.
 - Dashed lines are from scipy.optimize.curve_fit using the following fitting function.

def fitfunc(x, a,b,c):
    return a*np.sin(np.deg2rad(x-b))+c

Notebook: /opt/rtcds/caltech/c1/Git/40m/scripts/CAL/SensingMatrix/SensMatBHDvsLOPhase.ipynb

Next:
 - Lock MICH with BHDC_DIFF under LO phase locked to BH55_Q+MICHdither
 - Estimate LO phase noise contribution to MICH displacement sensitivity
 - Improve LO phase loop
 - Try audio+audio dither
 - Move on to FPMI
 - Move on to 44MHz
 - Estimate the amount of residual MICH offset and contrast defect from these plots

[Yehonathan, Yuta]

Here's a plot using same dataset from yesterday, but x-axis in raw BH55_Q data, not calibrated into degrees in LO phase.
This way you are free from calibration error in BH55_Q data to LO phase.
Elliptic fit is done using least squares.
dphi is calculated using the following equation where (ap, bp) are the semi-major and semi-minor axes, phi is the rotation of the semi-major axis from the x-axis.

beta=np.arctan(ap/bp/np.tan(phi))
dphi=-np.arctan(ap/bp*np.tan(phi))-beta

dphi gives you LO phase at zero-crossing.
For example, the top plot says that the sensitivity of BH55_Q to BS crosses zero at "-133.92 deg," which means BH55_Q+MICHdither can lock LO phase at -134 deg or 46 deg.
The top plot also says that the sensitivity of BHDC_DIFF to BS crosses zero at "127.45 deg," which means BHDC_DIFF sensitivity to MICH maximizes at 38 deg or 217 deg.
The middle plot says that the sensitivity of BH55_Q to LO1 crosses zero at "90.09 deg," which means BH55_Q+LO1dither can lock LO phase at 90 deg or -90 deg, and BH55_Q(no dither) can lock LO phase at 0deg or 180 deg (by definition).

Notebook: /opt/rtcds/caltech/c1/Git/40m/scripts/CAL/SensingMatrix/PlotSensMatBHDvsLOPhaseData.ipynb

Next
 - Use also BH55_Q+LO1/AS1dither to scan around 90 deg.

[Paco, Yuta]

We have recovered the IFO alignment (FPMI and BHD) for the first time after the power loss on Dec 15th 11:30 AM PST.

PMC and IMC transmission
~>cdsutils avg -s 10 C1:PSL-PMC_PMCTRANSPD C1:IOO-MC_TRANS_SUMFILT_OUT
C1:PSL-PMC_PMCTRANSPD 0.7332153499126435 0.0004376671844386436
C1:IOO-MC_TRANS_SUMFILT_OUT 14370.22451171875 38.87278766402092

BHDC PDs with LO beam only
~>cdsutils avg -s 10 C1:HPC-BHDC_A_OUT C1:HPC-BHDC_B_OUT
C1:HPC-BHDC_A_OUT 111.07633819580079 2.916326545853983
C1:HPC-BHDC_B_OUT 96.85788955688477 2.4419271045522506

Arm transmissions
~>cdsutils avg -s 10 C1:LSC-TRY_OUT C1:LSC-TRX_OUT
C1:LSC-TRY_OUT 1.0163812279701232 0.02008742269699207
C1:LSC-TRX_OUT 0.9274496018886567 0.027689954466601815

MICH fringe with ETM misaligned
See Attachment #1

LO-ITMX single bounce fringe with ITMY misaligned
See Attachment #2

[JC, Paco, Yuta]

It seems like PMC transmission (C1:PSL-PMC_PMCTRANSPD) dropped to ~0.68 from ~0.74 on Dec 27.
We tried to tweak PZT offset for PMC loop and input alignment to PMC, but PMC transmission didn't increased.
PSL laser temperature was also sweeped in the range 30.3 - 31.6 degC, but didn't help. The PSL temperature was reverted to original 30.61(1) degC.
Power measured at PSL output now is 893 mW (measured at our standard place shown in 40m/16672), which used to be 951 mW in June 2022 (40m/16886).
Power measured at PMC input (see attached photo) now is 1.18 W.

Next:
 - What was the previous PMC input power we had?
 - Sweep PSL laser temperature for larger range

[Paco, Yuta]

We recovered FPMI BHD, and sensitivity was estimated. High frequency sensitivity is improved by an order of magnitude compared with AS55 FPMI.
We also estimated the contribution from LO phase noise, and found that LO phase noise is not limiting the sensitivity.

Locking sequence:

1. Lock electronic FPMI

DARM:
 - 0.5 * POX11_I - 0.5 * POY11_I
 - DARM filter module, FM4,5 for acquisition, FM1,2,3,6,8 triggered, C1:LSC-DARM_GAIN = 0.015 (gain lowered from BHD FPMI in December (was 0.02) to have more gain margin)
 - Actuation on 0.5 * ETMX - 0.5 * ETMY
 - UGF ~ 150 Hz

CARM:
 - 0.5 * POX11_I + 0.5 * POY11_I
 - CARM filter module, FM4,5 for acquisition, FM1,2,3,6,8 triggered, C1:LSC-CARM_GAIN = 0.012
 - Actuation on -0.734 * MC2

2. Lock MICH

MICH:
 - 1.1 * REFL55_Q
 - MICH filter module, FM4,5,8 for acquisition, FM2,3,6 triggered, C1:LSC-MICH_GAIN = +10
 - Actuation on 0.5 * BS
 - UGF ~35 Hz (Attachment #1)

3. Hand over to real DARM/CARM

DARM:
 - 2.617 * AS55_Q

CARM:
 - 0.496 * REFL55_I
 - UGF ~200 Hz (Attachment #2)

4. Lock LO_PHASE

LO_PHASE:
 - 1 * BH55_Q (demod phase: -110 deg; this was chosen by hand to maximize fringe in BH55_Q when FPMI is locked. This seemed to make more robust LO_PHASE lock, compared with December)
 - FM5,8, C1:HPC-LO_PHASE_GAIN=-1
 - Acuation on 1 * LO1
 - UGF ~110 Hz (Attachment #3)

5. Hand over to BHD_DIFF

DARM:
 - -1.91 * BHD_DIFF (ratio between AS55_Q to BHD_DIFF was measured with a DARM line at 575.125 Hz, which was measured to be -0.455; BHD_DIFF was zeroed by balancing A and B before the measurement)
 - UGF ~150Hz (Attachment #4)

Measured sensing matrix:

/opt/rtcds/caltech/c1/Git/40m/scripts/CAL/SensingMatrix/ReadSensMat.ipynb

Sensing matrix with the following demodulation phases (counts/counts)
{'AS55': -168.5, 'REFL55': 92.32, 'BH55': -110.0}
Sensors       DARM @307.88 Hz           CARM @309.21 Hz            MICH @311.1 Hz        LO1 @315.17 Hz           
AS55_I       (-0.35+/-1.44)e-02        (+3.23+/-1.09)e-02        (-0.18+/-9.84)e-03    (+0.02+/-1.33)e-03    
AS55_Q       (-0.45+/-3.48)e-02        (-0.19+/-1.23)e-02        (+0.05+/-1.43)e-03    (-0.02+/-2.03)e-04    
REFL55_I     (+0.05+/-1.31)e-01        (+4.51+/-0.07)e-01        (-0.03+/-3.44)e-02    (+0.11+/-1.01)e-03    
REFL55_Q     (-0.39+/-4.74)e-04        (-1.36+/-0.50)e-03        (+0.16+/-2.80)e-04    (+0.10+/-4.92)e-05    
BH55_I       (-1.90+/-5.00)e-03        (+0.80+/-8.82)e-03        (-0.50+/-2.29)e-03    (-4.61+/-9.52)e-04    
BH55_Q       (-0.31+/-4.86)e-02        (-3.13+/-3.04)e-02        (-0.06+/-1.07)e-02    (-1.42+/-2.11)e-03    
BHDC_DIFF    (-5.56+/-0.21)e-02        (+0.10+/-1.64)e-02        (+0.10+/-1.40)e-03    (+0.77+/-2.27)e-04    
BHDC_SUM     (+1.75+/-3.90)e-03        (-1.22+/-3.22)e-03        (-0.35+/-1.52)e-03    (-0.36+/-6.42)e-04

Sensing matrix with the following demodulation phases (counts/m)
{'AS55': -168.5, 'REFL55': 92.32, 'BH55': -110.0}
Sensors       DARM @307.88 Hz           CARM @309.21 Hz           MICH @311.1 Hz        LO1 @315.17 Hz           
AS55_I       (-0.30+/-1.25)e+11        (+2.18+/-0.73)e+11       (-0.07+/-3.59)e+10    (+0.08+/-5.02)e+09    
AS55_Q       (-0.39+/-3.03)e+11        (-1.27+/-8.32)e+10       (+0.20+/-5.23)e+09    (-0.06+/-7.66)e+08    
REFL55_I     (+0.04+/-1.14)e+12        (+3.05+/-0.05)e+12       (-0.01+/-1.26)e+11    (+0.42+/-3.82)e+09    
REFL55_Q     (-0.34+/-4.12)e+09        (-9.20+/-3.37)e+09       (+0.06+/-1.02)e+09    (+0.04+/-1.86)e+08    
BH55_I       (-1.65+/-4.34)e+10        (+0.54+/-5.95)e+10       (-1.83+/-8.36)e+09    (-1.74+/-3.59)e+09    
BH55_Q       (-0.27+/-4.23)e+11        (-2.11+/-2.05)e+11       (-0.21+/-3.91)e+10    (-5.36+/-7.98)e+09    
BHDC_DIFF    (-4.83+/-0.18)e+11        (+0.07+/-1.11)e+11       (+0.35+/-5.09)e+09    (+2.90+/-8.56)e+08    
BHDC_SUM     (+1.52+/-3.39)e+10        (-0.82+/-2.17)e+10       (-1.29+/-5.56)e+09    (-0.14+/-2.42)e+09

NOTE that some of sensing matrix element (e.g. DARM to AS55_Q) is wrong, because of ill defined sign in C1:CAL-SENSMAT_XXXX_XXXX_AMPMON channels.

Locked GPS times:
 - 1357511737 to 1357513050
 - 1357513448 to 1357518188 (intentionally unlocked)

Sensitivity estimate:
 - See Attachment #5 and #6 (high frequency zoomed), dashed traces with darker colors are of AS55_Q FPMI from 40m/17369.
 - DARM_IN1 was calibrated using DARM content in BHD_DIFF, with 1 / (1.91 * 4.83e11 counts/m) = 1.086e-12 m/counts (which should be similar to DARM_IN1 calibration for AS55_Q, because we are balancing the error signals going to DARM_IN1, and it is as expected; see 40m/17369).
 - DARM_OUT calibration is the same as 40m/17369; ETM plant = 10.91e-9 / f^2 m/counts.
 - LO phase noise was estimated using BH55_Q, with the collowing calibration factor (BH55_Q calibrated into LO1 motion, into BHD_DIFF, and then into DARM).

1/5.36e9*2.9e8/4.83e11 = 1.15e-13 m/counts

 - Seismic noise was estimated using C1:IOO-MC_F_DQ, with the same calibration factor found in 40m/16975.
 - Dark noise was estimated using DARM_IN1 when the PSL shutter is closed.

Discussion:
 - Sensitivity below ~10 Hz is probably limited by seismic noise
 - Noise above 1 kHz might be limited by dark noise of BHD_DIFF, but not below 1 kHz. For AS55_Q FPMI, sensitivity above ~300 Hz was limited by dark noise.
 - LO phase noise is not limiting the sensitivity, from the estimated noise using BH55_Q.
 - Both BH55_Q and BHD_DIFF have funny structure like forest above ~100 Hz. This might be from suspensions in the AS path and LO path to BHD. It could also be that calibration lines for measuring sensing matrix were too much (BHD FPMI sensitivity was measured with calibration lines around ~310 Hz on).

Next:
 - Update the c1cal model to put correct signs to C1:CAL-SENSMAT_XXXX_XXXX_AMPMON channels.
 - Measure BHD FPMI sensitivity with calibration lines off.
 - Find better LO phase to improve the sensitivity.
 - Lock LO phase with RF+audio and RF44 and compare the sensitivity.
 - Move on to PRFPMI BHD.

[Paco, Yuta]

Gains in DARM are corrected to make it more calibration friendly.

DARM error signals
 - 0.19 * POX11_I - 0.19 * POY11_I
 - 1 * AS55_Q
 - -0.455 * BHD_DIFF (needs to be checked if LO phase is different)

DARM gain:
 - C1:LSC-DARM_GAIN = 0.04 (it was 0.015 to have UGF of ~150 Hz)

Online calibration:
 - FM2 for C1:CAL-DARM_CINV was turned on, which is a calibration for AS55_Q in FPMI; 1 / (3.64e11 counts/m) = 2.747e-12 m/counts (see sensing matrix below; consistent with 40m/17369).
 - FM2 for C1:CAL-DARM_A was updated to 10.91e-9 (40m/16977).
 - C1:CAL-DARM_W_OUT will be our calibrated FPMI displacement in meters. This is correct with BHD_DIFF locking, if the BHD_DIFF is balanced with AS55_Q before DARM_IN1.

FPMI sensitivity:
 - Attached plot shows the sensitivity of FPMI with AS55_Q and BHD_DIFF, plotted together with their dark noise.
 - The sensisitivty was measured with calibration lines off and notches off, which removed the forest of lines we saw on Jan 11 (40m/17392).

FPMI sensing matrix:
 - Attached is a screenshot of uncalibrated sensing matrix MEDM screen. Audio demodulation phase for DARM was tuned to have stable sign.
 - The following is calibrated sensing matrix measured today with FPMI locked with AS55_Q. BH55 and BHD_DIFF have large uncertainties because LO_PHASE locking is not stable.

Sensing matrix with the following demodulation phases (counts/m)
{'AS55': -168.5, 'REFL55': 92.32, 'BH55': -110.0}
Sensors       DARM @307.88 Hz           CARM @309.21 Hz           MICH @311.1 Hz           LO1 @315.17 Hz           
AS55_I       (-3.28+/-0.90)e+11 [90]    (-0.05+/-2.53)e+11 [0]    (+0.47+/-2.14)e+10 [0]    (-0.06+/-1.76)e+09 [0]    
AS55_Q       (-3.64+/-0.08)e+11 [90]    (-0.09+/-8.30)e+10 [0]    (-0.28+/-2.24)e+09 [0]    (+0.12+/-1.16)e+08 [0]    
REFL55_I       (+0.07+/-1.24)e+12 [90]    (+3.32+/-0.06)e+12 [0]    (-0.01+/-1.39)e+11 [0]    (+0.00+/-1.12)e+09 [0]    
REFL55_Q       (-0.98+/-2.46)e+09 [90]    (-4.68+/-2.04)e+09 [0]    (+0.08+/-1.00)e+09 [0]    (+1.73+/-5.69)e+07 [0]    
BH55_I       (-6.84+/-2.47)e+11 [90]    (+0.43+/-2.32)e+11 [0]    (-0.33+/-6.21)e+10 [0]    (-2.51+/-8.09)e+09 [0]    
BH55_Q       (+5.68+/-1.57)e+11 [90]    (-0.17+/-2.52)e+11 [0]    (-0.46+/-4.39)e+10 [0]    (+0.73+/-5.01)e+09 [0]    
BHDC_DIFF       (-2.48+/-3.47)e+11 [90]    (-0.09+/-1.99)e+11 [0]    (+1.56+/-4.11)e+10 [0]    (-0.49+/-6.78)e+09 [0]    
BHDC_SUM       (-2.12+/-1.84)e+10 [90]    (+0.35+/-1.44)e+10 [0]    (-1.30+/-4.18)e+09 [0]    (-0.03+/-8.17)e+08 [0]    

Current status:
 - After locking FPMI BHD to get the FPMI sensitivity today, we are struggling to re-lock again. LO_PHASE locking is glitchy today for some reason. To be investigated.
 

[Paco, Yuta]

After several hours of unattended IFO, we realigned the IFO and somehow BH55_Q error signal looked better, LO_PHASE locking was more robust, and we could lock FPMI BHD.

Sensing matrix:
 - Attached #1 is the sensing matrix with FPMI BHD locked, with LO_PHASE locked with BH55_Q using AS4, and below is the calibrated one.

Sensing matrix with the following demodulation phases (counts/m)
{'AS55': -168.5, 'REFL55': 92.32, 'BH55': -110.0}
Sensors       DARM @307.88 Hz           CARM @309.21 Hz           MICH @311.1 Hz           LO1 @315.17 Hz           
AS55_I       (+0.43+/-2.22)e+11 [90]    (+5.26+/-1.26)e+11 [0]    (-0.15+/-5.88)e+10 [0]    (+0.05+/-1.94)e+09 [0]    
AS55_Q       (-3.73+/-0.19)e+11 [90]    (-0.09+/-1.06)e+11 [0]    (+0.21+/-6.82)e+09 [0]    (-0.24+/-2.41)e+08 [0]    
REFL55_I       (+0.08+/-1.17)e+12 [90]    (+3.14+/-0.07)e+12 [0]    (-0.01+/-1.31)e+11 [0]    (+0.07+/-1.03)e+09 [0]    
REFL55_Q       (+0.51+/-1.09)e+10 [90]    (+1.80+/-2.35)e+10 [0]    (+0.01+/-1.44)e+09 [0]    (-0.79+/-5.41)e+07 [0]    
BH55_I       (+1.20+/-0.34)e+11 [90]    (+0.26+/-1.07)e+11 [0]    (-0.04+/-1.33)e+10 [0]    (-1.07+/-1.82)e+09 [0]    
BH55_Q       (+1.12+/-1.52)e+11 [90]    (-3.40+/-1.98)e+11 [0]    (-0.01+/-4.16)e+10 [0]    (-5.55+/-1.66)e+09 [0]    
BHDC_DIFF       (+6.95+/-0.25)e+11 [90]    (+0.01+/-1.60)e+11 [0]    (-0.70+/-4.49)e+09 [0]    (+4.24+/-4.14)e+08 [0]    
BHDC_SUM       (+7.06+/-2.64)e+10 [90]    (-0.32+/-3.24)e+10 [0]    (+0.05+/-4.26)e+09 [0]    (+1.30+/-6.33)e+08 [0]    

BH55_Q:
 - We are not sure why BH55_Q error signal got better. Peak to peak amptidue of BH55_Q when LO_PHASE is not locked is at around 800 (even if LO_PHASE locking cannot be stably locked), but when we couldn't lock LO_PHASE stably, it was noisier. This suggests that BHD alignment is not bad.
 - Attachment #2 shows the spectrum of BHDC_DIFF, BH55_Q (measured at LO_PHASE_IN1) and AS55_Q when FPMI is locked with AS55_Q, and LO_PHASE is locked with BH55_Q. Dashed darker curves are when we could not lock LO_PHASE stably. You can see broad noise in BH55_Q at around 60 Hz and its harmonics when it is noisy, but they are not present when LO_PHASE can be locked stably. Also, AS55_Q stays the same, which suggests that the cause is not in FPMI to AS55, but in BHD path to BH55.

Sensitivity comparison with different LO_PHASE locking actuator:
  - We compared FPMI sensitivity with LO_PHASE locked with LO1 (red) vs with AS4 (orange). Didn't change, as expected (Attachment #3).

Next:
 - Investigate what is causing noisy BH55. Is it FPMI alignment, BHD alignment, suspensions in LO/AS paths, or electronics?
 - Try locking LO_PHASE with an offset to BH55_Q to see if BHD_DIFF sensitivity to DARM changes, and to see if we can reduce BHD_DIFF dark noise contribution to FPMI sensitivity
 - Try BH55+audio and BH44 to lock LO_PHASE.
 - PRMI and PRFPMI

[Anchal, Paco, JC, Yuta]

We have started hardware installations for BH44 RF PD. 44 MHz LO generation and signal chain from IQ demodulator was checked successfully, but found that IQ demodulator is not orthogonal.

RF PD Interface:
 - We have unplugged Ch2 of the RFPD Interface (labeled "Special") to re-use it for BH44. Ch2 was used for "UNIDENTIFIED" RFPD. DB15 cable was routed to ITMY table and connected to BH44 RF PD (40m/17398) now sitting on the cover of ITMY table. See Attachment #1.
 - Finding a DB15 RF PD interface cable was easy because of the organization work!

44 MHz LO generation:
 - LO for BH44 was generate following the scheme proposed in 40m/17319.
 - 11 MHz LO from RF distributor labeled "+16 dBm" (measured to be 16.5 dBm) and 55 MHz LO labeled "SPARE 55" (measured to be 2.26 dBm) was mixed with a mixer ZFM-1H-S+ (using 11 MHz as LO, and 55 MHz highpass filtered with SHP-50+ as RF). The mixer output was lowpass filtered with SLP-50+, and amplified with ZFL-500LN+, which gave 8.07 dBm at 44 MHz. The second heighest peak was -11.53 dBm at 22 MHz, which seems low enough. See Attachment #2 for the photo of the setup.

IQ demodulation:
 - We have used unused IQ demodulator labeled "AS165" to use it for BH44. See Attachment #3.
 - We have quickly checked if the IQ demodulator is working or not with LO from BH55, PD input of 55 MHz generated using Moku to see I and Q outputs. The outputs are sine waves at frequency consistent with the difference between LO frequency and "PD input" frequency, and the phase was off as expected. Q output was ~4 dBm higher than I output.

Measured diff of BH44:
 - After CDS modifications where done, BH44 IQ demodulator was tested by using 44 MHz LO generated in a method mentioned above, and injecting 11.066195 * 4 MHz signal from Moku as PD input. This gave ~75 Hz signal in C1:LSC-BH44_I and C1:LSC-BH44_Q.
 - With 0dB whitening gain and whitening/unwhitening filters off, gain imbalance was measured to be Q/I=137.04/62.49=2.19, and measured phase difference to be PHASE_D=27.21 deg (see Attachment #4; gpstime=1358051213).
 - With 0dB whitening gain and whitening/unwhitening filters on, gain imbalance was measured to be Q/I=138.44/63.21=2.19, and measured phase difference to be PHASE_D=26.95 deg (see Attachment #5; gpstime=1358051325138). This is consistent with whitening/unwhitening off, and noise is smaller, which mean whitening/unwhitening filters are probably working fine.
 - IQ demodulator board might be not working properly, as I and Q signals are not quite orthogonal.

Model changes:

• We modified c1lsc, c1hpc and c1cal model for BH44.
• BH44 ADC pins were identified and connected for RFPD phase rotator.
• The signals are sent to c1hpc through IPC where BH44 is now available for feedback loops in single and dual demodulation.
• The whitening filter controls and anti-aliasing filter enable buttons were created in c1iscaux slow machine db files.
• MEDM screens are updated accordingly (see Attachment #6).

Next:
  - Use different IQ demodulator board that has better IQ orthogonality.
  - Connect BH44 RF PD and use 44 MHz test input to check the signal chain.
  - Install BH44 RF PD optical path.
  - Try locking LO_PHASE with BH44.

[JC, Paco, Anchal, Yuta]

- We installed the new BH44 beam path and RFPD.

- JC installed the new beam path for the LO/AS camera.

- We succeeded in locking LO Phase with BH44_Q_ERR, but didn't attempt FPMI BH44 because we noted large 60 Hz harmonics in most of our RF error signals.

BH44 RPFD/Camera installation:

- We picked off LO/AS beam path previously going to the camera, and installed a Y1 (45 deg, s-pol) mirror, a f=150mm lens and the RFPD (Attachment #1). We initially tested it using the incandescent light from a flashlight and then aligned the beam, we also made sure it's not saturating.

- Using the spurious transmission from the mirror mentioned above, we steered a new beam path for the camera and realigned it using another short focal length lens (f ~ 100 mm).

LO Phase control:

- We increased the whitening gain from 0 dB to 42 dB for both C1:LSC-BH44_I and C1:LSC-BH44_Q, and zeroed the offsets. Even before this step we could see a fringe from BH44, which is quite promising!

- After alignment was recovered on the LO/AS path, we succeeded in locking the single bounce (ITMX) LO phase using BH44_Q. Here the configuration was FM4, FM5 and a gain of ~ 5 * 0.5 = 2.5 (to match the typical BH55_Q error point).

- While BH44_Q was used to control the LO phase, we saw the BH55_Q was not zero but actually almost at max fringe value (see Attachment #2). This implies the BH44_Q is indeed orthogonal to BH55_Q with respect to the LO Phase!

FPMI lock:

- We locked electronic FPMI but noted a  large 60 Hz + harmonics component in the RF error signals including AS55, BH55, REFL55, and BH44 (see Attachment #3). We could hand off to FPMI and even locked the LO phase with BH44_Q, but we were not sure the BHD_DIFF error signal was fit for handoff to FPMI BHD. Therefore we stopped here.

60 Hz + harmonics:

- We did a quick investigation around the areas we have been working in the lab to see if we had introduced this noise in any obvious way. First we checked the new amplifier for the 44 MHz LO, we briefly removed its power but the 60 Hz noise remained. Then we checked the AP table, but nothing had really changed there. We also disconnected and removed the rolling cart with the marconi and other instruments from the LSC rack. Finally, we turned all the lights down. None of these quick fixes changed the amount of noise.

- We also tried looking at these error signals under different IFO alignment and feedback configurations. We always see the noise in the AS55 and REFL55 quadratures, but not in BH44, BH55 or BHD_DIFF unless MICH is locked.

Next steps:

- Investigate more into 60 Hz noise, why? where?

- Measure sensing matrix with LO Phase locked with BH44 and BH55 to make comparison.

- FPMI-BH44