Kiwamu, Rana Elli
Last night we comissioned Dhard WFS. We were engaging the DHard WFS (H1:ASC-DHARD_P and H1:ASC-DHARD_Y) after the DARM handover from DIFF to AS45 had taken place.
The WFS use signals from AS_B_RF45_PIT and YAW, and we set H1:ASC-AS_A_RF45_Q_YAW_GAIN and H1:ASC-AS_A_RF45_Q_PIT_GAIN to 1.
We were engaging the WFS with the following settings: (This is in the guardian but hasn't been tested yet.)
ezca['ASC-AS_B_RF45_WHITEN_GAIN'] = 3 (we changed this from 0 to 3dB)
ezca['ASC-INMATRIX_P_8_6'] = 1
ezca['ASC-INMATRIX_Y_8_6'] = 1
ezca['ASC-DHARD_P_GAIN'] = 0.003
ezca['ASC-DHARD_Y_GAIN'] = 0.003
ezca['ASC-OUTMATRIX_P_7_8'] = 1
ezca['ASC-OUTMATRIX_P_8_8'] = -1
ezca['ASC-OUTMATRIX_Y_7_8'] = 1
ezca['ASC-OUTMATRIX_Y_8_8'] = 1
ezca.switch('ASC-DHARD_P', 'FM1', 'FM2', 'FM3','FM6','FM7', 'ON')
ezca.switch('ASC-DHARD_Y', 'FM1', 'FM2', 'FM3','FM6','FM7', 'ON')
ezca.switch('ASC-DHARD_P', 'INPUT', 'ON')
ezca.switch('ASC-DHARD_Y', 'INPUT', 'ON')
Where FM1 is 10dB gain, FM2 is 20dB gain, FM3 is z1p0 integrator, FM6 is an approximate inverse of the suspension response with a sneaky 1e-5 gain thrown in (zpk([1],[50+i*86.6025;50-i*86.6025],1,"n")gain(1e-05)), FM7 is a 1Hz low pass filter.
The pressure sensor noise is not a function of the cable run or the satellite amplifier.
See attached--60 days of sensors in PSI. The channels are in order: 1) The low pressure return at BSC2, 2) the high pressure supply at BSC2,[these two signal are amplified at the BSC and then run ~120ft to the ADC,] 3) The 'CALC' difference of the previous two--this channel is what the servo controls, 4) the Main Supply line pressure after all the Pump Stations [In parallel,] 5 & 6) the final pressure on two Pump Stations before joining the Main Supply Line.
Again, the first two channels are the remote sensors traveling ~120ft to the ADC. Channel 3 is the difference of the first two. The last channels travel <20ft to the ADC.
I've zoomed the graphs to be all the same and the noise level on the far sensors is not more noisy than the close sensors. Channel1, the return pressure signal in fact appears to be the least noisy. I don't see the difference there to be too signiicant but maybe it is.
One thing that is clear is seen in channel 3. Before the large shift to 70psi on 14 January, that calc record was servoing on the Main Supply Pressure, Channel 4. Now that the Differential channel is actually a difference of two channels, its noise has increased.
Finally, the actual quieter nature of these channels is just seen at the beginning of the traces on Channel 3 4 & 5. Before this time, the remote sensors from BSC2 were not available. This shift in the noise occurred whilst I was connecting up these channels. Somehow at that time when I was shifting cables around to dress them, a cleaner environment was lost. I have tried to put everything back to as it was, that is, reconnecting the unused terminal to the servo computer etc but no luck. Maybe the ground is floating around or something but it isn't clear.
Further, see the second plot. Here is 65 days of the differential pressure and the output of the servo driving the motor speed. This VOUT would be flatlined when the servo is not on. My clear recollection of this 'quiet' period is that the VOUT might change 1 digit every few seconds, not several counts every second.
New code was loaded into h1ecatc1 which required a restart of the associated plcs.
Ed and Kiwamu,
In the morning of today and yesterday, we went through all the initial aligment steps. In the course of the process, we found that there were a few settings that were hidden and not set correctly. So we made a couplle of modifications on the guardians and filter in order to automate them:
The attached screenshot shows the current good settings in order to go through the alignment sequence. This configuration should be automatically achieved by the guardians.
Some additional notes:
We found two failure modes which we did not get a chance to fix the guardians. I will have a close look at them tomorrow morning or next week.
Kiwamu, Elli, Alexa, Evan, Rana, Daniel, Sheila
Today we were able to reach CARM offsets around 30 pm.
We transitioned DARM to AS45Q, at a CARM offset where sqrt(TRX+TRY) was -7, we then normalized the signal by sqrt(TRX) (with a factor of 0.23). One important step in getting there was to implement the ezca servo that adjusts the ALS DIFF offset to bring AS45Q into the linear range. We are now using that both at a CARM offset of 1 (in SQRT TRX+TRY), and after we transition to the QPDs. We then change the DARM loww pass filter from a 33Hz low pass to a 80 Hz low pass, to get better phase margin since we are no longer saturating the ESD with ALS noise. We did this transition several times sucessfully. As Rana mentioned in alog 16334, we installed an ND1 filter on ASAIR A. Since this we haven't transitioned to RF DARM again (for reasons that seem to be unrelated to the ND filter), so we will need to check the gain before we transition again.
After making this transition, Elli, Kiwamu and Rana worked on the DHARD WFS, which we turned on to reduce the fuctuations in AS DC. This allowed us to go to CARM offsets -25 in sqrt(TRX+TRY), which is about half of the total power we expect in the arms. If you assume that the recyling gain is 30 (we don't really know) it is something like 30 pm. In Refl DC we saw the power drop by about 20%. We saw that the linearized REFL 9 I signal had turned over, and that without linearization the signal had reached the peak. We made a few attempts to transition, and we were able to turn down the gain of the TR CARM signal to 50% of the nominal and turn the REFL 9 signal to what we think the nominal gain should be (-100 in the input matrix). We lost the lock when we turned off the TR CARM signal. Our next plan was to leave the TR CARM engaged with reduced gain, and keep reducing the CARM offset.
However, we have been having a hard time locking in the last few hours. We think that it might help to try transitioning to DARM to RF a little sooner, so that we could use WFS.
The sequence that was working earlier this afternoon is in guardian up to the RF DARM transition, although this might need to be re-worked. We have added but not tested a state for the DHARD WFS.
Attached is a screenshot of our striptool durring the sequence, which was all handled by the guardian this time.
As we were about to leave we had a nice stable lock. We were able to transition to RF DARM again at -7.0 cts CARM offset after having adjusted for the ND filters. We now have a +20dB filter in ASAIR_RF45Q, and the input matrix is now 25. We turned on FM4(z4^2:p1^2) in DARM loop which significantly helped the DARM noise. We then proceeded to adjust the CARM offset to -20 cnts. At this point we transitioned TR_RELF9 to 100% with TR CARM at 50%. We were able to reduce the CARM offset to zero, but this only lasted for about a second or so. We never fully turned off TR CARM, but we think it has a zero slope here since we are at zero offset. More tomorrow when we are awake ....
Lock loss time: 09:58:40 UTC Jan 29th
Great work!
Assuming I am looking at the right lock attempt, (data attached starting at 09:48:00 UTC) it seems that REFL DC is only ~30% less than at the beginning of the sequence when you transition to RF. There should be room to get closer. P.S: For comparison, trend of powers with "lossy" arms is here. The build up in the arms for same relative REFL DC power was about a factor of 3 lower (by eye numbers).
Daniel and Rana have mentioned that optical torques may become significant as we come in to resonance.
For 10 mm of miscentering and 46 kW of circulating arm power (at 0 pm of CARM offset), we get a torque of 3×10−6 N m. I estimate the stiffness constants of each pendulum to be 4.9 N m for pitch and 6.5 N m for yaw (a better estimate could be made using the actual suspension models). This means that the static misalignments induced by the radiation torque could be as large as 1 μrad. The attached code computes the torsional stiffnesses of the pendula.
As a next step, we might also consider the stiffness of the optical springs using, e.g., eqs. 31 in the paper by Sidles and Sigg. At 46 kW of circulating power, we get 15 N m for the major mode and −0.6 N m for the minor mode.
[Edit: Also, Kiwamu has pointed out an error in the expression for the moment of inertia for the test masses. This has been fixed in this entry and in the attached code.]
Here are some oplev trends from last night's final lock attempt.
The drop in the buildup of POP18 seems correlated with a drift in ETMX pitch (0.3 μrad), and to a lesser extent BS pitch (0.2 μrad), SR3 pitch (0.4 μrad) and ITMY pitch (0.2 μrad). There may also be some drift in PR3 yaw and pitch (≈0.1 μrad). All of these drifts happen on time scales much slower than the change in TRX buildup, which supports the idea that these are thermal drifts induced, e.g., by wire heating.
For the record here are some lock loss times from last night:
Early in the evening we were trying to transition DARM from ALS_DIFF to ASAIR_A_RF45_Q and the lock dropped at the following times:
Jan 28 20:46:40 UTC, Jan 28 21:06:32 UTC, Jan 28 21:06:17 UTC, Jan 28 22:55:00 UTC.
Speculating from the lock loss plots, we think that DARM noise causes a big spike in light leaking out of the AS port. This causes the power on the LSC-TR_X/Y_QPDs controlling CARM to fluctatue enough that CARM drops lock. Running the ASAIR centering servo should help minimise big spikes at ASAIR_A_LF. Once we were able to transition DARM to RF this type of lock loss stopped happening.
-----------------------
Here are some lock losses from after the transition DARM to RF. The cause of these lock losses remained unclear. MICH, PRCL and SRCL were ringing up signals at various frequencies (4Hz-20Hz) but this changed from lock loss to lock loss. Again there are big spikes in ASAIR_A_LF right before the lock loss. ETMy alignment needed frequent touching up.
Jan 29 00:28:31 UTC, Jan 29 00:49:41 UTC, Jan 29 05:34:23 UTC.
--------------------
Later in the evening we were having a hard time locking. Again we were loosing lock before DARM transition to RF. Again there are big spikes in ASAIR_A_LF, probably caused by DARM motion.
Jan 29 07:35:25 UTC, Jan 29 07:35:25 UTC, Jan 29 08:24:26 UTC, Jan 29 08:40:40 UTC
Here is a screen shot of the CARM offset reduction from earlier in the evening, when the alingment must hve been slightly better. Although we didn't reduce the CARM offset, and were locked on TR CARM, we had a recylcing gain of about 10.2. Also, some of the signal from POP18 and POP90 is rotated into the Q phase as CARM offset is reduced.
Peter, Matt, Lisa For the records, we had this theory that if the f_1 was tuned such as to make the 2f_1 resonant in the arms, the beat between the 2f_1 and the carrier in the recycling cavity could be responsible for the decay in POP18. Looking in the L1/H1 logs and MEDM screens, we arrived at the conclusion that in H1, given the arm length of 3994.4704 m and f_1 = 9.100230 MHz, the offset from resonance for the 2f_1 should be 380 Hz. In L1 the offset for the 2f_1 from resonance is 500 Hz (reported here), as nominal.
Alexa, Rana, Ellie
We noticed that one of our locklosses during the CARM offset reduction (after we're transitioned onto AS45 -> DARM) seems to be due to ADC sauration in ASAIR 45 so we put in an ND1.0 screwed onto ASAIR_A.
In the attached plot you see that ASAIR_A_45 goes to 0.0003. Since we have a -160 dB digital filter between the ADC and this point, this corresponds to ~32000 counts at the ADC.
OTOH, REFL9 seems to be pretty puny (~500 counts with half of the full buildup), so we're turning up the gain on there by ~10 and turning on the Whitening filter to see if it helps the CARM noise be low.
To get bigger plots in dataviewer (i.e. getting the plot window to auto-scale) for both playback and 'real-time' plots, you can now just set up an option in your GRACE environment:
1) In your user directory make a directory called .grace: mkdir .grace 2) Change to that directory: cd .grace 3) Copy my gracerc file into your directory: cp /ligo/home/rana.adhikari/.grace/gracerc.user .
Now, when you restart dataviewer, you will have autoscaling.
BTW, this is described somewhat in CDS Bugzilla #377 from Tobin Fricke, Jim Batch, & Keith Thorne.
Note that the Realtime display suffers a bit when the "PAGE LAYOUT FREE" option is used to autoscale the plot windows. As the attached screen shot shows, you no longer get the full plot in the window. It takes some inconvenient fiddling to get the full plot displayed.
true enough - this is usually fixed by hitting stop and then start after first resizing the window. It usually undoes the cutoff plots, but sometimes it just leaves it bad...
Completed the BSC file modifications today like the HAM-ISI yesterday and HEPIs last Thursday. Now all SEIs are monitored by SDF with Guardian doing the rest. There remain red issues (5 platforms) I did not clear to remind us of ongoing questions & issues to address.
Summary:
QE and transimpedance combined for Baffle PD1 and PD4 on all ITMs and ETMs agree with each other within 10%. In the table below, QE/transimpedance combined for each PD is normalized such that the average of all diodes becomes 1.
X | Y | |||
PD1 | PD4 | PD1 | PD4 | |
ETM | 1.01 | 1.00 | 1.09 | 1.03 |
ITM | 0.94 | 0.98 | 1.00 | 0.95 |
I haven't done any error analysis, but the biggest error should be the systematic in ITM-ETM comparison that was introduced by an assumption that IR mode matching is perfect.
PD1-PD4 comparison on the same baffle (e.g. ETMX PD1 and PD4) should be quite good.
Details:
This is yet another baffle health check.
The other day I used some short time available to misalign TMS and steer MMT3 to point IR beam on ETM baffle PD4, first on ETMX and then on ETMY. This could be used to make ETMX(PD4)/ETMY(PD4) comparison, assuming that the IR straight shot beam size is the same on ETMX and ETMY.
Then I proceeded to align MMT3, align ETMs to point the beam to PD4 on ITM baffles. Since ETM reflectivity is basically 1, and since the mode matching cannot be as bad as green, we can use these measurements combined with nominal IR beam size to make ETM(PD4)/ITM(PD4) comparison.
Also, when people do initial alignment they are basically doing PD1/PD4 comparison on a single baffle using green beam. (You cannot make comparison across different baffles using green because of the beam size difference described in my previous alog.)
So, stitching all these together, if we believe that the QE variation of diodes for green is the same as that for IR, we can map out the QE/transimpedance combined for all baffle diodes.
In the table below of the raw data, photo current in mA for straight shot measurements of green on PD1 and PD4 as well as IR on PD4 on all arm cavity baffles are shown.
PD1 green | PD4 green | PD4 IR | |
EX | 0.0453 | 0.0448 | 0.0127 |
EY | 0.0111 | 0.0104 | 0.0130 |
IX | 0.123 | 0.129 | 0.0176 |
IY | 0.0243 | 0.0230 | 0.0170 |
Green data can be used to make PD1/PD4 ratio on each of the optics.
IR ETM data is used to make PD4 ETMX/ETMY ratio.
IR ITM data, combined with the nominal IR beam radius on ITM and ETM (53.4 and 63.5 mm, respectively), is used to make PD4 ETMX/ITMX and ETMY/ITMY ratio.
After some math we get the table at the top of this entry.
06:45 Karen into the LVEA
07:00 Cris into the LVEA
08:00 set OIB to 'Commissioning'
08:05 unlocked DRMI to start practicing locking Green light to arms (ALS) No luck. Later learned the proper way to unlock DRMI. Thanks Kiwamu
08:30 Morning meeting
08:45 Started working with Kiwamu on getting my first arm locks. SUCCESS!
09:09 Betsy and Mitch to W bay
09:22 Jeff B to LVEA
09:29 Mike L into LVEA
09:30 Jeff B out of LVEA
09:49 Mike L out of LVEA
09:52 McCarthy to EX
10:54 Betsy and MItch out of LVEA
10:55 Caterers on site
11:03 McCarthy back from EX
11:07 Fil to MY
11:43 Fill back from MY
12:15 Vanessa to MY
12:24 S&K Electricians to EY
The temporary barrier originally installed for the DCS construction has been removed. Mechanical is 98% complete with remaining items being plenum grills and dampers, start up of the HVAC system by a factory rep, installation of cable trays and some touch up painting. Electrical is on the order of 70% complete.
Like the quad model scripts (LHO 16126), I updated the generate_Triple_Model_Production.m and generate_Double_Model_Production.m functions so that you can specify a foton filter for the top mass damping loops.
The generate scripts is in
.../SusSVN/sus/trunk/Common/MatlabTools/TripleModel_Production
and
.../SusSVN/sus/trunk/Common/MatlabTools/TripleModel_Production
If youhaven't done so in the past week, you will also need to svn up
.../SusSVN/sus/trunk/Common/MatlabTools
since this is where the foton file reading functions are located (imported from the SeiSVN)
You can still specify the usual .mat struct file as before. The generate script looks for the .txt extension to determine if you are sending it a foton file.
Here is an example of how to make a model with damping filters read from foton:
hsts_model = generate_Triple_Model_Production(frequency_vector_for_plots, 'hstsopt_metal', [], 0, 1, '/opt/rtcds/lho/h1/chans/H1SUSMC1.txt')
IMPORTANT: The foton file does not know which filter modules are engaged. This information is coded into the function with a variable called medm_engaged_modules, which must be updated at the time of running the script. The script will output to the command line which modules are being used for each filter, and where to change in case you forget. For a variety of reasons I thought it more convenient from a user-point-of-view to have this coded into the function rather than as another input to the function, however this could be modified. This is the same in the quad script.
These functions have more instructions commented into their headers.
Just have the single sus scripts to go.
Summary:
Green straight shot beam size on ITMX is a factor of 0.87 of nominal (i.e. too small) while ITMY a factor of 1.2 too large.
The beam radius ratio is (X/Y)=0.72, which makes the ITMX baffle PD output current to be a factor of 2 larger than that of Y when you hit these PDs using a straight shot beam.
There used to be a factor of 2 to 2.5 unexplained difference between the photocurrent in ITMX and ITMY baffle diodes during initial green alignment, but the beam size explains a factor of 2, so the remaining factor is about 1.2 or so, which I don't care for now.
Background:
This is one of those health check type stuff. Baffle PDs could be used for scattering measurement but the health of some of those PDs were in question.
One of the suspicions came from the fact that, when looking at green straight shot beam during full initial alignment, PD current from X arm baffle is much larger than Y arm. A part of it comes from the laser power (there's about a factor of 2 or so difference) but the baffle PD current is a factor of 4 to 5-ish different, so there always was a factor of 2 to 2.5 unexplained.
Details:
I first used dither align script to point the TMS beam to PD1 and then PD4 of the ITM baffle. Since the horizontal distance between PD1 and PD4 is 11.3", the script allows us to calibrate TMS alignment slider.
Immediately after the script finished, I pointed the beam to PD4 center, move TMS in YAW in one direction by more than the beam radius, and ran another script to move TMS in YAW, wait and measure PD4 current.
Attached is the YAW scan data (circles and crosses) as well as the fit to Gaussian profile assuming the same radius in PIT and YAW:
current = ofst+ A * 2/pi/w^2 * exp(-2*(X-X0)^2/w^2)
where w is the beam radius and A is an overall factor.
beam radius [m] | A [mA m^2] | QPD_SUM | |
IX PD4 | 0.033 (0.038 nominal) | 3.8E-7 | 46180 |
IY PD4 | 0.046 (0.038 nominal) | 2.7E-7 | 26124 |
X/Y | 0.715 | 1.4 | 1.77 |
The last column in the above table shows QPDA_SUM+QPDB_SUM. If everything makes sense, and if the beam was at the same height as the PD center during the scan, QPD_SUM ratio should agree with A ratio, but apparently it doesn't at 20%-ish level. This is good enough because the suspicion was that there was something grossly wrong about the baffle PDs.
BTW we can do the same thing for IR beam, scanning MMT3, to assess the IR matching to the arms if we want to (or better yet, do the spiral scan to see both PIT and YAW).
Looks like lock was from about 9:14 UTC to 16:09 UTC. For posterity.
Also for posterity: - This is a DRMI lock stretch. - IMC WFS DOF4 is OFF - No corner station WFS engaged. - ISS second loop is *OFF*. - 10 [W] input reqested. - SEI was in the most recent nominal configuration -- HPI Pump Servo ON, Sensor Correction to ISI XY and HEPI Z, (no ST0-1 Feed Forward yet), Using 01_28 blends on the HAM ISIs, "45 [mHz]" X&Y, LLO blends on BSC ISIs. Great for offline, data-mining studies of - Gain Problems with STS2B - HAM3 0.6 [Hz] feature - Coherence with HPI Pump Pressure. - DAC major-carry transition glitches. - Cavity performance with respect to SEI performance. - Lock Loss statistics Among other things...