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.
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...
http://earthquake.usgs.gov/earthquakes/eventpage/nc72387946
Ground velocity exceeds 20 μm/s from 0.03 Hz to 0.1 Hz. IMC MC2 transmission is noticeably wobbly.
Added 2 channels and removed 788.
Work Permit review.
SEI - continuing work on SDF system for BSC ISI. HEPI pump servo is ongoing tinvestigation. Maybe Jeff will run some TFs looking at 1Hz to 10Hz at some point today.
SUS - no big news. Thomas will ressurect the Drift Monitor this week sometime.
ISC/Commish - taking small steps in the positivedirection. Had a rather lengthy DRMI lock last night beginning ~01:30.
3IFO - no activities in LVEA. All activities across the street.
Facilities - Beam Tube cleaning is ongoing. Gutter work is don at LSB. Crew has moved to VPW.
CDS - variac installation at End station for Heater. Work in HAM6 area on Fast Shutter. R Abbott will be coming to work with this.
Gary will be moving ISS arrays and using the large airlock door in the LVEA
model restarts logged for Tue 27/Jan/2015
2015_01_27 01:43 h1fw0
2015_01_27 03:04 h1fw1
2015_01_27 08:57 h1iscex
2015_01_27 09:03 h1iscey
2015_01_27 09:18 h1isiham6
2015_01_27 16:32 h1fw1
2015_01_27 23:42 h1fw0
4 unexpected restarts. Light maintenance day, minor model changes for End ISCs and ISIHAM6 which did not require DAQ restart. Conlog frequently changing channels list attached.
I opened the gap between the bundle of copper pipes that were again rubbing on the ITMy optical lever TX pier causing the 10min-ish sawtooth signal. This is a temporary fix, but is more stable than the first time we saw this. It is in the planning to reroute the pipes and fix other similar issues with other oplvrs. Keita ran a new scan verifying the fix.
The OIB was set to commissioning upon my arrival for shift at 08:00. DRMI was still locked. I'm going to try my hand at locking the arm(s) as no one is here working yet.
I touched up the RM suspension controls today. RM1 had been tripping its watchdog often during the engagement of the scripts which use the RMs to center the beams on the REFL WFS.
Now, when the scripts turn these on the RMs do not trip and the beams are centered.
We have left the DRMI locked (on 1f), Guardian undisturbed starting around 1:13 AM for seismic performance evaluation.
Sheila, Evan, Rana, Eli, Kiwamu, Alexa
Today we were able to reduce the CARM offset from 800pm to 100pm. We made three major changes to acheive such an offset.
Transition Process from 800 pm to 400pm
To start,
With some TR buildup (approximately 1ct on LSC-X/Y_TR_A_LF_OUT),
With CARM on sqrt(TRX+TRY),
Transition Process from 400pm to 100pm
With all the settings as above, we bring the CARM offset to -3.2 (or ~150pm). At this point we transitioned the in-air TR PDs to the IR TR QPD_Bs. For reference, ASC-X/Y_TR_B_POW_NORM =1,1 and LSC-TR_X/Y_QPD_B_SUM_GAIN = 0.222, 0.175 respectively. Also, the QPD NSUM was about 400cts here. We continued to reduce the CARM offset to -4.5 cts (or ~80pm). At first,t the ASAIR_RF45_Q looked like a decent signal and responded well to a DIFF offset change. We turned on a servo for about 1.5min to ensure ASAIR was around zero (z servo -r LSC-ASAIR_A_RF45_Q_NORM_MON -g -11111 ALS-C_DIFF_PLL_CTRL_OFFSET). We measured ASAIR_RF45_Q/DARM_IN1 = 3e7 (at this point DARM_IN1 is just ALS DIFF signal). Note: we added a -160dB FM10 filter in LSC-ASAIR_A_RF45. Then the PD_DOF_MATRIX for ASAIR_RF45_Q is 3.3. So far, we were able to transition from ALS DIFF to RF DARM 50% of the way, and then we broke the lock. We started to reconsider whether we were actually in the linear range for RF DARM. Maybe we should be using AS45Q/SQRT(TRX) at this point? This requires a model change; however, this is already impleneted in the l1lsc.mdl that is in the SVN, so we could just copy that.
A little more ...
The farthest CARM offset we reached was -7 cts or ~60pm on sqrt(TRX+TRY) with the QPDs controlling CARM and ALS DIFF for DARM. Here we had 26 x single arm power, and REFL DC dropped from 80mW to ~75mW.
Some lock loss times for today:
Lock losses while transitioning DARM from ALS-DIFF to AS_AIR_45_Q:
06:01:40 UTC (28jan)
06:33:34
07:53:45
Almost there!
J. Kissel, B. Lantz, H. Radkins
Hugh and I, with the remote advice from Brian, have been investigating the coupling between HEPI differential pressure noise and HEPI platform motion (see preliminary studies in LHO aLOGs 16231 and 16309). In summary, the coupling is smattered between every chambers DOFs, limiting the performance between 0.02 and 2 [Hz]. There's no pattern as to how this coupling occurs except that the coherence for a given DOF happens in a similar frequency band in all chambers. All chambers show differential pressure moderate coherence to HP, but we're still working on estimating the transfer functions between this "pringle" mode and any translations / rotations to assess this impact. Further, we still have set up any successful long stretches of IFO / WFS data for the DRMI, so this pressure noise' impact can't yet be assessed. However, for now, I continue under the assumption that *any* low frequency improvement, especially to RZ / YAW, will improve the relative angular fluctuations of the IFO and stabilize optical gains, etc.
Hugh has already revealed the punchline: this morning's pump-servo OFF data shows that the problem is entirely because is of noise on differential pressure sensor. I'll attach a comment with complete set of plots for proof in a bit.
Further discussion on the noise as it stands with the pump servo (i.e. ON), based on the attached plots, which include
is below. I post all of these discussion points for future reference and design considerations to show how nastily and diversely imposing sensor noise on the pump servo imposes the platform motion. But mostly, because I sorted through a TON of data and plots to make these conclusions, before we knew the source of the problem.
** For the HEPI L4C ASDs, I plot the *raw* sensor noise because (a) I'm trying to cram as much information on as few plots as possible, and (b) the scale factors for both the IPS and L4Cs are the same for each DOF, and are all between 0.3 and 1.1, so the estimation estimation is off by at most a factor of three. Mostly it's there just to guide the eye, as a primitive noise budget.
(1) Krishna (before coherence between all DOFs and chambers was fully flushed out), suspected that because we now run Z sensor correction is on HEPI, the extra drive force meant excess pressure noise coupling. However, compared to 23e3 [ct_{rms}] range, the Vertical RMS of 5 [ct_{rms}] is peanuts. Further discussion of DC / mean requested drive is below.
(2) Brian suspected that the passive filtering system, the accumulating bladders might not be operating at their nominal air pressure. However, pump 8's (out of 12 for the corner) bladder pressure was checked last week, and was within spec. There're still lots lots to check, but Hugh suspects they're also OK.
(3) I suspect the PID controller for the servo is likely not well-tuned for new differential pressure sensing, suggesting we need to remeasure plant (which could mean significant down time for the IFO). There's been some discussion offline in the SEI group as to how we should systematically characterize the pump servo plant, given there's no "fast" excitation point in the pump servo make a tranditional frequency-response transfer function difficult (see SEI aLOGs 684 and 686). For now, we suspect we can do a sufficient characterization with simple steps using the existing EPICs infrastructure. Characterization to-date (though very early on in aLIGO) was also made by step response (see E1100508), but the assumptions about the plant don't seem sound in retrospect, and were done when the pressure was still on an absolute "single ended" readback.
- Brian's "Allowable Pump Servo Noise" defined in T1100198, but contains a LOT of assumptions. One of which is that the open loop gain transfer function of the feed back loops is "100 at low frequency." As such, Brian suggests, if the HEPI feedback loops are gain limited, try adding some boost at these frequencies to suppress the noise. However, I argue, with Hugh's recent redesign of the controller (see LHO aLOG 15308, the relevant attachment re-attached here), RZ Loops for example have *plenty* of gain, upwards of 3000 - 4000 at 0.1 [Hz], typically asymptoting to a few million a DC. This is far more than Brian assumed, so I thing we're alright in the gain department. Regrettably, the differential pump pressure is an EPICs channel, so I can't get an ASD of the performance above ~8 [Hz], but I'll post more on that in the future comment.
- Differential Pressure Noise supposedly coupling increases with actuator force request, but I've not found this to be true empirically. Check out the table below. Colors of the text help indicate which requested DC offsets are small, and which are large (dark green is smallest, then light green, gold, orange, with bright red is largest, in bins of 50 [micro-whatevers] and anything higher than 200 [micro-whatevers] is marked as bright red). Cells shaded in gray show coherence between the differential pressure and the HEPI L4Cs.
There is no discernable pattern:
X [um] | Y [um] | Z [um] | RX [urad] | RY [urad] | RZ [urad] | HP ["um"] | VP ["um"] | |
HAM2 |
24 |
-159 |
-142 |
27 |
160 |
-21 |
185 |
6 |
HAM3 |
3 |
109 |
-105 |
-24 |
-25 |
-124 |
131 |
93 |
HAM4 |
169 |
-124 |
-45 |
11 |
4 |
-78 |
141 |
-14 |
HAM5 |
47 |
-103 |
-17 |
97 |
-13 |
13 |
71 |
-19 |
HAM6 |
111 |
108 |
-283 |
1 |
115 |
1 |
320 |
-28 |
BS |
3 |
39 |
-105 |
-162 |
-177 |
-20 |
-140 |
-58 |
ITMX |
-118 |
13 |
127 |
32 |
-20 |
70 |
-1 |
11 |
ITMY |
600 |
-63 |
-189 |
15 |
60 |
-16 |
-87 |
27 |
Here're what patterns I was able to glean from the above table and attached plots:
- My best guess at this point is that it depends on the intricacy of how the piping is laid out between each chamber and the pumps. But it's only a guess, it's the only thing I know that significantly different between all chambers.
-------------------
DTT templates for this data live here:
/ligo/svncommon/SeiSVN/seismic/HEPI/H1/Common/
2015-01-23_H1HPI_ActuatorDrive.xml
2015-01-23_H1HPI_PumpControllerNoise.xml
J. Kissel, H. Radkins Here's the proof -- checking every DOF and every corner station HEPI -- that all coherence disappears when the pump servo is turned off. I also attach a comparison of the open and closed loop ASD of the differential pressure against the goals and requirements defined in T1100198. P.S. If you weren't excited about HEPI before now -- check this out. (Probably best to right-click and "save link as...")
In the attached trend of the OLPIT 3-10 Hz BLRMS, you can see the ringup start around 6AM local time this Sunday, the 25th.
Around noon on Monday (1800 UTC), Evan and Alexa start damping the mode using the bandpass filter in the OL PIT servo. They stop around 2100 UTC when the mode is small. It then rings up for the next 7 hours until we turn the stopband filter on.
Then it rings down with a 1/e time of ~4 hours (which implies a bounce mode Q of ~440,000 if nothing else is driving it).
We should make sure to install these on all of the test mass OL servos.
Added the Bounce + Roll bandstop filter to all the ETM/ITM OLDAMP filter banks. In some cases, I have overwritten existing 'notch' filters there for this purpose. Its not a good idea to use 'notch' for these mechanical modes which drift in frequency.
I made the Roll bandstop wider since there's no phase hit from this higher frequency filter. The vertical lines in the plot show the expected bounce frequency (9.775 Hz) and the roll frequency (13.81 Hz) which has some small natural spread between the optics, but I figure this should catch them all, even when the VEA temperatures drift.
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...