40m QIL Cryo_Lab CTN SUS_Lab TCS_Lab OMC_Lab CRIME_Lab FEA ENG_Labs OptContFac Mariner WBEEShop
  40m Log, Page 331 of 339  Not logged in ELOG logo
New entries since:Wed Dec 31 16:00:00 1969
Entry  Wed Nov 27 17:51:52 2019, gautam, Update, Wiener, MCL FF status MCL_FF.pdfPRC_FF.pdf
The old MCL filters are not completely useless - I find a factor of ~2 reduction in the MCL RMS when I turn the FF on. It'd be interesting to see
how effective the FF is during the periods of enhanced seismic activity we see. I also wonder if this means the old PRC angular FF filters are also working,
it'd help locking, tbc with PRMI carrrier...
Entry  Tue Mar 24 19:41:57 2020, gautam, Update, Wiener, Seismic feedforward for MCL IMCseisFF.pdffilterComp.pdfoldFilter_v_proposed.pdfMCL_ff_performance.pdf

I think the feedforward filters used for stabilizing MCL with vertex seismometers would benefit from a retraining (last trained in Sep 2015). 

    Reply  Wed Apr 1 00:51:41 2020, gautam, Update, Wiener, Slightly improved MCL FF MC2_act_calib.pdfIIR_fit_to_FIR.pdfFIRvIIR.pdf

Retraining the MCL filters resulted in a slight improvement in the performance. Compared to no FF, the RMS in the 0.5-5 Hz range
is reduced by approximately a factor of 3. 
Entry  Thu Dec 17 02:57:08 2009, Jenne, Update, WienerFiltering, L1 DARM Static Wiener Filtered data L1darmCompPlot_17Dec2009_4daysLong.png
This is perhaps best put in the LLO elog, but I'm not yet a 'person' there, so I can't write to their elog (yet another thing for the eternal to-do list). 
So for now, we're putting things here...
This isn't totally finalized, but I do want to get what I have posted before I hop on a plane in the morning.  Mostly it just needs more time
    Reply  Thu Dec 17 07:47:29 2009, Jenne, Update, WienerFiltering, L1 DARM Static Wiener Filtered data L1darmComp_17Dec2009_6day_residualsNormSurfacePlot.pngL1darmComp_17Dec2009_6day_ratioBLRMS.pngL1darmComp_17Dec2009_6day_rawBLRMS.pngL1darmComp_17Dec2009_6day_residualsBLRMS.png
This surface plot is the same as the previous one, with a little more data than I had previously.  
This time around, I also include the "BLRMS" plots for this data.  The first one takes each residual and normalizes it by the DARM_CTRL
signal at that time, separates the spectra into bands, and integrates underneath the spectra within that band.  The second one is the raw DARM_CTRL
Entry  Tue Jan 5 01:31:09 2010, Jenne, Update, WienerFiltering, New Wiener Filters installed in PEM IIR matrix on OAF screen MCwino-FFtest_4Jan10_sameFiltersAsinEPICS.pngOAF-FF_test_4Jan10.png
Using the techniques employed at LLO, and then by Rana here at the 40m a few weeks ago, Wiener filters have been installed on the inputs of all of the
PEM IIR channels which are hooked up to the 110B PEM ADCU.  Some slight modifications have taken place to the code, and it's all been checked in to
the 40m svn. 
    Reply  Tue Jan 5 17:32:59 2010, Jenne, Update, WienerFiltering, New Wiener Filters installed in PEM IIR matrix on OAF screen OAF-FF_test_5Jan10.png
EDIT 6 Jan 2010:  Shouldn't have
done this.  My bad.  The AA32 is on the other PEM matrix because the Adaptive code runs at 64Hz, so there's downsampling, calculating,
and upsampling which goes on.  The Feed Forward path all runs at 2kHz, the regular rate of the ASS/OAF machine.  All of these filters are turned
Entry  Sun Mar 7 07:01:21 2010, rana, Update, WienerFiltering, Guralp Huddle Test software huddlez.png
We need to do a new huddle test of the Guralps for the Wiener filtering paper. The last test had miserable results.
I tried to use recent data to do this, but it looks like we forgot to turn the Guralp box back on after the power outage or that they're far off
    Reply  Sat Mar 13 22:16:03 2010, rana, Update, WienerFiltering, Guralp Huddle Test software guralpznoise.png
I used some recent better data to try for better Z subtraction.
Dmass helped me understand that sqrt(1-Coherence) is a good estimate of the theoretical best noise subtraction residual. This should be added to
DTT. For reference the Jan statistic is the inverse of this.
Entry  Mon Aug 30 18:35:22 2010, Jenne, Update, WienerFiltering, Secret Hiding Place for Raw Data 
As it turns out, data seems to fall off the 16Tb drives after ~20 days.  Which makes it a good thing that I saved all of my raw data from my good
Mode Cleaner / seismic weekend for offline Wiener Filtering in the following secret place:
Entry  Sat Jan 29 03:13:44 2011, rana, Update, WienerFiltering, Improvement in H1 Wiener FF prediction by using weights and taps darmweight.pngdarm3000.png
(Jenne, Rana)
Tonight we noticed that there were significant improvements to be had in the predicted DARM Wiener filtering FF performance by using weighting
filters and more taps in the FIR filter.
    Reply  Mon Jan 31 10:31:30 2011, josephb, Update, WienerFiltering, Improvement in H1 Wiener FF prediction by using weights and taps 
Rossa is a rather beefy machine. It effectively has 8 Intel i7 Cores (2.67 Ghz each) and 12 Gigs of ram.  Megatron only has 8 Gigs of ram and just
8 Opterons (1 GHz each).  Rosalba has 4 Quad Core2  (2.4 GHz) with only 4 Gigs of ram. 

Entry  Thu Feb 24 10:56:04 2011, Jenne, Update, WienerFiltering, One month of H1 S5 data is now on Rossa 
Just in case anyone else wants to access it, we now have 30 days of H1 S5 DARM data sitting on Rossa's harddrive.  It's in 10min segments. 
This is handy because if you want to try anything, particularly Wiener Filtering, now we don't have to wait around for the data to be fetched from elsewhere.
Entry  Thu Jul 14 17:34:35 2011, Ishwita, Manuel, HowTo, WienerFiltering, Wiener-Hopf equations derivwf.pdf
Since we are using Wiener filtering in our project, we studied the derivation of Wiener-Hopf equations. Whatever we understood we have written it as
a pdf document which is attached below...
Entry  Mon Aug 1 23:29:24 2011, Manuel, Ishwita, Update, WienerFiltering, Getting Data by matlab 
We tried to acquire data from the seismometers and the mode cleaner using the Matlab function
datalist = NDS2_GetData({'C1:PEM-SEIS_GUR1_X_IN1_DQ'}, 996258376 , 10, CONFIG.nds.C)
and encountered the following error
Entry  Wed Aug 3 02:28:08 2011, Manuel, Ishwita, Update, WienerFiltering, Wiener Filtering in X-arm WFgur1X1Y2X2YN20000srate2048.pngWFgur2XN20000srate2048.pngWFgur2YN20000srate2048.pngSchermata-6.pngSchermata-7.png
Wiener Filtering was applied on the data collected from the X-arm during the time: GPS time-996380715 (Aug 02, 2011. 21:25:00. PDT) to GPS time-996382215 (Aug 02, 2011. 21:50:00. PDT) for a duration of 1500 seconds. During this time the X-arm was locked, we checked it by acquiring
Entry  Wed Aug 3 12:24:08 2011, Manuel, Ishwita, Update, WienerFiltering, Weekly summary 
Last Friday (Jul 29) we reinstalled the blue breakout box, and changed the names of the C1:PEM channels. Elog
We continued the work on the simulation ad applied wiener filter on the simulated ground motion, but the result is unsatisfactory, yet. We will
Entry  Tue Aug 9 16:21:59 2011, Manuel, Ishwita, Update, WienerFiltering, First results of our simulations simu1.pngsimu2w.pngTFwhitened2.png
We did the simulation of the stacks by defining a transfer function for one stack  (green plot) and another
similar transfer function for the other stack.
We simulated the ground motion by filtering a white noise with a low pass filter with a cutoff frequency at 10Hz. (
Entry  Fri Nov 25 12:46:13 2011, Mirko, Update, WienerFiltering, Wiener filtering tryout WienerFiltering.pdf
Tried the wiener filter with the TF from p.5900
Tried it out with the TFs from p.5900:
Entry  Mon Nov 28 18:53:35 2011, Den, Summary, WienerFiltering, seismic noise substraction 9x
There is still a problem why GUR, STS signals are poorly coherent to MC_L.  But at least we can see coherence at 2-5 Hz. It might be useful to do
something with adaptive filtering because it does not work at all for a long time. We start with Wiener filtering. I still doubt that static filtering
is useful. Adaptive filter output is linear to its coefficients, so why not to provide adaptive filter with a zero approximation equal to calculated Wiener
Entry  Tue Dec 6 00:58:34 2011, Den, Update, WienerFiltering, OAF current goal oaf_goal_psd.pngoaf_goal_coh.png
After reducing the digital noise I did offline Wiener filtering to see how good should be online filter. I looked at the MC_F and GUR1_X and GUR1_Y signals.
Here are the results of the filtering. The coherence is plotted for MC_F and GUR1_X signals.
Entry  Wed Aug 15 00:52:17 2012, Den, Update, WienerFiltering, filter calculation wiener.png
A Matlab script to calculate Wiener filter coefficients and convert fir to iir is ready. Input is a file with zero mean witness and desired signals,
output is a Foton zpk command to specify iir filter.
The plot shows comparison of offline fir , iir and online iir filtering. Spectrum below 4 Hz is still oscillating due to acoustic coupling, this
Entry  Thu Nov 1 20:20:33 2012, Den, Ayaka, Update, WienerFiltering, acoustic noise in PMC psd.pngcoh.png
We've subtracted acoustic noise from PMC using 1 EM 172 microphone. We applied a 10 Hz high-pass filter to PMC length signal and 100,200,300:30,30 to
whiten the signal.We used ~10 minutes of data at 2048 Hz as we did not see much coherence at higher frequencies.
We were able to subtract acoustic noise from PMC length in the frequency range 10-700 Hz. In the range 30-50 Hz error signal is less by a factor
Entry  Mon Nov 5 20:37:01 2012, Ayaka, Update, WienerFiltering, Play with wiener filtering estimate.pngWfilter_notch.png
 I am trying to find what limits the reduction rate with wiener filtering.
I did some calculations below:
Reduction rate estimation by microphone noise
    Reply  Thu Nov 8 10:11:58 2012, Ayaka, Update, WienerFiltering, Play with wiener filtering document.pdfestimate.pngR_model.png
I will attach a document which describes how the noise affect the wiener filter and the noise cancellation ratio.
And I re-estimate the SN ratio in the microphone (but still rough):
Entry  Tue Nov 13 21:05:22 2012, Ayaka, Update, WienerFiltering, Acoustic noise cancellation with MC 9x
 In order to perform acoustic noise cancellation with MCL signal, I am trying to find sweet spots for microphones.
I set microphones at various places around MC chambers, and see how coherent microphones and MC signals are.
I had checked the half part of MC.
Entry  Tue Nov 27 01:03:42 2012, Ayaka, Update, WienerFiltering, Sound on PSL gauss_psdcoh_mcf.pngtri40_psdcoh_mcf.pnggauss_filters_mcf.pngGauss_estimatefilters_mcf.png
 Last Thursday, I put the speaker and my laptop in the PSL table, and make triangular wave sound with the basic frequency of 40Hz, and Gaussian
distributed sound.
(I create the sounds from my laptop using the software 'NHC Tone Generator' because I could not find the connector from BNC to speaker plug.)
    Reply  Wed Nov 28 23:55:13 2012, Ayaka, Update, WienerFiltering, Sound on PSL tri40_psdcoh_mcf_varioust.pnggauss_psdcoh_mcf_varioust.png
 I have been searching for the way we can subtract signal better since I could see the acoustic coupling signal remains in the target signal even
though there are no coherence between them.
I changed the training time which is used to decide wiener filter.
Entry  Wed Dec 5 01:06:22 2012, Ayaka, Update, WienerFiltering, The microphones and the speaker on the AP table DSC_4956.JPGDSC_4961.JPGAPsig.pdf
 In order to see the acoustic coupling on arm signals, I set 6 microphones and the speaker on the AP table. The microphones are not seismically
isolated for now.
I have a signal generator under the AP table.
    Reply  Wed Dec 5 09:53:01 2012, rana, Update, WienerFiltering, The microphones and the speaker on the AP table 
  Don't try to re-invent the mic mount: just copy the LIGO mic mount for the first version.
Entry  Thu Dec 13 18:06:59 2012, Ayaka, Update, WienerFiltering, New microphone mounts IMG_0129.JPGIMG_0130.JPG
 Yesterday, I made new mounts for microphones.
Entry  Mon May 13 17:24:38 2013, Jenne, Update, WienerFiltering, PRM YAW Wiener filtering PRM_SUSYAW_subtractUsing_GUR1_xyz.png
I have done a quicky offline Wiener filter to check how much PRM yaw motion we can subtract using a seismometer in the corner station.  This work
may be redundant since Koji got the POP beam shadow sensor feedback loop working on Friday night.
Anyhow, for now, I used the GUR2 channels, since GUR2 was underneath the ITMX chamber (at the north edge of the POX table).  Note that Zach
Entry  Thu Aug 6 21:03:45 2015, Ignacio, HowTo, WienerFiltering, How to do online static IIR Wiener filtering iir.png
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
Entry  Tue Sep 1 06:58:11 2015, Ignacio, Update, WienerFiltering, MCL FF => WFS1 and WFS2 FF => ARMS FF 6x
I took some training data during Sunday night/Monday morning while the MCL MISO FF was turned on. We wanted to see how
much residual noise was left in the WFS1/WFS2 YAW and PITCH signals. 

The offline subtractions that can be achieved are:
Entry  Mon Dec 14 16:46:42 2015, ericq, Update, WienerFiltering, Noise Subtraction Puzzler subpuzz.pdf
Here's something to ponder.

Our online MCL feedforward uses perpendicular vertex T240 seismometer signals as input. When designing a feedforward filter, whether FIR Wiener
or otherwise, we posit that the PSD of the best linear subtraction one can theoretically achieve is given by the coherence, via Psub = P(1-C). 
    Reply  Thu Jan 7 15:29:54 2016, ericq, Update, WienerFiltering, Noise Subtraction Puzzler subpuzz2.pdfpuzzle.zip
The puzzle continues...

I found some reference for computing "multicoherence," which should properly estimate the potential MISO subtraction potential in situations
where the witness channels themselves have nontrivial coherence. Specifically, I followed the derivations in LIGO-P990002.
    Reply  Thu Jan 14 15:08:37 2016, ericq, Update, WienerFiltering, Noise Subtraction Puzzler subpuzz_resolved.pdf
The anticlimatic resolution to my subtraction confusion: Spectral leakage around 1Hz. Increasing the FFT length to 256 sec now shows that the FIR
WF pretty much achieves the ideal subtraction. 

    Reply  Thu Jan 14 18:41:36 2016, rana, Update, WienerFiltering, Noise Subtraction Puzzler 
Just not just pedagogical !  Freq domain MISO coherence based subtraction estimation is much faster than calculating MISO WF. And since each bin
is independent of each other, this gives us an estimate of how low the noise can go, whereas the Wiener filter is limited by Kramers-Kronig. We should
be able to use this on the L1 DARM channel to do the noise hunting as well as estimating the subtraction efficacy of the pseudo channels that you and Rory
Entry  Fri Apr 1 13:11:26 2016, ericq, Update, WienerFiltering, Bilinear Noise Testing bilinear_drive.pdf40m_bilin.pdf
I've been banging my head against bilinear noise subtraction, and figured I needed to test things on some real hardware to see if what I'm doing
makes sense.

I ran the ASS dither alignment on the Y arm, which ensures that the beam spots are centered on both mirrors. 
Entry  Tue Jun 10 11:54:36 2014, Zach, Configuration, Wiki, DokuWikis are still down 
Not sure if someone is already on this, but the nodus DokuWikis are still down (AIC, ATF, Cryo, etc.)
I get:

DokuWiki Setup Error
    Reply  Tue Jun 10 19:11:27 2014, Evan, Configuration, Wiki, DokuWikis are back up 
As of this writing, the DokuWiki appears to be
As you and I suspected, it looks like this was a clusterwhoops with the permissions for the NFS mount. Let's recap what happened in the past 24
    Reply  Wed Jun 11 05:16:14 2014, Zach, Configuration, Wiki, DokuWikis are back up 
It looks like auth is broken on the AIC wiki (though working fine on ATF and Cryo). I did some poking around
but can't see how anything we did could have broken it.

    Reply  Wed Jun 11 10:15:15 2014, Evan, Configuration, Wiki, DokuWikis are back up 


           It looks like auth is broken on the AIC wiki (though working fine on ATF and Cryo). I did some poking
    Reply  Wed Jun 11 20:09:37 2014, Zach, Configuration, Wiki, DokuWiki situation 
I have been looking into the DokuWiki situation, with no great progress so far.
The problem is with authentication. When you try to access the wiki, you get: "User authentication is temporarily unavailable. If this situation persists, please inform your Wiki
Entry  Fri Aug 8 16:45:40 2014, Nichin, Update, Wiki, PDFR wiki updated 
 The PDFR system has been documented in the 40m wiki and all the relevant information about making changes and keeping it updated have been mentioned.
This pretty much wraps up my SURF 2014 project at the 40m lab. 
    Reply  Mon Aug 11 23:50:38 2014, rana, Configuration, Wiki, DokuWikis are back up 


Entry  Thu Aug 14 20:52:36 2014, rana, Update, Wiki, Violin Mode table added to Wiki 
Mech Resonance Wiki
I've updated the wiki by trawling the elog for violin entries. Please keep it up to date so that we can make violin notches.
    Reply  Mon Apr 10 15:22:46 2017, rana, Configuration, Wiki, DokuWikis are back up 
AIC Wiki updated to latest stable version of DokuWiki: 2017-02-19b "Frusterick Manners" + CAPTCHA + Upgrade + Gallery PlugIns
Entry  Thu Apr 19 09:57:03 2018, gautam, Update, Wiki, AP and ETMX tables uploaded to wiki 
The most up to date pictures of the AP table and ETMX table that Steve took have been uploaded to the relevant
page on the wiki. It seems like the wiki doesn't display previews of jpg images - by using png, I was able to get the thumbnail of the attachment
to show up. It would be nice to add beam paths to these two images. The older versions of these photos were moved to the archive section on the same page.
    Reply  Fri Apr 20 16:25:08 2018, Steve, Update, Wiki, ETMX table layout uploaded to wiki 
ETMX table layout uploaded with beam paths to the wiki.   
ELOG V3.1.3-