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])
def ic86_new_season(year): return IceCubeSeason( season_name="IC86_{0}".format(year), sample_name=sample_name, exp_path=ps_data_dir + "IC86_{0}_exp.npy".format(year), mc_path=ps_data_dir + "IC86_2012_MC.npy", grl_path=grl_data_dir + "IC86_{0}_exp.npy".format(year), sin_dec_bins=get_ps_binning(f"IC86_{year}")[0], log_e_bins=get_ps_binning(f"IC86_{year}")[1])
old_seasons = ["IC40", "IC59", "IC79", "IC86_2011"] for season in old_seasons: ps_v003_p01.add_season(old_ic_season(season)) # Add in combined IC86 2012-2017 seasons new_years = ["2012", "2013", "2014", "2015", "2016", "2017"] ic86_234567 = IceCubeSeason( season_name="IC86_2012_17", sample_name=sample_name, exp_path=[ps_data_dir + "IC86_{0}_exp.npy".format(x) for x in new_years], mc_path=ps_data_dir + "IC86_2012_MC.npy", grl_path=[grl_data_dir + "IC86_{0}_exp.npy".format(x) for x in new_years], sin_dec_bins=get_ps_binning("IC86_2012")[0], log_e_bins=get_ps_binning("IC86_2012")[1]) ps_v003_p01.add_season(ic86_234567) # Add in each new season as an optional subseason def ic86_new_season(year): return IceCubeSeason( season_name="IC86_{0}".format(year), sample_name=sample_name, exp_path=ps_data_dir + "IC86_{0}_exp.npy".format(year), mc_path=ps_data_dir + "IC86_2012_MC.npy", grl_path=grl_data_dir + "IC86_{0}_exp.npy".format(year), sin_dec_bins=get_ps_binning(f"IC86_{year}")[0],
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) ic59 = IceCubeSeason( season_name="IC59", sample_name=sample_name, exp_path=ps_data_dir + "IC59_exp.npy", mc_path=ps_data_dir + "IC59_corrected_MC.npy", grl_path=ps_data_dir + "IC59_GRL.npy", sin_dec_bins=get_ps_binning("IC59")[0], log_e_bins=get_ps_binning("IC59")[1] )