![[PukiWiki]](http://www2.nao.ac.jp/~aste-www/aste_wiki/obs.jpg "[PukiWiki]")
[[Commissioning and Science Verification]] * Subref position offsets for single-beam receivers > Analysis [#x6970f30] :Revision history| -- 2016-08-01 - created by Kamazaki,T. :Table of contents| #contents ---- &size(16){&color(red){Since the analysis tools are still under preparation, it is UNNECESSARY to reduce obtained data.};}; ** Analysis [#b508adf1] + Login to an analysis computer (e.g. aste-ms1) with account=''asteXX'' (e.g. aste16) + Rsync CONT data from the ASTE site asteXX@aste-ms1: ~/analysis/bin/cont_rsync.sh <project> <obs. table> <time (YYYY-MM-DD)> + e.g.) asteXX@aste-ms1: ~/analysis/bin/cont_rsync.sh sbrf16 stsbd3 2016-08-24 ... CONT.stsbd3a.aste16.sbrf16.201608240044.LCK ... + Rsync antenna logs and also make a list of observation tables asteXX@aste-ms1: ~/analysis/bin/sbref_mkList.sh <project> <obs. list (e.g. stsbd3)> <start_time (YYYY-MM-DDThh:mm)> <end_time (YYYY-MM-DDThh:mm)> + e.g.) asteXX@aste-ms1: ~/analysis/bin/sbref_mkList.sh sbrf16 stsbd3 2016-08-24T00:00 2016-08-24T03:30 ... * ~/analysis/beam/bin/sbref_mkList.py 2016-08-23T23:55 2016-08-24T03:35 > 20160824.list --- 20160824.list --- stsbd3z CONT.stsbd3z.aste16.sbrf16.201608240019 real20160824001819 ... stsbd3z CONT.stsbd3z.aste16.sbrf16.201608240321 real20160824031956 -- "20160824.list" will be created in this example. + Edit the created list of observation tables so as to include only subref position-offset measurements + e.g.) --- 20160824.list --- stsbd3z CONT.stsbd3z.aste16.sbrf16.201608240019 real20160824001819 <--- remove stsbd3z CONT.stsbd3z.aste16.sbrf16.201608240041 real20160824004012 <--- 1st loop start stsbd3a CONT.stsbd3a.aste16.sbrf16.201608240044 real20160824004322 stsbd3b CONT.stsbd3b.aste16.sbrf16.201608240051 real20160824004957 stsbd3c CONT.stsbd3c.aste16.sbrf16.201608240058 real20160824005658 stsbd3d CONT.stsbd3d.aste16.sbrf16.201608240105 real20160824010357 stsbd3e CONT.stsbd3e.aste16.sbrf16.201608240112 real20160824011056 stsbd3a CONT.stsbd3a.aste16.sbrf16.201608240119 real20160824011756 stsbd3f CONT.stsbd3f.aste16.sbrf16.201608240126 real20160824012458 stsbd3g CONT.stsbd3g.aste16.sbrf16.201608240133 real20160824013157 stsbd3h CONT.stsbd3h.aste16.sbrf16.201608240140 real20160824013857 stsbd3i CONT.stsbd3i.aste16.sbrf16.201608240147 real20160824014556 stsbd3a CONT.stsbd3a.aste16.sbrf16.201608240154 real20160824015257 <--- 1st loop end stsbd3z CONT.stsbd3z.aste16.sbrf16.201608240201 real20160824015956 <--- 2nd loop start stsbd3a CONT.stsbd3a.aste16.sbrf16.201608240204 real20160824020323 stsbd3b CONT.stsbd3b.aste16.sbrf16.201608240211 real20160824020956 stsbd3c CONT.stsbd3c.aste16.sbrf16.201608240218 real20160824021657 stsbd3d CONT.stsbd3d.aste16.sbrf16.201608240225 real20160824022357 stsbd3e CONT.stsbd3e.aste16.sbrf16.201608240232 real20160824023056 stsbd3a CONT.stsbd3a.aste16.sbrf16.201608240239 real20160824023756 stsbd3f CONT.stsbd3f.aste16.sbrf16.201608240246 real20160824024457 stsbd3g CONT.stsbd3g.aste16.sbrf16.201608240253 real20160824025157 stsbd3h CONT.stsbd3h.aste16.sbrf16.201608240300 real20160824025856 stsbd3i CONT.stsbd3i.aste16.sbrf16.201608240307 real20160824030557 stsbd3a CONT.stsbd3a.aste16.sbrf16.201608240314 real20160824031257 <--- 2nd loop end stsbd3z CONT.stsbd3z.aste16.sbrf16.201608240321 real20160824031956 <--- remove + Run ''subref_redux.sh'' to do Gaussian fit for each scan asteXX@aste-ms1: ~/analysis/bin/subref_redux.sh <project> <list> -- This will create .sbrefx, .sbrefy and .sbrefz files for x-axis, y-axis and z-axis / .strefx and .strefy files for tx-axis, ty-axis, respectively. The files contain Gaussian fit results of each scan in the format described in Appendix. + e.g.) aste16@aste-ms1: ~/analysis/bin/subref_redux.sh sbrf16 20160824.list ... -- Many plots will be displayed, which plots show fit results of each scan toward a planet + Run ''sbref_xyfit.py'' to estimate SubRef positions of the peak amplitude in each measurement loop + Run ''sbref_elfit.py'' to estimate SubRef position parameters ** Appendix [#ic9b3a0e] - ''lck.list'' e.g.) 2016-08-16 00:08:53.551 aste-1c astcnt DEB START 160816001000.000 aste16.sbrf16.stsbb8z< 2016-08-16 00:12:25.150 aste-1c astcnt DEB START 160816001300.000 aste16.sbrf16.stsbb8a< 2016-08-16 00:19:00.050 aste-1c astcnt DEB START 160816002000.000 aste16.sbrf16.stsbb8b< 2016-08-16 00:25:59.651 aste-1c astcnt DEB START 160816002700.000 aste16.sbrf16.stsbb8c< 2016-08-16 00:32:59.350 aste-1c astcnt DEB START 160816003400.000 aste16.sbrf16.stsbb8d< 2016-08-16 00:39:58.850 aste-1c astcnt DEB START 160816004100.000 aste16.sbrf16.stsbb8e< 2016-08-16 00:47:00.250 aste-1c astcnt DEB START 160816004800.000 aste16.sbrf16.stsbb8a< 2016-08-16 00:53:59.650 aste-1c astcnt DEB START 160816005500.000 aste16.sbrf16.stsbb8f< 2016-08-16 01:00:59.450 aste-1c astcnt DEB START 160816010200.000 aste16.sbrf16.stsbb8g< 2016-08-16 01:07:59.450 aste-1c astcnt DEB START 160816010900.000 aste16.sbrf16.stsbb8h< 2016-08-16 01:14:59.150 aste-1c astcnt DEB START 160816011600.000 aste16.sbrf16.stsbb8i< 2016-08-16 01:22:00.150 aste-1c astcnt DEB START 160816012300.000 aste16.sbrf16.stsbb8a< - ''antlog.list'' - Data format of .sbrfx and .sbrfy files x y z az el file peak_amp[0] FWHM[0] peak_amp[1] FWHM[1] daz[0] del[0] daz[1] del[1] fit_error[0,0-3] ... ^^^^^^^^^^^^^^^^^^^ ^^^^^^^^^^^^^^^^^^^ ^^^^^^^^^^^^^ ^^^^^^^^^^^^^ ^^^^^^^^^^^^^^^^ 1st scan 2nd scan 1st scan 2nd scan 1st scan ... - Data format of .sbrfz files [K] [mm] [mm] [K] [mm] [mm] date time az el peak center width d(peak) d(center) d(width) - Data format of .strfx and .strfy files tx ty z az el file peak_amp[0] FWHM[0] peak_amp[1] FWHM[1] ... fit_error[0,0-3] ... ^^^^^^^^^^^^^^^^^^^ ^^^^^^^^^^^^^^^^^^^ ^^^^^^^^^^^^^^^^ 1st scan 2nd scan ... 1st scan ...
