IceCube HESE 12-year data release

This data release describes a revisited analysis of the IceCube high-energy starting events (HESE) sample with an additional 4.5 years of data that employs a full-resimulation reconstruction approach, DirectFit, which exactly incorporates recent updates in glacial ice modeling. These updates include a microscopic description of the observed ice anisotropy due to birefringence and a full description of ice layer undulations in three dimensions. Each event is reconstructed with the Approximate Bayesian Computation method, with the posterior sample points then fitted to an asymmetric directional probability density function. For further details, refer to the IceCube publication PoS-ICRC23-1030 / arXiv:2307.13878 Note on interpretation Reconstructed quantities for subsets of these events have been published previously. IceCube utilizes several event reconstructions, each with different strengths and employed in different contexts. DirectFit, the reconstruction in this data release, performs accurate event reconstruction, but is too computationally expensive to run on large scale Monte Carlo (MC). Therefore this data release should not be employed in a forward-folding likelihood fit for MC-based spectral analyses. Instead, as the reconstruction utilizes a complete description of the ice, the dataset is meant to provide accurate directional constraints for HESE events and for use in neutrino source searches. We highlight the differences between the reconstructions listed here, and those available elsewhere for subsets of this sample, below. For overlapping events already present in the previous 7.5 year dataset (arxiv:2011.03545), the reconstructions provided here may be compared to those in Appendix J. The reconstructions in Appendix J were obtained with DirectFit, but using an older ice model. Differences primarily arise from the updates to ice optical modeling and/or changes in the reconstruction routine (shower vs track) applied for select events. As the new ice model leads to a generally improved data-MC agreement in both calibration and physics data, the resulting fits are as a result preferred. The per-event differences between the fits may allow to conservatively gauge the impact of systematic detector uncertainties. For overlapping events already present in the previous 7.5 year dataset (arxiv:2011.03545), the reconstructions provided here utilize a different algorithm from those in Appendix F and listed as part of that data release. Specifically, the energies and directions published here should not be combined with the spectral fitting tools accompanying the 7.5 year dataset. Diffuse spectral analysis rely on likelihood-based reconstructions (arXiv:1311.4767) that leverage tabulated light yield approximations. Such algorithms are substantially more computationally efficient than DirectFit and are fast enough to be applied to large scale Monte Carlo, as required for forward-folding fits, at a cost of reduced accuracy. For overlapping events already present in the IceCat-1 catalog of alert tracks and released in real time via GCN/AMON, the visible energy reported here should not be compared to arxiv:2304.01174 (Tab. 3) and published in the accompanying data release . The IceCat-1 data release reports the most likely neutrino energy assuming a power-law spectral index of -2.19 at a given declination. The energies reported here correspond to the visible energy detected by IceCube. Note on event times The Modified Julian Date (MJD) is provided here for each event. In previously published HESE datasets that include the MJD, events that occurred between 2012-06-30 and 2015-05-17 were one second behind UTC time due to a master clock firmware issue. This offset has been corrected for in this data release; the overlapping events that fall within those dates will have a one second difference when compared to previous data releases. For those cases the MJDs provided here should be preferred. Data release Provided in this release is data.npy, which is a numpy structured array that contains the 164 reconstructed events that pass the HESE selection. The structured array contains the following fields: ('id', '< i4') ('mjd', '< f8') ('ra', '< f4') ('dec', '< f4') ('params', [('f0', '< f8'), ('f1', '< f8'), ('f2', '< f8'), ('f3', '< f8'), ('f4', '< f8'), ('f5', '< f8'), ('f6', '< f8'), ('f7', '< f8')]) ('reconstruction', 'S10') ('energy', '< f4') ('drlogl', '< f4') For each event, id corresponds to the event-ID consistent with prior HESE data releases, mjd is the Modified Julian Date when the event was observed, ra and dec correspond to the most-probable arrival direction, params corresponds to the eight parameters that describe the directional probability density function (PDF) using a generalized Fisher-Bingham distribution, reconstruction corresponds to the algorithm used by DirectFit, and energy the most-probable visible energy. For more details on the directional PDF construction see this paper....