# 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,
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)
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,
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)