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 Jul 13 16:10:34 2017, Kaustubh, Summary, Computer Scripts / Programs, Cavity Scan Simulation Code 
    Reply  Mon Jul 31 23:39:02 2017, rana, Update, COC, Cavity Scan Simulation Code manyCavityScans-SIS.pdfmanyCavityScans-SIS.pdfMonteCarlo_CavityScans.mat
Message ID: 13155     Entry time: Mon Jul 31 23:39:02 2017     In reply to: 13116
Author: rana 
Type: Update 
Category: COC 
Subject: Cavity Scan Simulation Code 

Hiro Yamamoto has updated SIS (Static Interferometer Simulation) to allow us to do the MCMC based inference of the 40m arm cavity mirror maps. 

The latest version is in git.ligo.org: IFOsim/SIS/

In the examples directory I have put 3 files:

  1. mcmcCavityScans.m - runs many cavity scans using parfor and saves the data
  2. plotCavityScans.m - loads the .mat file with the data and plots it
  3. plotCavityScans.py - python file which also loads & plots, but nicer since python has a transparency option for the traces.

Attached is the plots and the data. The first attached plot is a low resolution one: 200 scans of 100 frequency points each. Second plot is 200 scans of 300 points each.

The run was done assuming perfect LIGO arm params with a random set of Zernike perturbations for each run. The amplitude of each Zernike was chosen from a Normal distribution with a standard deviation of 10 nm.

We need to come up with a better guess for the initial distribution from which to sample, and also to use the more smart sampling that one does using the MCMC Hammer.

Attachment 1: manyCavityScans-SIS.pdf  134 kB  Uploaded Tue Aug 1 01:03:00 2017  | Hide | Hide all
Attachment 2: manyCavityScans-SIS.pdf  197 kB  Uploaded Tue Aug 1 01:23:28 2017  | Hide | Hide all
Attachment 3: MonteCarlo_CavityScans.mat  564 kB  Uploaded Tue Aug 1 01:23:28 2017
ELOG V3.1.3-