Dima's hole-ice model - Githubissues

fiedl commented 6 years ago

https://wiki.icecube.wisc.edu/index.php/MSU_Forward_Hole_Ice

Another advantage of using Dima's hole ice instead of H2 to extend the parametrization is that Dima's hole ice has a somewhat simple formula with a single parameter p as systematic, while the H2 model was made via fitting the angular acceptance curve of different simulated hole ices with different scattering lengths which means each of the fits turned out somewhat different as seen in the figure above and one parametrization would not be as easily transferable to the hole ice at different scattering parameters.

https://docushare.icecube.wisc.edu/dsweb/Get/Document-75847/20151116_jpamdeandre_DC_holeice.pdf

https://github.com/fiedl/hole-ice-study/issues/80#issuecomment-410478799

import matplotlib.pyplot as plt
import numpy as np

theta = np.arange(0, np.pi, 0.01)

# Dima's model has a parameter p
# typically between 0.2 and 0.4
p = 0.3

A_dima = 0.34 * (1 + 1.5 * np.cos(theta) - np.cos(theta)**3/2) + p * np.cos(theta) * (np.cos(theta)**2 - 1)**3

f, ax = plt.subplots()
ax.plot(costheta, A_dima)
ax.set_xlabel("cos theta")
ax.set_ylabel("relative sensitivity")
plt.show()

p = 0.4
A_dima = 0.34 * (1 + 1.5 * np.cos(theta) - np.cos(theta)**3/2) + p * np.cos(theta) * (np.cos(theta)**2 - 1)**3

f, ax = plt.subplots()
ax.plot(costheta, A_dima)
ax.set_xlabel("cos theta")
ax.set_ylabel("relative sensitivity")
ax.set_yscale("log")
plt.ylim(ymin=1e-3)
plt.show()

Compare this to the new simulations.

fiedl commented 6 years ago

Plot the simulation with parameters from the H2 hole-ice model (#80) in comparison to Dima's model.

[2018-08-05 18:24:33] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./plot_angular_acceptance.py --no-reference-curve --dima-reference-curve ../../results/angular_acceptance_karle_h2

bildschirmfoto 2018-08-05 um 18 25 03

fiedl commented 6 years ago

Comparison of Dima's model to the a priori angular-acceptance curve. The curve called H2 in this plot is the reference plot (#10).

h2-dima-baselines-jpa14

Image source: https://wiki.icecube.wisc.edu/index.php/File:H2_Dima_baselines.png

The H1, H2, H3 curves are fitted to polynomials. The parameters are here:

fiedl commented 6 years ago

[2018-08-07 18:20:45] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./plot_angular_acceptance.py --dima-reference-curve

Dima's function is described above.

fiedl commented 6 years ago

[2018-08-07 18:21:23] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./plot_angular_acceptance.py --dima-reference --h0-reference --h2-reference

fiedl / hole-ice-study

Dima's hole-ice model #102