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]
)
