OnOff观测

方法一 和 方法二 主要区别在 是否考虑驻波相位的漂移。

方法一

1. 用 hifast.pos_swi 分离噪音管on、off 和 源（src）和参考点（ref）观测，然后噪音管定标，src减去ref后合并成一条谱线（保留两个偏振）。

M1-S1.sh

#!/bin/bash

# enable **/
shopt -s globstar

# only need input chunk *0001.fits
files="$(ls RAW_data/*/*/*M01*W_0001.fits)"
# check files
printf "${files}"
echo
# commands
# python -m hifast sep -p 3 同时执行5个 python -m hifast.sep
python -m hifast pos_swi -p 5 \
    $files \
   -d 4 -m 4 -n 36 \
   --frange 1378  1410  \
   --smooth gaussian --s_sigma 5 \
   --noise_mode high \
   --t_src 120 --t_ref 120 --t_change 30 --n_repeat 1 \
   --only-off True \
   --outdir 'output/%(project)s/%(date)s'

2. 手动去基线，sin拟合驻波。

下载: M1-S2.ipynb

3. 用 hifast.flux 流量定标 和 hifast.multi --fc True 坐标系修正。

方法二

1. 先不区分 源（src）和参考点（ref）观测。用跟mapping一样的方法来定标。

   M2-S1.sh

   #!/bin/bash
   
   # enable **/
   shopt -s globstar
   # 1. Temperature calibration
   echo "Temperature calibration:"
   # only need input chunk *0001.fits
   files="$(ls RAW_data/*/*/*M01*W_0001.fits)"
   # check files
   printf "${files}"
   echo
   # commands
   # python -m hifast sep -p 3 同时执行5个 python -m hifast.sep
   python -m hifast sep -p 5 \
       $files \
      -d 4 -m 4 -n 36 \
      --frange 1369  1419  \
      --smooth gaussian --s_sigma 2 \
      --check_cal A --pcal_vary_lim_bin 0.02 \
      --merge_pcals True --method_merge median --method_interp nearest \
      --save_pcals True \
      --outdir 'output2/%(project)s/%(date)s'
   
   # RA-DEC
   # # only need beam 01 specs_T.hdf5 file
   python -m hifast radec **/*-M01*-specs_T.hdf5
   

2. 用fft方法对每条谱线去驻波

   M2-S2.sh S2-sw.ini

   #!/bin/bash
   
   # enable **/
   shopt -s globstar
   
   files="$(ls output2/*/202*/*M01*-specs_T.hdf5)"
   # check fpatten
   printf "${files}"
   echo
   
   # using hifast.sh to run multi files
   # commands
   commands=$(cat <<EOF
   # limit frange again
   python -m hifast.bld  | --nproc 5 \
                           --frange 1369 1419 \
                           --method PLS-asym2 --lam 1e9 \
                           --njoin 40 \
                           --s_method_freq gaussian --s_sigma_freq 3 \
                           --exclude_type auto2 \
                           --post_method poly-asym2 --post_deg 3 \
                           --post_s_method_freq gaussian --post_s_sigma_freq 3 \
                           --post_exclude_type auto2
   python -m hifast.sw   | -c conf/S2-sw.ini --nobld True
   python -m hifast.flux  |  
   EOF
   )
   # check commands
   echo "commands: ""$commands"
   # run multiple files
   hifast.sh "$files" -c "$commands" -n 4
   

3. 用 hifast.pos_swi_2 分离源（src）和参考点（ref）观测，然后src减去ref后合并成一条谱线（保留两个偏振）。

   M2-S3.sh

   #!/bin/bash
   
   # enable **/
   shopt -s globstar
   
   # only need input chunk *0001.fits
   files="$(ls -v output2/*/202*/*M01*-sw_nobld-flux.hdf5)"
   # check files
   printf "${files}"
   echo
   # commands
   python -m hifast pos_swi_2 -p 5 \
       $files \
      --frange 1378  1400  \
      --t_src 120 --t_ref 120 --t_change 30 --n_repeat 1 \
      --only-off True
      
   

4. 手动去基线

   下载: M2-S4.ipynb

5. hifast.multi --fc True 坐标系修正。