I've made some improvements, and done a bit of cleanup, to the SUS guardian code.
Some of this was to clarify the SUS/ISC interface, after discussion with the ISC team.
To the main SUS.py (H1 is using a separate copy of this in sus/h1/guardian/SUS.py, so we should eventually merge these changes into the common SUS.py):
I then cleaned up the various ALIGNED_TO_PD* states that have been added to align optics to various baffle PDs.
All code has been committed, and SUS guardians have been restarted to pull in these changes.
The ETM HWS computers currently have outdated .bashrc files that specify aliases and paths. I'm updating these.
With all SUSes in their aligned states, I have updated the Driftmon threshold values. This was mostly to address OM and RM values that have been ignored in the past. I trended their alignments in Dataviewer and found them to be in their nominal state this morning. The script to update these thresholds updates ALL SUSes, so the new values will reflect the latest in the constantly changing 'good' alignment of the IFO.
All SEI GS13s are now set in High Analog Gain without Whitening except HAM6 and BS (Stage2.) This standard High Gain configuration is with both FM4 & 5 Off in the GS13 INF module, H1:ISI-HAM5_GS13INF_H1 for example.
For HAM6, because of the shutter hit to the ISI, the GS13s are run with Low Analog Gain without Whitening. This configuration is with FM4 On and FM5 Off.
For BS, because of the hit to the ISI when the MICH acquires, we run the GS13s in the standard low Gain configuration of Analog Low Gain and Whitening. This has both FM4 & 5 On.
These configurations are set in SDF.
Here are the PSL DBB/ISS scans for this week. No significant change from last week.
One of the ISCT6 enclosure handle could not be locked and it was guarded with an yellow caution tape.
Corey brought a new handle, but when we assessed the situation it turns out that the new one and the old one were identical, it was just that the old one was installed at a wrong angle. We removed the old one and re-installed after rotating 90 degrees, and voila, it can be locked when the door is closed and the handle becomes horizontal. Only caveat is that the handle tip points up, not down, when the door is open.
We locked the doors and removed the yellow caution tape.
Here are the past 10 days trends:
The computer that generates web-accessible MEDM screen shots and model simlink captures is not running well since about March 11, 17:59 PDT. It will hopefully be repaired soon.
MEDM screen shots are working now, simlink capture will run at 12:00 PDT.
Today we have spent most of the day trying to lock. We have locked, but only after some diifictulties.
As Sheila mentioned we had problems engaging the ASC in DRMI. In particular, with today's initial alignment, for some reason, both POB_A and AS_C see pointing numbers between 0.2 and 0.46. When we engage the PRC1 loop (which uses POB_A), the control signal very quickly saturates the bottom stage of PRM, leading to a lock loss. The top mass relief loops are not fast enough this engagement. To a lesser extent the same is true for the SRC2 loop. For now I implemented a Guardian section that copies the error signal into the offset (with opposite sign), engages all loops at nominal bandwidth, and starts a 5min slow offset ramp back to the nominal point. We should do several things to mitigate this, probably in this order: - 1) Find an initial alignment procedure that puts us closer to where we want to be, making the WFS turn-on transient less violent. - 2) Commission a much higher bandwidth relief servo to the top mass. - 2) Implement a saturable integrator in the ASC control filter bank - like we have in the tidal servo. This would make engaging the WFS trivial, and guarantee to avoid saturation.
I also added limiters to the PRC1_P, PRC1_Y, SRC2_P and SRC2_Y filter modules: - PRC1: 1e4 cts limit. This corresponds to a drive of - 53168cts to PRM M3 coils - 8741cts to PR2 M3 coils - 24221cts to SRM M3 coils - 3457cts to SR2 M3 coils - 809cts to IM4 coils - SCR2: 1e3cts limit. This corresponds to a drive of - 40368cts to SRM M3 coils - 5762cts to SM2 M3 coils This should under all circumstances leave us enough range for length control on the M3 coils. In addition I added limiters to all involved to top stages angular controls: H1:SUS-BS_M1_LOCK_P_LIMIT = 500 H1:SUS-BS_M1_LOCK_Y_LIMIT = 500 H1:SUS-IM4_M1_LOCK_P_LIMIT = 1000 H1:SUS-IM4_M1_LOCK_Y_LIMIT = 1000 H1:SUS-PR2_M1_LOCK_P_LIMIT = 1000 H1:SUS-PR2_M1_LOCK_Y_LIMIT = 1000 H1:SUS-PR3_M1_LOCK_P_LIMIT = 500 H1:SUS-PR3_M1_LOCK_Y_LIMIT = 500 H1:SUS-PRM_M1_LOCK_P_LIMIT = 1000 H1:SUS-PRM_M1_LOCK_Y_LIMIT = 1000 H1:SUS-SR2_M1_LOCK_P_LIMIT = 1000 H1:SUS-SR2_M1_LOCK_Y_LIMIT = 1000 H1:SUS-SR3_M1_LOCK_P_LIMIT = 500 H1:SUS-SR3_M1_LOCK_Y_LIMIT = 500 H1:SUS-SRM_M1_LOCK_P_LIMIT = 1000 H1:SUS-SRM_M1_LOCK_Y_LIMIT = 1000 They are large enough that we shoul never run into them under normal operations, but small enough to avoid "astronomical" kicks on lock-loss.
Chris, Alexa
The first two attachments show the OLTFs of the various length DOFs. In these plots, the model is the blue trace and the red trace is measured data. The DARM TF agrees pretty well. The PRCL and SRCL loops appear to have a gain missing. The CARM and MICH model loops are not correct.
The fourth attachment depicts the noise budget produced by the model for the 8W lock on March 04, 2015 15:30 UTC. This roughly agrees with the conclusions Evan has made with his NB (LHO#17101). At high frequencies we are dominated by shot noise. At around 300 Hz we suspect intensity noise (or maybe beam jitter), as Gabriele noticed. From 30-100 Hz we are limited by DAC noise. And below 10 Hz we see some angular fluctuations.
The modeled MICH loop gain is looking better now. There were some aggressive BSFM plant inversion filters that weren't being linearized accurately in Simulink. This was fixed by using numerical TFs for them. The shape is still off above 60 Hz, but that's probably due to the filter change from alog 17148.
I put in sign flips as needed to ensure an even number of them are present in each OLTF. There remains some 10 dB of mystery gain to be hunted for in the DRMI loops.
SudarshanK, RickS
Today we had a quick look at our calibration data from Yend, taken after installing a polarizing beamsplitter downstream of the AOM (see aLOG 17145).
We will post a more detailed log after we have fully digested the data, but the preliminary result is that the new calibration of the Tx PD signal is within 6% of the previously reported number (see aLOG 16718).
This is about what we expected from our recent measurements of the signal drifts resulting from depolarization in the AOM.
Added 689 channels. Removed 123 channels.
The X direction at end X and the Y direction at the Y end are blended at 90mHz. Sensor correction is in the normal conifugration (no BRS). This was done at around 22:58 because we have had several lock losses due to large motions in the arms that ALS cannot handle.
This was at about -30 minutes on this screenshot, as you can see the arm control signals are reduced with the higher blends.
There is nothing that look like high ground motion on the seismic FOMs, although the winds are a little elevated (gusty 10-25 mph)
Here's a screen shot of the wind trend. Sheila switched the blends roughly 1 hour in the past on this trend.
Attached is a screenshot of the PEM's seismic FOM. So now we want to look at ISI T240s - 5 hours ago (~18:45 UTC -- after the Earthquake at -6 on the SEI FOM), - 3 hours ago (~20:45 UTC -- when the wind picks up, before the blends were switched) - and now (~22:45 UTC -- when the wind is still up, and when the blends have switched) Take a look at the time series of these few IFO channels: H1:ALS-(X/Y)_REFL_CTRL_OUT(MON or _DQ, whichever) this is what Sheila typically uses to assess badness. (and they're calibrated into [um]'s )
Alexa, Sheila, Evan, Chris, Stefan We implemented Peter's DARM 'SUScomp' from alog 16728. Since we can't lock with that filter directly (see alog 16840), deleted the old 'LLO' filter, but instead loaded a difference filter called 'acqLP' that makes the 'SUScomp' look like the old 'LLO' filter: zpk([80;500-800i;500+800i],[50;70;200],1,"n") Guardian was updated to turn off acqLP' in FM8 instead of turning on a lead. A note on the previous filter is also found in alog 16381.
I have attached a plot comparing our various configurations:
1. Red trace: RF DARM sus compenstation as designed by Peter to obtain more phase margin. The LSC DARM configuration is: FM1(suscomp), FM2(2:0), FM3(resG), FM4(4^2:1^2), FM5(2:0), Gain 800 (LHO#16728, 16840)
2. Blue trace: Our old RF DARM sus compenstation where we used the LLO control filer and a 200Hz lead filter. (LHO#16381)
3. Green trace: ALS DIFF sus compensation. This is the configuration we use to lock ALS DIFF. The LSC DARM configuration is: FM1(suscomp), FM2(2:0), FM3(resG), FM7(SB60), FM8(acqLP),FM10(RLP33), Gain 400. As Stefan mentioned, FM1+FM8 returns our old LLO control filter.
I have also attached the RF DARM OLTF model with the new (red) and old (blue) configuration as described above, along with the respective measured data. The RF DARM UGF is now 55 Hz with a phase margin of ~45deg.
Measurement of new DARM loop on dc readout is attached.
I've saved Evan's .xml to the calibration repository here: /ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/PreER7/H1/Measurements/DARMOLGTFs/2015-03-09_DARM_OLGTF_LHOaLOG17153.xml and exported text files of the transfer function and coherence, /ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/PreER7/H1/Measurements/DARMOLGTFs/ 2015-03-09_H1_DARM_OLGTF_LHOaLOG17153_coh.txt 2015-03-09_H1_DARM_OLGTF_LHOaLOG17153_tf.txt Transfer function contains the following columns (i.e. I exported IN1 / IN2 "as is"): Frequency [Hz] Real Part [ ] Imaginary Part [ ] We'll use later for calibration / noisebudget model verification!
Also including the CARM OLTF that we took at 9 W.
I checked the measured DARM open loop transfer function posted by Evan against my DARM open loop model. Even though I did not do a fitting or any fancy analysis yet, it seems that the optical gain was consistent -- the measurement matched the model with an optical gain of 1.1x106 cnts/m or 9.09x10-7 m/cnts which we have been using since Feb. 21st (alog 16843) for the CAL-CS front end model.
Here is a plot showing the model and measured one:
The matlab script to generate this plot is archived in calSVN:
/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/PreER7/H1/Scripts/DARM_OLTFGTF_LHOaLOG17153.m
Note for myself:
optical gain in the model = 1.1e6
ESD strength in the model = 2.8e-10 [N/V^2] (see alog 16843)
LSC_DARM_GAIN = 800 (instead of 400)
Only ETMX was actuated
Some specifics:
H1HWSEY
The .bashrc file is the same as that on H1HWSMSR with the following modifications:
.bashrcto.bashrc.old.20150312OPTIC_0toETMY_HWSOPTIC_1=ITMY_HWSenvironment variablePREFIX1=$IFO:TCS-$OPTIC1environment variablestream_image_EY' alias:'/opt/Hartmann_Sensor_SVN/release/bin/stream_image_Y/distrib/stream_image_Y'(this file is called 'stream_image_Y' but that's just a name, the important thing is that it points to the first port on the capture cardRun_HWS_ITMYandRun_HWS_ITMXaliasesRun_HWS_H1ETMY='/opt/HWS/Run_HWS_Y/distrib/Run_HWS_Y'The libraries for the MATLAB Runtime Environment are now properly identified in the path. I confirmed that
stream_image_EYruns and displays the image from the camera. See the second attached image.H1HWSEX
The .bashrc file is the same as that on H1HWSMSR with the following modifications:
.bashrcto.bashrc.old.20150312OPTIC_0toETMX_HWSOPTIC_1=ITMY_HWSenvironment variablePREFIX1=$IFO:TCS-$OPTIC1environment variablestream_image_EX' alias:'/opt/Hartmann_Sensor_SVN/release/bin/stream_image_Y/distrib/stream_image_Y'(this file is called 'stream_image_Y' but that's just a name, the important thing is that it points to the first port on the capture card)Run_HWS_ITMYandRun_HWS_ITMXaliasesRun_HWS_H1ETMX='/opt/HWS/Run_HWS_Y/distrib/Run_HWS_Y' (this is not a mistake - the Run_HWS_Y code points to the first port on the capture card)The libraries for the MATLAB Runtime Environment are now properly identified in the path. I confirmed that
stream_image_EXruns and displays the image from the camera. See the first attached image.Additionally on both:
/home/controls/archived_images/opt/HWS -> /opt/Hartmann_Sensor_SVN/release/HWS_v1.1.8//home/controls/framearchive/ -> /data/ETM/Note that when the HWS code stores data, it does so in the environment variable named FRAME_ARCHIVE_DIR which is defined by:
export FRAME_ARCHIVE_DIR=/home/controls/framearchive/$HOSTNAMESo, the files will be stored, ultimately in the following way:
/data/ITM/h1hwsmsr/ITMX//data/ITM/h1hwsmsr/ITMY//data/ETM/h1hwsex/ETMX//data/ETM/h1hwsey/ETMY/