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

Compare cable-shadow simulation to 7-string flasher data #60

Open fiedl opened 6 years ago

fiedl commented 6 years ago

Basic idea

Scenario: DOM 63_30 (in the middle of string 63) emits a flasher pulse. Count receiving DOM hits for all DOMs in the neighbouring strings 62, 54, 55, 64, 71, and 70.

Flasher data brightess: 127, pulse width: 127.

In the simulation, a photon-absorbing cable is placed at -134.849 degrees from the grid north according to Dima's flasher study results (https://icecube.wisc.edu/~dima/work/IceCube-ftp/leds/). Also add cables for the receiving DOMs.

As hole-ice parameters, use Martin's best values: r = 0.5 r_DOM esca = 0.05m abs = 100m

39135197-b1f4c5c2-4718-11e8-94c5-f2ce4ba927e4

Image source: https://wiki.icecube.wisc.edu/index.php/File:Distances.i86.jpg

fiedl commented 6 years ago

Cable positions

For this simulation, use cable positions from: https://icecube.wisc.edu/~dima/work/IceCube-ftp/leds/

If these results are stored locally at ~/icecube/cable-data/icecube.wisc.edu/~dima/work/IceCube-ftp/leds/results.txt, read out from ruby like this:

# ruby

require_relative 'cable_positions'
CablePosition.find_by(string: 63, dom: 30).angle_from_grid_north_in_degrees
# => -134.849
fiedl commented 6 years ago

Have a look at simulated cables in steamshovel:

[2018-05-03 16:37:08] fiedl@fiedl-mbp ~/hole-ice-study/scripts/FlasherSimulation master ⚡
▶ ./run.rb --cables
▶ steamshovel tmp/gcd_with_hole_ice.i3 
fiedl commented 6 years ago

Could the different distances be responsible for the asymmetry?

Receiving string   distance [m]   max peak
62                 124.8          358.98994325
54                 125.2          361.805595814
55                 125.4          649.15679599
64                 122.0          603.13827089
71                 124.6          521.9574258
70                 125.2          705.07347013

mean distance = 124.5m
mean peak = 533.35

mean distance - max distance = 0.9m = 0.7 percent of the mean distance
abstandsquadratgesetzt -> expected intensity reduction for this distance:
  1/(1.007)^2 = 0.986
0.986 * mean peak = 525.88

which is much larger than the min observed peak of 358.99.

No.