NRAO Home > CASA > CASA Toolkit Reference Manual
imager.apparentsens - Function

2.3.1 Calculate rms sensitivity directly from weights


Description

This function calculates the point source sensitivity for the data selected by im.selectvis(...), and according to the imaging weighting parameters specified in im.weight(...) and im.defineimage(...). The calculation is performed solely using the weight information stored in the MS WEIGHT column (WEIGHT_SPECTRUM tbd), and as adjusted by the net imaging weighting function (natural, uniform, robust, taper, etc.). Therefore, it is assumed that the MS WEIGHTs have been properly initialized and calibrated along with the visibility data. As long as the WEIGHTs are in the inverse square units of the visibilities (i.e., inverse variance weights), the calculation should yield the real theoretical imaging sensitivity for data at any stage of the calibration (though data at early and intermediate stages of calibration may not be sufficiently coherent for imaging at high–or even modest–fidelity).

Two values are reported in the logger and returned (see example below). First, the apparent sensitivity (in the units implied by the WEIGHTs’ units), for the specified imaging weighting scheme. Second, a unitless factor describing the ratio of the apparent sensitivity to that obtained with pure ’natural’ weighting (the nominal peak sensitivity). When ’natural’ weighting is selected, this ratio factor will be 1.0; all other weighting choices will yield an apparent sensitivity ratio greater than 1.0.

Currently, this function reports only the continuum sensitivity for the selected data, and in particular, for the aggregate bandwidth indicated by the spectral window selection. The calculation further assumes that the visibility samples are each entirely independent (i.e., no redundant samples such as would occur for overlapping spectral windows).

A future version of this function will support reporting a sensitivity spectrum for the spectral line case (including support for WEIGHT_SPECTRUM). For now, spectral line sensitivity may be reasonably estimated by dividing the reported sensitivity by the square root of the fractional bandwidth of a single image channel, or by selecting a bandwidth matching the width of a single image channel.

Arguments





Outputs

pointsource

Calculated apparent point source sensitivity (in units implied by the MS weights)

allowed:

double

Default:

relative

Ratio of apparent sensitivity relative to natural weighting

allowed:

double

Default:

Inputs

async

Run asynchronously in the background

allowed:

bool

Default:

false

Returns
bool

Example

 
 
# open and set up selection and image plane parameters  
im.open(’mydata.ms’)  
im.selectvis(field=’2’,spw=’0’)  
im.defineimage(mode=’mfs’,spw=0,stokes=’I’,cellx=’15arcsec’,celly=’15arcsec’,nx=256,ny=256)  
 
# report natural weighting sensitivity  
im.weight(type=’natural’)  
nat=im.apparentsens();  
print ’Natural Sensitivity =’, nat[1];  
print ’Relative to Natural Weighting = ’, nat[2];  
 
# switch to uniform weighting  
im.weight(type=’uniform’)  
uni=im.apparentsens();  
print ’Uniform Sensitivity =’, uni[1];  
print ’Relative to Natural Weighting = ’, uni[2];  
 
# switch to briggs weighting  
im.weight(type=’briggs’,robust=0.0)  
rob=im.apparentsens();  
print ’Briggs Sensitivity =’, rob[1];  
print ’Relative to Natural Weighting = ’, rob[2];  
 
im.close()  
 

__________________________________________________________________


More information about CASA may be found at the CASA web page

Copyright © 2016 Associated Universities Inc., Washington, D.C.

This code is available under the terms of the GNU General Public Lincense


Home | Contact Us | Directories | Site Map | Help | Privacy Policy | Search