A Wavelet-Based Algorithm for the Spatial Analysis of Poisson Data
Wavelets are scalable, oscillatory functions that deviate from zero only within a limited spatial regime and have average value zero, and thus may be used to simultaneously characterize the shape, location, and strength of astronomical sources. But in addition to their use as source characterizers, wavelet functions are rapidly gaining currency within the source detection field. Wavelet-based source detection involves the correlation of scaled wavelet functions with binned, two-dimensional image data. If the chosen wavelet function exhibits the property of vanishing moments, significantly nonzero correlation coefficients will be observed only where there are high-order variations in the data; e.g., they will be observed in the vicinity of sources. Source pixels are identified by comparing each correlation coefficient with its probability sampling distribution, which is a function of the (estimated or a priori known) background amplitude. In this paper, we describe the mission-independent, wavelet-based source detection algorithm "WAVDETECT," part of the freely available Chandra Interactive Analysis of Observations (CIAO) software package. Our algorithm uses the Marr, or "Mexican Hat" wavelet function, but may be adapted for use with other wavelet functions. Aspects of our algorithm include: (1) the computation of local, exposure-corrected normalized (i.e., flat-fielded) background maps; (2) the correction for exposure variations within the field of view (due to, e.g., telescope support ribs or the edge of the field); (3) its applicability within the low-counts regime, as it does not require a minimum number of background counts per pixel for the accurate computation of source detection thresholds; (4) the generation of a source list in a manner that does not depend upon a detailed knowledge of the point spread function (PSF) shape; and (5) error analysis. These features make our algorithm considerably more general than previous methods developed for the analysis of X-ray image data, especially in the low count regime. We demonstrate the robustness of WAVDETECT by applying it to an image from an idealized detector with a spatially invariant Gaussian PSF and an exposure map similar to that of the Einstein IPC; to Pleiades Cluster data collected by the ROSAT PSPC; and to simulated Chandra ACIS-I image of the Lockman Hole region.

title: Evidence for reionization at z ∼ 6: Detection of a gunn-peterson trough in a z = 6.28 quasar
abstract: We present moderate-resolution Keck spectroscopy of quasars at z = 5.82, 5.99, and 6.28, discovered by the Sloan Digital Sky Survey (SDSS). We find that the Ly? absorption in the spectra of these quasars evolves strongly with redshift. To z ~ 5.7, the Ly? absorption evolves as expected from an extrapolation from lower redshifts. However, in the highest-redshift object, SDSSp J103027.10+052455.0 (z = 6.28), the average transmitted flux is 0.0038 ? 0.0026 times that of the continuum level over 8450 ? 20, on the optical depth to Ly? absorption at z = 6. This is a clear detection of a complete Gunn-Peterson trough, caused by neutral hydrogen in the intergalactic medium. Even a small neutral hydrogen fraction in the intergalactic medium would result in an undetectable flux in the Ly? forest region. Therefore, the existence of the Gunn-Peterson trough by itself does not indicate that the quasar is observed prior to the reionization epoch. However, the fast evolution of the mean absorption in these high-redshift quasars suggests that the mean ionizing background along the line of sight to this quasar has declined significantly from z ~ 5 to 6, and the universe is approaching the reionization epoch at z ~ 6.

title: Likelihood methods and classical burster repetition
abstract: We develop a likelihood methodology which can be used to search for evidence of burst repetition in the BATSE catalog, and to study the properties of the repetition signal. We use a simplified model of burst repetition in which a number Nr of sources which repeat a fixed number of times Nrep are superposed upon a number Nnr of non‐repeating sources. The instrument exposure is explicitly taken into account. By computing the likelihood for the data, we construct a probability distribution in parameter space that may be used to infer the probability that a repetition signal is present, and to estimate the values of the repetition parameters. The likelihood function contains contributions from all the bursts, irrespective of the size of their positional errors—the more uncertain a burst’s position is, the less constraining is its contribution. Thus this approach makes maximal use of the data, and avoids the ambiguities of sample selection associated with data cuts on error circle size. We present the results ...

title: Evidence that gamma-ray bursts repeat and are Galactic in origin
abstract: We describe work which has provided evidence that there are two distinct classes of γ-ray bursts, that the bursts repeat, and that they are Galactic in origin.

title: The Extinction toward the GRB 970228 Field
abstract: We determine the local Galactic extinction toward the field of gamma-ray burst GRB 970228 using a variety of methods. We develop a maximum-likelihood method for measuring the extinction that compares galaxy counts in different bands in the field of interest to those in a field of known extinction, and we apply this method to the GRB 970228 field. We also measure the extinction by comparing the observed broadband colors of stars in the GRB 970228 field to the colors of library spectra of the same spectral type. Finally, we estimate the extinction using the neutral hydrogen column density and the amount of infrared dust emission toward this field. Combining the results of these methods, we find a best-fit Galactic extinction in the optical of AV = 1.09, which implies that the intrinsic afterglow of GRB 970228 is significantly brighter and bluer than it appears.

title: High-energy observations of XRF 030723: Evidence for an off-axis gamma-ray burst?
abstract: We report High Energy Transient Explorer 2 (HETE-2) Wide Field X-ray Monitor (WXM) and French Gamma Telescope observations of XRF 030723 along with observations of the XRF afterglow made using the 6.5 m Magellan Clay telescope and Chandra. The observed peak energy E of the νFν burst spectrum is found to lie within (or below) the WXM 2-25 keV passband at 98.5% confidence, and no counts are detected above 30 keV. Our best-fit value is E = 8.4 keV. The ratio of X-ray to γ-ray flux for the burst follows a correlation found for GRBs observed with HETE-2, and the duration of the burst is similar to that typical of long-duration GRBs. If we require that the burst isotropic equivalent energy Eiso and Epk satisfy the relation discovered by Amati et al. (2002), a redshift of z = 0.38 can be determined, in agreement with constraints determined from optical observations. We are able to fit the X-ray afterglow spectrum and to measure its temporal fade. Although the best-fit fade is shallower than the concurrent fade in the optical, the spectral similarity between the two bands indicates that the X-ray fade may actually trace the optical fade. If this is the case, the late-time rebrightening observed in the optical cannot be due to a supernova bump. We interpret the prompt and afterglow X-ray emission as arising from a jetted GRB observed off-axis and possibly viewed through a complex circumburst medium that is due to a progenitor wind.

title: Measurement of di-jet cross-sections in photoproduction and photon structure
abstract: The production of hard di-jet events in photoproduction at HERA is dominated by resolved photon processes in which a parton in the photon with momentum fraction x γ is scattered from a parton in the proton. These processes are sensitive to the quark and gluon content of the photon. The differential di-jet cross-section d σ /dlog( x γ ) is presented here, measured in tagged photoproduction at HERA using data taken with the H1 detector, corresponding to an integrated luminosity of 7.2 pb −1 . Using a restricted data sample at high transverse jet energy, E T,jet >6 GeV, the effective parton density f γ ,eff (x γ )= q(x γ )+ q (x γ )+9/4 g(x γ ) in the photon in leading order QCD is measured down to x γ = 0.05 from which the gluon density in the photon is derived.

title: Search for leptoquark bosons in ep collisions at HERA
abstract: A search for scalar and vector leptoquarks coupling to first generation fermions is performed using the e(+) p and e(-) p scattering data collected by the H1 experiment between 1994 and 2000. The data correspond to a total integrated luminosity of 117 pb(-1). No evidence for the direct or indirect production of such particles is found in data samples with a large transverse momentum final state electron or with large missing transverse momentum. Constraints on leptoquark models are established. For leptoquark couplings of electromagnetic strength, leptoquarks with masses up to 275-325 GeV are ruled out at 95% confidence level. These limits improve and supercede earlier H1 limits based on subsamples of the data used here.

title: Search for compositeness, leptoquarks and large extra dimensions in eq contact interactions at HERA
abstract: The reaction e + p → e + X is studied with the H1 detector at Hera . The data cover momentum transfers Q 2 between 200 GeV 2 and 30 000 GeV 2 and correspond to an integrated luminosity of 35.6 pb −1 . The differential cross section d σ /d Q 2 is compared to the Standard Model expectation for neutral current scattering and analysed to search for ( e e)( q q) contact interactions. No evidence for new phenomena is observed. The results are used to set limits on scales within models of electron–quark compositeness, quark form factors and the exchange of virtual heavy leptoquarks. A search for gravitational effects mediated through the exchange of virtual gravitons which propagate into large extra dimensions is presented.

title: Deep inelastic inclusive e p scattering at low x and a determination of alpha(s)
abstract: A precise measurement of the inclusive deep-inelastic \(e^+ p\) scattering cross section is reported in the kinematic range \(1.5 \le Q^2 \le 150\) GeV\(^2\) and \( 3 \cdot 10^{-5} \le x \le ). The data were recorded with the H1 detector at HERA in 1996 and 1997, and correspond to an integrated luminosity of 20 pb\(^{-1}\). The double differential cross section, from which the proton structure function \(\) and the longitudinal structure function \(F_L(x,Q^2)\) are extracted, is measured with typically 1% statistical and 3% systematic uncertainties. The measured derivative \((\partial F_2(x,Q^2) / \partial \ln Q^2)_x \) is observed to rise continuously towards small x for fixed \(Q^2\). The cross section data are combined with published H1 measurements at high \(Q^2\) for a next-to-leading order DGLAP QCD analysis. The H1 data determine the gluon momentum distribution in the range \(3 \cdot 10^{-4} \le x \le 0.1\) to within an experimental accuracy of about 3% for \(Q^2 =20\) GeV\(^2\). A fit of the H1 measurements and the \(\mu p \) data of the BCDMS collaboration allows the strong coupling constant \(\alpha_s\) and the gluon distribution to be simultaneously determined. A value of \(\alpha_s (M_Z^2) =0.1150 \pm 0.0017 (exp)^{+0.0009}_{-0.0005} (model)\) is obtained in NLO, with an additional theoretical uncertainty of about \(\pm 0.005\), mainly due to the uncertainty of the renormalisation scale.

title: Measurement and QCD Analysis of Neutral and Charged Current Cross Sections at HERA
abstract: The inclusive e(+)p single and double differential cross sections for neutral and charged current processes axe measured with the H1 detector at HERA. The data were taken in 1999 and 2000 at a centre-of-mass energy of roots = 319 GeV and correspond to an integrated luminosity of 65.2 pb(-1). The cross sections are measured in the range of four-momentum transfer squared Q(2) between 100 and 30 000 GeV2 and Bjorken x between 0.0013 and 0.65. The neutral current analysis for the new e(+)p data and the earlier e-p data taken in 1998 and 1999 is extended to small energies of the scattered electron and therefore to higher values of inelasticity y, allowing a determination of the longitudinal structure function F-L at high Q(2) (110 - 700 GeV2). A new measurement of the structure function x (F) over tilde (3) is obtained using the new e(+)p and previously published e p neutral current cross section data at high Q(2). These data together with H1 low Q(2) precision data are further used to perform new next-to-leading order QCD analyses in the framework of the Standard Model to extract flavour separated parton distributions in the proton.