Transient reduction pipeline

From CASA Guides
Jump to navigationJump to search


Transient Reduction Pipeline

2010 May 11

There are a class of observations for which only relatively simple data reduction steps are needed. One such example is that of transient observations, which are typically conducted in continuum mode and for which one is merely trying to determine a flux density of an essentially unresolved or only partially resolved source. This guide describes the steps in such a pipeline.

Warning: This guide was written at the time when the EVLA was still in its commissioning phase. As the instrument matures, specific steps taken here may need to be adjusted. Caveat emptor.

In order to process these data in a semi-automatic fashion, certain assumptions are made

  • Data stored as an ALMA Science Data Model (ASDM) file on disk.
  • Sources, listed in order of appearance, are structured as Option A or Option B
    • Option A: phase calibrator, target source, amplitude calibrator
    • Option B: amplitude calibrator, phase calibrator, target source

User Input

1. Obtain the name of the ASDM file, use it to construct the Measurement Set name.

# In CASA
ASDM_name=raw_input('ASDM file name:')
msname=ASDM_name+'.ms'

2. Determine the order in which sources appear in the observation.

# In CASA
obs_structure=raw_input('Structure of observations [A|B]:')
if obs_structure = 'A':
    sources = {'acal':2, 'pcal':0, 'target':1}
elif obs_structure = 'B':
    sources = {'acal':0, 'pcal':1, 'target':2}


Processing Steps

1. Read the data from the ASDM file converting to a Measurement Set with importevla. Apply basic flagging operations (zeros, shadowing) here.

importevla(asdm=ASDM_name, vis=msname,
           flagzero=True, flagpol=True, shadow=True)

2. Extract the amplitude calibrator name from the data, for later calibration.

vishead(vis=msname, mode='get',
        hdkey='field'; hdindex=sources['acal'])
acal_name=vishead()

3. Flag first (dummy) scan.

flagdata(vis=msname,
         mode='manualflag',
         selectdata=T; scan='1')

4. Flag first 10 seconds of each scan.

flagdata(vis=msname,
         mode='quack',quackinterval=10,quackmode='beg',
         selectdata=F)


5. Compress data in time.

   (You've checked that time average smearing isn't an issue ...?)
outmsname=ASDM_name+'_split.ms'
split(vis=msname,outputvis=outmsname,
      datacolumn='data',
      timebin='10s')