N891 simdata (CASA 3.4): Difference between revisions

From CASA Guides
Jump to navigationJump to search
No edit summary
Line 15: Line 15:


==Simulation==
==Simulation==
To begin, reset all simobserve parameters to their default values, set '''project''' to define the prefix for all output files, and set '''skymodel''' to the above FITS file.


<source lang="python">
<source lang="python">
Line 25: Line 27:
</source>
</source>


* Units: K - first convert to flux surface brightness: Jy/Sr = 2x10<sup>23</sup> k T / &lambda;<sup>2</sup>, <!-- <math>\frac{Jy}{Sr} = \frac{2\times 10^{23} k T}{\pi D^2 \Omega}</math>, where <math>\Omega</math> is the beam solid angle --> = 4x10<sup>8</sup>T at 110GHz.
''Units: K - first convert to flux surface brightness: Jy/Sr = 2x10<sup>23</sup> k T / &lambda;<sup>2</sup>, <!-- <math>\frac{Jy}{Sr} = \frac{2\times 10^{23} k T}{\pi D^2 \Omega}</math>, where <math>\Omega</math> is the beam solid angle --> = 4x10<sup>8</sup>T at 110GHz.'''


<source lang="python">
<source lang="python">
Line 33: Line 35:
</source>
</source>


* Now we need to decide if this model data will work at the desired pixel scale
* The GRS resolution of 40" at ~10kpc is 0.04" at 10Mpc, so we should be able to do a simulation of observing at ~0.1-0.2".  The resolution plot (See Figure 1) indicates that for ALMA at 100GHz, configuration 20 is appropriate.
[[Image:Beamsummary.png|thumb|Figure 1: Resolution plot.]]
[[Image:Beamsummary.png|thumb|Figure 1: Resolution plot.]]
* If we intend to set <tt>incell=0.2arcsec</tt> in <tt>simobserve</tt>, then the cube needs to be multiplied by 4x10<sup>8</sup> * (.04/206265)<sup>2</sup> = 1.4x10<sup>-5</sup> to obtain Jy/pixel.  The cube peaks at ~20K, so we can perform the simulation with <tt>inbright=3e-4</tt>, which should yield a peak of ~1mJy/bm.
 
We want to simulate an observation of this model data as though it were 10-Mpc distant.  The Galactic Ring Survey resolution is 40" at ~10 kpc, which translates to 0.04" at 10 Mpc.  The plot at right shows how the synthesized beam size changes with the full science array configurations provided with CASA.  For a 100-GHz observation, configuration 20 will give us the appropriate synthesized beam size.
 
''If we intend to set <tt>incell=0.2arcsec</tt> in <tt>simobserve</tt>, then the cube needs to be multiplied by 4x10<sup>8</sup> * (.04/206265)<sup>2</sup> = 1.4x10<sup>-5</sup> to obtain Jy/pixel.  The cube peaks at ~20K, so we can perform the simulation with <tt>inbright=3e-4</tt>, which should yield a peak of ~1mJy/bm.''


* Will we be dominated by the noise in the input model? Input noise ~150mK or S/N~20, so at our scaled intensity, ~0.05 mJy/bm. The [http://almascience.eso.org/call-for-proposals/sensitivity-calculator ALMA Sensitivity Calculator] says that ALMA will achieve 2.5mJy/bm in 2 hours for the input 212m/s channel width (0.075MHz), so the noise in the input model should not affect our results.   
* Will we be dominated by the noise in the input model? Input noise ~150mK or S/N~20, so at our scaled intensity, ~0.05 mJy/bm. The [http://almascience.eso.org/call-for-proposals/sensitivity-calculator ALMA Sensitivity Calculator] says that ALMA will achieve 2.5mJy/bm in 2 hours for the input 212m/s channel width (0.075MHz), so the noise in the input model should not affect our results.   

Revision as of 12:16, 12 June 2012

Simulating Observations in CASA

This tutorial simulates ALMA interferometric data for a nearby edge-on spiral galaxy.

Roughly modeled after NGC891

Data

For our model, we will use a data cube of 13CO from the Galactic Ring Survey on the 14-meter FCRAO. To speed the simulation, this FITS file contains the data we want binned to a coarser velocity resolution.

Simulation

To begin, reset all simobserve parameters to their default values, set project to define the prefix for all output files, and set skymodel to the above FITS file.

# In CASA
# Initializing simobserve
# Laying down some basic ground rules
default 'simobserve'
project = 'n891d'
skymodel = 'grs-12kms.fits'

Units: K - first convert to flux surface brightness: Jy/Sr = 2x1023 k T / λ2, = 4x108T at 110GHz.'

# In CASA
# Setting the new frequency of the central channel
incenter = '110.1777GHz'
Figure 1: Resolution plot.

We want to simulate an observation of this model data as though it were 10-Mpc distant. The Galactic Ring Survey resolution is 40" at ~10 kpc, which translates to 0.04" at 10 Mpc. The plot at right shows how the synthesized beam size changes with the full science array configurations provided with CASA. For a 100-GHz observation, configuration 20 will give us the appropriate synthesized beam size.

If we intend to set incell=0.2arcsec in simobserve, then the cube needs to be multiplied by 4x108 * (.04/206265)2 = 1.4x10-5 to obtain Jy/pixel. The cube peaks at ~20K, so we can perform the simulation with inbright=3e-4, which should yield a peak of ~1mJy/bm.

  • Will we be dominated by the noise in the input model? Input noise ~150mK or S/N~20, so at our scaled intensity, ~0.05 mJy/bm. The ALMA Sensitivity Calculator says that ALMA will achieve 2.5mJy/bm in 2 hours for the input 212m/s channel width (0.075MHz), so the noise in the input model should not affect our results.
# In CASA
# Setting the new channel width
inwidth = '0.075MHz'
  • We do have a sensitivity issue though - if we decrease the spectral resolution by a factor of 6 (bin the input channels in some other program - simobserve will know how to do that in the future but not yet), and plan for 3 8-hr tracks, then the sensitivity calculator suggests that we'll get <0.25mJy rms, or S/N>10 per beam. Rather than simulate 3 days of observing, I'll increase inbright by sqrt(3) and simulate one 8 hour track.
Figure 2: here's the cube with the simobserve's scaling and World Coordinate System


# In CASA
# Scaling the surface brightness
inbright = '1.4e-4'
  • the ALMA 12m primary beam is 50" so we'd space a mosaic by 25", but the model cube has 326x357 pixels, or 13 arcsec with our small pixels. That's a lot smaller than the primary beam, so it doesn't matter much what output image size we ask for.
# In CASA
# Finish up the image model, and setting up the pointing
indirection = 'J2000 7h00m34 -23d03m00'
incell = '0.2arcsec'
setpointings = True
integration = '300s'
pointingspacing = '25arcsec'
mapsize = '60arcsec'

There are 659 channels in the input cube, but as noted above we want to bin those to 109 channels of 1.2 km/s each.

# In CASA
# Finish up the rest of the settings for this run of simobserve
graphics = 'both'
verbose = True
overwrite = True
obsmode = 'int'
antennalist = 'alma;0.5arcsec'
totaltime = '3600s'
simobserve()  # Run simobserve to create the simulated data we need
default 'simanalyze'  
project = 'n891d'
image=T
imsize = [336,360]
vis = project+'.alma_0.5arcsec.ms'
simanalyze()  # All other default settings are OK in simanalyze


N891.grs-24-cube.coord18-59-59.976-40d00m01.972.png
Figure 3: a spectral profile in the box marked in green

Input:
N891d.alma 0.5arcsec.skymodel.png
Predict:
N891d.alma 0.5arcsec.observe.png
Image:
N891d.alma 0.5arcsec.image.png
Analyze:
N891d.analysis.png

Figure 4: Sample results

Last checked on CASA Version 3.4.0.