https://casaguides.nrao.edu/index.php?title=ACA_Simulation_(CASA_4.1)&feed=atom&action=historyACA Simulation (CASA 4.1) - Revision history2024-03-19T11:59:41ZRevision history for this page on the wikiMediaWiki 1.38.6https://casaguides.nrao.edu/index.php?title=ACA_Simulation_(CASA_4.1)&diff=14816&oldid=prevRindebet: /* Deconvolve the visibilities back into an image */2013-07-01T20:50:55Z<p><span dir="auto"><span class="autocomment">Deconvolve the visibilities back into an image</span></span></p>
<table style="background-color: #fff; color: #202122;" data-mw="interface">
<col class="diff-marker" />
<col class="diff-content" />
<col class="diff-marker" />
<col class="diff-content" />
<tr class="diff-title" lang="en">
<td colspan="2" style="background-color: #fff; color: #202122; text-align: center;">← Older revision</td>
<td colspan="2" style="background-color: #fff; color: #202122; text-align: center;">Revision as of 16:50, 1 July 2013</td>
</tr><tr><td colspan="2" class="diff-lineno" id="mw-diff-left-l150">Line 150:</td>
<td colspan="2" class="diff-lineno">Line 150:</td></tr>
<tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td></tr>
<tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div><font color="red"></div></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div><font color="red"></div></td></tr>
<tr><td class="diff-marker" data-marker="−"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div>If you decide to concatenate the two measurement sets and image all visibilities simultaneously, it is critical that the relative weights be set between the two different interferometric arrays. Simulated data has weights=1, since the thermal noise is generated uniformly per baseline. However, in reality the 7m baselines have lower sensitivity than the 12m baselines, and their weights must be decreased by that sensitivity ratio. <tt>simalma</tt> uses the <tt>visweightscale</tt> parameter of <tt>concat</tt> to apply that lower weight of (7./12)**2 to the 7m visibilities (See [Simalma] guide)</div></td><td class="diff-marker" data-marker="+"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div>If you decide to concatenate the two measurement sets and image all visibilities simultaneously, it is critical that the relative weights be set between the two different interferometric arrays. Simulated data has weights=1, since the thermal noise is generated uniformly per baseline. However, in reality the 7m baselines have lower sensitivity than the 12m baselines, and their weights must be decreased by that sensitivity ratio. <tt>simalma</tt> uses the <tt>visweightscale</tt> parameter of <tt>concat</tt> to apply that lower weight of (7./12)**2 to the 7m visibilities (See <ins style="font-weight: bold; text-decoration: none;">[</ins>[Simalma<ins style="font-weight: bold; text-decoration: none;">]</ins>] guide)</div></td></tr>
<tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div></font></div></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div></font></div></td></tr>
<tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td></tr>
</table>Rindebethttps://casaguides.nrao.edu/index.php?title=ACA_Simulation_(CASA_4.1)&diff=14815&oldid=prevRindebet: /* Deconvolve the visibilities back into an image */2013-07-01T20:50:25Z<p><span dir="auto"><span class="autocomment">Deconvolve the visibilities back into an image</span></span></p>
<table style="background-color: #fff; color: #202122;" data-mw="interface">
<col class="diff-marker" />
<col class="diff-content" />
<col class="diff-marker" />
<col class="diff-content" />
<tr class="diff-title" lang="en">
<td colspan="2" style="background-color: #fff; color: #202122; text-align: center;">← Older revision</td>
<td colspan="2" style="background-color: #fff; color: #202122; text-align: center;">Revision as of 16:50, 1 July 2013</td>
</tr><tr><td colspan="2" class="diff-lineno" id="mw-diff-left-l150">Line 150:</td>
<td colspan="2" class="diff-lineno">Line 150:</td></tr>
<tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td></tr>
<tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div><font color="red"></div></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div><font color="red"></div></td></tr>
<tr><td class="diff-marker" data-marker="−"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div>If you decide to concatenate the two measurement sets and image all visibilities simultaneously, it is critical that the relative weights be set between the two different interferometric arrays. Simulated data has weights=1, since the thermal noise is generated uniformly per baseline. However, in reality the 7m baselines have lower sensitivity than the 12m baselines, and their weights must be decreased by that sensitivity ratio. <tt>simalma<tt> uses the <tt>visweightscale</tt> parameter of <tt>concat<tt> to apply that lower weight of (7./12)**2 to the 7m visibilities (See [Simalma] guide)</div></td><td class="diff-marker" data-marker="+"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div>If you decide to concatenate the two measurement sets and image all visibilities simultaneously, it is critical that the relative weights be set between the two different interferometric arrays. Simulated data has weights=1, since the thermal noise is generated uniformly per baseline. However, in reality the 7m baselines have lower sensitivity than the 12m baselines, and their weights must be decreased by that sensitivity ratio. <tt>simalma<<ins style="font-weight: bold; text-decoration: none;">/</ins>tt> uses the <tt>visweightscale</tt> parameter of <tt>concat<<ins style="font-weight: bold; text-decoration: none;">/</ins>tt> to apply that lower weight of (7./12)**2 to the 7m visibilities (See [Simalma] guide)</div></td></tr>
<tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div></font></div></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div></font></div></td></tr>
<tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td></tr>
</table>Rindebethttps://casaguides.nrao.edu/index.php?title=ACA_Simulation_(CASA_4.1)&diff=14814&oldid=prevRindebet: /* Deconvolve the visibilities back into an image */2013-07-01T20:50:05Z<p><span dir="auto"><span class="autocomment">Deconvolve the visibilities back into an image</span></span></p>
<table style="background-color: #fff; color: #202122;" data-mw="interface">
<col class="diff-marker" />
<col class="diff-content" />
<col class="diff-marker" />
<col class="diff-content" />
<tr class="diff-title" lang="en">
<td colspan="2" style="background-color: #fff; color: #202122; text-align: center;">← Older revision</td>
<td colspan="2" style="background-color: #fff; color: #202122; text-align: center;">Revision as of 16:50, 1 July 2013</td>
</tr><tr><td colspan="2" class="diff-lineno" id="mw-diff-left-l147">Line 147:</td>
<td colspan="2" class="diff-lineno">Line 147:</td></tr>
<tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td></tr>
<tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>If given a total power and interferometric measurement set, '''simanalyze''' will automatically create the total power image, </div></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>If given a total power and interferometric measurement set, '''simanalyze''' will automatically create the total power image, </div></td></tr>
<tr><td class="diff-marker" data-marker="−"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div>then use it as a model and deconvolve the interferometric image. <del style="font-weight: bold; text-decoration: none;"> It </del>is <del style="font-weight: bold; text-decoration: none;">not recommended </del>to <del style="font-weight: bold; text-decoration: none;">do both interferometric images simultaneously</del>.</div></td><td class="diff-marker" data-marker="+"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div>then use it as a model and deconvolve the interferometric image.</div></td></tr>
<tr><td colspan="2" class="diff-side-deleted"></td><td class="diff-marker" data-marker="+"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div> </div></td></tr>
<tr><td colspan="2" class="diff-side-deleted"></td><td class="diff-marker" data-marker="+"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div><ins style="font-weight: bold; text-decoration: none;"><font color="red"></ins></div></td></tr>
<tr><td colspan="2" class="diff-side-deleted"></td><td class="diff-marker" data-marker="+"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div><ins style="font-weight: bold; text-decoration: none;">If you decide to concatenate the two measurement sets and image all visibilities simultaneously, it is critical that the relative weights be set between the two different interferometric arrays. Simulated data has weights=1, since the thermal noise </ins>is <ins style="font-weight: bold; text-decoration: none;">generated uniformly per baseline. However, in reality the 7m baselines have lower sensitivity than the 12m baselines, and their weights must be decreased by that sensitivity ratio. <tt>simalma<tt> uses the <tt>visweightscale</tt> parameter of <tt>concat<tt> </ins>to <ins style="font-weight: bold; text-decoration: none;">apply that lower weight of (7</ins>.<ins style="font-weight: bold; text-decoration: none;">/12)**2 to the 7m visibilities (See [Simalma] guide)</ins></div></td></tr>
<tr><td colspan="2" class="diff-side-deleted"></td><td class="diff-marker" data-marker="+"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div><ins style="font-weight: bold; text-decoration: none;"></font></ins></div></td></tr>
<tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td></tr>
<tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>====First image total power and ACA with total power as a model====</div></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>====First image total power and ACA with total power as a model====</div></td></tr>
</table>Rindebethttps://casaguides.nrao.edu/index.php?title=ACA_Simulation_(CASA_4.1)&diff=14156&oldid=prevJbraatz at 14:38, 10 May 20132013-05-10T14:38:50Z<p></p>
<table style="background-color: #fff; color: #202122;" data-mw="interface">
<col class="diff-marker" />
<col class="diff-content" />
<col class="diff-marker" />
<col class="diff-content" />
<tr class="diff-title" lang="en">
<td colspan="2" style="background-color: #fff; color: #202122; text-align: center;">← Older revision</td>
<td colspan="2" style="background-color: #fff; color: #202122; text-align: center;">Revision as of 10:38, 10 May 2013</td>
</tr><tr><td colspan="2" class="diff-lineno" id="mw-diff-left-l1">Line 1:</td>
<td colspan="2" class="diff-lineno">Line 1:</td></tr>
<tr><td class="diff-marker" data-marker="−"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div><del style="font-weight: bold; text-decoration: none;">{{Simulations Intro 4.0}}</del></div></td><td colspan="2" class="diff-side-added"></td></tr>
<tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>[[Category: Simulations]]</div></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>[[Category: Simulations]]</div></td></tr>
<tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td></tr>
<tr><td colspan="2" class="diff-lineno" id="mw-diff-left-l184">Line 184:</td>
<td colspan="2" class="diff-lineno">Line 183:</td></tr>
<tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div></source></div></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div></source></div></td></tr>
<tr><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td><td class="diff-marker"></td><td style="background-color: #f8f9fa; color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><br/></td></tr>
<tr><td class="diff-marker" data-marker="−"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div><del style="font-weight: bold; text-decoration: none;">{{Simulations Intro 4.0}}</del></div></td><td class="diff-marker" data-marker="+"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div>{{Checked 4.<ins style="font-weight: bold; text-decoration: none;">1</ins>.0}}</div></td></tr>
<tr><td class="diff-marker" data-marker="−"></td><td style="color: #202122; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div>{{Checked 4.<del style="font-weight: bold; text-decoration: none;">0</del>.0}}</div></td><td colspan="2" class="diff-side-added"></td></tr>
</table>Jbraatzhttps://casaguides.nrao.edu/index.php?title=ACA_Simulation_(CASA_4.1)&diff=14001&oldid=prevJbraatz: Created page with "{{Simulations Intro 4.0}} Category: Simulations To create a script of the Python code on this page see Extracting scripts from these tutorials. __NOTOC__ == ALMA 12m ..."2013-04-12T21:27:22Z<p>Created page with "{{Simulations Intro 4.0}} <a href="/index.php?title=Category:Simulations" title="Category:Simulations">Category: Simulations</a> To create a script of the Python code on this page see <a href="/index.php?title=Extracting_scripts_from_these_tutorials" title="Extracting scripts from these tutorials">Extracting scripts from these tutorials</a>. __NOTOC__ == ALMA 12m ..."</p>
<p><b>New page</b></p><div>{{Simulations Intro 4.0}}<br />
[[Category: Simulations]]<br />
<br />
To create a script of the Python code on this page see [[Extracting scripts from these tutorials]].<br />
__NOTOC__<br />
== ALMA 12m + ACA + Total Power ==<br />
<br />
We can simulate ALMA observations that use the main array of 12 m antennas augmented by the ACA and Total Power antennas by generating simulated observations for each component separately, and then combining the Measurement Sets in '''simanalyze'''. This technique is general and can be used to simulate observations using multiple 12 m array configurations, as well. Total power observations can be simulated either in an independent run of simobserve, or along with an interferometric simulation. Note that if you simulate total power and an interferometric observation simultaneously with simobserve, they must have the same set of pointing centers and the same integration and total time, which is probably not realistic. (For example it is generally recommended to observe a larger area by 1/2 primary beam in total power mode to combine with a 12 m ALMA mosaic).<br />
<br />
=====Set simobserve as current task=====<br />
Reset all parameters to default, and then set the project name to ''m51c''<br />
<source lang="python"><br />
# Set simobserve to default parameters<br />
default("simobserve")<br />
# Our project name will be m51c, and all simulation products will be placed in a subdirectory m51c/<br />
project="m51c"<br />
</source><br />
<br />
=====Specify sky model image=====<br />
We'll use an Halpha image of M51 as a model of the sky, for this example. The ''curl'' command will copy the file from a URL and rename it.<br />
<br />
<source lang="python"><br />
# Model sky = Halpha image of M51 <br />
os.system('curl http://casaguides.nrao.edu/images/3/3f/M51ha.fits.txt -f -o M51ha.fits')<br />
skymodel = "M51ha.fits"<br />
</source> <br />
<br />
Note that '''simobserve''' will not modify your original input image. Rather, it will make a copy ''m51c/m51c.skymodel''.<br />
<br />
We will override most of the parameters in the Halpha FITS image to make the image more suitable to a sub-millimeter ALMA observation. We will:<br />
* place the source in the southern hemisphere with the ''indirection'' parameter, <br />
* set the pixel size to 0.1arcsec, to simulate an observation of a galaxy that is smaller in angular size. (M51 itself would require a quite large mosaic, and in any case we'd like the input model pixels to be significantly smaller than the synthesized beam.)<br />
* set the peak brightness to 0.004 Jy/pixel<br />
* set the frequency to 330GHz, and since it's a 2D image we'll set the single "channel" width to be 50MHz, and peak brightness of 0.004 Jy/pixel. These parameters are plausible for observing a sub-mm emission line in a galaxy.<br />
<source lang="python"><br />
# Set model image parameters:<br />
indirection="B1950 23h59m59.96s -34d59m59.50s"<br />
incell="0.1arcsec"<br />
inbright="0.004"<br />
incenter="330.076GHz"<br />
inwidth="50MHz"<br />
</source><br />
<br />
----<br />
<br />
====Simulate 12m interferometric observation====<br />
[[Image:M51c.ALMA_0.5arcsec.skymodel.png|thumb|hexagonal mosaic overplotted on sky model]]<br />
<br />
We'll begin by simulating the observation as seen by the main 12 m ALMA array. We'll have '''simobserve''' calculate a hexagonal mosaic of pointings.<br />
<br />
The interface for '''simobserve''' provides an ''integration'' parameter, which is the dwell time at each pointing in the mosaic. We'll set that to 10 s. A real observation might dwell on each pointing for ~ 20 s and record 3 sets of visibilities (i.e. 6 s integration time) for each pointing. '''simobserve''' does not currently offer a capability to reproduce this behavior exactly, but the difference in UV coverage should be very minor between the real observation and the simulated one. You could try setting ''integration'' to ~ 20 s, which would more closely match the dwell time for some real observations.<br />
<br />
<source lang="python"><br />
# have simobserve calculate mosaic pointing locations:<br />
setpointings = True<br />
integration = "10s"<br />
mapsize = "1arcmin"<br />
maptype = "hex"<br />
pointingspacing = "9arcsec" # this could also be specified in units of the primary beam e.g. "0.5PB"<br />
</source> <br />
<br />
=====Calculate Visibilities=====<br />
When you provide a desired resolution in the ''antennalist'' parameter, '''simobserve''' will determine what array configuration to use. <br />
<source lang="python"><br />
obsmode = "int"<br />
antennalist = "ALMA;0.5arcsec"<br />
totaltime = "3600s"<br />
</source><br />
<br />
Run '''simobserve''', displaying graphics to the screen and also writing graphics output to files that can be found in the project directory, ''m51c''.<br />
<br />
<source lang="python"><br />
graphics = "both"<br />
simobserve()<br />
</source><br />
<br />
----<br />
<br />
====Simulate 12m total power observation====<br />
[[Image:M51c.aca.tp.skymodel.png|thumb|rectangular map overplotted on sky model]]<br />
<br />
Next we'll simulate a total power raster map of the same area, on a square grid. CASA simulation tools can not simulate <br />
true on-the-fly mapping (with smearing on timescales smaller than an integration time), but a square grid with a short integration time will provide a very accurate approximation.<br />
<br />
By virtue of CASA's global parameters, we already have project and image world coordinate system parameters set correctly. <br />
<br />
We need to change the mapping parameters to specify a square region a bit larger than the interferometric map.<br />
<source lang="python"><br />
integration = "10s"<br />
mapsize = "1.3arcmin"<br />
maptype = "square"<br />
</source> <br />
<br />
We'll observe on a different day. This doesn't really matter, but if you choose to simulate two different 12m ALMA configurations and combine them, you will have issues concatenating the datasets if they are simulated on the same day with the same antenna names. So, it's a good habit to change the day.<br />
<br />
<source lang="python"><br />
obsmode = "sd"<br />
sdantlist = "aca.tp.cfg"<br />
sdant = 0<br />
refdate = "2012/12/01"<br />
totaltime = "2h"<br />
</source><br />
<br />
Run '''simobserve''', displaying graphics to screen and to files.<br />
<br />
<source lang="python"><br />
simobserve()<br />
</source><br />
<br />
----<br />
<br />
====Simulate 7m ACA observation====<br />
[[Image:M51c.aca.i.skymodel.png|thumb|hexagonal map overplotted on sky model]]<br />
<br />
Next we'll add an ACA mosaic, with its larger primary beam.<br />
<br />
<source lang="python"><br />
integration = "10s"<br />
mapsize = "1arcmin"<br />
maptype = "hex"<br />
pointingspacing = "15arcsec"<br />
</source> <br />
<br />
We can specify an integral number of times to repeat the mosaic by setting ''totaltime'' to an integer string without units.<br />
<br />
<source lang="python"><br />
obsmode = "int"<br />
refdate = "2012/12/02"<br />
antennalist = "aca.i.cfg"<br />
totaltime = "3"<br />
</source><br />
<br />
Run '''simobserve''', displaying graphics to screen and to files.<br />
<br />
<source lang="python"><br />
simobserve()<br />
</source><br />
<br />
----<br />
<br />
====Deconvolve the visibilities back into an image====<br />
<br />
Next we use '''simanalyze''' to combine the three measurement sets and create a single image. <br />
<br />
There are many ways to do this, and you may wish to discuss options with scientists at your ARC. In this example we will use the total power image as a model when deconvolving the ACA image, and then use the result as a model when deconvolving the 12m interferometric image. This method tends to give low weight to the large spatial scales, but is simple to illustrate.<br />
<br />
It's possible to get better results if one used multiscale clean in the clean task (again using the lower resolution image as a model when deconvolving the higher resolution one). An alternative would be to create an image independently from each dataset, and then use the CASA feather task to combine them entirely in the image plane.<br />
<br />
If given a total power and interferometric measurement set, '''simanalyze''' will automatically create the total power image, <br />
then use it as a model and deconvolve the interferometric image. It is not recommended to do both interferometric images simultaneously.<br />
<br />
====First image total power and ACA with total power as a model====<br />
[[Image:M51c.aca.i.analysis.png|thumb|Total power combined (with relatively low weight) with ACA]]<br />
<br />
<source lang="python"><br />
default("simanalyze")<br />
project = "m51c"<br />
vis = '$project.aca.i.ms,$project.aca.tp.sd.ms' <br />
imsize = [512,512]<br />
cell = '0.2arcsec'<br />
analyze = True<br />
showpsf = False<br />
showresidual = False<br />
showconvolved = True<br />
simanalyze()<br />
</source><br />
<br />
====Next add the 12m interferometric data====<br />
[[Image:M51c.ALMA_0.5arcsec.analysis.3.4.png|thumb|TP+ACA combined (with relatively low weight) with 12m ALMA]]<br />
Here we explicitly have to set the ''modelimage'' to the result from the previous run.<br />
<br />
<source lang="python"><br />
default("simanalyze")<br />
project = "m51c"<br />
vis = '$project.ALMA_0.5arcsec.ms'<br />
imsize = [512,512]<br />
cell = '0.1arcsec'<br />
modelimage = "$project.aca.i.image"<br />
analyze = True<br />
showpsf = False<br />
showresidual = False<br />
showconvolved = True<br />
simanalyze()<br />
</source><br />
<br />
{{Simulations Intro 4.0}}<br />
{{Checked 4.0.0}}</div>Jbraatz