Change logs: Version 2018/07/25 * initial version """ from flarestack.data.icecube.ic_season import IceCubeDataset, get_dataset_dir from flarestack.data.icecube.northern_tracks import NTSeason, get_diffuse_binning icecube_dataset_dir = get_dataset_dir() nt_data_dir = icecube_dataset_dir + "northern_tracks/version-002-p05/" nt_v002_p05 = IceCubeDataset() sample_name = "northern_tracks_v002_p05" def generate_diffuse_season(name): season = NTSeason( season_name=name, sample_name=sample_name, exp_path=nt_data_dir + "dataset_8yr_fit_{0}_exp_compressed.npy".format(name), mc_path=nt_data_dir + "dataset_8yr_fit_{0}_MC_compressed.npy".format(name), grl_path=nt_data_dir + "GRL/dataset_8yr_fit_{0}_exp_compressed.npy".format(name), sin_dec_bins=get_diffuse_binning(name)[0],
"""PS Tracks v003_p01, as used by Tessa in the 10 year PS analysis. """ from flarestack.data.icecube.ic_season import IceCubeSeason, \ IceCubeDataset, icecube_dataset_dir from flarestack.data.icecube.ps_tracks import get_ps_binning import numpy as np ps_data_dir = icecube_dataset_dir + "ps_tracks/version-003-p01/" grl_data_dir = ps_data_dir + "GRL/" ps_v003_p01 = IceCubeDataset() sample_name = "ps_tracks_v003_p01" # Add in old seasons before full detector operation, and IC86_1 def old_ic_season(season): return IceCubeSeason(season_name=season, sample_name=sample_name, exp_path=ps_data_dir + "{0}_exp.npy".format(season), mc_path=ps_data_dir + "{0}_MC.npy".format(season), grl_path=grl_data_dir + "{0}_exp.npy".format(season), sin_dec_bins=get_ps_binning(season)[0], log_e_bins=get_ps_binning(season)[1]) old_seasons = ["IC40", "IC59", "IC79", "IC86_2011"] for season in old_seasons: ps_v003_p01.add_season(old_ic_season(season))
gfu_data_dir = icecube_dataset_dir + "gfu/version-002-p04/" grl_dir = gfu_data_dir + "GRL/" def get_exp_path(year): return gfu_data_dir + "IC86_{0}_data.npy".format(year) def get_grl_path(year): return grl_dir + "IC86_{0}_data.npy".format(year) all_years = ["2011", "2012", "2013", "2014", "2015", "2016", "2017", "2018"] gfu_v002_p04 = IceCubeDataset() sample_name = "gfu_v002_p04" gfu_8_year = IceCubeSeason( season_name="gfu_8_year", sample_name=sample_name, exp_path=[get_exp_path(x) for x in all_years], mc_path=gfu_data_dir + "IC86_2011_MC.npy", grl_path=[get_grl_path(x) for x in all_years], sin_dec_bins=gfu_binning[0], log_e_bins=gfu_binning[1], ) gfu_v002_p04.add_season(gfu_8_year)
/data/ana .. /current, with the following readme: * better pull-correction using splines, instead of polynomial fits This is the sample that was used for the original TXS analysis in 2017 (NOT THE ONE THAT INCLUDED ADDITIONAL DATA UP TO SPRING 2018!) """ from flarestack.data.icecube.ps_tracks.ps_v002_p01 import ps_v002_p01 from flarestack.data.icecube.ic_season import IceCubeDataset, \ IceCubeSeason, icecube_dataset_dir from flarestack.data.icecube.gfu import gfu_binning import numpy as np gfu_data_dir = icecube_dataset_dir + "gfu/version-002-p01/" gfu_v002_p01 = IceCubeDataset() sample_name = "gfu_v002_p01" gfu_season = IceCubeSeason( season_name=sample_name, sample_name=sample_name, exp_path=[ gfu_data_dir + "SplineMPEmax.MuEx.IC86-2015.npy", gfu_data_dir + "SplineMPEmax.MuEx.IC86-2016.npy", gfu_data_dir + "SplineMPEmax.MuEx.IC86-2017.npy" ], mc_path=gfu_data_dir + "SplineMPEmax.MuEx.MC.npy", grl_path=gfu_data_dir + "SplineMPEmax.MuEx.GRL.npy", sin_dec_bins=gfu_binning[0], log_e_bins=gfu_binning[1]
MC File IC86-2013_corrected_MC_v2.npy IC86, 2014: Data File IC86-2014_exp_v2.npy MC File IC86-2014_corrected_MC_v2.npy """ from flarestack.data.icecube.ic_season import IceCubeDataset, \ IceCubeSeason, icecube_dataset_dir from flarestack.data.icecube.ps_tracks import get_ps_binning import numpy as np import copy ps_data_dir = icecube_dataset_dir + "ps_tracks/version-002-p01/" ps_v002_p01 = IceCubeDataset() sample_name = "ps_tracks_v002_p01" ic40 = IceCubeSeason( season_name="IC40", sample_name=sample_name, exp_path=ps_data_dir + "IC40_exp.npy", mc_path=ps_data_dir + "IC40_corrected_MC.npy", grl_path=ps_data_dir + "IC40_GRL.npy", sin_dec_bins=get_ps_binning("IC40")[0], log_e_bins=get_ps_binning("IC40")[1] ) ps_v002_p01.add_season(ic40)
""" from flarestack.data.icecube.ps_tracks.ps_v002_p01 import ps_v002_p01 from flarestack.data.icecube.ic_season import ( IceCubeDataset, IceCubeSeason, get_dataset_dir, ) from flarestack.data.icecube.gfu import gfu_binning import numpy as np icecube_dataset_dir = get_dataset_dir() gfu_data_dir = icecube_dataset_dir + "gfu/version-002-p02/" gfu_v002_p02 = IceCubeDataset() sample_name = "gfu_v002_p02" gfu_season = IceCubeSeason( season_name=sample_name, sample_name=sample_name, exp_path=[ gfu_data_dir + "SplineMPEmax.MuEx.IC86-2015.npy", gfu_data_dir + "SplineMPEmax.MuEx.IC86-2016.npy", gfu_data_dir + "SplineMPEmax.MuEx.IC86-2017.npy", ], mc_path=gfu_data_dir + "SplineMPEmax.MuEx.MC.npy", grl_path=gfu_data_dir + "SplineMPEmax.MuEx.GRL.npy", sin_dec_bins=gfu_binning[0], log_e_bins=gfu_binning[1],
from flarestack.data.icecube.ic_season import IceCubeDataset, get_dataset_dir from flarestack.data.icecube.northern_tracks import ( NTSeasonNewStyle, get_diffuse_binning, ) icecube_dataset_dir = get_dataset_dir() nt_data_dir = icecube_dataset_dir + "northern_tracks/version-005-p00/" nt_v005_p00 = IceCubeDataset() sample_name = "northern_tracks_v005_p00" def generate_diffuse_season(name): season = NTSeasonNewStyle( season_name=name, sample_name=sample_name, exp_path=nt_data_dir + f"{name}_exp.npy", mc_path=nt_data_dir + "IC86_pass2_MC.npy", grl_path=nt_data_dir + f"GRL/{name}_exp.npy", sin_dec_bins=get_diffuse_binning(name)[0], log_e_bins=get_diffuse_binning(name)[1], ) nt_v005_p00.add_season(season) seasons = [f"IC86_201{i}" for i in [1, 2, 3, 4, 5, 6, 7, 8, 9]] for season in seasons:
dataset folder and the presentations here: https://events.icecube.wisc.edu/event/125/contributions/7204/attachments/5554/6427/200915_pstracks_updates.pdf https://drive.google.com/file/d/1E1I8kgRmlWdLdXX_p7UweXn7gJIYx9Mg/view?usp=sharing """ import numpy as np import copy from flarestack.data.icecube.ic_season import IceCubeDataset, \ IceCubeSeason, icecube_dataset_dir from flarestack.data.icecube.ps_tracks import get_ps_binning ps_data_dir = icecube_dataset_dir + "ps_tracks/version-004-p00/" grl_data_dir = ps_data_dir + "GRL/" ps_v004_p00 = IceCubeDataset() sample_name = "ps_tracks_v004_p00" ########################## # START ADDING SEASONS # ########################## ############ # - IC40 - # ############ sinDec_bins = np.unique( np.concatenate([ np.linspace(-1., -0.25, 10 + 1), np.linspace(-0.25, 0.0, 10 + 1),
from flarestack.data.icecube.ic_season import ( IceCubeDataset, IceCubeSeason, get_dataset_dir, ) from flarestack.data.icecube.ps_tracks import get_ps_binning import numpy as np import logging logger = logging.getLogger(__name__) icecube_dataset_dir = get_dataset_dir() ps_data_dir = icecube_dataset_dir + "ps_tracks/version-002-p03/" grl_data_dir = ps_data_dir + "GRL/" ps_v002_p03 = IceCubeDataset() sample_name = "ps_tracks_v002_p03" logger.debug(f"building {sample_name}") logger.debug(f"adding IC40") ic40 = IceCubeSeason( season_name="IC40", sample_name=sample_name, exp_path=ps_data_dir + "IC40_exp.npy", mc_path=ps_data_dir + "IC40_MC.npy", grl_path=grl_data_dir + "IC40_exp.npy", sin_dec_bins=get_ps_binning("IC40")[0], log_e_bins=get_ps_binning("IC40")[1], ) ps_v002_p03.add_season(ic40)