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],
log_e_bins=get_ps_binning(f"IC86_{year}")[1])
for year in new_years:
ps_v003_p01.add_subseason(ic86_new_season(year))
ps_v002_p01.add_season(ic86_1)
# Add optional subseasons for IC86 2, 3, and 4, that can be called instead of
# the combined season
for i in range(2, 5):
ic86_i = IceCubeSeason(
season_name="IC86_{0}".format(i),
sample_name=sample_name,
exp_path=ps_data_dir + "IC86-201{0}_exp_v2.npy".format(i),
mc_path=ps_data_dir + "IC86-2012_corrected_MC_v2.npy",
grl_path=ps_data_dir + "IC86-201{0}_GRL.npy".format(i),
sin_dec_bins=get_ps_binning("IC86_2012")[0],
log_e_bins=get_ps_binning("IC86_2012")[1]
)
ps_v002_p01.add_subseason(ic86_i)
ic86_234 = IceCubeSeason(
season_name="IC86_234",
sample_name=sample_name,
exp_path=[
ps_data_dir + "IC86-2012_exp_v2.npy",
ps_data_dir + "IC86-2013_exp_v2.npy",
ps_data_dir + "IC86-2014_exp_v2.npy"
],
mc_path=ps_data_dir + "IC86-2012_corrected_MC_v2.npy",
grl_path=[
ps_data_dir + "IC86-2012_GRL.npy",
ps_data_dir + "IC86-2013_GRL.npy",
ps_data_dir + "IC86-2014_GRL.npy"
],
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)
for year in all_years:
gfu_year = IceCubeSeason(
season_name="IC86_{0}".format(year),
sample_name=sample_name,
exp_path=get_exp_path(year),
mc_path=gfu_data_dir + "IC86_2011_MC.npy",
grl_path=get_grl_path(year),
sin_dec_bins=gfu_binning[0],
log_e_bins=gfu_binning[1],
)
gfu_v002_p04.add_subseason(gfu_year)
energy_bins = np.arange(1., 9.5 + 0.01, 0.125)
mc_path = ps_data_dir + "IC86_2016_MC.npy"
# full detector configurations have a unified processing now
# so it is not necessary anymore to treat IC86, 2011 differently
for i in range(1, 10):
ic86_i = IceCubeSeason(
season_name="IC86_{0}".format(i),
sample_name=sample_name,
exp_path=ps_data_dir + "IC86_201{0}_exp.npy".format(i),
mc_path=mc_path,
grl_path=grl_data_dir + "IC86_201{0}_exp.npy".format(i),
sin_dec_bins=sinDec_bins,
log_e_bins=energy_bins)
ps_v004_p00.add_subseason(ic86_i)
# add the combined season
ic86_123456789 = IceCubeSeason(
season_name="IC86_1-9",
sample_name=sample_name,
exp_path=[f"{ps_data_dir}IC86_201{i}_exp.npy" for i in range(1, 10)],
mc_path=mc_path,
grl_path=[f"{grl_data_dir}IC86_201{i}_exp.npy" for i in range(1, 10)],
sin_dec_bins=sinDec_bins,
log_e_bins=energy_bins)
ps_v004_p00.add_season(ic86_123456789)
########################