.gitignore

-*.png
-*.pdf
-*.ipynb
-*.py
+./*.png
+./*.pdf
+./*.ipynb
+./*.py
+./*.pkl
+./*.hdf5

NuRadioMC_UL.py

+import numpy as np
+from NuRadioReco.utilities import units
+from NuRadioMC.utilities import fluxes
+import pandas as pd
+detectors = ["baseline", "hex_shallow", "hex_shallowheavy", "hex_hybrid_only"]
+
+sensitivities = {}
+for detector in detectors:
+    dfs = [pd.read_json(f"param/effective_volume/Veff__{detector}_array__config_config_ARZ2020_noise__D01detector_sim__{flav}.json") for flav in ["e","mu","tau"]]
+    df = pd.concat(dfs)
+    df.reset_index(inplace=True)
+    print(df['veff'][0])
+    df["veff_energy"] = np.array([item['veff']['combined_4channelPA'][0] for ii, item in df.iterrows()])
+    veffmean = df.groupby("energy").mean()
+
+    livetime=10*units.year
+    sensitivity = fluxes.get_limit_e2_flux(energy=np.array(veffmean.index), veff_sr=np.array(veffmean.veff_energy) * 4 * np.pi,
+                          livetime=livetime)
+
+    sensitivities[detector] = sensitivity
+    sensitivities["energy"] = np.array(veffmean.index)
+
+out_df = pd.DataFrame(sensitivities)
+out_df.to_csv("sensitivities.csv")