Robert G. Radocinski

Multiparameter linear least-squares fitting to Poisson data one count at a time

A standard problem in gamma-ray astronomy data analysis is the decomposition of a set of observed counts, described by Poisson statistics, according to a given multicomponent linear model, with underlying physical count rates or fluxes which are to be estimated from the data. Despite its conceptual simplicity, the linear least-squares (LLSQ) method for solving this problem has generally been limited to situations in which the number n(sub i) of counts in each bin i is not too small, conventionally more than 5-30. It seems to be widely believed that the failure of the LLSQ method for small counts is due to the failure of the Poisson distribution to be even approximately normal for small numbers. The cause is more accurately the strong anticorrelation between the data and the wieghts w(sub i) in the weighted LLSQ method when square root of n(sub i) instead of square root of bar-n(sub i) is used to approximate the uncertainties, sigma(sub i), in the data, where bar-n(sub i) = E(n(sub i)), the expected value of N(sub i). We show in an appendix that, avoiding this approximation, the correct equations for the Poisson LLSQ (PLLSQ) problems are actually identical to those for the maximum likelihood estimate using the exact Poisson distribution. We apply the method to solve a problem in high-resolution gamma-ray spectroscopy for the JPL High-Resolution Gamma-Ray Spectrometer flown on HEAO 3. Systematic error in subtracting the strong, highly variable background encountered in the low-energy gamma-ray region can be significantly reduced by closely pairing source and background data in short segments. Significant results can be built up by weighted averaging of the net fluxes obtained from the subtraction of many individual source/background pairs. Extension of the approach to complex situations, with multiple cosmic sources and realistic background parameterizations, requires a means of efficiently fitting to data from single scans in the narrow (approximately = 1.2 keV, HEAO 3) energy channels of a Ge spectrometer, where the expected number of counts obtained per scan may be very low. Such an analysis system is discussed and compared to the method previously used.

A BATSE Earth Occultation Catalog of 0.03-1.8 MeV Gamma-Ray Source Spectra and Light Curves for Phases 1-3 (1991-1994)

Using the powerful Earth-occultation technique, long-term, nearly continuous monitoring of the entire low-energy gamma-ray sky is now possible with the advent of BATSE, the Burst and Transient Source Experiment on board the Compton Gamma Ray Observatory (CGRO). In this paper, we present a catalog of 34 moderately strong gamma-ray sources measured by BATSE. It consists of 0.03-1.8 MeV photon spectra averaged over weeks and months, and light curves of the 35-200 keV flux, with 1 day resolution, covering the first three phases of the CGRO mission (1991 May through 1994 October). These results have been obtained using the JPL Enhanced BATSE Occultation Package (EBOP) developed under the CGRO Guest Investigator Program. The EBOP concept and approach are also described in some depth. This paper presents highlights extracted from a large EBOP database which has now been archived at the Compton Observatory Science Support Center (COSSC). This database contains a complete record of ~1200 daily source count rates in 14 energy channels along with the corresponding Poisson and systematic errors for 64 sources, including 30 not described here. An interface to XSPEC is included in the archive allowing conversion from count rates to photon fluxes. This paper therefore serves also as a reference and entree into the archive and provides an index and guide for those investigators using the EBOP system and database for their respective scientific investigations.

Gamma-ray observations of Co-56 in SN 1987A

During observations of SN 1987A 286 days after the supernova explosion using a JPL high-resolution gamma-ray spectrometer, a line feature was found in the net supernova spectrum at an energy of 1240.8 + or - 1.7 keV, with an intrinsic width of 8.2 + or - 3.4 keV FWHM, and a net flux of 0.0021 + or - 0.0007 photons/sq cm per s. This feature is interpreted as the 1238 keV line from the decay of Co-56 in the supernova remnant. The blueshift and intrinsic broadening of the 1238 keV line are consistent with models incorporating mixing of the radioactive cobalt in the expanding eject. 26 references.