Bayesian posterior for IceCube 7year pointsource search with neutrinocount statistics (released 4 Sep 2019)
Click anywhere to dismiss this message.
Bayesian posterior for IceCube 7year pointsource search with neutrinocount statistics (preview)
Download Data
We ask for users who wish to download our data to fill out the information below. We will use this information to contact you in the event that a revised data file is released due to errors or significant improvements. Your personal information will never be released for any other purposes.
Introduction
The presence of a population of point sources in a dataset modifies the underlying neutrinocount statistics from the Poisson distribution. This deviation can be exactly quantified using the nonPoissonian template fitting technique, and in this work we present the first application of this approach to the IceCube highenergy neutrino dataset. Using this method, we search in 7 years of IceCube data for pointsource populations correlated with the disk of the Milky Way, the Fermi bubbles, the Schlegel, Finkbeiner, and Davis dust map, or with the isotropic extragalactic sky. In addition, a northernskyonly spatial template is also considered. In order to enhance the flexibility of the results, we publish the full posterior from our analysis, which can be used to establish limits on specific population models that would contribute to the observed IceCube neutrino flux.
A full description of the method can be found in http://arxiv.org/abs/1909.08623.
Data release
Suggested citation for this dataset:
IceCube Collaboration (2019): Bayesian posterior for IceCube 7year pointsource search with neutrinocount statistics. DOI:10.21234/dtrse557
Click here to download (.zip, 12MB)
Included in the download are the following files:
Data files
This release contains the posterior and a short example, written in Python, on how to load and use the posterior. We emphasize that our analysis and hence these posteriors are constructed with the assumption that the astrophysical population produces neutrinos with an E^(2) energy spectrum.

posterior.h5: Within this HDF5 formatted file, five tables named Isotropic, Galactic_disk, Fermi_bubble, SFD_dust, and Northern_sky contain the posterior samples for their respective templates.
Each table describes equally weighted samples using five columns. Four columns, labeled ln_A, ln_Fb, n1, and n2, contain the coordinates for the sample in natural logarithmic parameter space for the differential sourcecount function normalization A and break F_b, while the power indices n_1 and n_2 are in linear space. The fifth column — labeled loglikelihood — gives the natural logarithm of the likelihood function at the location of the corresponding sample. In addition, each table has two attributes named P_log_evidence and NP_log_evidence that contain the natural logarithm of the evidence integral for the Poissonian model and nonPoissonian model, respectively.
The root node of the HDF5 file also contains a series of attributes named units_ln_A, units_ln_Fb, units_n1, and units_n2 that specify the units that the posterior sample coordinates are given in. Another series of root attributes named prior_ln_A, prior_ln_Fb, prior_n1, and prior_n2 give the probability density of the uniform priors for each of the model parameters.  example.py: An example program that loads the posterior samples, computes the Bayes factor at specific model locations in parameter space, and generates a reproduction of Figure 10 in the accompanying paper. This program requires the numpy, scipy, matplotlib and pytables (named tables in pip) python libraries to be installed.
For any questions about this data release, please write to data@icecube.wisc.edu