Robert, Sheila, Evan, Gabriele
I tried to look at one of Robert's injections from yesterday, and we noticed a dangerous bug, which had previously been reported by Annamaria and Robert 20410. This is also the subject of https://bugzilla.ligo-wa.caltech.edu/bugzilla3/show_bug.cgi?id=804
When we changed the Stop frequency on the template, without changing anything else, the noise in DARM changes.
This means we can't look at ISI, ASC, PEM, or SUS channels at the same time as DARM channels and get a proper representation of the DARM noise, which is what we need to be doing right now to improve our low frequency noise. Can we trust coherence measurements between channels that have different sampling rates?
This is not the same problem as reported by Robert and Keita alog 22094
people have looked at the DTT manual and speculate that this could be because of the aggressive whitening on this channel, and the fact that DTT downsmaples before taking the spectrum.
If there is no near term prospect for fixing the problem in DTT, then we would want to have less aggressive whitening for CAL_DELTA_L_EXTERNAL
I spent a little time looking into this and added some details to the bug report. As you said, it seems to be an issue of high frequency noise leaking through the downsampling filter in DTT.
Until this gets fixed, any reason you can't use DARM_IN1 instead of DELTAL_EXTERNAL as your DARM channel? It's better whitened, so it doesn't suffer from this problem.
The dynamic range issue in the whitened channel can be improved by switching to five zeros at 0.3 Hz and five poles at 30 Hz.
The current whitening settings (five zeros at 1 Hz, five poles at 100 Hz) produce more than 70 dB of variation from 10 Hz to 8 kHz, and 130 dB of variation from 0.05 Hz to 10 Hz.
The new whitening settings can give less than 30 dB of variation from 10 Hz to 8 kHz, and 90 dB of variation from 0.05 Hz to 10 Hz.
We could also use 6 zeros at 0.3 Hz and 6 poles at 30 Hz, which would give 30 dB of variation from 10 Hz to 8 kHz, and 66 dB of variation from 0.05 Hz to 10 Hz.
The 6x p/z solution was implemented: LHO#25778
At low frequency, DARM is coherent with many ASC and LSC channels. However, ASC channels are not uncorrelated to each other, and the same is true for DRMI LSC channels. So if one wants to use coherence to estimate how much of DARM is explanable with coupling from a set of auxliary channels, it is important to avoid over counting due to the cross coheence between the aux signals.
To cope with this, I wrote a multi-channel cohrence code, which takes into account the cohrence between the aux channels and estimates the total contribution of a set of correlated aux channels to a target channel. The attached script multicoherence.m does this in MATLAB.
I computed this multi-channel coherence (and the projection of noise into DARM) for DRMI LSC signals (plot 1) and all ASC signals (plot 2). In each of those plots, the top panel shows the total coherence of DARM with the set of channels being considered. This coherence is defined as (PSD_DARM - PSD_LINEAR_COMBINATION_AUX) / PSD_DARM and so gives an idea of how well the combination of aux channels can explain the noise we see in DARM. In each plot, the bottom panel shows a coherence-baed projection of the aux channels into DARM.
As expected, we see that the low frequency region of DARM (<20 Hz) is dominated by the sum of all ASC noises. Similarly, the LSC noises are not far from the measured sensitivity (as already known, we need a MICH FF retuning).
I attached a brief not describing how the multi-channel coherence is computed. This kind of computation is not new material, it's been around for a long time.
Reset counters in attached screenshot.
TITLE: 01/19 [DAY Shift]: 16:00-00:00 UTC (08:00-16:00 PDT), all times posted in UTC STATE OF H1: Left IFO at DC readout for Jenne to run tests. Jim's ISI blend filter changes are on all of the BSC chambers. SHIFT SUMMARY: Received SNEWS test alert at ~ 09:00 AM PDT (see attached). The verbal alarms script reported it as a GRB. SUPPORT: Jim W., Jenne, Sheila INCOMING OPERATOR: Ed ACTIVITY LOG: 15:45 UTC Richard set observing bit to commissioning 15:45 UTC Ken to end X air handler room to look at heater wiring 15:50 UTC Chris to X2 for beam tube sealing 16:11 UTC Chistina and Karen to end stations for cleaning 16:28 UTC Joe D. to LVEA to charge forklift batteries 16:30 UTC Ryan B. patching alog 16:47 UTC Filiberto pulling cables in LVEA (WP 5681) 16:56 UTC Kyle to mid Y (alog 25024) 17:05 UTC Karen and Christina from end Y to end X 17:12 UTC LN2 delivery 17:20 UTC Hugh checking HEPI fluid pump levels at end stations 17:26 UTC Joe D. done 17:42 UTC Karen and Christina done 17:45 UTC Jim W. loaded filter for ETMY ISI, filter not on 17:52 UTC Paradise water delivery 17:59 UTC Hugh done at end stations, going to mechanical room 18:04 UTC Ryan B. patching alog, bugzilla, svn 18:08 UTC Hugh done in mechanical room 18:11 UTC US Linen on site 18:14 UTC Ryan B. done 18:14 UTC Jeff B. to end X, end Y and mechanical room to attach fittings to dust monitor vacuum pumps 18:14 UTC Jason resetting PSL WD (alog 25021) 18:19 UTC Lock loss (reason unknown) 18:23 UTC Jason done, going to LVEA to retrieve equipment 18:26 UTC Jeff K. starting charge measurements on ETMX and ETMY 18:29 UTC Jason back from LVEA 18:41 UTC Kyle back from mid Y 18:41 UTC Karen to LVEA to clean 18:54 UTC LN2 truck leaving 18:56 UTC Joe D. to LVEA to unplug forklift batteries 19:08 UTC Jeff B. done 19:08 UTC US Linen leaving site 19:16 UTC Joe D. done 19:25 UTC Coca Cola truck on site 19:26 UTC Dan M. replacing faulty controller on gds-h1 in MSR (WP 5682) 19:42 UTC Jason to LVEA west bay to retrieve part 19:42 UTC Jim W. setting new filters on all BSCs 19:55 UTC Filiberto done for lunch 19:55 UTC Jim W. setting ETMX ISI to damped 19:57 UTC Dan M. done 19:57 UTC Jason done 20:29 UTC Starting lock acquisition to check Jim W.'s change to blend filters Lock loss on engaging BS ISI stage 2 21:13 UTC Locked at DC readout with Jim W.'s filter change. Jenne starting tests. 22:04 UTC Kyle to mid Y 22:07 UTC Filiberto and Ed to LVEA to pull cables Lock loss 22:33 UTC Kyle back 23:09 UTC Filiberto done 23:11 UTC Ken to midY to check grounding 23:33 UTC Ken done
JeffreyK, DarkhanT,
Summary
After taking another set of DARM OLG TF and PCALY to DARM TF measurements on Jan 7, 2016 (see LHO alog 24756) we compared the DARM actuation, sensing and OLG TF measurements against their models (the most recent similar analysis report was posted in LHO alog 24569).
As in the previous similar analysis we see that residuals of the actuation, the sensing and the OLG TF measurements vs. κC and κtst corrected DARM model are smaller compared to residuals against not-kappa-corrected DARM model.
This time we also produced residuals against the full kappa corrected DARM models (κpu, κtst, κC and fc). Sensing function residuals at higher frequencies and actuation function residualsat frequencies close to UGF,as we exptected, became smaller (closer to 1.0).
Details
We produced a comparison plots of DARM OLG TF and PCALY to DARM TF measurements taken during O1. For each of the TF measurements we also took "kappas" (DARM temporal variations) calculated from calibration lines close to times of the TF measurements.
For measurements taken on Jan 7, 2016 the kappa values used in the analysis are 30 minute average values (starting at GPS 1136215829) recorded about two hours before the TF measurements in H1_HOFT_C00 frames. Values were taken only at times when 0:HOFT_OK_BIT in online H1:GDS-CALIB_STATE_VECTOR was set.
par.kappa_tst = 1.059039;
par.kappa_pu = 1.019854;
par.kappa_A = 1.041371;
par.kappa_C = 1.002553;
par.fc = 335.351749; % [Hz]
Coherence of DARM OLG TF measurements at frequencies below crossover, ~30 Hz, is above 0.98, but an accuary of the DARM acutation function model at these frequencies is not very good, it can be seen on the DARM OLG TF measurement vs. all-kappa-corrected model residuals plot (page 9 of the attachment).
The parameter files for these measurements were committed to CalSVN at
CalSVN/aligocalibration/trunk/Runs/PostO1/H1/Scripts/DARMOLGTFs/H1DARMparams_1136225488.m
CalSVN/aligocalibration/trunk/Runs/PostO1/H1/Scripts/DARMOLGTFs/H1DARMparams_1136225488_kappa_corr.m
CalSVN/aligocalibration/trunk/Runs/PostO1/H1/Scripts/DARMOLGTFs/H1DARMparams_1136225488_kappa_corr_full.m
The comparison script and plots were committed to SVN at
CalSVN/aligocalibration/trunk/Runs/PostO1/H1/Scripts/DARMOLGTFs/CompareDARMOLGTFs_O1andPostO1.m
CalSVN/aligocalibration/trunk/Runs/PostO1/H1/Scripts/DARMOLGTFs/CompareDARMOLGTFs_O1andPostO1_full.m
CalSVN/aligocalibration/trunk/Runs/PostO1/H1/Results/DARMOLGTFs/2016-01-07_CompareDARMOLGTFs_O1andPostO1.pdf
Note that DCS/LDAS continues to archive raw (H1_R and L1_R) and aggregated hoft (H1_HOFT_C00 and L1_HOFT_C00) data 24/7. However, DCS/LDAS switched the directory location of the data from O1 to postO1 today at this time: 1137258496 == Jan 19 2016 09:07:59 PST == Jan 19 2016 11:07:59 CST == Jan 19 2016 17:07:59 UTC This switch will be transparent to users using gw_data_find or NDS2 to find the data, but in the directory path to the data, O1->postO1 at the above time.
Hugh worked on a script looking at the ISI T240 centering in alog 24902. I noticed on that list that ETMX had one T240 that had all 3 legs mis-centered. Eager for any straw to grasp on the ETMX ISI ring up, I took a couple of minutes to recenter all of ETMX's T240s. First I repeated Hugh's measurement, using his script from 24902, then I took the ISI down to Damped and opened all three of the T240 filter banks from the overview. I then clicked on FM6 (labeled AutoZ) on the first (X) bank for each T240 for several seconds. This triggers the T240's Auto-zero servo, which centers up the masses. I watched the readouts on the ISI overview and waited for them to settle down to below ~1000 counts, retook the centering measurement, then re-isolated the ISI. This fixed the centering of the out of range T240, according the to voltages Hugh's script read out.
My before measurement of the corner 2 T240 was:
ETMX T240 2 DOF X/U = -1.988 [V]
ETMX T240 2 DOF Y/V = -2.227 [V]
ETMX T240 2 DOF Z/W = -1.135 [V]
After, from Hugh's script is:
ETMX T240 2 DOF X/U = 0.084 [V]
ETMX T240 2 DOF Y/V = 0.065 [V]
ETMX T240 2 DOF Z/W = 0.097 [V]
Much better. It's entirely likely this doesn't matter that much, we'll have to wait and see. I've been watching ETMX all day and it hasn't really rung up, so far.
Today we had a Norco LN2 delivery -> The increase in dewar LN2 column height following the delivery was enough to over fill CP3 by itself (recall that the LLCV is kept 15% open in "MANUAL" mode) -> I went ahead and opened the exhaust check valve bypass and the LLCV bypass to confirm the CP3 is over filled -> closed valves -> OK Next over fill of CP3 is scheduled for Thursday, Jan 21st before 4:00 pm
Today I perfomed several broad-band noise injections to check how large the coupling of MICH, PRCL and SRCL noises are to DARM, and to check if they are stationary.
In summary:
Some more details follow.
While injecting noise in an AUX channel, if the coherence with DARM is good, we can estimate the linear transfer functon TF. Then we can check if the PSD of DARM during the noise injection matches the PSD of AXU times the transfer function: the two should be equal if the coupling is dominated by a linear term. If there is strong non-linearity or strong non-stationarity, the two can be different. This is the case for SRCL (fig. 6).
The plan for the next days is to repeat this same kind of tests with ASC degrees of freedom.
As a side note, while injecting PRCL and SRCL I caused two lock losses. In both cases I'm quite confident that nothing was saturating. However, after the aborted SRCL injection, we relocked and found a mode at 41.02 Hz largely excited (SNR ~ 1000). This mode is unidentified, but seems to roughly match the roll mode of the triple suspensions.
For future reference, attached the logfile of the injections.
Maybe the ring-up hints at the origin of one of the two lines near 41 Hz that we suspect is a roll mode of a triple suspension. This line is discussed in alog 21696 and its comments. We should check that the mode just below 41 Hz is not also rung up. I wasn't immediately able to find which times to compare - could someone post good times to look at?
The line at 40.X Hz was not excited by the injection. Only the line at 41.02 Hz was.
The line at 41.02 Hz was very large during the first lock right after 137193090 Mon Jan 18 14:57:53 PST 2016
TITLE: Jan 19 EVE Shift 20:00-04:00UTC (12:00-20:00 PDT), all times posted in UTC
STATE Of H1: Commissioning
DAY OPERATOR: Patrick
QUICK SUMMARY:
Arrived in time for Journal Club Exodus from the control room but managed to get some preliminary commissioning plans from a few folks.
JEFF K: Reports that he has nothing planned for himself for this evening
JENNE: To take measurements at 2. 10 and 20 W LASER power.; Also, sensing matrix measurement(s)
EVAN:
SHEILA: Darm/actuator/sensor Upconversion studies
ROBERT SCHOFIELD: Reports that he has nothing further to do. PEM injections/measurements are done and he will be going over his data. Any trouble that he would find with his measurements would be the only cause for a re-do.
IFO unlocked for ≈2 hrs prior to my arrival
0915 -0945 hrs local (ran QDP80 and purge air compressors) Recently, Bubba G. had cut out and welded a section of the turbo header in order to create pallet jack access beneath the BT in the Y-mid VEA. -> Now iems can be stored on either side of the BT without the need to crane heavy items over the BT -> Today I pumped purge air through the header for a while and then verified that there weren't any leaks -> Achieved 9.2 x 10-3 torr after 20 minutes -> OK
J. Kissel I got charge measurements today, data files live in /ligo/svncommon/SusSVN/sus/trunk/QUAD/H1/ETMX/SAGL3/Data/ data_2016-01-19-18-26-30 data_2016-01-19-18-37-10 data_2016-01-19-18-48-53 data_2016-01-19-19-00-36 data_2016-01-19-19-12-29 /ligo/svncommon/SusSVN/sus/trunk/QUAD/H1/ETMY/SAGL3/Data/ data_2016-01-19-18-26-45 data_2016-01-19-18-38-30 data_2016-01-19-18-50-13 data_2016-01-19-19-02-04 data_2016-01-19-19-13-57 Data will be processed later.
Attached are 7 day pitch, yaw, and sum trends for all active H1 optical levers.
I reset the 35W FE watchdog at 18:17 UTC (10:17 PST).
I also noticed the crystal chiller water reservoir was very low; I added 200 mL of water.
Related: 25010
I wanted to see how well the LSC feedforward performs from 10 to 600 Hz.
I injected noise into MICH (and then SRCL) twice: once with the FF on and once with the FF off.
In both cases, the FF improves the coupling in the bucket, but makes it worse elsewhere. We already knew this was true around a few hertz (22586), but it seems to also be true above 500 Hz or so.
The intrinsic couplings (no FF) seem to have an uptick approaching 1000 Hz.
Looking both at BruCo reports and my yesterday noise injections, it looks like it's possible and necessary to improve the feedforward at low frequency, while the high frequency part is not very important.
Locked in NLN and observing bit set to preventive maintenance upon arrival. Tuesday maintenance has begun per Richard's entry.
IFO has been locked steadily last 12 hours.80 MPC. 0745 took IFO out of Observing for start of Maintenance. People were ready to go and LLO was down. 0745 S and K electric to EX air handler room to look at heater wiring. 0750 Chris S. will be resuming beam tube sealing on X2.