search

fiedl / hole-ice-study

This project aims to incorporate the effects of hole ice into the clsim photon propagation simulation of the icecube neutrino observatory.
Other
4 stars 2 forks source link

Where do the reference plots (nominal and hole ice) come from? #10

Open fiedl opened 6 years ago

fiedl commented 6 years ago

Reference plot

As reference for comparing the angular acceptance plots, I'm using figure 7 from the "Measurement of South Pole ice transparency with the IceCube LED calibration system" paper (2013).

bildschirmfoto 2018-02-14 um 12 57 33

I've taken the parameters for this polynomial from clsim/python/GetIceCubeDOMAngularSensitivity.py. See also: https://github.com/claudiok/clsim/commit/340637e6367d53abe34f6903d953a70eda813a32.

Where do these plots come from?

fiedl commented 6 years ago

Assumed hole-ice parameters

The "nominal" curve corresponds to the H0 model, the "hole ice" curve to the H2 model in a hole-ice study by Karle using a YAG laser in 1998: http://icecube.berkeley.edu/kurt/interstring/hole-ice/yak.html

Hole ice models Model H0: hole ice identical to bulk ice Model H1: L_scatt = 100 cm Model H2: L_scatt = 50 cm Model H3: L_scatt = 30 cm Model H4: L_scatt = 10 cm Scattering length L_scatt = geometric scattering length on air bubbles. Hole diameter = 60 cm.

Using the photonics software, effective angular-acceptance curves have been created accounting for the hole-ice effects of the respective models.

Curve parameters

The angular-acceptance curves have been fitted to polynomials, for example for the H0 model:

bildschirmfoto 2018-08-07 um 21 19 53

Coefficients for the models: http://code.icecube.wisc.edu/projects/icecube/browser/IceCube/projects/ice-models/trunk/resources/models/angsens?order=name

Coefficients for H0: 0.26266 0.47659 0.15480 -0.14588 0.17316 1.3070 0.44441 -2.3538 -1.3564 1.2098 0.81569

H1: 0.29958 0.52578 0.18281 -0.11287 0.48188 1.0777 -0.96676 -2.2425 0.70419E-01 0.95132 0.13426

H2: 0.32813 0.63899 0.20049 -1.2250 -0.14470 4.1695 0.76898 -5.8690 -2.0939 2.3834 1.0435

H3: 0.35191 0.68928 0.35927 -1.7176 -1.1212 5.3179 2.2695 -7.2229 -3.2619 2.9822 1.4572

fiedl commented 6 years ago

Plot of H1, H2, H3 against Dima's hole ice model (#102):

h2-dima-baselines-jpa14

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