View source for Reference Gain Calibration

Reference complex gain calibration refers to determining the "system" amplitude and phase with high precision as a function of IF band (i.e. integrated over each 500 MHz IF band). A typical reference calibration would cover 34 bands, 2 polarizations, and 13 antennas. The routines to retrieve and analyze the data are in the module refcal_anal.py. Start from a time range during which the reference calibration observation occurred (check the [http://ovsa.njit.edu/phasecal/ phasecal results] to see which 1-hr scan had stable phases)

<pre>
import refcal_anal as ra
from util import Time
#example time range of the refcal scan that took place BEFORE the observation
trange=Time(['2017-04-08T14:00','2017-04-08T15:00'])
</pre>

<!--BIN'S PREVIOUS TEXTS<pre>
import refcal_anal as ra
from util import Time
trange=Time(['2017-04-08T05:00','2017-04-08T15:30'])
</pre>-->

Now read the data from that timerange
<pre>
out=ra.rd_refcal(trange)
</pre>
Note that there are optional arguments that can be given, such as projid='PHASECAL' (the default) and srcid='1229+060', but these are generally not needed.

To correct for differential feed rotation (for refcal data after 05:00 UT on 2017-Jul-01), you can use:
<pre>
out_corr = ra.unrot_refcal(out)
</pre>

Take a look at the phases (averaged over each band): 
<pre>
ra.graph(out_corr)
</pre>
By default it displays only bands 5, 11, 17, 23 and polarization XX. To change the bands and polarization to display, use the "bandplt" and 'pol' parameters (better to check more bands).
<pre>
ra.graph(out_corr, bandplt=[5,13,19,25], pol=1)
</pre>

<!--BIN'S PREVIOUS TEXTS: From the plot, one can determine which scans are the good ones to use and record the indices of useful scans as "scanidx". In this example, we were alternating between the low- and high-frequency receivers. The low-frequency receiver observations did not have delays adjusted. So we use scans with indices 1, 3, 5, 7, 9, 11, 13 which are associated with the high-frequency receiver scans. 
<pre>
scanidx=[1, 3, 5, 7, 9, 11, 13]
</pre>,

Then, we can further determine a finer selection of time range which has clean and stable phases (better to check more bands using ra.graph()). Observations within this time range (after the scanidx selection, if provided) will be used for averaging to get the reference phases/amplitudes. Data will be flagged if the SNR is less than minsnr, which can be supplied. The default minsnr is set to 1.0.

<pre>
timerange=Time(['2017-04-08T05:00','2017-04-08T07:00'])

refcal=ra.refcal_anal(out_corr,timerange=timerange,scanidx=scanidx)
</pre>-->

Observations within this time range will be used for averaging to get the reference phases/amps. Data will be flagged if the SNR is less than minsnr, which can be supplied. The default minsnr is set to 1.0.

<pre>
refcal=ra.refcal_anal(out_corr,timerange=trange,scanidx=[0])
</pre>

It will generate four plots. First and second are phases and amplitudes of the data similar as those produced by ra.graph(), but with the selected time range for averaging highlighted (one can use "bandplt" and "pol" to choose bands and polarization to display as well). The 3rd plot is the averaged phase vs. frequency and the 4th is the averaged amplitude vs. frequency. Those data will be our "analyzed" reference calibration data saved in the returned dictionary (named "refcal" here). The dictionary contains the actual complex values (refcal['refcal']), a flag array with the same shape (refcal['flag']), and some timestamps, e.g., refcal['t_mid'] is the middle of the time range used for averaging. Both refcal['refcal'] and refcal['flag'] have a shape of (15, 2, 34) -- which corresponds to # of antennas, # of polarizations, and # of bands.

If you determine some antenna, band, polarization need to be flagged/unflagged, they have to be manually changed in the returned refcal['flag'] array. In this example, Antenna 12 was not tracking, and the algorithm did not pick up all of them, so we have to do this manually -- setting refcal['flag'][11, :, :]=1

<pre>
refcal['flag'][11]=1
</pre>

If satisfied with the results (by looking at the plots), we can go ahead and send them to the SQL database
<pre>
import cal_header as ch
ch.refcal2sql(refcal)
</pre>

To retrieve the refcal data back from SQL database based on, e.g., a solar data one wish to calibrate. We can provide a time: 
<pre>
t=Time('2017-04-07T19:40')
refcal=ra.sql2refcal(t)
</pre>
The resulted refcal have three keys: 'amp' (15 x 2 x 34 array), 'pha' (15 x 2 x 34 array), and 'timestamp' (in jd). The refcal data found will be the closest one PRIOR to the supplied observation time.

== Updating the Reference Calibrations Table ==
The table below lists the known good reference calibrations.  To facilitate adding entries to this table, do the following, where the time is after the desired refcal, but before any later refcal:
<pre>
t=Time('2017-06-23T19:40')  # retrieves the refcal at 13:07 UT on that day
refcal=ra.sql2refcal(t)
ra.graph_results(refcal, savefigs=True)
</pre>
This will create summary plots of the amplitude and phase, will write the two corresponding *.png files into the /common/webplots/refcal folder, and will print as a text string the lines needed to add an entry. Just cut from the terminal and paste into the wiki editor at the appropriate place to preserve time order.

== List of Analyzed Reference Calibrations ==

{| class="wikitable" style="text-align: center; width: 800px; height: 200px;"
|-
| Date || Timestamp || Source || Obs Trange || Scan Idx || Avg Trange || Bands || Phase Plot || Amp Plot || Comments
|-
| 2017/04/02  || 06:36:50 || 3c273 || 04:20~11:00 || All || 05:47~07:26 || 5, 7, 9, 11, 13, 15, 17, 19, 21, 23 || [[:File:20170402_refcal_pha.png|Phase]] || [[:File:20170402_refcal_amp.png|Amp]] || Ant 12 was not tracking. Delay center change at 11:38:17.
|-
| 2017/04/04  || 09:32:50 || 3c273 || 04:43~10:57 || All || 08:46~10:19 || 5, 7, 9, 11, 13, 15, 17, 19, 21, 23 || [[:File:20170404_refcal_pha.png|Phase]] || [[:File:20170404_refcal_amp.png|Amp]] || Ant 12 was not tracking. Delay center change at 16:43:51.
|-
| 2017/04/05  ||06:44:10 || 3c273 || 04:17~10:30 || All || 04:27~09:34 || 5, 7, 9, 11, 13, 15, 17, 19, 21, 23 || [[:File:20170405_refcal_pha.png|Phase]] || [[:File:20170405_refcal_amp.png|Amp]] || Ant 12 was not tracking.
|-
| 2017/04/06  || 09:54:50|| 3c273 || 04:14~10:49 || All || 05:39~08:11 || 5, 7, 9, 11, 13, 15, 17, 19, 21, 23 || [[:File:20170406_refcal_pha.png|Phase]] || [[:File:20170406_refcal_amp.png|Amp]] || Ant 12 was not tracking.
|-
| 2017/04/07  || 08:30:51|| 3c273 || 04:00~10:30 || 1, 3, 5, 7, 9, 11, 13, 15, 17 || 07:53~09:10 || 5~34 || [[:File:20170407_refcal_pha.png|Phase]] || [[:File:20170407_refcal_amp.png|Amp]] || Ant 12 was not tracking.
|-
| 2017/04/08  || 06:14:51 || 3c273 || 05:00~10:30 || 1, 3, 5, 7, 9, 11, 13 || 05:00~07:00 || 5~34 || [[:File:20170408_refcal_pha.png|Phase]] || [[:File:20170408_refcal_amp.png|Amp]] || Ant 12 was not tracking. Ant 13 not working. Delay center change at 2017-04-08T03:12:26. Another refcal record is added at 2017-04-09T06:14:45.
|-
| 2017/04/10  || 06:57:10 || 3c273 || 04:30~10:30 || 2, 6, 8 ||  || 5~34 || [[:File:20170410_refcal_pha.png|Phase]] || [[:File:20170410_refcal_amp.png|Amp]] || Ant 12 was not tracking. Delay center change at 2017-04-08 13:20:39. 
|-
| 2017/04/16  || 06:36:48 || 3c273 || 03:00~11:00 ||  || 05:55~07:20 || 5~34 || [[:File:20170416_refcal_pha.png|Phase]] || [[:File:20170416_refcal_amp.png|Amp]] || Ant 12 was not tracking. 
|-
| 2017/04/17  || 06:41:51 || 3c273 || 03:29~10:07 || All || 05:55~07:29 || 5~34 || [[:File:20170417_refcal_pha.png|Phase]] || [[:File:20170417_refcal_amp.png|Amp]] || Ant 12 was not tracking. Delay center change at 19:12:58.
|-
| 2017/04/26  ||  08:06:52 || 3c273 || 02:50~09:28 || All || 07:21~08:57 || 5~34 || [[:File:20170426_refcal_pha.png|Phase]] || [[:File:20170426_refcal_amp.png|Amp]] || Ant 12 was not tracking. Delay center change at 04:37:28 and 11:51:26.
|-
| 2017/06/14 || 02:33:19 ||  1229+020  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170614_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170614_refcal_amp.png Amp] || No calibration for: Ant 8 Ant 12
|-
| 2017/06/14 || 13:01:19 ||  2253+161  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170614_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170614_refcal_amp.png Amp] || No calibration for: Ant 12
|-
| 2017/06/15 || 13:05:48 ||  ||   || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170615_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170615_refcal_amp.png Amp] || No calibration for: Ant 12
|-
| 2017/06/16 || 13:05:48 ||  ||  || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170616_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170616_refcal_amp.png Amp] || No calibration for: Ant 12
|-
| 2017/06/20 || 13:06:26 ||  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170620_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170620_refcal_amp.png Amp] || No calibration for: Ant 1 Ant 12
|-
| 2017/06/21 || 13:06:38 ||  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170621_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170621_refcal_amp.png Amp] || No calibration for: Ant 12
|-
| 2017/06/22 || 13:06:53 ||  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170622_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170622_refcal_amp.png Amp] || No calibration for: Ant 12
|-
| 2017/06/23 || 13:07:09 ||  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170623_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170623_refcal_amp.png Amp] || No calibration for: Ant 12
|-
| 2017/06/24 || 02:25:49 ||  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170624_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170624_refcal_amp.png Amp] || No calibration for: Ant 12
|-
| 2017/06/24 || 13:07:25 ||  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170624_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170624_refcal_amp.png Amp] || No calibration for: Ant 12
|-
| 2017/06/25 || 13:07:43 ||  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170625_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170625_refcal_amp.png Amp] || No calibration for: Ant 10 Ant 12
|-
| 2017/06/26 || 13:15:47 ||  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170626_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170626_refcal_amp.png Amp] || No calibration for: Ant 10 Ant 12
|-
| 2017/06/27 || 13:02:49 ||  2253+161  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170627_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170627_refcal_amp.png Amp] || No calibration for: Ant 12
|-
| 2017/06/30 || 02:47:50 ||  1229+020  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170630_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170630_refcal_amp.png Amp] || No calibration for: Ant 12
|-
| 2017/07/01 || 02:44:49 ||  1229+020  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170701_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170701_refcal_amp.png Amp] || No calibration for: Ant 12
|-
| 2017/07/02 || 03:01:20 ||  1229+020  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170702_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170702_refcal_amp.png Amp] ||
|-
| 2017/07/03 || 02:47:49 ||  1229+020  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170703_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170703_refcal_amp.png Amp] || No calibration for: Ant 10
|-
| 2017/07/04 || 02:47:49 ||  1229+020  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170704_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170704_refcal_amp.png Amp] ||
|-
| 2017/07/05 || 02:55:19 ||  1229+020  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170705_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170705_refcal_amp.png Amp] ||
|-
| 2017/07/06 || 02:47:49 ||  1229+020  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170706_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170706_refcal_amp.png Amp] || 
|-
| 2017/07/07 || 13:12:49 ||  2253+161  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170707_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170707_refcal_amp.png Amp] ||
|-
| 2017/07/09 || 02:35:49 ||  1229+020  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170709_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170709_refcal_amp.png Amp] ||
|-
| 2017/07/09 || 13:16:49 ||  2253+161  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170709_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170709_refcal_amp.png Amp] ||
|-
| 2017/07/10 || 13:14:17 ||  0319+415  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170710_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170710_refcal_amp.png Amp] ||
|-
| 2017/07/11 || 13:14:18 ||  0319+415  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170711_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170711_refcal_amp.png Amp] ||
|-
| 2017/07/12 || 13:14:24 ||  0319+415  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170712_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170712_refcal_amp.png Amp] ||
|-
| 2017/07/13 || 13:14:59 ||  0319+415  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170713_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170713_refcal_amp.png Amp] ||
|-
| 2017/07/14 || 13:15:35 ||  0319+415  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170714_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170714_refcal_amp.png Amp] ||
|-
| 2017/07/16 || 02:48:33 ||  1229+020  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170716_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170716_refcal_amp.png Amp] ||
|-
| 2017/07/16 || 13:16:48 ||  0319+415  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170716_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170716_refcal_amp.png Amp] ||
|-
| 2017/07/17 || 13:14:54 ||  0319+415  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170717_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170717_refcal_amp.png Amp] ||
|-
| 2017/07/18 || 13:15:31 ||  0319+415  || || 0 ||  || 5~34 || [http://ovsa.njit.edu/refcal/20170718_refcal_pha.png Phase] || [http://ovsa.njit.edu/refcal/20170718_refcal_amp.png Amp] ||
|-