40m QIL Cryo_Lab CTN SUS_Lab TCS_Lab OMC_Lab CRIME_Lab FEA ENG_Labs OptContFac Mariner WBEEShop
 40m Log, Page 229 of 341 Not logged in
ID Date Author Type Category Subject
11474   Sat Aug 1 17:04:29 2015 EveUpdateSummary PagesStates and Triggers in SPs

I've added states to the summary pages to only show data for times at which one certain channel is above a specified threshold. So far, I've incorporated states for the IOO tab to show when the mode cleaner is locked.

You can see these changes implemented in the IOO tab of my personal summary pages for June 30: https://ldas-jobs.ligo.caltech.edu/~eve.chase/summary/day/20150630/ioo/.

I've written a description of how to add states to summary pages here: https://wiki-40m.ligo.caltech.edu/DailySummaryHelp#How_to_Define_and_Implement_States.

11475   Sat Aug 1 20:46:29 2015 KojiUpdatePEMX seismo station short cable removed

## OMG

11476   Mon Aug 3 08:16:19 2015 SteveUpdateSUSETMX damping restored
Attachment 1: ETMXrestored.png
11477   Mon Aug 3 18:19:09 2015 JessicaUpdateGeneralAnodization of front panels accounted for

Previously, I had gotten the same results for the conductive and the isolated front panels. Today, I sanded off the anodized part on the back of the conductive front panel. I checked afterwards with a mulitmeter to ensure that it was indeed conductive through all the SMA connectors.

I drove a frequency of 29.359 Hz through the X Arm cable and 29.3592 Hz through the Y Arm cable, giving a difference of 200 Hz. Previously, there would only be a spike in the Y Arm at the difference, while the X Arm did not change if the Y arm was on or off. Now that the panel is fully conductive, a spike can also be seen in the X arm, indicating that crosstalk may possibly be happening with this panel, now that the spike corresponds to both the X arm and Y arm. These results are only after one set of data. Tomorrow I'll take two more sets of data with this panel and do a more in depth comparison of these results to what had been previously seen.

Attachment 1: redo_conduct1X.png
Attachment 2: redo_conduct1Y.png
11478   Tue Aug 4 03:02:30 2015 ericqUpdateLSCBeat note Alignment fluctuation effects measured

Notes from tonight's work:

• PMC alignment tweaked. Not much gained
• WFS/MC2 offsets tweaked after recentering beams on WFS and some hand alignment.
• Vertex oplevs realigned for the first time in forever
• With an RF coupler, measured the X green beatnote to be +5dBm into the splitter. This resulted in -33dBm at the control room analyzer.
• Switched the ALS demod board inputs, from piping the delayed signal to the RF input, to sendingit to the LO input. This was motivated by wanting the mixer closer to compression, hopefully to reduce beatnote amplitude fluctuation sensitivity. This won some noise >100Hz.
• This led to record ALS noise levels - X:217Hz, Y:203Hz
• +2dBm into the board still leaves us some headroom for futher amplification. Board schematic lists +10dBm LO as "nominal," but maybe this isn't worth it...
• PRFPMI locking is still stalled at bringing in the RF signals. Debugging continues.
• Some beatnote amplitude fluctuation investigations (see below)
• Note to self: demod board schematics include an unspecified RF lowpass. Check out what got stuffed in there.

I've explored the beatnote fluctuations a bit further.

First, I realized that we already had a channel than functions much like an RF level monitor: the phase tracker Q output. I verified that indeed, the Q signal agrees with the RF monitor signals from the demod board within the phase tracker bandwidth. This simplifies things a little.

I also found that the Y beat suffers a fair bit less from these effects; which isn't too surprising given our experience with the alignment stability.

One possible caveat to my earlier conclusions is that the beatnote amplitude could be fluctuating due to real RIN of the green light transmitted through the cavity. In fact, this effect is indeed present, but can't explain all of the coherence. If it did, we would expect the DC green PDs (ALS-TR[X/Y]) to show the same coherence profile as the RF monitors, which they don't.

The next thing I was interested was whether the noise level predicted via coherence was realistic.

To this end, I implemented a least-squares subtraction of the RF level signal from the phase tracker output. I included a quadratic term of the RF power, but this turned out to be insiginficant.

Indeed, using the right gain, it is possible to subtract some noise, reproducing nearly the same spectrum as the coherence based estimate. The discrepency at 1Hz is possible from 1Hz cavity RIN, as suggested by the presence of some coherence with TRX.

However, this is actually kind of weird. In reality, I would've expected the coupling of RF level fluctuations to be more like a bilinear coupling; changing the gain of the mixer, rather than directly introducing a linearly added noise component. Maybe I just discovered the linear part, and the bilinear coupling is the left over low frequency noise... I need to think this over a little more.

Attachment 1: coherences.png
Attachment 2: linX.png
11479   Wed Aug 5 10:56:07 2015 ericqUpdateCDSMany models crashed

Last night around 1AM, many of the the frontend models crashed due to an ADC timeout. (But none of the IOPs, and all the c1lsc models were fine.)

First, on c1sus (Wed Aug  5 00:56:46 PDT 2015)
[1502036.695639] c1rfm: ADC TIMEOUT 0 46281 9 46153 [1502036.945259] c1pem: ADC TIMEOUT 0 56631 55 56695 [1502036.965969] c1mcs: ADC TIMEOUT 1 56706 2 56770 [1502036.965971] c1sus: ADC TIMEOUT 1 56706 2 56770

Then, simultaneously on c1ioo, c1iscex, and c1iscey. (Wed Aug  5 01:10:53 PDT 2015)

[1509007.391124] c1ioo: ADC TIMEOUT 0 46329 57 46201 [1509007.702792] c1als: ADC TIMEOUT 1 63128 24 63192

[2448096.252002] c1scx: ADC TIMEOUT 0 46293 21 46165 [2448096.258001] c1asx: ADC TIMEOUT 0 46669 13 46541

[1674945.583003] c1scy: ADC TIMEOUT 0 46297 25 46169 [1674945.685002] c1tst: ADC TIMEOUT 0 52993 1 52865

I'm still working on getting things back up and running. Just restarting models wasn't working, so I'm trying some soft reboots...

UPDATE: A soft reboot of all frontends seems to have worked,

Attachment 1: crashes.png
11480   Wed Aug 5 17:15:08 2015 EveUpdateSummary PagesFixed ASC Tab

I've fixed the ASC tab on the summary pages to populate the graphs with data without causing an error.

Motivation: The ASC tab was showing no data. It resulted in a name error when generated.

What I did:

A name error indicates a bad channel name in the plot definition. I identified two errors in the code:

1. I said C1:SUS_MC1_ASCPIT_OUT16.mean instead of C1:SUS-MC1_ASCPIT_OUT16.mean (underscore should be dash)
2. The channel C1:ASX-XARM_M1_PUT_OSC_CLKGAIN was resulting in a name error. I removed it.

Results:

The plots are not processing without error. However, no titles or axis labels are present on the plots- I'll work on adding these.

11481   Thu Aug 6 01:38:19 2015 ericqUpdateComputer Scripts / ProgramsChiara gets new Ethernet card

Since Chiara's onboard ethernet card has a reputation to be flaky in Linux, Koji suggested we could just buy a new ethernet card and throw it in there, since they're cheap.

I've installed a Intel EXPI9301CT ethernet card in Chiara, which detected it without problems. I changed over the network settings in /etc/networking/interfaces to use eth1 instead of eth0, restarted nfs and bind9, and everything looked fine.

Sadly, EPICS/network slowdowns are still happening. :(

11482   Thu Aug 6 04:36:41 2015 ericqUpdateASCReviving PRC angular feedforward

Tonight, I've taken a bunch of data where the PRC is carrier locked and the ITM oplevs have the DC coupling FM turned on, as we use during locking. This is to inform new feedforward filters to stabilize the PRC angular motion, by using Wiener filtering with the POP QPD as the target, and local seismometers/accelerometers as witnesses. So far I've looked at the 1800 seconds leading up to GPS time 1122885600, but there has been plenty of locked time tonight if I need to retrieve more.

I've also measured the PRM ASC output torque -> POP QPD spot motion with high (>0.95) coherence from 0.1Hz to 10Hz.

Prefiltering so far consists of a 4th order elliptic LP at 5 Hz, with the target subtraction band being the 1-3Hz range.

With offline FIR filtering, the RMS pitch motion is reduced by a factor of 3 just with the STS1_X data. IIR fitting remains to be done.

The PRC yaw motion, which is marginally noisier, is a little more tangled up across X and Y.

Plots / filters forthcoming pending more analysis.

11483   Thu Aug 6 09:16:57 2015 SteveUpdateVACChiara gets new Ethernet card
1. Yesterday we did put the vacuum system into  for the ethernet card swap.
2. Atm1, IFO pressure
3. Atm2, as prepared valve configuration where " moving " means closed and disconnected

Attachment 1: ChiaraNetCard.png
Attachment 2: 23.png
11484   Thu Aug 6 11:45:01 2015 JessicaUpdateGeneralALS Delay Line front panel testing

Koji had suggested that I sync up the two function generators to ensure that they have the same base frequency and so that crosstalk will actually appear at the expected frequency. After syncing up the two function generators, I drove the following frequencies through each cable:

Conductive SMAs:

X: 29.537 MHz           Y: 29.5372 MHz

Isolated SMAs:

X: 29.545 MHz           Y: 29.5452 MHz

Each time, the difference between the frequencies was 200 Hz, so if there was crosstalk, a spike should appear in the PSDs at 200 Hz when frequencies are being driven through both cables simulataneously, but not when just one is on. We very clearly see a spike at 200 Hz in both the X arm and the Y arm with the conductive SMAs, indicating crosstalk. For the front panel with isolated SMAs, we see a spike at 200 Hz when both frequencies are on, but it is much less pronounced than with the conductive SMAs. It seems as though there will be crosstalk using either panel, just less with the isolated SMAs.

Attachment 1: conductive_X.png
Attachment 2: conductive_Y.png
Attachment 3: isolated_X.png
Attachment 4: isolated_Y.png
11485   Thu Aug 6 21:03:45 2015 IgnacioHowToWienerFilteringHow to do online static IIR Wiener filtering

In order to do online static IIR Wiener filtering one needs to do the following,

1) Get data for FIR Wiener filter witnesses and target.

2) Measure the transfer function (needs to be highly coherent, ~ 0.95 for all points) from the actuactor to the control signal of interest (ie. from MC2_OUT to MC_L).

3) Invert the actuator transfer function.

4) Use Vectfit or (LISO) to find a ZPK model for the actuator transfer and inverse transfer functions.

5) Prefilter your witness data with the actuator transfer function, to take into account the actuator to control transfer function when designing the offline Wiener FIR filter.

6) Calculate the frequency response for each witness from the FIR coefficients.

7) Vectfit the frequency reponses to a ZPK model, this is the FIR to IIR Wiener conversion step.

8) Now, either, divide the IIR transfer function by the actuator transfer function or more preferably, multiply by the inverse transfer function.

9) Use Quack to make SOS model of the IIR Wiener filter and inverse transfer function product that goes into foton for online implementation.

10) Load it into the system.

The block diagram below summarizes the steps above.

Attachment 1: iir.png
11486   Mon Aug 10 11:57:45 2015 ericqUpdateGeneralIMC work
 Quote: Often when I come to manually lock the mode cleaner due to a long unlocked period, I find that the sliders are not in the state specified by the mcdown script. Furthermore, it's not the same channels every time; sometimes the servo gain is left high, sometimes the boosts are left on. I fear that some of the caput commands are failing to execute. Ugh.

This continues to happen. I believe the network latency boogeyman is to blame.

There was a long unlocked period because the enable switch for the MC servo fast path (FASTSW) was left off. Running the mcdown script fixed this, but included the error message:

Channel connect timed out: 'C1:IOO-MC_REFL_GAIN' not found. CA Client Library: Ignored duplicate create channel response from CA server?

which means the IN1 gain didn't get touched. A second pass of the script produced no errors.

I'm thinking of adding some logic that if the autolocker has failed to lock for some period (5 minutes?), it should rerun mcdown.

11487   Mon Aug 10 15:25:05 2015 SteveUpdatePEMwasp nest

The wasp nest will be removed tomorrow from from the out side of the east arm window.

The resonant frequency of the newly arrived gravity bee detector is not known.

Attachment 1: waspN.jpg
11488   Mon Aug 10 22:18:19 2015 IgnacioUpdateIOOReady to do some online mode cleaner subtraction

I'm attaching a SISO IIR Wiener filter here for reference purposes that will go online either tonight or tomorrow evening. This is a first test to convince myself that I can get this to work, MISO IIR filters are close to being ready and will soon be employed.

This Wiener filter uses the STS-X channel as a witness and MCL as target. The bode plot for the filter is shown below,

The performance of the FIR and IIR Wiener filters and the ammount of subtraction achive for MCL is shown below,

Output from quack to be loaded with foton: filter.zip

K bye.

Attachment 1: stsx.png
Attachment 2: performance.png
Attachment 3: filter.zip
11489   Tue Aug 11 02:26:46 2015 ericqUpdateASCPRC Angular FF Lives!

# PRC Angular FF is back in action!

Short and sweet of it:

• Took witness (T240 channels) and target (POP QPD) with DC coupled oplevs on. About 25 minutes of nice stationary data.
• Downsampled everything to 32Hz, since coherence suggests subtraction only really possible from 1-5Hz.
• Prefiltering done by detrending and ellip(3,3,40,5Hz)
• 4 second FIR impulse time was enough
• Filtered target by inverse actuator TF before sending to wiener code. The only difference between this and filtering the witnesses with the actuator TF directly is an effective RMS cost function, i.e. prefiltering.
• Spending time tweaking IIR fitting pays off. Divided out zpk(0, [p3, p3*],1), where p3 is some well fit stack/suspension resonance, so that vectfit fits remaining portion with equal numbers of poles and zeros, guaranteeing AC coupling and 1/f rolloff to prevent noise injection
• Quack->foton->OAF all worked fine
• ## All in all, seems to work well. POPDC RMS goes down by a factor of 2

• Code used lives in /users/ericq/2015-08-PRCFF and the NoiseCancellation github repo

## Subtraction prediction vs. reality (positive dB is good)

Attachment 1: fitExample.png
Attachment 2: FFspectra.png
Attachment 3: PITsub.png
Attachment 4: YAWsub.png
11490   Tue Aug 11 02:40:29 2015 ericqUpdateLSC50m delay lines - Rough calibrations

Jessica will soon ELOG about some measurements suggesting that the conductive connector-ized ALS delay line enclosure is the way to go, when considering crosstalk between the delay lines. It is currently mounted and hooked up on the LSC rack, though I need to make a bunch of new SMA cables now that I think a semi-permanent arrangement has been reached.

I did a rough re-calibration of the phase tracker output, since the increased cable delay changes the degree/Hz gain. This was done by fitting a line to a slow sawtooth FM of the SRS DS345's (1Hz rate, 10kHz deviation, 30MHz carrier). This resulted in the following calibration updates

• ## ALSY: 19425 -> 12585 Hz/count, 3.8db more sensitive

Again, this is a rough calibration. Nevertheless, it is not so surprising we don't get the 50m/30m = 4.4dB increase we would expect just from the lengths; the (I presume) increased cable loss matters. Also, the loss' frequency dependance is an additional reason that the phase tracker calibration is not constant over all frequencies.

I took spectra with the arms in IR lock, but didn't see any real improvement beyond a possible dip in the floor from 100-200Hz. This doesn't surprise me too much, however, since I don't believe that we are currently dominated by electronic noises that this gain increase would help overcome.

Last week, Koji mentioned the ALS phase noise added due to the post-cavity table motion the arm-transmitted green beams experience before hitting the beat PD. I should estimate the size of this effect for our situation.

11491   Tue Aug 11 10:13:32 2015 JessicaUpdateGeneralConductive SMAs seem to work best

After testing both the Conductive and Isolated front panels on the ALS delay line box using the actual beatbox and comparing this to the previous setup, I found that the conductive SMAs improved crosstalk the most. Also, as the old cables were 30m and the new ones are 50m, Eric gave me a conversion factor to apply to the new cables to normalize the comparison.

I used an amplitude of 1.41 Vpp and drove the following frequencies through each cable:
X: 30.019 MHz          Y: 30.019203 MHz

which gave a difference of 203 Hz.

In the first figure, it can be seen that, for the old setup with the 30m cables, in both cables there is a spike at 203 Hz with an amplitude of above 4 m/s^2/sqrt(Hz). When the 50m cables were measured in the box with the conductive front panel, the amplitude drops at 203 Hz by a factor of around 3. I also compared the isolated front panel with the old setup, and found that the isolated front panel worse by a factor of just over 2 than the old setup. Therefore, I think that using the conductive front panel for the ALS Delay Line box will reduce noise and crosstalk between the cables the most.

Attachment 1: best4.png
Attachment 2: isolated4.png
11492   Tue Aug 11 11:30:19 2015 IgnacioUpdateIOOSISO (T240-X) FF of MCL

Last night we finally got some online subtraction going. The filter used is described in the post this eLOG is @eLOG 11488

The results were as follow:

The filter worked as expected when subtracting noise out of MCL,

There is about a factor of 6 subtraction at the ~3Hz resonant peak. The static IIR filter predicted a factor of 6-7 subtraction of this peak as well.

The 1.2 Hz resenonant feature improved by a factor of 3. This should improve quite drastically when I implement the y-channel of the T240 seismo.

There is some high frequency noise being injected, not very noticeable, but present.

We then took a look at the power in the MC when the filter was on,

The power being transmitted in the cavity was not as stable as with the feedforward on. We believe that the filter is not at fault for this as Eric mentioned to me that the MC2 actuator lacked some sort of compensation that I need to understand a bit better.

YARM was then locked when the filter was on and we took a look at how it was doing. There was stationary sound arising from the locking of the YARM, leading us to believe that the filter might have injected some noise in the signal. IT DID.

The filter injected nasty high frequency noise at YARM from 11 Hz and on. This is to be expected since the filter did not roll off to zero at high frequencies. Implementing a 1/f rolloff should mitigate some of the injected noise.

Also, as one can see above, subtraction by around a factor of 2 or so, was induced by the mode cleaner feedforward subtraction.

Attachment 1: MCL.png
Attachment 2: MCTRANS.png
Attachment 3: YARM.png
11493   Tue Aug 11 11:56:36 2015 Ignacio, JessicaUpdatePEMWasps obliterated maybe...

The wasp terminator came in today. He obliterated the known wasp nest.

We discovered a second wasp nest, right next to the previous one...

Jessica wasn't too happy the wasps weren't gone!

11494   Tue Aug 11 16:13:28 2015 EveUpdateGeneralGaussianity tests

I’m working on a code to determine the Gaussianity of a PSD.

Motivation:

It can be difficult to distinguish between GW events and non-Gaussian noise, especially in burst searches. By characterizing noise Gaussianity, we can better recognize noise patterns and distinguish between GW events and noise.

What I did:

I analyzed an hour of S5 L1 data. First, I plotted a timeseries, just to see what I was working with. Then, I produced a PSD (technically, an ASD) for the timeseries using Welch’s method in Python.

I split the data segment into smaller time-chunks and then produced a PSD for each chunk. All PSDs were superimposed in one plot. Here’s a plot for 201 time-chunks of equal length:

For a specific frequency, I can view the spread in PSD value through the production of a histogram.

Results:

I’ve made histograms displaying varying PSD values for the 201 PSD plot at 100 Hz, 500 Hz, and 1kHz.

For Gaussian noise, an exponential decay plot is expected. I will continue this analysis by following the statistical method in Ando et al. 2003 to calculate specific values indicative of the Gaussianity of various distributions. I’ll then look at different periods of time in the S5 L1 data to find periods of time suggesting non-Gaussian behavior.

11495   Tue Aug 11 18:43:42 2015 JessicaUpdateIOOMCL Online Subtraction

Today I finished fitting the transfer function to a vectfit model for seismometers T240_X and T240_Y, and then used these to filter noise online from the mode cleaner.

The Bode plot for T240_X is in figure 1, and T240_Y is in figure 2. I made sure to weight the edges of the fit so that no DC coupling or excessive injection of high frequency noise occurs at the edges of the fit.

I used C1:IOO-MC_L_DQ as the first channel I filtered, with C1:IOO-MC_L_DQ(RMS) for RMS data. I took reference data first, without my filter on. I then turned the filter on and took data from the same channel again. The filtered data, plotted in red, subtracted from the reference and did not inject noise anywhere in the mode cleaner.

I also looked at C1:LSC-YARM_OUT_DQ and C1:LSC-YARM_OUT_DQ(RMS) for its RMS to see if noise was being injected into the Y-Arm when my filter was implemented. I took reference data here also, shown in blue, and compared it to data taken with the filter on. My filter, in pink, subtracted from the Y-Arm and injected no noise in the region up to 10 Hz, and only minimal noise at frequencies ~80 Hz. Frequencies this high are noisy and difficult to filter anyways, so the noise injection was minimal in the Y-Arm.

Attachment 1: SeisX_bode.png
Attachment 2: SeisY_bode.png
Attachment 3: MCL_first.png
Attachment 4: Yarm_first.png
11496   Wed Aug 12 01:32:18 2015 IgnacioUpdateIOOImproved SISO (T240-X) FF of MCL

In my previous elog:11492, I stated that in order to improve the subtraction and reduce the injection of high frequency noise we want the filter's magnitude to have a 1/f rolloff.

I implemented this scheme on the filter SISO filter previously analyzed. The results are shown below.

The filters bode plot:

The nice 1/f rollof is the main change here. Everything else remained pretty much the same.

The predicted FIR and IIR subtractions:

Everything looks right but that hump at 8 Hz. I used 8 pairs of poles/zeros to get this subtraction.

The online MCL subtraction:

This looks better than I expected. One has to keep in mind that I ran this at 1 AM. I wonder how well this filter will do during the noisier hours of the day. The RMS at high frequencies doesn't look great, there will definitely be noise being injected into the YARM signal at high frequencies.

Measuring the YARM signal:

There is still noise being injected on YARM but it is definitely much better than the previous filter. I'm thinking about doing some IIR subtraction on the arms now to see if I can get rid of the noise that is being injected that way, but before I embark on that quest I will rething my prefiltering.

The plot below shows the ratio of the unfiltered versus filtered ASDs for the FIR and IIR subtraction predictions as well as for the measured online IIR subtraction. Positive dB means better subtraction.

Attachment 1: filter.png
Attachment 2: stsx.png
Attachment 3: mclonline.png
Attachment 4: yarmonline.png
Attachment 5: sub.png
11497   Wed Aug 12 11:44:13 2015 ericqUpdatePEMGur2 Channels still wonky

In previous elogs, we saw that the X and Y spectra out of GUR2 (X end Guralp seismometer) looked strange (i.e. inconsistent with the GUR1 spectra).

This morning, Steve and I brought the handheld control unit to the Guralp to center the test mass, by adjusting the centering potentiometers inside the unit while monitoring the voltage readout corresponding to the DC mass position (manual has instructions).

At first glance, this seemed like the likely culprit, as the offsets for the horizontal directions were much larger than the vertical one. We zeroed all three to within a mV or two. Unfortuntately, the spectra look the same as they did 10 hours ago.

Since we already had the kit out, we checked the offsets for GUR1. Only the "East/West" had an offset over 50mV, so we zeroed that one, but left the others alone.

Attachment 1: gur2_centered.png
11498   Wed Aug 12 14:35:46 2015 ericqUpdateComputer Scripts / ProgramsPDFs in ELOG

I've tweaked the ELOG code to allow uploading of PDFs by drag-and-drop into the main editor window. Once again we can bask in the glory of

(You may have to clear your browser's cache to load the new javascript)

Attachment 1: smooth.pdf
11499   Wed Aug 12 16:39:46 2015 IgnacioUpdateIOOMISO WIener (T240-X and T240-Y) FF of MCL

Last night I performed some MISO FF on MCL using the T240-X and T240-Y as witnesses. Here are the results:

Filter:

T240-X

T240-Y

Training data + Predicted FIR and IIR subtraction:

Online subtraction results:

MCL
YARM

Subtraction performace:

Attachment 1: stsx.png
Attachment 2: stsy.png
Attachment 3: performance.png
Attachment 4: sub.png
Attachment 5: mcliir.png
Attachment 6: yarmiir.png
11500   Wed Aug 12 16:48:26 2015 IgnacioUpdateIOOBetter? Nope. MISO WIener (T240-X and T240-Y) FF of MCL

Last night, I also worked on a "better" (an improvement, I thought) of the MISO Wiener filter (T240-X and T240-Y witnesses) IIR filter. The FF results are shown below:

Filter:

T240-X

T240-Y

Training data + Predicted FIR and IIR subtraction:

Online subtraction results:

MCL
YARM

Subtraction performace:

Although the predicted FIR and IIR results are "better" than the previous MISO filter, the subtraction performance for this filter is marginally better if not worse (both peak at 15 dB, in slightly different regions).

Attachment 1: stsx.png
Attachment 2: stsy.png
Attachment 3: performance.png
Attachment 4: mcliir.png
Attachment 5: yarmiir.png
Attachment 6: sub.png
11501   Wed Aug 12 22:33:36 2015 IgnacioUpdateIOORe-measured MC2 -> MCL TF

Since I will need to do transfer function measurements in order to implement FF for the arms and the MC2's yaw and pitch channels, I decided to practice this by replicating the transfer function measurement Eric did for MC2 to MCL. I followed his procedure and the data that I aquired for the TF looked as shown below,

About five minutes of data were taken (0.05 Hz resolution, 25 averages) by injecting noise from 1 to 100 Hz. The TF coherence looked as below,

Attachment 1: bode_TF.png
Attachment 2: Coherence.png
11502   Thu Aug 13 12:06:39 2015 Jessica SummaryIOOBetter predicted subtraction did not work as well Online

Yesterday I adjusted the preweighting of my IIR fit to the transfer function of MC2, and also managed to reduce the number of poles and zeros from 8 to 6, giving a smoother rolloff. The bode plots are pictured here:

The predicted IIR subtraction was very close to the predicted FIR subtraction, so I thought these coefficients would lead to a better online filter.

However, the actual subtraction of the MCL was not as good and noise was injected into the Y arm.

The final comparison of the subtraction factors between the online and offline data showed that the preweighting, while it improved the offline subtraction, needs more work to improve the online subtraction also.

Attachment 1: newBodeX.png
Attachment 2: newBodeY.png
Attachment 3: pred_Sub.png
Attachment 4: MCLSub.png
Attachment 5: YarmSub.png
Attachment 6: comparison.png
11503   Thu Aug 13 20:32:07 2015 IgnacioUpdateLSCWorking towards YARM FF

The mode cleaner FF static filtering is by no means done. More work has to be done in order to succefuly implement it, by the means of fine tuning the IIR fit and finding better MISO Wiener filters.

I have begun to look at implementing FF to the YARM cavity for several reasons.

1) Even if the mode cleaner FF is set up as best as we can, there will still be seismic noise coupling into the arm cavities.

2) YARM is in the way of the beam path. When locking the IFO, one locks YARM first, then XARM. This means that it makes sense to look at YARM FF first rather than XARM.

3) XARM FF can't be done now since GUR2 is sketchy.

I'm planning on using this eLOG entry to document my Journey and Adventures (Chapter 2: YARM) to the OPTIMAL land of zero-seismic-noise (ZSN) at the 40m telescope.

11504   Thu Aug 13 23:57:33 2015 IgnacioUpdateLSCYARM coherence plots

I took data from 1123495750 to 1123498750 GPS time (Aug 13 at 3AM, 50 mins of data) for  C1:LSC-YARM_OUT_DQ, and all T240 and GUR1 channels.

Here is the PSD of the YARM_OUT, showing the data that I will use to train the FIR filter:

Coherence plots for YARM and all channels of T240 and GUR1 sesimometers are shown below. This will help determine what regions to preweight the best before computing FIR filter. They also show how GUR1 is back to work compared to those of elog:11457.

Attachment 1: YARM_psd.png
Attachment 2: YARM_GUR1_COH.png
Attachment 3: YARM_STS_COH.png
Attachment 4: YARM_GUR1_COH.png
11505   Fri Aug 14 09:07:45 2015 SteveUpdatePEMGur interface box is wonky

Atm1, Before cable swap

Atm2, The long cables were swapped at the input of the interface box.

We can conclude that the problem is in the interface box

I wonder if interface box input 3 is wired?

Attachment 1: Gur1&2.png
Attachment 2: longCablesSwapped.png
11506   Fri Aug 14 12:10:08 2015 KojiUpdatePEMGur interface box is wonky

Let's dismantle the I/F unit from the rack and connect the cable with the lid open.
We need to trace the signal.

11507   Fri Aug 14 17:20:01 2015 JenneUpdatePEMGur interface box is wonky

IIRC, the Guralp box's 3rd set of channels do not have all of the modifications that were made on channels 1 and 2.

11508   Fri Aug 14 21:40:26 2015 IgnacioUpdateLSCQuick static offline subtractions of YARM

Plotte below are the resultant subtractions for YARM using different witness configurations,

The best subtraction happens with all the channels of both the GUR1 and T240 seismometers, but one gets just as good subtraction without using the z channels as witnesses.

Also, why is the T240 seismometer better at subtracting noise for YARM compared to what GUR1 alone can acomplish? Using only the X and Y channels for the T240 gave the third best subtraction(purple trace).

My plan for now is as follows:

1) Measure the transfer function from the ETMY actuator to the YARM control signal

2) Collect data for YARM when FF for MCL is on in order to see what kind of subtractions can be done.

Attachment 1: arms_wiener.png
11509   Fri Aug 14 23:49:34 2015 KojiSummaryGeneralB&K Shaker fixed

I fixed a shaker that was claimed to be broken. I had to cut the rubber membrane to open the head.

Once it was opened, the cause of the trouble was obvious. The soldering joint could not put up with the motion of the head.

It is interesting to see that the spring has the damping layer between the metal sheets.

After the repair the DC resistance was measured. It was 1.9Ohm. The side of the shaker chassis said "3.5Ohm, Max 15VA". So it can take more than 4A (wow).

I gave 2A DC from the bench top supply and turn the current on and off. I could confirm the head was moving.

I'll claim the use of this shaker for the seismometer development.

Attachment 1: IMG_1947.JPG
11510   Sat Aug 15 02:10:35 2015 IgnacioUpdateLSCMCL FF => YARM FF

In my last post I calculated the different subtractions (offline) that could be done to YARM alone just to get a sense of what seismometers were better witnesses for the Wiener filter calculation.

In this eLOG I show what subtractions can be done when the MCL has FF on (as well as Eric's PRC FF), with the SISO filter described on elog:11496.

The plot below shows what can be done offline,

What is great about this results is that the T240-X and T240-Y channels are plenty enough to mitigate any remaining YARM seismic noise but also to get rid of that nasty peak at 55 Hz induced by the MCL FF filter.

The caviat, I haven't measured the TF for the ETMY actuator to YARM control signal. I need to do this and recompute the FIR filters with the prefiltered witnesses in order to move on to the IIR converions and online FF!

Attachment 1: YARM_LIVES.png
11511   Sun Aug 16 23:26:40 2015 EveUpdateGeneralGaussianity tests

I've continued to work on my Gaussianity tests for S5 L1 data.

Following the statistical measure in Ando et al. (2003), I've calculated the Laguerre coefficient, c2, for all frequencies present in my S5 L1 PSD as a metric of Gaussianity. When c2 is zero, the distribution is Gaussian. A positive c2 corresponds to glitchy noise, while a negative c2 suggests stationary noise.

Below is a plot displaying variation in c2 for this PSD:

By observing the c2 value and histogram of distribution of various PSD values at a given frequency, we can elucidate statistical differences in the Gaussian nature of noise at that frequency which are unclear in the standard PSD.

Attachment 1: Gaussianity_noc1.png
11512   Mon Aug 17 17:48:12 2015 KojiUpdatePEMWasps obliterated maybe...

We found the same wasp in the 40m. Megan found it walking behind Steve desk!

11513   Tue Aug 18 03:56:09 2015 ericqUpdateLSClocking efforts

Now that the updated ALS is stable, and the PRC angular FF is revived, I've been working on relocking PRFPMI. While the RMS arm fluctuations are surely smaller than they used to be, there is no noticible difference to the ears when buzzing around resonance, but this doesn't really mean much.

Frustratingly, I am not able to stably blend in any RF CARM error signal into the slow length control path (i.e. CARM_B). Bringing AS55 Q into DARM with the 20:0 integrator is working fine, but we really need to supress CARM to get anywhere. I'm not sure why this isn't working; poking around into the settings that were used when we were regularly locking didn't turn up any differences as far as I could tell. Investigations continue...

Some minor changes to the locking script were made, to account for the increased ALS displacement sensitivity from the longer delay line.

Since the ALS is now in a fairly stable state, I've updated the calibrated PSD template at /users/Templates/ALS/ALS_outOfLoop_Ref.xml, and added some coherence plots for some commonly coupled quantities (beat signal amplitude, IR error signal, green PDH error signal and green transmission).

Attachment 1: newALSref.pdf
Attachment 2: xCoh.pdf
Attachment 3: yCoh.pdf
11514   Tue Aug 18 11:16:17 2015 SteveUpdatePEMGur interface box is wonky

The Guralp ADC  interface box  D060506 is ready for inspection. It is in front of 1X1 with open top and running.

Obviously c7 as 1 miroF cap is missing.

 Quote: Let's dismantle the I/F unit from the rack and connect the cable with the lid open. We need to trace the signal.

Attachment 1: IMG_0009.JPG
11515   Wed Aug 19 00:55:35 2015 IgnacioUpdateLSCLSC-YARM-EXC to LSC-YARM-IN1 TF measurement + error analysis

Yesterday, Rana, Jessica and I measured the Transfer function from LSC-YARM-EXC to LSC-YARM-IN1.

The plot below shows the magnitude and the phase of the measured transfer function. It also shows the normalized standard error in the estimated transfer function magnitude; the same quantity can be applied to the phase, only in this case it is interpreted as its standard deviation (not normalized). It is given by

$\frac{[1-\gamma_{xy}^2(f)]^{1/2}}{|\gamma_{xy}(f)|\sqrt{2n_{d}}}$

where $\gamma_{xy}^2(f)$ is the ordinary coherence function and $n_{d}$ is the number of averages used at each point of the estimate, in the case here we used 9 averages. This quantity is of interest to us in order to understand how the accuracy of transfer function measurement affects the ammount of subtraction that can be achieved online.

Since this transfer function is flat from 1-10 Hz (out of phase by 180 deg), this means that we can apply our IIR wiener filters direclty into YARM without taking into account the TF by prefiltering our witnesses with it. Of course this is not the case if we care about subtractions at frequencies higher than 10 Hz, but since we are dealing with seismic noise this is not a concern.

The coherence for this transfer function measurement is shown below,

11516   Wed Aug 19 01:45:10 2015 IgnacioUpdateIOODoubly Improved SISO (T240-X) FF of MCL

Today I tried and doubly-improved SISO FF filter on MCL. This filter has a stronger rolloff than the previous SISO filters I have tried. The rolloff most definelty helped towards reducing the ammount of noise being injected into YARM. Below is the usual stuff:

Filter:

T240-X (SISO)

Training data + Predicted FIR and IIR subtraction:

Online subtraction results:

MCL
YARM
MCL TRANS

Subtraction performace:

11517   Wed Aug 19 07:58:25 2015 SteveUpdatePEMGur interface box

Koji and Steve,

We took transferfunctions of each channel yesterday. They were identical ?. I will check the cables from ADC to DAQ next.

11518   Thu Aug 20 02:31:09 2015 ericqUpdateLSCPRFPMI is back

# PRFPMI locking has been revived.

I've had 6 5min+ locks so far; arm powers usually hit ~125 for a recycling gain of about 7; visibility is about 75%

The locking script takes a little under 4 minutes to take you from POX/POY lock to PRFPMI if you don't have to stop and adjust anything.

At Koji's suggestion, I used digital REFL11 instead of CM_SLOW, which got me to a semistable lock with some RF, at which time I could check the CM_SLOW situtation. It seemed like the whitening Binary IO switch got out of sync with the digital FM status somehow...

I've been making the neccesary changes to the carm_cm_up script. I also added a small script which uses the magnitude of the I and Q signals to set the phase tracker gain automatically based on some algebra Koji posted in an ELOG some years ago.

The RF transition seems much smoother now, most likely due to the improved PRC and ALS stability. In fact, it is possible to hold at arm powers of >100 solely on the digital servos; I don't think we were able to do this before until the AO had kicked in.

Right now I'm losing lock when trying to engage the CARM super boost. I also haven't switched the PRMI over to 1F signals yet. Would be good to hook the SR785 back up for a loop TF, but I'll stop here for tonight since our SURFs are presenting bright and early tomorrow morning.

Attachment 1: lock.pdf
11519   Thu Aug 20 11:09:10 2015 ranaUpdateIOOsome points about seismic FF
• When plotting the subtraction performance, we mainly care about the 0.5 - 10 Hz band, so we care about the RMS in this band. Don't integrate over the whole band.
• When calculating the Wiener filter, you must use the pre-weighting so as to not let the Wiener residual be dominated by the out of band signals. We don't want the filter to try to do anything outside of the 0.5 - 10 Hz band.
• Somehow, we want to assign a penalty for the filter to have high frequency gain. We do NOT want to slap on an ad-hoc low pass filter. The point of the Wiener filtering is to make the optimum.
• What is the reason for the poor filter performance from 0.5 - 2 Hz ? If we use the frequency domain (Dmass) subtraction technique, we can do better, so there's some inefficiency in this process.
• we're getting too much of the 3 Hz stack mode coupling into MCL. I think this means that our damping filters should be using RG around the suspension eigenmodes rather than just simple velocity damping. We had this years ago, but it caused some weird interaction with the angular loops...to be puzzled out.
11520   Thu Aug 20 11:31:55 2015 KojiUpdatePEMGur interface box

As reported previously, the transfer functions of the channels look fine. (i.e. All channels almost identical.)
I checked the chain from the unit input to the DAQ BNC connectors. They were all OK.

Today I have been checking the signals on the unit with the long DB37 cables connected.
I could not see anything on the Gur2 channels on the board. I looked at the DB37 for Gur2 and felt something is wrong.

I opened the housing of the cable and realized that all the pins are not fully inserted.
The wires were crimped improperly and prevents them to be fully inserted.

=> We need to redo the crimping to insert them.
=> We need to check the other side too.

Attachment 1: IMG_1958.JPG
11521   Thu Aug 20 18:08:28 2015 IgnacioFrogs40m upgradingFatality. Something broke.

So I made coffee at 1547 and was astonished to find the above. Its a sad, very sad day.

At first I thought that something (a gravity wave?) or someone, accidentally hit the thing and it fell and broke. But Koji told me that the janitor was cleaning around the thing and it did indeed fell accidentally.

11522   Fri Aug 21 08:38:37 2015 SteveUpdateIOOPMC locked

PMC needed to be locked manually.

Attachment 1: PMClocked.png
11523   Fri Aug 21 17:15:13 2015 SteveUpdatePEMGur interface box

 Quote: As reported previously, the transfer functions of the channels look fine. (i.e. All channels almost identical.) I checked the chain from the unit input to the DAQ BNC connectors. They were all OK. Today I have been checking the signals on the unit with the long DB37 cables connected. I could not see anything on the Gur2 channels on the board. I looked at the DB37 for Gur2 and felt something is wrong. I opened the housing of the cable and realized that all the pins are not fully inserted. The wires were crimped improperly and prevents them to be fully inserted. => We need to redo the crimping to insert them. => We need to check the other side too.

I'm making a new long cable. Both connector ends of this X arm long cable were terrible. It was removed from the cable tray yesterday.

The Y arm Guralp is running fine. The interface box is open in front of 1X1 on a cart and it is alive!    Please be aware of it !

ELOG V3.1.3-