Daily Gain Calibration: Difference between revisions
Line 29: | Line 29: | ||
filename = '/data1/eovsa/fits/UDB/2017/UDB20170821155<tab>' | filename = '/data1/eovsa/fits/UDB/2017/UDB20170821155<tab>' | ||
pc.udb_corr(filename) | pc.udb_corr(filename) | ||
</pre> | |||
== 20:00 UT == | == 20:00 UT == | ||
Run the phasecal analysis, answering y to questions about plots and accepting: | Run the phasecal analysis, answering y to questions about plots and accepting: |
Revision as of 15:03, 20 August 2017
Similar as reference complex gain calibration, daily gain calibration refers to tuning the "system" amplitude and phase as a function of IF band at a finer level. Typically, daily gain calibration would be carried out three times between two reference gain calibrations every day. We have a routine to retrieve and analyze the daily gain calibration data, named phacal_anal.py. Start from a timestamp (Note: The time of day is important! Make sure to set it to somewhere in the middle of the day, e.g. 15:00 UT)
import phacal_anal as pa from util import Time t = Time('2017-07-03 15:00') pa.phacal_anal(t)
phacal_anal will find all daily gain calibrations occurred during that natural day (CA local time).
Commands for calibrating the eclipse [temporary]
15:00 UT
First run the refcal analysis for the day (can be done by 15:00 UT):
from util import Time import refcal_anal as ra import cal_header as ch out = ra.rd_refcal(Time(['2017-08-21 11:00','2017-08-21 14:00')) out_corr = ra.unrot_refcal(out) refcal = ra.refcal_anal(out) ch.refcal2sql(refcal)
19:20 UT
Next, as soon as the solar (eclipse) UDB file appears (around 19:20 UT) run the pipeline processing. Change directory to the one where the final eclipse files will be, /data1/dgary/solar/eclipse. Note: the <tab> will complete the filename, but delete the trailing /.:
from util import Time import pipeline_cal as pc filename = '/data1/eovsa/fits/UDB/2017/UDB20170821155<tab>' pc.udb_corr(filename)
20:00 UT
Run the phasecal analysis, answering y to questions about plots and accepting:
from util import Time import phacal_anal as pa pa.phacal_anal(Time('2017-08-21 19:00'))
Immediately after that, exit ipython, cd to /data1/dgary/solar/eclipse and get into CASA, then do
default importeovsa idbfiles = 'UDB20170821<tab>' go from suncasa.tasks import task_calibeovsa as tc tc.calibeovsa('UDB2017081915<tab>',doimage=True,stokes='I')
This will generate the first eclipse map, but it will be a single image of the entire eclipse! There may be some visible indication of the Moon in the image, but perhaps not. To make better images, it will be necessary to run clean on restricted time ranges.
default clean <set parameters> go