# Corrupt: Difference between revisions

No edit summary |
No edit summary |
||

Line 1: | Line 1: | ||

{{Simulations Intro}} | {{Simulations Intro 4.0}} | ||

[[Category: Simulations]] | [[Category: Simulations]] | ||

## Revision as of 10:18, 2 November 2012

↵ **Simulating Observations in CASA**

The descriptions on this page are rather anemic at the moment, so please read this pdf explanation of how this is all done

## Thermal noise

The task allows you to add thermal noise from the atmosphere and from the
ALMA receivers according to current ALMA specifications according to observing band.
At time of writing, `simobserve` knows the receiver temperature specifications for ALMA and EVLA, and a little about the altitude and pressure at the ALMA site.

For interferometric simulation, one can select **thermalnoise** = *tsys-manual* or *tsys-atm*. For single dish simulation, only *tsys-manual* is currently available. The manual method requires specifying **t_ground**, i.e. spillover temperature, **t_sky** the physical temperature in the atmosphere, and **tau0** the zenith opacity in the center of the simulated band. *tsys-atm* constructs an ATM atmospheric model from the site characteristics and the user specified Precipitable Water Vapor *user_pwv*, and the user must still specify **t_ground**.

NOTE: adding noise can be slow (you have to generate a random number for every visibility). For purposes of learning the software, and doing zeroeth order parameter studies, consider increasing the integration time (to what would be an unreasonably large value for a real observation, but may still work to explore simulated parameter space). See Etime study for details.

For more flexibility and control, one can use the `sm` tool. run simobserve *without* noise, and niter=0 to not produce a clean image. Optionally make a copy of the noise-free measurement set, and then do the following:

```
sm.openfromms("my.ms")
sm.setnoise( ... options ... )
sm.corrupt()
sm.done()
```

`sm.setnoise` takes as arguments the atmospheric and ground/ambient temperatures, antenna efficiencies, site characteristics, etc. There are three modes:

**simple**: specify**simplenoise**="1Jy" to get random Gaussian noise with 1Jy RMS**tsys-atm**: use environment temperatures, antenna parameters, and the aatm library to create a model of the troposphere and add random noise of the appropriate magnitude to the visibilities**tsys-manual**: specify atmospheric brightness temperature and optical depth yourself (rather than let aatm calculate it for you) and apply noise of the corresponding magnitude.

- see the pdf mentioned above for the precise equations used

For any of the `sm` corruption methods, an actual calibration table can be written out, and plotted with `plotcal` to make sure you're introducing the magnitude and type of corruption that you desire.

## Atmospheric Phase Noise

The troposphere is refractive, so introduces a phase delay, and is not uniform, so that phase delay differs between antennas. Excessive phase noise decorrelates the signal. In CASA, you can use ` sm.settrop` to corrupt your measurement set with phase noise in two different ways:

Full corruption will make your data look worse than real ALMA data should look. ALMA will use Water Vapor Radiometers (WVRs) to correct for the atmospheric phase delay in real time during the observation. Once the full hardware and software system is operational, the residual phase errors should be small (a few microns of equivalent PWV variation), roughly Gaussian, and independent on each antenna, so a more realistic portrayal of that may be to use `sm.setgain` with the real part of the gain equal to zero, or to use **mode**="individual" below.
The full phase screen treatment is interesting to explore the full uncorrected corruption, to test other telescopes, and to test heuristics for WVR calculations. For more information, please see ALMA memos 587, 582,
573, 568, 535, 523, 515, 496, 495, 491, 490, 451, 415, 404, 361, 332, 252, 210, 209, 176 ...

For more information about the atmosphere at Chajnantor, see memos 542, 529, 521, 517, 512, 500, 497, 459, 384, 365, 363, 345, 334, 238, 237, 187, 186 ...

**mode**="individual": the delay varies with time according to fractional Brownian Motion (fBM), or generalized 1/f noise. This mathematical description is very good for many natural stochastic processes, and used very widely for e.g. CGI simulations of fractal landscapes, turbulent motions and structures in clouds and water, electronic gain drift, etc.**mode**="screen": a 2-dimensional phase screen is generated with the same fBM mathematical model, and blows across the array at the specified wind speed. This creates much more realistic correlations in time and space between antennas. In practice, what is created is a fluctuation of precipitable water vapor (PWV) screen, and then the aatm library and tropospheric model is used to calculate the corresponding phase delay as a function of frequency.

old method, unsupported, YMMV: Rob Reid wrote a python script to add phase noise,
noisify_alma.py. noisify_alma gets its parameters from the ALMA specifications where it can, and an approximate Komolgorov model of the atmosphere such that the phase difference
between two antennas increases with baseline length.

- note that noisify_alma.py also does thermal noise. don't add twice as much as you intended.

## Gain

`sm.setgain`**mode**="random" generates random Gaussian complex gains with the given real and imaginary RMS**mode**="fbm" generates a fractional Brownian Motion complex gain drift or time-dependent variation, independently for each antenna

## Leakage a.k.a. Cross-polarization

`sm.setleakage`- random Gaussian (time-independent) complex cross-polarization leakage terms, with specified real and imaginary amplitudes, and an optional systematic offset from zero
- at time of writing, sm and simdata only calculate cross-polarization in the linear approximation (i.e. your cross-hands will change but Stokes I should not). Full treatment will be added in the future.

see also alma memo 288