Create a Component List for Selfcal: Difference between revisions

From CASA Guides
Jump to navigationJump to search
(Created page with " If you want to selfcal on a point source (or a list of point sources) with known position, you may create your own model component list. To do so, you can use the tools like thi...")
 
No edit summary
 
(6 intermediate revisions by the same user not shown)
Line 1: Line 1:


If you want to selfcal on a point source (or a list of point sources) with known position, you may create your own model component list. To do so, you can use the tools like this:
If you want to selfcal on one or more point sources with known position (or a more complicated shapes such Gaussians), you may create your own model component list. To do so, you can use the tools like this:


<source lang="python">


# for a point source with no spectral index
cl.addcomponent(flux=0.39, fluxunit='Jy',shape='point', dir='J2000 19h33m09s 15d01m20s')
# for a Gaussian with a spectral index
cl.addcomponent(flux=1.25, fluxunit='mJy', polarization='Stokes',
cl.addcomponent(flux=1.25, fluxunit='mJy', polarization='Stokes',
dir='J2000 19h30m00 15d00m00', shape='gaussian', majoraxis='10arcsec',
dir='J2000 19h30m00s 15d00m00s', shape='gaussian', majoraxis='10arcsec',
minoraxis='6arcsec', positionangle='0deg', freq='1.25GHz',
minoraxis='6arcsec', positionangle='0deg', freq='1.25GHz',
spectrumtype='spectral index', index=-0.8)
spectrumtype='spectral index', index=-0.8)
###you can add more components if you wish by calling addcomponent
###you can add more components if you wish by calling addcomponent repeatedly with different params
repeatedly with different params


##save it to disk
##save it to disk
cl.rename('my_1_component.cl')
cl.rename('my_component.cl')
cl.close()
cl.close()


 
## write the model into the MODEL_DATA column of your data ('myms')
ft(vis='myms', complist='my_1_component.cl')
ft(vis='myms', complist='my_component.cl')
</source>


Then proceed with the gaincal as in normal selfcal routines.
Then proceed with the gaincal as in normal selfcal routines.

Latest revision as of 15:12, 29 February 2012

If you want to selfcal on one or more point sources with known position (or a more complicated shapes such Gaussians), you may create your own model component list. To do so, you can use the tools like this:

# for a point source with no spectral index
cl.addcomponent(flux=0.39, fluxunit='Jy',shape='point', dir='J2000 19h33m09s 15d01m20s')

# for a Gaussian with a spectral index
cl.addcomponent(flux=1.25, fluxunit='mJy', polarization='Stokes',
dir='J2000 19h30m00s 15d00m00s', shape='gaussian', majoraxis='10arcsec',
minoraxis='6arcsec', positionangle='0deg', freq='1.25GHz',
spectrumtype='spectral index', index=-0.8)
###you can add more components if you wish by calling addcomponent repeatedly with different params

##save it to disk
cl.rename('my_component.cl')
cl.close()

## write the model into the MODEL_DATA column of your data ('myms')
ft(vis='myms', complist='my_component.cl')

Then proceed with the gaincal as in normal selfcal routines.