Local to Equatorial Coordinates

Where do my neutrinos come from?

__author__ = "moritz"

import numpy as np
import pandas as pd

from km3astro.coord import local_event, Sun, neutrino_to_source_direction

Detector Coordinates

Let’s define some random events.

theta = np.array([10, 45, 70, 23, 20, 11, 24, 54]) * np.pi / 180
phi = np.array([4, 23, 200, 320, 10, 45, 29, 140]) * np.pi / 180
time = pd.to_datetime(
    [
        "2015-01-12T15:10:12",
        "2015-06-12T13:48:56",
        "2015-03-09T21:57:52",
        "2015-03-15T14:24:01",
        "2015-01-12T15:10:12",
        "2015-06-12T13:48:56",
        "2015-03-09T21:57:52",
        "2015-03-15T14:24:01",
    ]
).values

Note:

Phi, theta: Where the neutrino is pointing to

Zenith, azimuth: where the neutrino is coming from

Create event in local coordinates (aka AltAz or Horizontal Coordinates)

This returns an astropy.SkyCoord instance.

evt_local = local_event(theta, phi, time, location="orca")

sun = Sun(time)

sep = evt_local.separation(sun)

print(sep)

Out:

[91d55m41.63665769s 101d11m43.36401381s 87d52m39.96833153s
 115d30m09.9465609s 82d41m44.34697762s 135d50m15.92986853s
 46d42m34.01756926s 116d48m28.00919023s]