40m QIL Cryo_Lab CTN SUS_Lab TCS_Lab OMC_Lab CRIME_Lab FEA ENG_Labs OptContFac Mariner WBEEShop
  40m Log  Not logged in ELOG logo
Entry  Sun Apr 29 22:53:06 2018, gautam, Update, General, DARM actuation estimate darmEst.pdf
    Reply  Tue May 1 19:37:50 2018, gautam, Update, General, DARM actuation estimate darmEst.pdfdarmEst_time.pdf
       Reply  Mon May 7 16:23:06 2018, gautam, Update, General, DARM actuation estimate darmProj.pdf37.pngMCalign_20180507.png
          Reply  Mon May 7 20:01:14 2018, Rorpheus, Update, General, Use anti-dewhitening + show CARMA/DARMA 
             Reply  Thu May 10 14:13:22 2018, gautam, Update, General, More refinement of DARM control signal projection darmProj.pdfdarmOLTF.pdfDWcomparison.pdf
                Reply  Fri May 11 13:58:42 2018, rana, Update, General, More refinement of DARM control signal projection 
                   Reply  Fri May 11 19:02:52 2018, gautam, Update, General, More refinement of DARM control signal projection darmProj.pdfdarmOLTF.pdfDWcomparison.pdfDARMsuppression.pdfControlLoop.pdf
                      Reply  Sat May 12 10:02:03 2018, rana, Update, General, More refinement of DARM control signal projection 
                         Reply  Sun May 20 17:43:01 2018, rana, Update, Electronics, How to choose resistors 
                      Reply  Mon May 14 18:58:32 2018, Kevin, Update, PonderSqueeze, Squeezing with no SRM homodyne_heatmap.pdfdisplacement_noise.pdfnoise_budget.pdf
                         Reply  Tue Oct 2 23:57:16 2018, gautam, Update, PonderSqueeze, Squeezing scenarios PonderSqueeze_NB_LORIN.pdf
                   Reply  Tue May 15 21:56:57 2018, gautam, Update, General, Stack measurement setup decommissioned 
                      Reply  Thu May 17 09:14:38 2018, Steve, Update, General, Stack measurement setup decommissioned Load_Cell_Measurement_Set_Up.jpgETMX_stack_up_down.png
Message ID: 13805     Entry time: Tue May 1 19:37:50 2018     In reply to: 13799     Reply to this: 13822
Author: gautam 
Type: Update 
Category: General 
Subject: DARM actuation estimate 

Here is an updated plot - the main difference is that I have added a trace that is the frequency domain wiener filter subtraction of the coherent power between the L_X and L_Y time series. I tried reproducing the calculation with the time domain wiener filter subtraction as well, using half of the time series (i.e. 5 mins) to train the wiener filter (with L_X as target and L_Y as witness), but I don't get any subtraction above 5 Hz on the half of the data that is a test data set. Probably I am not doing the pre-filtering correctly - I downsampled the signal to 1 kHz, weighted it by low passing the signal above 40 Hz and trained the Wiener filter on the resulting time series. But this frequency domain Wiener filter subtraction should be at least a lower bound on DARM - indeed, it is slightly lower everywhere than simply taking the time domain subtraction of the two data streams.

To do:

  • Re-measure calibration numbers used.
  • Redo calculation once the numbers have been verified.

Putting a slightly cleaned up version of this plot in now - I'm only including the coherence-inferred DARM estimate now instead of the straight up time domain subtraction. So this is likely to be an underestimate. At low (<10 Hz) frequencies, the time domain computation lines up fairly well, but I suspect that I am getting huge amounts of spectral leakage (see Attachment #2) in the way I compute the spectrum using scipy's filtering routine (once the Wiener filter has been computed). Note that Attachment #2, I didn't break up the data into a training/testing set as in this case, we just care about the one-off offline performance in order to get an estimate of DARM.

The python version of the wiener filter generating code only supports [b,a] output of the digital filter, an sos filter might give better results. Need to figure out the least painful way of implementing the low-noise digital filtering in python...

Attachment 1: darmEst.pdf  332 kB  Uploaded Wed May 2 00:16:19 2018  | Hide | Hide all
Attachment 2: darmEst_time.pdf  358 kB  Uploaded Wed May 2 00:16:43 2018  | Hide | Hide all
ELOG V3.1.3-