Coordinates on the sky and the chip¶
In every marx simulation, one or more sources are placed at some sky position. marx simulates photons coming from that position, traces them through the mirror and gratings and finally places them on the chip. With a known aspect solution, chip coordinates can then be transformed back to sky coordinates. In general, this will not recover the exact sky position where a photon started out. A big part of that is scatter in the mirrors, which blurs the image (see Point Spread Function (PSF) for tests of the PSF). However, with a large number of photons, we can fit the average position which should be close to the real sky position.
In real observations, other factors contribute, such as the finite
resolution of the detectors (marx usually takes that into account, but it can
be switched of through the
--pixadj="EXACT" switch in
and the uncertainty of the aspect solution.
Within a single observation, positions will be less certain for fainter sources (due to Poisson statistics) and for sources at a larger off-axis angles (due to the larger PSF).
Chandra Orion Ultradeep project¶
The Orion Nebula Cluster (ONC) is a dense star forming region with about 1600 X-ray sources observed in the COUP survey by Getman et al (2005) . We simulate this field with marx and then run a source detection to check how well we recover the input coordinates. This will depend on the number of counts detected and the position in the field. To simplify the simulation input, we assume that all sources have flat lightcurves and are monoenergetic at the observed mean energy (the energy matters because the effective area is energy dependent and so is the PSF). We write a short C code that reads an input coordiante list and generates the photons in this manner. We compile the code, and call it as a USER Source.
For this field, we know the true input coordinates so we can check how well marx reproduces those. In the center of the field (about one armin) the coordiante error is less than the size of an ACIS pixel for all sources and the average error never grows much beyond 1 ACIS pixel even for far off-axis source. The upper envelope of the distribution of errors is approximate linear and reaches 1 arcsec at a distance of 200 arcsec. No strong correlation of coordiante error and count rate of the source is apparent, indicating that the dominant error is not just due to Poisson counting statistics.
Regular Grid (ACIS)¶
In this example we place a radial grid of sources on the sky. Each source emits an equal number of photons (exactly, no Poisson statistics) so that we can compare the accuracy of the position we recover. Note that the detected number of photons will be smaller for off-axis photons because of vignetting!
We write a short C code that generates the photons in this manner, compile it, and call is as a USER Source.
The input position is typically recovered to much better than one pixel for sources with a few hundred counts. There is a small systematic trend that needs to be studied further.
Regular grid (HRC)¶
Same as above, but with HRC-I as a detector.
The field-of-view for the HRC-I is larger for than for ACIS-I, but the PSF becomes very large at large off-axis angles and thus the positional uncertainty will be so large that a comparison to marx is no longer helpful to test the accuracy of the marx simulations.
In the central few arcmin the input position is typically recovered to better than 0.2 pixels for sources with a few hundred counts.