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  Thu Dec 1 20:20:15 2016, gautam, Update, IMC, IMC loss measurement plan IMG_3471.JPG
    Reply  Fri Dec 2 16:40:29 2016, gautam, Update, IMC, 24V fuse pulled out 
       Reply  Mon Dec 5 01:58:16 2016, gautam, Update, IMC, IMC ringdowns IMCringdown.pdfIMCringdown_2.pdf
          Reply  Mon Dec 5 15:55:25 2016, gautam, Update, IMC, IMC ringdowns 7x
             Reply  Mon Dec 5 19:29:52 2016, gautam, Update, IMC, IMC ringdowns REFLanomaly.pdf
                Reply  Tue Dec 6 00:43:41 2016, gautam, Update, IMC, more IMC ringdowns 
                   Reply  Wed Dec 7 11:52:48 2016, ericq, Update, IMC, Partial IMC ringdowns IMCpartial.pngIMCpartial.zip
                      Reply  Thu Dec 8 19:01:21 2016, rana, Update, IMC, Partial IMC ringdowns 
                   Reply  Sun Dec 2 17:26:58 2018, gautam, Update, IMC, IMC ringdown fitting tauTheoretical.pdfringdownFit.pdfringdownScatter.pdfcavPole.pdf
                      Reply  Fri Dec 7 12:51:06 2018, gautam, Update, IMC, IMC ringdown fitting time_reflAndTrans.pdfcorner_transOnly.pdfcorner_reflAndTrans.pdf
Message ID: 14334     Entry time: Fri Dec 7 12:51:06 2018     In reply to: 14328
Author: gautam 
Type: Update 
Category: IMC 
Subject: IMC ringdown fitting 

I started putting together some code to implement some ideas we discussed at the Tuesday meeting here. Pipeline isn't setup yet, but i think it's commented okay so if people want to play around with it, the code lives on the 40m gitlab

Model parameters:

  • T+ --- average transmission of MC1 and MC3.
  • T- --- difference in transmission between MC1 and MC3 (this basis is used rather than T1 and T3, because the assumption is that since they were coated in the same coating run, the difference in transmission should be small, even if there is considerable uncertainty in the actual average transmission number.
  • T2 --- MC2 transmission.
  • Lrt --- Round trip loss in the cavity.
  • "sigma" --- a nuisance parameter quantifying the error in the time domain ringdown data.


  • Using these model parameters, calculate some simulated time-domain ringdowns. Optionally, add some noise (assumed Gaussian).
  • Try and back out the true values of the model parameters using emcee - priors were assumed to be uniformly distributed, with a +/- 20% uncertainty around the central value.
  • For a first test, see if there is any improvement in the parameter estimation uncertainty using only transmission ringdown vs both transmission and reflection.

Initial results and conclusions:

  • Attachment #1 - Simulated time series used for this study. The "fit" trace is computed using the median values from the monte-carlo.
  • Attachment #2 - Corner plots showing the distribution of the estimated parameter values, using only transmission ringdown. The "true" values are indicated using the thick blue lines.
  • Attachment #3 - Corner plots showing the distribution of the estimated parameter values, using both transmission and reflection ringdowns.
  • The overall approach seems to work okay. There seems to be only marginal improvement in the uncertainty in estimated parameters using both ringdown signals, at least in the simulation.
  • However, everything seems pretty sensitive to the way the likelihood and priors are coded up - need to explore this a bit more.

Next steps:

  • Add more simulated measurements, see if we can constrain these parameters more tightly. 
  • Use linear error analysis to see if that tells us which measurements we should do, without having to go through the emcee.

There still seems to be some data quality issues with the ringdown data I have, so I don't think we really gain anything from running this analysis on the data I have already collected - but in the future, we can do the ringdown with complete extinguishing of the input light, and repeat the analysis.

As for whether we should clean the IMC mirrors - I'm going to see how much power comes out at the REFL port (with PRM aligned) this afternoon, and compare to the input power. This technique suffers from uncertainty in the Faraday insertion loss, isolation and IMC parameters, but I am hoping we can at least set a bound on what the IMC loss is.

Attachment 1: time_reflAndTrans.pdf  22 kB  Uploaded Fri Dec 7 14:01:00 2018  | Hide | Hide all
Attachment 2: corner_transOnly.pdf  919 kB  Uploaded Fri Dec 7 14:01:08 2018  | Hide | Hide all
Attachment 3: corner_reflAndTrans.pdf  890 kB  Uploaded Fri Dec 7 14:01:21 2018  | Hide | Hide all
ELOG V3.1.3-