Absolute efficiency estimation of photon-number-resolving detectors using twin beams

A nonclassical light source is used to demonstrate experimentally the absolute efficiency calibration of a photon-number-resolving detector. The photon-pair detector calibration method developed by Klyshko for single-photon detectors is generalized to take advantage of the higher dynamic range and additional information provided by photon-number-resolving detectors. This enables the use of brighter twin-beam sources including amplified pulse pumped sources, which increases the relevant signal and provides measurement redundancy, making the calibration more robust

A Lyman-alpha blob in the GOODS South field: evidence for cold accretion onto a dark matter halo

We report on the discovery of a z = 3.16 Lyman-alpha emitting blob in the GOODS South field. The blob has a total Ly-alpha luminosity of ~ 10^(43) erg s^(-1) and a diameter larger than 60 kpc. The available multi-wavelength data in the GOODS field consists of 13 bands from X-rays (Chandra) to infrared (Spitzer). Unlike other discovered Ly-alpha blobs, this blob shows no obvious continuum counter-part in any of the broad-bands. In particular, no optical counter-parts are found in the deep HST/ACS imaging available. For previously published blobs, AGN (Active Galactic Nuclei) or 'superwind' models have been found to provide the best match with the data. We here argue that the most probable origin of the extended Ly-alpha emission from the blob in the GOODS South field is cold accretion onto a dark matter halo.Comment: 4 pages, 2 tables, 2 figures, Accepted to A&A Letters, minor changes to tex

Correlator Bank Detection of GW chirps. False-Alarm Probability, Template Density and Thresholds: Behind and Beyond the Minimal-Match Issue

The general problem of computing the false-alarm rate vs. detection-threshold relationship for a bank of correlators is addressed, in the context of maximum-likelihood detection of gravitational waves, with specific reference to chirps from coalescing binary systems. Accurate (lower-bound) approximants for the cumulative distribution of the whole-bank supremum are deduced from a class of Bonferroni-type inequalities. The asymptotic properties of the cumulative distribution are obtained, in the limit where the number of correlators goes to infinity. The validity of numerical simulations made on small-size banks is extended to banks of any size, via a gaussian-correlation inequality. The result is used to estimate the optimum template density, yielding the best tradeoff between computational cost and detection efficiency, in terms of undetected potentially observable sources at a prescribed false-alarm level, for the simplest case of Newtonian chirps.Comment: submitted to Phys. Rev.

Risk profiling and efficacy of albendazole against the hookworms Necator americanus and Ancylostoma ceylanicum in Cambodia to support control programs in Southeast Asia and the Western Pacific

Background: Hookworm disease is endemic throughout many parts of the Asia Pacific, despite targeted control programs of at-risk populations. The success of these programs has been hindered by the limited efficacy of widely-used mebendazole, rapid re-infection rates linked to persistent reservoirs of untreated people and dogs, and the low sensitivity of conventional coprodiagnostic techniques employed. Methods: Here, we used standard faecal flotation (SFF) and a multiplex qPCR (mqPCR) assay to calculate and compare species-specific cure and egg reduction rates of single dose albendazole (400 mg) against hookworm infections at community level. Data from a cross-sectional survey in 1,232 people from Cambodia were used to inform a generalised linear mixed model to identify risk factors linked to hookworm infection(s) at baseline. Furthermore, we calculated risk factors associated to the probability of being cured after albendazole administration. Findings: Overall, 13.5% of all 1,232 people tested by SFF were positive for hookworm infection(s). Most (80.1%) infected people were >12 years of age, hence above the age targeted by the WHO control program. We estimate that as age increases, the odds of being infected increases at a faster rate for females than for males. We revealed a substantial difference in cure rate of hookworm infection(s) following albendazole treatment using the SFF (81.5%) and mqPCR (46.4%) assays, and provide the first data on the efficacy of this drug against the zoonotic hookworm Ancylostoma ceylanicum. We estimated that as age increases by one year, the odds of being cured decreases by 0.4%-3.7%. Similarly, the odds of being cured for people who boiled drinking water was estimated to be between 1.02 and 6.82. Interpretation: These findings show that the adoption of refined diagnostic techniques is central to monitoring hookworm infection(s) and the success of control strategies, which can ultimately aid in reducing associated morbidity in human populations. The approach taken is likely to be directly applicable to other parts of Southeast Asia and the Western Pacific, where specific epidemiological conditions might hamper the success of targeted treatment programs. Funding: Faculty of Veterinary and Agricultural Sciences Strategic Research Funds, The University of Melbourne

Bayesian Blocks, A New Method to Analyze Structure in Photon Counting Data

I describe a new time-domain algorithm for detecting localized structures (bursts), revealing pulse shapes, and generally characterizing intensity variations. The input is raw counting data, in any of three forms: time-tagged photon events (TTE), binned counts, or time-to-spill (TTS) data. The output is the most likely segmentation of the observation into time intervals during which the photon arrival rate is perceptibly constant -- i.e. has a fixed intensity without statistically significant variations. Since the analysis is based on Bayesian statistics, I call the resulting structures Bayesian Blocks. Unlike most, this method does not stipulate time bins -- instead the data themselves determine a piecewise constant representation. Therefore the analysis procedure itself does not impose a lower limit to the time scale on which variability can be detected. Locations, amplitudes, and rise and decay times of pulses within a time series can be estimated, independent of any pulse-shape model -- but only if they do not overlap too much, as deconvolution is not incorporated. The Bayesian Blocks method is demonstrated by analyzing pulse structure in BATSE $\gamma$-ray data. The MatLab scripts and sample data can be found on the WWW at: http://george.arc.nasa.gov/~scargle/papers.htmlComment: 42 pages, 2 figures; revision correcting mathematical errors; clarifications; removed Cyg X-1 sectio

The evolution of the terrestrial-terminating Irish Sea glacier during the last glaciation

Here we reconstruct the last advance to maximum limits and retreat of the Irish Sea Glacier (ISG), the only land‐terminating ice lobe of the western British Irish Ice Sheet. A series of reverse bedrock slopes rendered proglacial lakes endemic, forming time‐transgressive moraine‐ and bedrock‐dammed basins that evolved with ice marginal retreat. Combining, for the first time on glacial sediments, optically stimulated luminescence (OSL) bleaching profiles for cobbles with single grain and small aliquot OSL measurements on sands, has produced a coherent chronology from these heterogeneously bleached samples. This chronology constrains what is globally an early build‐up of ice during late Marine Isotope Stage 3 and Greenland Stadial (GS) 5, with ice margins reaching south Lancashire by 30 ± 1.2 ka, followed by a 120‐km advance at 28.3 ± 1.4 ka reaching its 26.5 ± 1.1 ka maximum extent during GS‐3. Early retreat during GS‐3 reflects piracy of ice sources shared with the Irish‐Sea Ice Stream (ISIS), starving the ISG. With ISG retreat, an opportunistic readvance of Welsh ice during GS‐2 rode over the ISG moraines occupying the space vacated, with ice margins oscillating within a substantial glacial over‐deepening. Our geomorphological chronosequence shows a glacial system forced by climate but mediated by piracy of ice sources shared with the ISIS, changing flow regimes and fronting environments

A procedure for the change point problem in parametric models based on phi-divergence test-statistics

This paper studies the change point problem for a general parametric, univariate or multivariate family of distributions. An information theoretic procedure is developed which is based on general divergence measures for testing the hypothesis of the existence of a change. For comparing the accuracy of the new test-statistic a simulation study is performed for the special case of a univariate discrete model. Finally, the procedure proposed in this paper is illustrated through a classical change-point example

Tensor Regression with Applications in Neuroimaging Data Analysis

Classical regression methods treat covariates as a vector and estimate a corresponding vector of regression coefficients. Modern applications in medical imaging generate covariates of more complex form such as multidimensional arrays (tensors). Traditional statistical and computational methods are proving insufficient for analysis of these high-throughput data due to their ultrahigh dimensionality as well as complex structure. In this article, we propose a new family of tensor regression models that efficiently exploit the special structure of tensor covariates. Under this framework, ultrahigh dimensionality is reduced to a manageable level, resulting in efficient estimation and prediction. A fast and highly scalable estimation algorithm is proposed for maximum likelihood estimation and its associated asymptotic properties are studied. Effectiveness of the new methods is demonstrated on both synthetic and real MRI imaging data.Comment: 27 pages, 4 figure

Data analysis of gravitational-wave signals from spinning neutron stars. III. Detection statistics and computational requirements

We develop the analytic and numerical tools for data analysis of the gravitational-wave signals from spinning neutron stars for ground-based laser interferometric detectors. We study in detail the statistical properties of the optimum functional that need to be calculated in order to detect the gravitational-wave signal from a spinning neutron star and estimate its parameters. We derive formulae for false alarm and detection probabilities both for the optimal and the suboptimal filters. We assess the computational requirements needed to do the signal search. We compare a number of criteria to build sufficiently accurate templates for our data analysis scheme. We verify the validity of our concepts and formulae by means of the Monte Carlo simulations. We present algorithms by which one can estimate the parameters of the continuous signals accurately.Comment: LaTeX, 45 pages, 13 figures, submitted to Phys. Rev.