A Compound Poisson Generator approach to Point-Source Inference in Astrophysics

Abstract

The identification and description of point sources is one of the oldest problems in astronomy; yet, even today the correct statistical treatment for point sources remains one of the field's hardest problems. For dim or crowded sources, likelihood based inference methods are required to estimate the uncertainty on the characteristics of the source population. In this work, a new parametric likelihood is constructed for this problem using Compound Poisson Generator (CPG) functionals which incorporate instrumental effects from first principles. We demonstrate that the CPG approach exhibits a number of advantages over Non-Poissonian Template Fitting (NPTF) - an existing method - in a series of test scenarios in the context of X-ray astronomy. These demonstrations show that the effect of the point-spread function, effective area, and choice of point-source spatial distribution cannot, generally, be factorised as they are in NPTF, while the new CPG construction is validated in these scenarios. Separately, an examination of the diffuse-flux emission limit is used to show that most simple choices of priors on the standard parameterisation of the population model can result in unexpected biases: when a model comprising both a point-source population and diffuse component is applied to this limit, nearly all observed flux will be assigned to either the population or to the diffuse component. A new parametrisation is presented for these priors which properly estimates the uncertainties in this limit. In this choice of priors, CPG correctly identifies that the fraction of flux assigned to the population model cannot be constrained by the data.

arXiv comment30+10 pages, 15 figures, comments welcome; v2 updated to published version. Accompanying code available at https://github.com/ghcollin/cpg_likelihood

References in corpus (46)

• emcee: The MCMC Hammer
• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• The Nuclear Spectroscopic Telescope Array (NuSTAR) Mission
• Dark Matter Annihilation in The Galactic Center As Seen by the Fermi Gamma Ray Space Telescope
• On The Origin Of The Gamma Rays From The Galactic Center
• The Characterization of the Gamma-Ray Signal from the Central Milky Way: A Compelling Case for Annihilating Dark Matter
• Detection of a Gamma-Ray Source in the Galactic Center Consistent with Extended Emission from Dark Matter Annihilation and Concentrated Astrophysical Emission
• Background model systematics for the Fermi GeV excess
• Dark Matter and Pulsar Model Constraints from Galactic Center Fermi-LAT Gamma Ray Observations
• Fermi-LAT Observations of High-Energy Gamma-Ray Emission Toward the Galactic Center
• Astrophysical and Dark Matter Interpretations of Extended Gamma-Ray Emission from the Galactic Center
• Strong support for the millisecond pulsar origin of the Galactic center GeV excess
• X-ray Emissions from Accreting White Dwarfs: a Review
• Evidence for Unresolved Gamma-Ray Point Sources in the Inner Galaxy
• Calibration of the NuSTAR High Energy Focusing X-ray Telescope
• NuSTAR Observations of the Bullet Cluster: Constraints on Inverse Compton Emission
• Galactic Bulge Preferred Over Dark Matter for the Galactic Center Gamma-Ray Excess
• Dark Matter Strikes Back at the Galactic Center
• Statistics of gamma-ray point sources below the Fermi detection limit
• Foreground Mismodeling and the Point Source Explanation of the Fermi Galactic Center Excess
• The High-Energy Tail of the Galactic Center Gamma-Ray Excess
• Spurious Point Source Signals in the Galactic Center Excess
• NuSTAR Hard X-ray Survey of the Galactic Center Region II: X-ray Point Sources
• Faint Source Counts from Off-source Fluctuation Analysis on Chandra Observations of the Hubble Deep Field-North
• The Gamma-Ray-Flux Probability Distribution Function from Galactic Halo Substructure
• Deciphering Contributions to the Extragalactic Gamma-Ray Background from 2 GeV to 2 TeV
• The Enigmatic Galactic Center Excess: Spurious Point Sources and Signal Mismodeling
• A New Mask for An Old Suspect: Testing the Sensitivity of the Galactic Center Excess to the Point Source Mask
• Unveiling the Gamma-ray Source Count Distribution Below the Fermi Detection Limit with Photon Statistics
• Distinguishing Dark Matter from Unresolved Point Sources in the Inner Galaxy with Photon Statistics
• Characterizing the Nature of the Unresolved Point Sources in the Galactic Center
• The GCE in a New Light: Disentangling the $\gamma$-ray Sky with Bayesian Graph Convolutional Neural Networks
• Statistical Measurement of the Gamma-ray Source-count Distribution as a Function of Energy
• Discovery of a New Galactic Center Excess Consistent with Upscattered Starlight
• Dissecting the inner Galaxy with $\gamma$-ray pixel count statistics
• NPTFit: A code package for Non-Poissonian Template Fitting
• Diffusive Synchrotron Radiation from Pulsar Wind Nebulae
• NuSTAR measurement of the cosmic X-ray background in the 3-20 keV energy band
• One-point fluctuation analysis of the high-energy neutrino sky
• Testing gamma-ray models of blazars in the extragalactic sky
• Improved Point Source Detection in Crowded Fields using Probabilistic Cataloging
• Constraining Galactic dark matter with gamma-ray pixel counts statistics
• Inference of Unresolved Point Sources At High Galactic Latitudes Using Probabilistic Catalogs
• Modelling the flux distribution function of the extragalactic gamma-ray background from dark matter annihilation
• Dark matter substructure cannot explain properties of the Fermi Galactic Centre excess
• A Search for Neutrino Point-Source Populations in 7 Years of IceCube Data with Neutrino-count Statistics

Cited by in corpus (6)

• The Return of the Templates: Revisiting the Galactic Center Excess with Multi-Messenger Observations
• Millisecond Pulsars from Accretion Induced Collapse as the Origin of the Galactic Centre Gamma-ray Excess Signal
• Extracting the gamma-ray source-count distribution below the Fermi-LAT detection limit with deep learning
• PCAT-DE: Reconstructing point-like and diffuse signals in astronomical images using spatial and spectral information
• Characterizing the Expected Behavior of Non-Poissonian Template Fitting
• On the Energy Distribution of the Galactic Center Excess' Sources