40m QIL Cryo_Lab CTN SUS_Lab TCS_Lab OMC_Lab CRIME_Lab FEA ENG_Labs OptContFac Mariner WBEEShop
  PSL  Not logged in ELOG logo
Entry  Sat Feb 17 18:04:53 2018, Craig, DailyProgress, Mode matching, Set up beam profiler for North path mode matching NorthPathGeneralSetup.jpgSecondNorthPathSteeringMirror.jpgBeamProfilerInstalled.jpgFirstBeamProfilerDataPoint.jpg
    Reply  Fri Feb 23 07:51:51 2018, rana, DailyProgress, Mode matching, emcee Hammer my mode 
       Reply  Fri Feb 23 20:44:12 2018, Craig, DailyProgress, Mode matching, emcee Hammer my mode MCMCAlamodeOptimizedModeMatchingSolutionXAxis.pngNorthPathModeMatchingMCMC.tar
          Reply  Sat Feb 24 21:25:08 2018, rana, DailyProgress, Mode matching, emcee Hammer my mode 
             Reply  Tue Feb 27 03:25:34 2018, rana, DailyProgress, Mode matching, emcee Hammer my mode 
Message ID: 2104     Entry time: Fri Feb 23 20:44:12 2018     In reply to: 2102     Reply to this: 2105
Author: Craig 
Type: DailyProgress 
Category: Mode matching 
Subject: emcee Hammer my mode 

Today I altered beamPath.m in alamode to support gwmcmc.m, which is pretty much emcee for matlab.  (The gw in gwmcmc stands for Goodman-Weare MCMC algorithm, not gravitational waves). 

It's capabilities are pretty limited, all it can do is move around optics.  It cannot set hard limits on where an optic can be placed yet, and it cannot vary a lens' focal length.  It would be easy to add these capabilities for the future.

My log likelihood for the MCMC algorithm is -0.5 * [ weight1 * (1 - overlap)^2 + weight2 * positionSensitivity^2 ], where weight1 = 100 and weight2 = 1.  The user can choose how much to weight each metric which quantifies the quality of the mode match.  Overlap and positionSensitivity were already calculated for me by alamode.

I ran the MCMC on our North path mode matching, allowing two of our lens to move around by about 20 cm for our old lens focal lengths.  Initially, Lens1 was 0.9625 m from the PMC and Lens2 was 1.3717 m.
I get two peaks of different mode matching, which I thought was pretty cool.  My final results for the Lens2 position ended up pretty far away from where I started, which is not physically possible in our setup.  Need to add hard limits to this. 

Not sure how to distribute this.  It's in the tar in this ELOG.  Will also copy this to Git/labutils/alm.  Once it's cleaned up, our future person could make a pull request on github/alm.


This fragility of MM is one of those things we ought to develop a MCMC for. Jenne did something like this for the 40m input mode matching a long time ago, but that was before I understood what corner plots were for.

I wonder if JamMT or alaMode have the ability to be run to generate corner plots. If so, we could find out the sensitivity to lens positions and lens focal length errors in a visually useful way for each mode matching solution (generated from the discrete set of available focal lengths in the CVI catalog).

Might be a project to give to the CTN SURF this summer.


Attachment 1: MCMCAlamodeOptimizedModeMatchingSolutionXAxis.png  181 kB  | Hide | Hide all
Attachment 2: NorthPathModeMatchingMCMC.tar  288 kB
ELOG V3.1.3-