18 September 2018 Process IF data to create new KL Basis Set with binning=4

17:42 Start KL analysis
  • aoidl
  • @startup
IDL> TN=['20180912_022930','20180912_023204','20180912_023410','20180912_023548','20180912_023823','20180912_024037','20180912_024314','20180912_024423','20180912_024545','20180912_024659','20180912_024818','20180912_024937','20180912_025043','20180912_025330','20180912_025504','20180912_025606','20180912_025705','20180912_025808','20180912_025916','20180912_030026','20180912_030432','20180912_030539','20180912_030642','20180912_030742','20180912_030843','20180912_031001','20180912_031108','20180912_031419','20180912_031517','20180912_031616','20180912_031718','20180912_031816','20180912_031941','20180912_032148','20180912_032304','20180912_032359','20180912_032504','20180912_032600','20180912_032706']
  • IDL> mm=obj_new('mirmodes', TN, binning=4)
IDL> mm->summary_opds
tracknum 		 nn 	 amp [nm] 	 first 	  last 	 missing    sampling
20180912_022930 	 125 	  500.0 	 0 	 4 	 0 	 990 x 1002
20180912_023204 	 125 	  350.0 	 5 	 9 	 0 	 990 x 1002
20180912_023410 	 51 	  250.0 	 10 	 19 	 0 	 990 x 1002
20180912_023548 	 35 	  150.0 	 20 	 37 	 0 	 990 x 1002
20180912_023823 	 35 	  150.0 	 38 	 55 	 0 	 990 x 1002
20180912_024037 	 35 	  150.0 	 56 	 73 	 0 	 990 x 1002
20180912_024314 	 35 	  150.0 	 74 	 91 	 0 	 990 x 1002
20180912_024423 	 35 	  100.0 	 92 	 109 	 0 	 990 x 1002
20180912_024545 	 35 	  100.0 	 110 	 127 	 0 	 990 x 1002
20180912_024659 	 35 	   80.0 	 128 	 145 	 0 	 990 x 1002
20180912_024818 	 35 	   80.0 	 146 	 163 	 0 	 990 x 1002
20180912_024937 	 35 	   80.0 	 164 	 181 	 0 	 990 x 1002
20180912_025043 	 35 	   80.0 	 182 	 199 	 0 	 990 x 1002
20180912_025330 	 35 	   70.0 	 200 	 217 	 0 	 990 x 1002
20180912_025504 	 35 	   50.0 	 218 	 235 	 0 	 990 x 1002
20180912_025606 	 35 	   50.0 	 236 	 253 	 0 	 990 x 1002
20180912_025705 	 35 	   50.0 	 254 	 271 	 0 	 990 x 1002
20180912_025808 	 35 	   50.0 	 272 	 289 	 0 	 990 x 1002
20180912_025916 	 35 	   40.0 	 290 	 307 	 0 	 990 x 1002
20180912_030026 	 35 	   40.0 	 308 	 325 	 0 	 990 x 1002
20180912_030432 	 35 	   40.0 	 326 	 343 	 0 	 990 x 1002
20180912_030539 	 35 	   40.0 	 344 	 361 	 0 	 990 x 1002
20180912_030642 	 35 	   40.0 	 362 	 370 	 0 	 990 x 1002
20180912_030742 	 35 	   40.0 	 371 	 379 	 0 	 990 x 1002
20180912_030843 	 35 	   40.0 	 380 	 397 	 0 	 990 x 1002
20180912_031001 	 35 	   40.0 	 398 	 415 	 0 	 990 x 1002
20180912_031108 	 35 	   40.0 	 416 	 433 	 0 	 990 x 1002
20180912_031419 	 35 	   40.0 	 434 	 451 	 0 	 990 x 1002
20180912_031517 	 35 	   40.0 	 452 	 469 	 0 	 990 x 1002
20180912_031616 	 35 	   40.0 	 470 	 487 	 0 	 990 x 1002
20180912_031718 	 35 	   40.0 	 488 	 505 	 0 	 990 x 1002
20180912_031816 	 35 	   40.0 	 506 	 523 	 0 	 990 x 1002
20180912_031941 	 35 	   40.0 	 524 	 541 	 0 	 990 x 1002
20180912_032148 	 35 	   40.0 	 542 	 559 	 0 	 990 x 1002
20180912_032304 	 35 	   40.0 	 560 	 577 	 0 	 990 x 1002
20180912_032359 	 35 	   40.0 	 578 	 595 	 0 	 990 x 1002
20180912_032504 	 35 	   40.0 	 596 	 613 	 0 	 990 x 1002
20180912_032600 	 35 	   40.0 	 614 	 631 	 0 	 990 x 1002
20180912_032706 	 35 	   40.0 	 632 	 645 	 0 	 990 x 1002
  • IDL> mm->collect_opds
  • IDL> mm->mask_calibration, /SHOW
    • images
  • IDL> mm->produce_ifmat
  • IDL> mm->check_noise
    • images
  • IDL> mm->inv_ifmat
  • IDL> mm->interactive_tweak_mask
    • tag "good"
    • images
  • IDL> mm->tweak_mask, 'better'
  • IDL> mm->tweak_ifmat, 'better'
  • IDL> mm->showme_the_modes, nrows=16, ncol=20, zoom=0.22
  • IDL> list=indgen(327)+320
  • IDL> mm->showme_the_modes, mode_num_idx=list, nrows=16, ncol=20, zoom=0.22
    • images
  • IDL> klf = obj_new('klfit', mm, 'better')
  • IDL> klf->create_theo_kls
  • IDL> save, klf, file='/local/aomeas/adsec_calib/if_functions/20180912_mirmodes/klf_1.sav'
  • IDL> klf->create_ortho_fitted_kls
  • IDL> save, klf, file='/local/aomeas/adsec_calib/if_functions/20180912_mirmodes/klf_2.sav'
  • IDL> klf->check_fitting
    • 405 modes used for fitting
    • images

20180918

  • IDL> restore, '/local/aomeas/adsec_calib/if_functions/20180912_mirmodes/KLdata_fit_better_/force_data.sav', /v
  • IDL> for i=0,12 do begin display, forcemat[i,*], adsec.act_w_cl, title=strtrim(string(i),2), /sh, /no_num, min_v=-1.e5, max_v=2.e5, /no_back, /smooth & wait, 3.0 & endfor
  • 20180918_12_modes_bin_4.mov: 20180918 12 Mode Force Plot with Bin=4

-- %USERSIG{DougMiller - 2018-09-18}%

Comments

 

  • KL Modes fitting error and peak force:
    peak_force.jpg

  • KL Modes fitting error and peak force:
    fitting_error.jpg
I Attachment Action Size Date Who Comment
20180918_12_modes_bin_4.movmov 20180918_12_modes_bin_4.mov manage 10 MB 18 Sep 2018 - 19:02 DougMiller 20180918 12 Mode Force Plot with Bin=4
fitting_error.jpgjpg fitting_error.jpg manage 84 K 18 Sep 2018 - 19:03 DougMiller KL Modes fitting error and peak force
peak_force.jpgjpg peak_force.jpg manage 89 K 18 Sep 2018 - 19:03 DougMiller KL Modes fitting error and peak force
Topic revision: r1 - 18 Sep 2018, DougMiller
This site is powered by FoswikiCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding Foswiki? Send feedback