casa::refim::SynthesisImager Class Reference

#include <SynthesisImager.h>

Public Member Functions
	SynthesisImager ()
	Default constructor.
virtual	~SynthesisImager ()
virtual Bool	selectData (const SynthesisParamsSelect &selpars)
	Copy constructor and assignment operator Imager(const Imager&); Imager& operator=(const Imager&);.
virtual Bool	selectData (const String &msname, const String &spw="", const String &freqBeg="", const String &freqEnd="", const MFrequency::Types freqFrame=MFrequency::LSRK, const String &field="", const String &antenna="", const String &timestr="", const String &scan="", const String &obs="", const String &state="", const String &uvdist="", const String &taql="", const Bool usescratch=False, const Bool readonly=False, const Bool incrementModel=False)
	make all pure-inputs const
virtual Bool	defineImage (SynthesisParamsImage &impars, const SynthesisParamsGrid &gridpars)
virtual Bool	defineImage (const String &imagename, const Int nx, const Int ny, const Quantity &cellx, const Quantity &celly, const String &stokes, const MDirection &phaseCenter, const Int nchan, const Quantity &freqStart, const Quantity &freqStep, const Vector< Quantity > &restFreq, const Int facets=1, const String ftmachine="GridFT", const Int nTaylorTerms=1, const Quantity &refFreq=Quantity(0,"Hz"), const Projection &projection=Projection::SIN, const Quantity &distance=Quantity(0,"m"), const MFrequency::Types &freqFrame=MFrequency::LSRK, const Bool trackSource=False, const MDirection &trackDir=MDirection(Quantity(0.0,"deg"), Quantity(90.0,"deg")), const Bool overwrite=False, const Float padding=1.0, const Bool useAutocorr=False, const bool useDoublePrec=True, const Int wprojplanes=1, const String convFunc="SF", const String startmodel="", const Bool aTermOn=True, const Bool psTermOn=True, const Bool mTermOn=False, const Bool wbAWP=True, const String cfCache="", const Bool doPointing=False, const Bool doPBCorr=True, const Bool conjBeams=True, const Float computePAStep=360.0, const Float rotatePAStep=5.0)
	When having a facetted image ...call with (facets > 1) first and once only .
virtual Bool	defineImage (CountedPtr< SIImageStore > imstor, const String &ftmachine)
	Define image via a predefine SIImageStore object.
virtual void	setComponentList (const ComponentList &cl, Bool sdgrid=False)
	Defining componentlist to use while degriding This should be called once...if multiple lists are used..they can be merged in one if sdgrid=True then image plane degridding is done.
Bool	weight (const String &type="natural", const String &rmode="norm", const Quantity &noise=Quantity(0.0,"Jy"), const Double robust=0.0, const Quantity &fieldofview=Quantity(0.0,"arcsec"), const Int npixels=0, const Bool multiField=False, const String &filtertype=String("Gaussian"), const Quantity &filterbmaj=Quantity(0.0,"deg"), const Quantity &filterbmin=Quantity(0.0,"deg"), const Quantity &filterbpa=Quantity(0.0,"deg"))
void	resetMappers ()
	the following get rid of the mappers in this object
CountedPtr< SIImageStore >	imageStore (const Int id=0)
void	executeMajorCycle (Record &controls)
	Record getMajorCycleControls();.
void	makePSF ()
	make the psf images i.e grid weight rather than data
void	predictModel ()
Protected Member Functions
Choose between different types of FTMs *void	createFTMachine (CountedPtr< refim::FTMachine > &theFT, CountedPtr< refim::FTMachine > &theIFT, const String &ftname, const uInt nTaylorTerms=1, const String mType="default", const Int facets=1, const Int wprojplane=1, const Float padding=1.0, const Bool useAutocorr=False, const Bool useDoublePrec=True, const String gridFunction=String("SF"), const Bool aTermOn=True, const Bool psTermOn=True, const Bool mTermOn=False, const Bool wbAWP=True, const String cfCache="", const Bool doPointing=False, const Bool doPBCorr=True, const Bool conjBeams=True, const Float computePAStep=360.0, const Float rotatePAStep=5.0, const String interpolation=String("linear"), const Bool freqFrameValid=True, const Int cache=1000000000, const Int tile=16)
	SIIterBot& getLoopControls();.
CountedPtr< SIImageStore >	createIMStore (String imageName, CoordinateSystem &cSys, IPosition imShape, const Bool overwrite, String mappertype="default", uInt ntaylorterms=1, Quantity distance=Quantity(0.0,"m"), uInt facets=1, Bool useweightimage=False)
	Choose between different types of ImageStore types (single term, multiterm, faceted).
CountedPtr< SIMapper >	createSIMapper (String mappertype, CountedPtr< SIImageStore > imagestore, CountedPtr< refim::FTMachine > ftmachine, CountedPtr< refim::FTMachine > iftmachine, uInt ntaylorterms=1)
	Choose between different types of Mappers (single term, multiterm, imagemosaic, faceted).
Block< CountedPtr< SIImageStore > >	createFacetImageStoreList (CountedPtr< SIImageStore > imagestore, Int facets)
void	setPsfFromOneFacet ()
Bool	toUseWeightImage (CountedPtr< refim::FTMachine > &ftm, String mappertype)
void	createVisSet (const Bool writeaccess=False)
void	createAWPFTMachine (CountedPtr< FTMachine > &theFT, CountedPtr< FTMachine > &theIFT, const String &ftmName, const Int facets, const Int wprojPlane, const Float padding, const Bool useAutocorr, const Bool useDoublePrec, const String gridFunction, const Bool aTermOn, const Bool psTermOn, const Bool mTermOn, const Bool wbAWP, const String cfCache, const Bool doPointing, const Bool doPBCorr, const Bool conjBeams, const Float computePAStep, const Float rotatePAStep, const Int cache, const Int tile)
refim::ATerm *	createTelescopeATerm (const MeasurementSet &ms, const Bool &isATermOn)
void	runMajorCycle (const Bool dopsf=False, const Bool savemodel=False)
void	appendToMapperList (String imagename, CoordinateSystem &csys, String ftmachine, Quantity distance=Quantity(0.0,"m"), Int facets=1, const Bool overwrite=False)
	It associated the ftmachine with a given field.
*void	appendToMapperList (String imagename, CoordinateSystem &csys, IPosition imshape, CountedPtr< refim::FTMachine > &ftm, CountedPtr< refim::FTMachine > &iftm, Quantity distance=Quantity(0.0,"m"), Int facets=1, const Bool overwrite=False, String mappertype=String("default"), uInt ntaylorterms=1)
Protected Attributes
*SIMapperCollection	itsMappers
CountedPtr< SIImageStore >	unFacettedImStore_p
	if facetting this storage will keep the unsliced version
CoordinateSystem	itsMaxCoordSys
IPosition	itsMaxShape
Data Selection Image Definition Imaging Gridding Other Options Block< const MeasurementSet >	mss_p
vi::FrequencySelections	fselections_p
CountedPtr < vi::VisibilityIterator2 >	vi_p
Bool	writeAccess_p
MPosition	mLocation_p
MDirection	phaseCenter_p
Int	facetsStore_p
VisImagingWeight	imwgt_p
Bool	imageDefined_p
Bool	useScratch_p
Bool	readOnly_p
refim::FTMachine::Type	datacol_p
	Bool freqFrameValid_p;.

Detailed Description

Definition at line 56 of file SynthesisImager.h.

Constructor & Destructor Documentation

casa::refim::SynthesisImager::SynthesisImager ( )

Default constructor.

virtual casa::refim::SynthesisImager::~SynthesisImager ( ) [virtual]

Member Function Documentation

* void casa::refim::SynthesisImager::appendToMapperList	(	String	imagename,
		CoordinateSystem &	csys,
		IPosition	imshape,
		CountedPtr< refim::FTMachine > &	ftm,
		CountedPtr< refim::FTMachine > &	iftm,
		Quantity	distance = `Quantity(0.0,"m")`,
		Int	facets = `1`,
		const Bool	overwrite = `False`,
		String	mappertype = `String("default")`,
		uInt	ntaylorterms = `1`
	)			`[protected]`

void casa::refim::SynthesisImager::appendToMapperList	(	String	imagename,
		CoordinateSystem &	csys,
		String	ftmachine,
		Quantity	distance = `Quantity(0.0,"m")`,
		Int	facets = `1`,
		const Bool	overwrite = `False`
	)			`[protected]`

It associated the ftmachine with a given field.

void casa::refim::SynthesisImager::createAWPFTMachine	(	CountedPtr< FTMachine > &	theFT,
		CountedPtr< FTMachine > &	theIFT,
		const String &	ftmName,
		const Int	facets,
		const Int	wprojPlane,
		const Float	padding,
		const Bool	useAutocorr,
		const Bool	useDoublePrec,
		const String	gridFunction,
		const Bool	aTermOn,
		const Bool	psTermOn,
		const Bool	mTermOn,
		const Bool	wbAWP,
		const String	cfCache,
		const Bool	doPointing,
		const Bool	doPBCorr,
		const Bool	conjBeams,
		const Float	computePAStep,
		const Float	rotatePAStep,
		const Int	cache,
		const Int	tile
	)			`[protected]`

Parameters:

	wprojPlane	----------------------------
	aTermOn	---------------------------

Block<CountedPtr<SIImageStore> > casa::refim::SynthesisImager::createFacetImageStoreList	(	CountedPtr< SIImageStore >	imagestore,
		Int	facets
	)			`[protected]`

Choose between different types of FTMs* void casa::refim::SynthesisImager::createFTMachine	(	CountedPtr< refim::FTMachine > &	theFT,
		CountedPtr< refim::FTMachine > &	theIFT,
		const String &	ftname,
		const uInt	nTaylorTerms = `1`,
		const String	mType = `"default"`,
		const Int	facets = `1`,
		const Int	wprojplane = `1`,
		const Float	padding = `1.0`,
		const Bool	useAutocorr = `False`,
		const Bool	useDoublePrec = `True`,
		const String	gridFunction = `String("SF")`,
		const Bool	aTermOn = `True`,
		const Bool	psTermOn = `True`,
		const Bool	mTermOn = `False`,
		const Bool	wbAWP = `True`,
		const String	cfCache = `""`,
		const Bool	doPointing = `False`,
		const Bool	doPBCorr = `True`,
		const Bool	conjBeams = `True`,
		const Float	computePAStep = `360.0`,
		const Float	rotatePAStep = `5.0`,
		const String	interpolation = `String("linear")`,
		const Bool	freqFrameValid = `True`,
		const Int	cache = `1000000000`,
		const Int	tile = `16`
	)			`[protected]`

SIIterBot& getLoopControls();.

Parameters:

	wprojplane	------------------------------
	aTermOn	------------------------------

CountedPtr<SIImageStore> casa::refim::SynthesisImager::createIMStore	(	String	imageName,
		CoordinateSystem &	cSys,
		IPosition	imShape,
		const Bool	overwrite,
		String	mappertype = `"default"`,
		uInt	ntaylorterms = `1`,
		Quantity	distance = `Quantity(0.0,"m")`,
		uInt	facets = `1`,
		Bool	useweightimage = `False`
	)			`[protected]`

Choose between different types of ImageStore types (single term, multiterm, faceted).

CountedPtr<SIMapper> casa::refim::SynthesisImager::createSIMapper	(	String	mappertype,
		CountedPtr< SIImageStore >	imagestore,
		CountedPtr< refim::FTMachine >	ftmachine,
		CountedPtr< refim::FTMachine >	iftmachine,
		uInt	ntaylorterms = `1`
	)			`[protected]`

Choose between different types of Mappers (single term, multiterm, imagemosaic, faceted).

Parameters:

ftmachine

/ make this inside !!!!!

refim::ATerm* casa::refim::SynthesisImager::createTelescopeATerm	(	const MeasurementSet &	ms,
		const Bool &	isATermOn
	)			`[protected]`

void casa::refim::SynthesisImager::createVisSet ( const Bool writeaccess = False ) [protected]

virtual Bool casa::refim::SynthesisImager::defineImage	(	CountedPtr< SIImageStore >	imstor,
		const String &	ftmachine
	)			`[virtual]`

Define image via a predefine SIImageStore object.

virtual Bool casa::refim::SynthesisImager::defineImage	(	const String &	imagename,
		const Int	nx,
		const Int	ny,
		const Quantity &	cellx,
		const Quantity &	celly,
		const String &	stokes,
		const MDirection &	phaseCenter,
		const Int	nchan,
		const Quantity &	freqStart,
		const Quantity &	freqStep,
		const Vector< Quantity > &	restFreq,
		const Int	facets = `1`,
		const String	ftmachine = `"GridFT"`,
		const Int	nTaylorTerms = `1`,
		const Quantity &	refFreq = `Quantity(0,"Hz")`,
		const Projection &	projection = `Projection::SIN`,
		const Quantity &	distance = `Quantity(0,"m")`,
		const MFrequency::Types &	freqFrame = `MFrequency::LSRK`,
		const Bool	trackSource = `False`,
		const MDirection &	trackDir = `MDirection(Quantity(0.0,"deg"), Quantity(90.0,"deg"))`,
		const Bool	overwrite = `False`,
		const Float	padding = `1.0`,
		const Bool	useAutocorr = `False`,
		const bool	useDoublePrec = `True`,
		const Int	wprojplanes = `1`,
		const String	convFunc = `"SF"`,
		const String	startmodel = `""`,
		const Bool	aTermOn = `True`,
		const Bool	psTermOn = `True`,
		const Bool	mTermOn = `False`,
		const Bool	wbAWP = `True`,
		const String	cfCache = `""`,
		const Bool	doPointing = `False`,
		const Bool	doPBCorr = `True`,
		const Bool	conjBeams = `True`,
		const Float	computePAStep = `360.0`,
		const Float	rotatePAStep = `5.0`
	)			`[virtual]`

When having a facetted image ...call with (facets > 1) first and once only .

\. Easier to keep track of the imstores that way CAREFUL: make sure you donot overwrite if you want to predict the model or subtract it to make residual

Parameters:

aTermOn

The extra params for WB-AWP

virtual Bool casa::refim::SynthesisImager::defineImage	(	SynthesisParamsImage &	impars,
		const SynthesisParamsGrid &	gridpars
	)			`[virtual]`

void casa::refim::SynthesisImager::executeMajorCycle ( Record & controls )

Record getMajorCycleControls();.

CountedPtr<SIImageStore> casa::refim::SynthesisImager::imageStore ( const Int id = 0 )

void casa::refim::SynthesisImager::makePSF ( )

make the psf images i.e grid weight rather than data

void casa::refim::SynthesisImager::predictModel ( )

void casa::refim::SynthesisImager::resetMappers ( )

the following get rid of the mappers in this object

void casa::refim::SynthesisImager::runMajorCycle	(	const Bool	dopsf = `False`,
		const Bool	savemodel = `False`
	)			`[protected]`

virtual Bool casa::refim::SynthesisImager::selectData	(	const String &	msname,
		const String &	spw = `"*"`,
		const String &	freqBeg = `""`,
		const String &	freqEnd = `""`,
		const MFrequency::Types	freqFrame = `MFrequency::LSRK`,
		const String &	field = `"*"`,
		const String &	antenna = `""`,
		const String &	timestr = `""`,
		const String &	scan = `""`,
		const String &	obs = `""`,
		const String &	state = `""`,
		const String &	uvdist = `""`,
		const String &	taql = `""`,
		const Bool	usescratch = `False`,
		const Bool	readonly = `False`,
		const Bool	incrementModel = `False`
	)			`[virtual]`

make all pure-inputs const

virtual Bool casa::refim::SynthesisImager::selectData ( const SynthesisParamsSelect & selpars ) [virtual]

Copy constructor and assignment operator Imager(const Imager&); Imager& operator=(const Imager&);.

virtual void casa::refim::SynthesisImager::setComponentList	(	const ComponentList &	cl,
		Bool	sdgrid = `False`
	)			`[virtual]`

Defining componentlist to use while degriding This should be called once...if multiple lists are used..they can be merged in one if sdgrid=True then image plane degridding is done.

void casa::refim::SynthesisImager::setPsfFromOneFacet ( ) [protected]

Bool casa::refim::SynthesisImager::toUseWeightImage	(	CountedPtr< refim::FTMachine > &	ftm,
		String	mappertype
	)			`[protected]`

Bool casa::refim::SynthesisImager::weight	(	const String &	type = `"natural"`,
		const String &	rmode = `"norm"`,
		const Quantity &	noise = `Quantity(0.0,"Jy")`,
		const Double	robust = `0.0`,
		const Quantity &	fieldofview = `Quantity(0.0,"arcsec")`,
		const Int	npixels = `0`,
		const Bool	multiField = `False`,
		const String &	filtertype = `String("Gaussian")`,
		const Quantity &	filterbmaj = `Quantity(0.0,"deg")`,
		const Quantity &	filterbmin = `Quantity(0.0,"deg")`,
		const Quantity &	filterbpa = `Quantity(0.0,"deg")`
	)