Wunderer, Cornelia

Imaging with the Test Setup for the Coded-Mask INTEGRAL Spectrometer SPI

Performance of a Coded Aperture Gamma-Ray Telescope at 60 keV - 8 MeV

Abbildung mit dem Testaufbau für das INTEGRAL-Spektrometer SPI
Verhalten eines Gamma-Teleskops mit kodierter Maske bei 60 keV bis 8 MeV

Thesis

Filetyp: PDF (.pdf)
Size: 14570 Kb

Schlüsselwörter:

kodierte Maske, INTEGRAL, SPI, Astronomie, Gamma-Strahlung, Teleskop, Entropiemaximierung, spiros, spiskymax, Punktquellenlokalisierung, Winkelauflösung

coded mask, coded aperture, INTEGRAL, SPI, astronomy, gamma-rays, telescope, maximum entropy, iterative removal of sources, spiros, spiskymax, point source localization accuracy, angular resolution

TUB Systematik
PHY 923d; PHY 909d
Schlagwortnormdatei
Gammastrahlenteleskop / Maske (Photographie) / Abbildung (Physik)
Sachgruppe der DNB
29 Physik, Astronomie


Dissertation eingereicht bei: Technische Universität München , Fakultät für Physik, 2002-12-20
Tag der mündlichen Prüfung: 2003-02-07


Abstract in English

The European Space Agency's INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) was launched on October 17, 2002. It carries two gamma-ray telescopes, the imager IBIS and the spectrometer SPI, together with X-ray and optical monitors. The spectrometer SPI is capable of measuring 20 keV - 8 MeV photons with an energy resolution of 2.3 keV at 1.3 MeV using 19 high-purity Germanium detectors. A tungsten-alloy coded aperture provides spatial information with an angular resolution of ~2.5° and a 16° fully coded field of view.
In order to investigate the imaging properties of SPI in more detail than was possible during the flight instrument calibration - which had to fit in the tight satellite integration schedule -, I designed and built the SPI Imaging Test Setup (SPITS) at the Max-Planck-Institute for extraterrestrial Physics. SPITS consists of a coded aperture and two Ge detectors. The coded mask is based on a SPI mask development model, has the same Tungsten-alloy HURA mask coding as SPI, and is made of SPI flight model materials. The two hexagonal Ge detectors in their Al capsules, taken from the SPI flight instrument manufacturing line, are housed in a common vacuum cap. Mounted on an XY-table, they can be moved to cover the 19 SPI Ge detector positions. Data from 11 XY-table positions are combined to one SPITS measurement of 19 Ge spectra. This approach is feasible in a laboratory environment using radioactive sources since neither source activity nor background vary significantly during one measurement.
SPITS was exposed to three types of sources: Laboratory point-like calibration sources with strong lines from 60 keV to 1836 keV, an extended source of 511 keV emission built from an 88Y source and a lead disk, and proton-irradiated accelerator targets emitting photons up to 9 MeV. All sources were placed 9 m from the SPITS detection plane. The laboratory point source measurements covered many incident directions not included in the SPI imaging calibrations, noteably seven incident directions which could be combined to the smaller of the standard SPI dithering patterns and a series of measurements that allowed the synthetization of dithered observations of two sources separated by 0.5°, 1.0°, 1.5°, etc. Neither a hexagonal dither pattern nor a dithered observation of sources separated by less than 1° can be synthesized from the SPI calibration data. In addition, the imaging of extended emission and an imaging observation of emission above 3 MeV have not been attempted with SPI during the instrument calibration.
To obtain the image of the "sky" from a deconvolution of the measured line counts in the 19 Ge detectors, the instrument response to photons from all possible incident directions must be known. The required Instrument Response Functions (IRFs) for SPITS with sources at 9 m distance have been generated using the ray-tracing coded aperture modelling platform CAPTIF. I discuss the CAPTIF mass model and compare it to a GEANT mass model and measured properties of both SPI and SPITS components. A comparison of IRF-predicted with SPITS-measured detector count rates for two incident directions and several photon energies shows that, after introducing a correction factor, IRF and measurement generally agree to within 2 sigma.
Both standard SPI image reconstruction algorithms, spiskymax and spiros, have been used to analyze the SPITS data. spiskymax, developed primarily for imaging extended emission, uses entropy maximization to determine the flux distribution. spiros ("SPI iterative removal of sources") is primarily geared towards point sources. In a cross-correlation image of the sky, the (remaining) brightest source is located and subsequently removed in dataspace; this process is repeated until all significant sources are found.
Since there is no easy and clear stopping criterion for maximum entropy analysis, I studied the performance of spiskymax with simulated datasets to determine the best number of iterations for the SPITS point source analysis as a function of source intensity.
The capability of the SPITS instrument to locate a point source depends strongly on the significance of the source and on the method of observation used ("staring" or "dithering"), as expected, and slightly on the source's angular distance from the instrument axis. For the observation of laboratory point sources, spiskymax and spiros location accuracies are routinely better than 10' even in non-dithered observations. Using both spiskymax and spiros, I find a correlation of the reconstructed positions with the spatial sky grid on which the IRFs were calculated.
By adding spectra obtained with a source at different positions, observations of two sources could be synthesized. From these, the capability of the imaging system (instrument, IRFs, and analysis software) to disentangle and correctly locate two close sources can be determined. Of course, this capability again depends on the significance of the sources. For non-dithered observations, spiros source localizations are not reliable (roughly correct locations in far fewer than 80% of the cases), and the two brightest local maxima in the spiskymax image are often more than 0.2° away from their true positions even for angular separations of 2°. With dithering, results improve dramatically: all spiros and almost all spiskymax reconstructions of sources separated by 1° or more result in both sources being localized to better than 0.2°. With dithering, for example, spiros location accuracies in the case of two 22Na sources are better than 8' for source separations of 1° and better than 5' for larger separation angles.
spiskymax reconstructions of the proton-irradiated targets in their respective gamma-ray lines also yielded reasonable results for non-dithered observations of one and two sources - the 19F target, for example, is positioned to better than 0.1° in several lines and for several different target positions. With spiros, however, the targets can often not be reasonably reconstructed from the same datasets. This is most likely due to the variability of "source activity" which makes the later combination of measurements at the 11 SPITS XY-table positions to one observation problematic. That these problems appear only with spiros might be an indication that it is more sensitive to discrepancies between the IRFs and the true instrument response than spiskymax.
The fluxes of the laboratory sources are reconstructed to within the combined statistical and systematic error bars by both spiskymax and spiros. The largest contribution to the overall error stems from the IRF and IRF correction factor uncertainties, followed by the uncertainty of the source activity. Source fluxes are underestimated more often than overestimated. This fits both with the tendency of IRF-predicted count rates to overpredict measured photopeak counts and with a property of the spiskymax algorithm seen in the simulated observations: source fluxes are initially underestimated and approach true fluxes asymptotically at higher iterations.
Finally, SPITS also proved the capability of such an imaging system (with the spiskymax algorithm) to reconstruct both extent and flux from an extended 511 keV source reasonably well.
Many of the lessons learned with SPITS, namely concerning properties of the imaging algorithms, the simulation results, and the (unintended) recovery of the IRF spatial grid in the distribution of point source locations, can be transferred directly to SPI. For others, the properties of the SPITS dataset and the IRFs used have to be kept in mind: strong sources were imaged in a low-background laboratory environment, the sources were placed 9 m from the Ge detection plane, and the IRFs used are cruder than those available for SPI.

Abstract in Deutsch

Das "INTERnational Gamma-Ray Astrophysics Laboratory" (INTEGRAL) der Europäischen Raumfahrtbehörde ESA wurde am 17. Oktober 2002 gestartet. Es hat zwei Gamma-Teleskope, den Imager IBIS und das Spektrometer SPI, sowie Röntgen- und optische Monitore an Bord. Das Spektrometer SPI kann 20 keV - 8 MeV Photonen mit einer Auflösung von 2,3 keV bei 1,3 MeV messen; es verwendet hierzu 19 hochreine Ge-Detektoren. Eine kodierte Maske aus einer Wolframlegierung liefert räumliche Information mit einer Auflösung von ~2,5° in einem voll kodierten Blickfeld von 16°.
Um die Abbildungseigenschaften von SPI detaillierter zu untersuchen, als dies während der Kalibration des Fluginstruments möglich war - die in den engen Zeitplan der Satellitenintegration eingefügt war -, erstellte ich das "SPI Imaging Test Setup" (SPITS) am Max-Planck-Institut für extraterrestrische Physik. SPITS besteht aus einer kodierten Maske und zwei Ge-Detektoren. Die kodierte Maske basiert auf einem SPI-Entwicklungsmodell und ist aus denselben Materialien wie die SPI-Maske. Die zwei hexagonalen Ge-Detektoren in ihren Al-Kappen stammen aus der Fluginstrumentfertigung und befinden sich in einer gemeinsamen Vakuumkappe. Sie können auf einem XY-Tisch bewegt werden, um die 19 Positionen der SPI Ge-Detektoren abzudecken. Daten von 11 Positionen des XY-Tisches werden zu einer SPITS-Messung von 19 Ge-Spektren kombiniert. Dies ist in einer Laborumgebung unter Verwendung von radioaktiven Quellen machbar, da weder Quellaktivität noch Hintergrund während einer Messung signifikant schwanken.
Mit SPITS wurden drei verschiedene Quelltypen gemessen: punktförmige Laboreichquellen mit starken Linien zwischen 60 keV und 1836 keV, eine ausgedehnte 511 keV-Quelle, bestehend aus einer 88Y-Quelle und einer Bleischeibe, und mit Protonen bestrahlte Beschleunigertargets, die Photonen bis 9 MeV emittierten. Alle Quellen wurden 9 m von der SPITS-Detektorebene plaziert. Die Messungen mit punktförmigen Laborquellen deckten viele Einfallsrichtungen ab, die in den SPI-Eichmessungen nicht vorkamen - insbesondere sieben Einfallsrichtungen, die zum kleineren der standardisierten SPI-Pointierungsrastermuster zusammengesetzt werden können, sowie eine Reihe von Messungen, die das Synthetisieren von gerasterten Beobachtungen von zwei Quellen ermöglichen, die 0,5°, 1,0°, 1,5°, usw. voneinander entfernt sind. Weder ein hexagonales Rastermuster noch eine gerasterte Beobachtung von Quellen mit Entfernungen von weniger als 1° können aus den SPI-Kalibrationsdaten generiert werden. Außerdem wurde weder eine Abbildung von ausgedehnter Emission noch eine abbildende Messung von Quellen oberhalb 3 MeV während der SPI-Instrument-Kalibration versucht.
Um das Bild des "Himmels" aus einer Dekonvolution der gemessenen Photonenanzahl in einer Linie in den 19 Ge-Detektoren zu erhalten, muß die Instrumentantwort auf Photonen aus allen möglichen Einfallsrichtungen bekannt sein. Die erforderlichen "Instrument Response Functions" (IRFs) für SPITS mit Quellen in 9 m Abstand wurden mit der Simulationsplattform für Teleskope mit kodierten Masken CAPTIF im "ray-tracing"-Verfahren erstellt. Ich diskutiere das CAPTIF-Massenmodell und vergleiche es mit einem GEANT-Massenmodell und gemessenen Eigenschaften von SPI- und SPITS-Komponenten. Nach Einführung eines Korrekturfaktors stimmen von den IRFs vorhergesagte und gemessene Detektorzählraten für zwei Einfallsrichtungen und mehrere Photonenenergien generell innerhalb 2 sigma überein.
Beide SPI-Standardbildanalyseverfahren, spiskymax und spiros, wurden zur Analyse von SPITS-Daten herangezogen. spiskymax, primär für die Abbildung ausgedehnter Emission entwickelt, verwendet Entropiemaximierung, um die Flußverteilung zu bestimmen. spiros ("SPI Iterative Removal of Sources") ist in erster Linie für Punktquellen gedacht. In einem Cross-Korrelationsbild des Himmels wird die hellste (verbleibende) Quelle lokalisiert und anschließend im Datenraum subtrahiert; dieser Prozeß wird wiederholt, bis alle signifikanten Quellen gefunden sind.
Da es kein klares, einfaches Kriterium für die Beendigung der Entropiemaximierungsanalyse gibt, habe ich das Verhalten von spiskymax mit simulierten Datensätzen studiert, um die optimale Iterationsanzahl für die SPITS-Punktquellenanalyse als Funktion der Quellstärke zu bestimmen.
Die Fähigkeit von SPITS, eine Punktquelle zu lokalisieren, ist, wie erwartet, stark von der Signifikanz der Quelle und der Beobachtungsstrategie (mit oder ohne Rastern) abhängig; außerdem besteht eine leichte Abängigkeit vom Winkelabstand der Quelle zur Instrumentachse. Bei der Beobachtung von Labor-Punktquellen erreichen spiskymax und spiros auch ohne Rastern regelmäßig Lokalisierungsgenauigkeiten unterhalb 10'. Sowohl mit spiskymax als auch mit spiros finde ich eine Korrelation der rekonstruierten Quellpositionen mit den räumlichen Gitterpunkten am Himmel, für die die IRFs berechnet wurden.
Beobachtungen von zwei Quellen konnten durch Addieren von Spektren, die mit einer Quelle an verschiedenen Positionen aufgenommen wurden, erzeugt werden. Auf dieser Basis kann die Fähigkeit des Abbildungssystems (Instrument, IRFs und Analysesoftware) bestimmt werden, zwei benachbarte Quellen zu trennen und korrekt zu lokalisieren. Selbstverständlich ist diese Fähigkeit wieder von der Signifikanz der Quellen abhängig. Bei nicht gerasterten Beobachtungen erweist sich spiros als nicht sehr verläßlich (grob korrekte Positionierung in weniger als 80% der Fälle), und die zwei hellsten Maxima im spiskymax-Bild sind selbst für Quellabstände über 2° oft mehr als 0,2° von den wahren Positionen entfernt. Rastern verbessert die Ergebnisse enorm: in allen spiros- und fast allen spiskymax-Rekonstruktionen von Quellen, die mehr als 1° voneinander entfernt sind, sind beide Quellen innerhalb von 0,2° korrekt lokalisiert. Zum Beispiel werden zwei 1° voneinander entfernte 22Na-Quellen von spiros in einer gerasterten Beobachtung auf 8' genau lokalisiert (auf 5' genau für größere Quellabstände).
Auch die spiskymax-Rekonstruktionen der Beschleunigertargets in den jeweiligen Gammalinien ergaben für nicht gerasterte Beobachtungen von ein und zwei Quellen vernünftige Resultate - das 19F-Target zum Beispiel wird in mehreren Linien und an mehreren Positionen auf besser als 0,1° lokalisiert. spiros kann die Targets jedoch aus denselben Datensätzen oft nicht sinnvoll rekonstruieren. Das liegt wahrscheinlich an der Variabilität der "Quellaktivität", die die spätere Kombination von Messungen an den 11 XY-Tischpositionen zu einer Beobachtung problematisch macht. Diese Probleme treten nur mit spiros auf - vielleicht ein Anzeichen, daß dieser Algorithmus sensibler auf Diskrepanzen zwischen IRFs und wahrer Instrumentantwort reagiert.
Die Flüsse der Laborquellen werden innerhalb der kombinierten statistischen und systematischen Fehler von spiskymax und spiros korrekt rekonstruiert. Der größte Beitrag zum Gesamtfehler ist die Unsicherheit in den IRFs und dem IRF-Korrekturfaktor, gefolgt von der Unsicherheit der Quellaktivität. Flüsse werden öfter unter- als überschätzt. Das passt sowohl zu der Tendenz der in den IRF vorhergesagten Zählraten, zu hohe Werte zu liefern, als auch zu einer Eigenschaft von spiskymax, die in den simulierten Beobachtungen sichtbar wurde: Quellflüsse werden initial unterschätzt und nähern sich bei höheren Iterationen asymptotisch den wahren Werten.
Schließlich hat SPITS auch die Fähigkeit eines solchen Abbildungssystems (mit spiskymax) unter Beweis gestellt, Ausdehnung und Fluß einer flächigen 511 keV-Quelle recht gut zu rekonstruieren.
Viele Ergebnisse der SPITS-Studie, so z.B. Eigenschaften der Abbildungsalgorithmen, Simulationsergebnisse und das (unbeabsichtigte) Wiederfinden der IRF-Gitter in der Verteilung der Punktquellenpositionen, können direkt auf SPI übertragen werden. Bei anderen muß man die Eigenschaften des SPITS-Datensatzes und der IRFs im Auge behalten: starke Quellen wurden in einer Umgebung mit niedrigem Hintergrund aufgenommen, die Quellen waren nur 9 m von der Detektorebene entfernt und die IRFs sind grober als die für SPI verfügbaren.

Betreuer Schönfelder, V.; apl. Prof. Dr.
Gutachter Böni, P. K.; Univ.-Prof. Dr.

Upload: 2003-02-19
URL of Theses: http://tumb1.biblio.tu-muenchen.de/publ/diss/ph/2003/wunderer.pdf

Technische Universität München, Universitätsbibliothek
Arcisstr. 21, D-80333 München

Unversehrtheit der Publikation
Unversehrtheit der Publikation