MPE/OGIP Calibration Memo CAL/ROS/93-015

ROSAT PSPC

The Off-Axis Point Spread Function

ROSAT PSPC

The Off-Axis Point Spread Function

Günther Hasinger, Günter Boese & Peter Predehl

MPE
Garching
Germany

&

Codes 668 and 664,
NASA/GSFC,
Greenbelt, MD20771

Version: 1995 May 08

SUMMARY

This memo details the extension of the description of the point spread function of the ROSAT PSPC to include off-axis effects. We present the off-axis psf algorithm and compare the model to in-flight data. We discuss the limitations of this parameterization with respect to data analysis.

LOG OF SIGNIFICANT CHANGES

Release Sections Changed Brief Notes
Date
1993 Sep 18 Published in Legacy, 4, 40
1994 Apr 22 All Updates and corrections
1995 Jan 31 All Made compatible with LaTeX2HTML software
1995 Feb 07 Section 4 Added Figures 4a - 4e as EPSFs

1  THE DATA & MODEL

As detailed in (Hasinger et al 1992), the point spread function (psf) of the ROSAT X-ray mirror assembly (XMA) + PSPC is a convolution of 5 components:

  1. The XMA scattering profile
  2. The off-axis blur of the XMA
  3. The intrinsic spatial resolution of the PSPC
  4. Focus and detector penetration effects
  5. A widening of the PSPC resolution due to the existence of `ghost images' (for PHA channels less than 20 only)

There is also some residual ellipsoidal blur due to residual attitude motion, but for the PSPC, this component is usually negligable compared to the other effects (see for more details). The XMA scattering profile, the intrinsic spatial resolution of the PSPC, and focus and detector penetration effects have already been described in and (Hasinger et al. 1993). The psf was parameterized for on-axis observations from PANTER ground calibration measurements. In it was shown that the derived analytical functions satisfactorily described 5 moderate signal-to-noise datasets in all but the softest band where `Ghost Imaging' was already known to be a problem. Here we extend the treatment to include a parameterization of the change in shape of the psf with off-axis angle.

Examination of a bright point source observed off-axis reveals asymmetry in the distribution of source counts which becomes easily noticable once the source is observed outside of the central ring of the PSPC (i.e. greater than 20 arcmin off-axis). The asymmetry is produced because for off-axis irradiation half of the mirror sees a different grazing angle to the other. The two relatively sharp outer boundaries of the blur area, which have different radius of curvature, are produced by the front and rear aperture planes of the mirrors, respectively. The hole in the psf is due to the central hole in the mirror aperture (the mirror consists of three shells) defocussed at large angles as the detector is flat but the focal plane is curved. The radial struts also interfere with the photon distribution as they shadow part of the detector. This combination produces the diffraction-like radial pattern seen on bright sources observed at large offset angles. The images of a source at an offset angle of 40 arcmin for inflight data and for the PANTER calibration data are shown in Figures 1a and 1b (respectively), where the center of the detector is towards the upper left hand side of the image in each case.

While early attempts to parameterize the change in shape of the psf concentrated on the ground calibration data taken at the PANTER facility, a comparison of those data with in-flight data showed significant differences. These differences are due to the finite source distance in the PANTER facility, which results in a beam divergence of about 10 arcmin. The consequence of this is that the PANTER data have the hole on the inner side of the peak of counts for off-axis sources, while the in-flight data show the opposite effect (compare Figures 1a and 1b). In addition, the PANTER image shows greater interference between the psf and the PSPC wire grid (effects which are reduced by the wobble of in-flight observations). As a consequence of these differences, in flight data, collected from a variety of bright serendipitous sources detected in pointed observations, were finally used to parameterize the off-axis dependence of the psf, while the components which are independant of off-axis angle were carried over from the higher-quality on-axis data.

The off-axis blur of the telescope, although highly structured and asymmetric, has been modeled by a simple Gaussian for comparison to a radially integrated profile. This Gaussian term for the off-axis blur is added in quadrature to the Gaussian term describing the intrinsic PSPC resolution. Since the model approximates the behaviour of the detector by the addition of Gaussian, exponential and Lorentzian terms, instead of a mathematically correct convolution of terms, then we have to adjust the contribution of the other terms relative to the on-axis case, to allow for the off-axis behaviour of the Gaussian term. The new term for the exponential fraction describes the relative diminishing of this term with increasing offset angle, due to the increase of the Gaussian term.

2  The Off-Axis Parameterization

Throughout,

q is the off-axis angle (in arcmin) from the center of the PSPC field of view,
r is the radius (in arcmin) from the PSF(E,q) centroid,
PSF(E,r, q) is the psf (in normalized  cts  arcmin-2),
and E is the photon energy (in keV).

2.1  The mirror scattering term

The analytical form derived for the fraction of photons within this component as a function of energy is given by:

Fscatt = 0.075 E+1.43
(1)
The increase in Fscatt with E is slower than the E2 powerlaw expected due to the detailed structure (micro-roughness) of the mirror surface.

As in the on-axis case, this steepens to a powerlaw at large r, ie

PSFscatt(E,r)
= Ascatt(E)
( 1 + (2 r/rscatt(E))2)
 
 for  r rb
= PSF(E,rb)  (r/rb)-a
 for  r > rb
(2)
where
rb(E) = 861.9/60  E-1  arcmin
(3)
and
a = 2.119  +  0.212 E
(4)

The width of the lorentzian, rscatt, as a function of energy was found to be:

rscatt(E) = 79.9/60  E-1   arcmin
(5)

The normalization of the lorentzian, Ascatt,is given by

Ascatt = Fscatt(E) /Iscatt normalized  cts  arcmin-2
(6)
where1
Iscatt = p(rscatt/2)2 \ [ ln(1+ (2/rscatt)2 rb2) + 2 rb2 ( (rscatt/2)2 + rb2 )-1 (a- 2)-1 ]
(7)
is the renormalization factor such that2



0 
PSFscatt(E,r)  2 pr  dr = Fscatt(E)
(8)

2.2  The exponential focus and penetration depth term

The analytical form derived for the fraction of photons within this component as a function of energy is given by whichever is the smaller of the following two terms:

 Fexp(E,q) = 10-1.635  +  0.639 E  +  0.052 E2 e-0.5 (q/12)2
(9)
OR
 Fexp(E,q) = 1.0 - Fscatt
(10)

Equation 10 is necessary to avoid the ëxponential artifact" of which made the psf parameterization invalid above 2 keV.

The combined contribution of focussing and the finite penetration of photons into the counter on the psf can still be modelled as an exponential function where r is the radial distance from the centroid position of the radial profile, in arcmin.

PSFexp(E,r) = Aexp(E)  exp

- r
rt(E)


 normalized  cts  arcmin-2
(11)

The e-folding angle, rt, was found to be the same as for the on-axis case

rt(E) = 1
60
 ( 50.61 E-1.472 + 6.80 E+5.62 )1/2  arcmin
(12)

The normalization, Aexp(E), is given by

Aexp(E) = Fexp(E) / ( 2 p(rt(E))2 ) normalized  cts  arcmin-2
(13)
and as for the other components has been normalized such that the integral of PSFexp(E,r) 2 pr dr over the range 0 r equals Fexp(E).

2.3  The Gaussian intrinsic resolution and mirror blur term

All photons not in the other 2 components are assumed to be in this component. Thus the fraction of photons within this component as a function of energy is given by:

Fint = 1  -  Fscatt  -  Fexp
(14)

Note that this becomes zero above 2 keV from equation 10.

PSFint(E,r) = Aint(E)  exp

- 1
2
( r
Gs(E)
)2

 normalized  cts  arcmin-2
(15)
The derived analytical form for the width of the detector gaussian, Gs(E, q), as a function of energy is given by:
rs(E) = 1/60 \ ( 108.7 E-0.888 + 1.121 E+6 )1/2 arcmin
(16)

The analytical form for the Gaussian mirror blur is

ms(q) = 1/60 \ ( 0.219 q2.848 )1/2 arcmin
(17)

These two Gaussian widths are combined as follows:

Gs(E, q) = ( rs2 + ms2)1/2
(18)

The normalization, Aint(E), is given by

Aint(E) = Fint(E) / ( 2 p(Gs(E))2 ) normalized  cts  arcmin-2
(19)
where the denominator is the renormalization factor such that the integral of PSFint(E,r) 2 pr dr over the range 0 r equals Fint(E).

Figure 2 compares the mirror Gaussian width (ms) with the detector Gaussian width (rs), versus angle (q). Figure 3 shows the total Gaussian width (Gs), versus angle. From these figures we see that above ~ 14 arcminutes, the mirror Gaussian term becomes dominant over the detector Gaussian term. Figures 4,a 4b, 4c, 4d & 4e show the change of shape of the components of the psf with increasing off-axis angle, for 5 illustrative energies.

2.4  Test Results

The data shown here are a representative subset of the data used for testing the off-axis psf algorithm. Radial profiles were extracted centred on the source. When comparing the model and data at large angles (i.e. in the outer region of the PSPC) it was crucial to define the center of weight of the counts (which lies between the peak and hole) such that we could produce a radial average of these features in each bin. We used the IRAF imcntr routine (proto.imcntr). Given the approximate coordinates of the center of an object, imcntr will compute a more accurate center using the algorithms described in the Kitt Peak publication "Stellar Magnitudes from Digital Images" under the Mountain Photometry Code section. Briefly, this algorithm computes the sum of all the rows and the sum of all the columns in the extraction box. Because the centers are computed independently for x and y, the result may be considered inferior to a true two-dimensional centering algorithm, but in practice the results appear to be very usable.

The lowest 11 PI channels were rejected altogether to exclude problems due to the variable lower limit discriminator for valid events due to the variable instrument gain which is folded into these data. As described in we do not expect an acceptable fit to be obtained in the B band, due to the "ghost imaging" effect. Additional sources falling within the specified annuli were masked out of the analysis. No background subtraction was carried out. Background rates were measured from the events file and were later added to the predicted profile template for each band. As in a spectrum was also extracted for each source and the predicted model could then be appropriately weighted with source spectrum, for the energy band selected in each case. In addition, the model was weighted with the amount of time the source spent at each position on the detector, as the detector was wobbled across the sky during the observations.

3  RESULTS

These predicted psf templates are overlaid on some sample datasets in Figures 5a, 5b, 5c, 5d & 5e. We stress that the normalization of the model psf in each case was also calculated using the equations presented here, thus NO FITTING was performed. The normalization of the predicted psf was calculated such that the integral under the predicted template is equal to the integral under the observed psf.

It can be seen that even at this high signal-to-noise ratio, the psf model provides an acceptable parameterization of the source profile for all but the B-band (as expected). This model can be used to search for extended emission (the limits of accuracy are evident by comparison with the plots in Figures 5a, 5b, 5c, 5d & 5e) and to determine the extraction cell required for spectral and timing analysis.

4  CONCLUSIONS

An acceptable model for the off-axis dependance of the ROSAT PSPC psf has been produced which can be used to estimate the correct extraction radii required for data analysis of off-axis sources. In addition, this can be used to determine whether off-axis sources are significantly extended. The plots provided illustrate the limits of the accuracy of this parameterization (particularly important for determination of significant extent) and PIs are invited to contact the GOF team if they wish to discuss any aspects of this parameterization. The comparison between model and data is most difficult for sources observed in the outer region of the PSPC detector, where radial profiles must be carefully centred to obtain meaningful results.

The contact person at GSFC is Jane Turner (LHEAVX::TURNER)

ACKNOWLEDGMENTS

We thank Dave Davis for his help in the original extraction of some of the off-axis data.

REFERENCES

Hasinger, G., Turner, T.J., George, I.M. & Boese, G., 1992. Legacy, 2, 77 Hasinger, G., Turner, T.J., George, I.M. & Boese, G., 1993. Legacy, 3, 46

FIGURES

Figure

Figure 1: (a) The in-flight psf for a point source located 40 arcmin off-axis. The center of the detector is towards the upper left of the figure.

Figure

Figure 1: (b) The PANTER psf for a point source located 40 arcmin off-axis. The center of the detector is towards the upper left of the figure.

Figure

Figure 2: The widths of the mirror (ms) and detector (rs) components to the Gaussian term. This illustrates how the mirror effect begins to dominate outside of ~ 14 arcmins. The lines of constant detector component represent (from top to bottom) 0.28, 0.50, 2.0 and 1.0 keV (the detector Gaussian has a minimum between 1 and 2 keV, as detailed in Hasinger et al 1992, thus the inversion).

Figure

Figure 3: The width of the mirror plus detector Gaussian term versus off-axis angle. The lines represent (from top to bottom) 0.28, 0.50, 2.0 and 1.0 keV (the detector Gaussian has a minimum between 1 and 2 keV, as detailed in Hasinger et al 1992, thus the inversion).

Figure

Figure 4: (a) The components of the PSPC psf versus off-axis angle at 0.188 keV.

Figure

Figure 4: (b) The components of the PSPC psf versus off-axis angle at 0.284 keV.

Figure

Figure 4: (c) The components of the PSPC psf versus off-axis angle at 0.5 keV.

Figure

Figure 4: (d) The components of the PSPC psf versus off-axis angle at 1.0 keV.

Figure

Figure 4: (e) The components of the PSPC psf versus off-axis angle at 1.7 keV.

Figure 5: (a) The observed psf of the ROSAT PSPC in the B-band (0.1-0.188 keV) compared to that predicted from the ground calibration data, for sources at off-axis angles 11, 13 and 38 arcmin. The misfit in the B-band is due to the ghost imaging problem.

Figure 5: (b) As for Figure 5 but for the C-band (0.188-0.284 keV). The fits for sources within the ring ( < 20 arcmin) are very good. The problem with the 38 arcmin offset case is the large degree of asymmetry in the data. We consider these fits to be adequate for most purposes.

Figure 5: (c) As for Figure 5 but for the R1-band (0.284-0.5 keV).

Figure 5: (d) As for Figure 5 but for the R2-band (0.5-1.0 keV).

Figure 5: (e) As for Figure 5 but for the R3-band (1.0-2.0 keV).

USEFUL LINKS TO OTHER HTML PAGES

The following useful links are available (in the HTML version of this document only):


Footnotes:

1 Versions prior to 1994 April 22 contained an additional erroneous term

2 Note that due to a typographical error there was a factor 2 missing from the following equation in all versions of prior to 1992 Oct 05 (the software used to generate the figures was however correct)


File translated from TEX by TTH, version 2.25.
On 22 Aug 2001, 22:50.