Python custom_sources Examples

Example #1

# t_end = 56474.00 #true end hard flare

# Ra and dec of source
ra = 160.134167
dec = 6.29

# Distance to source -> z=0.048 according to TeVCat. With lambdaCDM, this gives:
z = 0.73
lumdist = Distance(z=z).to("Mpc").value

# Creates the .npy source catalogue
catalogue = custom_sources(
    name="GB6_J1040_0617",
    ra=ra,
    dec=dec,
    weight=1.0,
    distance=lumdist,
    start_time=t_start,
    end_time=t_end,
)

cat_path = transients_dir + "GB6_J1040_0617.npy"
np.save(cat_path, catalogue)

search_window = float(t_end - t_start)
max_window = 150.0
# Initialise Injectors/LLHs

injection_energy = {
    "energy_pdf_name": "PowerLaw",
    "gamma": 2.0,

Example #2

File: make_txs_catalogue.py Project: afedynitch/flarestack

import numpy as np
import os
from flarestack.shared import transients_dir

# Start and end time of neutrino flare, taken from box fit in
# https://arxiv.org/abs/1807.08794.

t_start = 56937.81
t_end = 57096.21

# Ra and dec of source, from Science paper (https://arxiv.org/abs/1807.08794)
ra = 77.3582
dec = 5.69314

# Distance to source, according to https://arxiv.org/abs/1802.01939, is 0.3365
z = 0.3365
lumdist = Distance(z=z).to("Mpc").value

# Creates the .npy source catalogue
txs_catalogue = custom_sources(name="TXS_0506+056",
                               ra=ra,
                               dec=dec,
                               weight=1.,
                               distance=lumdist,
                               start_time=t_start,
                               end_time=t_end,
                               ref_time=t_start)

txs_cat_path = transients_dir + "/TXS_0506+056.npy"
np.save(txs_cat_path, txs_catalogue)

Example #3

File: benchmark_blazar_flare_search.py Project: icecube/flarestack

t_end = 57595.00

# Ra and dec of source
ra = 300.00
dec = 65.15

# Distance to source -> z=0.048 according to TeVCat. With lambdaCDM, this gives:
z = 0.048
lumdist = Distance(z=z).to("Mpc").value

# Creates the .npy source catalogue
catalogue = custom_sources(
    name="1ES_1959+650",
    ra=ra,
    dec=dec,
    weight=1.0,
    distance=lumdist,
    start_time=t_start,
    end_time=t_end,
)

cat_path = transients_dir + "1ES_1959+650.npy"
np.save(cat_path, catalogue)

max_window = float(t_end - t_start)

# Initialise Injectors/LLHs

injection_energy = {
    "Name": "Power Law",
    "Gamma": 2.0,

Example #4

File: make_AT2018cow_catalogue.py Project: afedynitch/flarestack

t_start = ref_time - 30.
t_end = ref_time + 100.

# Ra and dec of source, from TNS (https://wis-tns.weizmann.ac.il/object/2018cow)
ra = 244.000927647
dec = 22.2680094118

# Distance to source,from http://www.astronomerstelegram.org/?read=11727
z = 0.014145
lumdist = Distance(z=z).to("Mpc").value

# Creates the .npy source catalogue
at2018cow_catalogue = custom_sources(name="AT2018cow",
                                     ra=ra,
                                     dec=dec,
                                     weight=1.,
                                     distance=lumdist,
                                     start_time=t_start,
                                     end_time=t_end,
                                     ref_time=t_start)

at2018cow_cat_path = individual_tde_cat("AT2018cow")

try:
    os.makedirs(os.path.dirname(at2018cow_cat_path))
except OSError:
    pass

np.save(at2018cow_cat_path, at2018cow_catalogue)