Tutorial

Guide to UHECR Public Data

The era of open science has transformed cosmic ray research. Major observatories now release their data publicly, enabling anyone — students, citizen scientists, independent researchers — to explore the universe's most energetic particles. Here's your complete guide to finding and using this data.

Pierre Auger Open Data

The Pierre Auger Observatory released its first public dataset in 2021, marking a milestone for open science in astroparticle physics.

What's Available

Website: opendata.auger.org

Data Format

Cosmic ray events are provided in JSON format with fields including:

📄 Auger Event Fields

  • timestamp: GPS time of event (seconds since Jan 6, 1980)
  • energy: Reconstructed energy in EeV
  • theta, phi: Zenith and azimuth angles (degrees)
  • ra, dec: Right ascension and declination (degrees)
  • gal_lat, gal_lon: Galactic coordinates (degrees)
  • sd_stations: Number of triggered surface detectors

Getting Started

1. Visit opendata.auger.org
2. Navigate to "Cosmic ray data" → "10% data release"
3. Download JSON files directly or use the provided Jupyter notebooks
4. The site includes Python tutorials for basic analysis

Analysis Ideas

GWOSC: Gravitational Wave Data

The Gravitational Wave Open Science Center provides all LIGO/Virgo/KAGRA data for multi-messenger studies.

What's Available

Website: gwosc.org

Key Catalogs

🌊 GWTC Event Catalogs

  • GWTC-1: 11 events from O1/O2 (2015-2017)
  • GWTC-2: 39 events from O3a (2019)
  • GWTC-2.1: 8 additional O3a events
  • GWTC-3: 35 events from O3b (2019-2020)
  • GWTC-4.0: 80+ events from O4 (ongoing)

For UHECR Correlation Studies

The most useful GWOSC products for UHECR research:

Download parameter estimation files (posterior samples) for detailed source properties.

Fermi GBM Burst Catalog

The Fermi Gamma-ray Burst Monitor has detected over 3,500 gamma-ray bursts since 2008.

What's Available

Website: HEASARC Fermi GBM Catalog

Accessing the Data

Use the HEASARC Browse interface to query the catalog. You can filter by date range, duration, position, and other parameters. Results can be downloaded in various formats (FITS, CSV, VOTable).

IceCube Data Releases

IceCube has released several datasets of high-energy neutrino events.

What's Available

Website: icecube.wisc.edu/science/data-releases

Multi-Messenger Potential

Neutrinos travel in straight lines and arrive without delay (like light). Finding neutrino-UHECR correlations would strongly constrain common sources.

KASCADE Cosmic Ray Data Centre

KASCADE and KASCADE-Grande operated in Germany from 1996-2013, detecting cosmic rays from 10¹⁵ to 10¹⁸ eV.

What's Available

Website: kcdc.iap.kit.edu

Energy Range

KASCADE covers lower energies than Auger — the "knee" region. Useful for understanding the galactic-extragalactic transition and calibrating composition methods.

Zenodo: General Repository

Many research groups deposit supplementary data on Zenodo, CERN's open data repository.

Website: zenodo.org

Search for: "cosmic rays," "UHECR," "Pierre Auger," "Telescope Array," etc. You'll find simulation datasets, analysis code, and supplementary materials from published papers.

Tools for Analysis

Python Libraries

Coordinate Systems

UHECR analysis requires fluency with multiple coordinate systems:

Astropy's SkyCoord class handles conversions elegantly.

Example Code

import json
import numpy as np
from astropy.coordinates import SkyCoord
import astropy.units as u

# Load Auger data
with open('auger_10percent.json') as f:
    events = json.load(f)

# Convert to SkyCoord
coords = SkyCoord(
    ra=[e['ra'] for e in events] * u.deg,
    dec=[e['dec'] for e in events] * u.deg,
    frame='icrs'
)

# Transform to Galactic
gal = coords.galactic
print(f"Galactic l: {gal.l}")
print(f"Galactic b: {gal.b}")

# Filter by energy
high_e = [e for e in events if e['energy'] > 50]
print(f"Events above 50 EeV: {len(high_e)}")
            

Tips for Success

Start Simple

Begin by reproducing known results — the energy spectrum, sky maps, basic anisotropy searches. This validates your analysis pipeline before attempting novel studies.

Understand Systematics

Public datasets have limitations:

The Auger open data documentation discusses these carefully — read it!

Use the Community

The Power of Open Data

Public cosmic ray data enables:

The universe's secrets are hidden in this data. The tools to find them are available to everyone.

Continue Reading