Plot arrival-time histograms for different hole-ice parameters

[fiedl]

Can one see a difference between hole-ice parameters in these histograms?

Other cross checks: #63

[fiedl]

Use data from the flasher run and flasher simulation (#59).

Flashing DOM: (63,30) Receiving DOM: (70,30)

[fiedl]

[2018-07-16 19:26:46] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ python create_time_distribution_histogram.py ../../results/flasher_without_hole_ice ../../results/flasher_parameter_scan_full_brightness_with_thinning/esca0.8_r0.5rdom_abs100_width127_bright127 ../../results/flasher_parameter_scan_full_brightness_with_thinning/esca0.1_r1.0rdom_abs100_width127_bright127 ../../results/flasher_standard_clsim_hole_ice_approximation

In the above plot, all hits from DOMs of the receiving strings are summed:

receiving_strings = [70, 71, 64, 55, 62 ,54]
eceiver_data = flasher_data[flasher_data.string_number.isin(receiving_strings)]

[fiedl]

General expectations if there is more scattering:

• Less photons should arrive as they are scattered away or effectively reflected.
• The photons should arrive later as scattering costs time.

[fiedl]

I'm producing a simulation with a large hole-ice cylinder to see if there is any time delay.

Only place a cylinder around string 70.

strings_with_hole_ice = [70]

[2018-07-18 16:09:48] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/FlasherSimulation.tmp
▶ ./run.rb --width=127 --brightness=127 --thinning-factor=0.1 --effective-scattering-length=0.1 --absorption-length=100 --hole-ice-radius-in-dom-radii=40
▶ mv results/Run-2018-aje3ook7 $SCRATCH/hole-ice-study/results/flasher_arrival_time_large_cylinder

[fiedl]

Place cylinders around all sending and receiving strings.

strings_with_hole_ice = [63] + [62, 54, 55, 64, 71, 70] + [61, 53, 44, 45, 46, 56, 65, 72, 78, 77, 76, 69]

[2018-07-18 17:18:15] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/FlasherSimulation.tmp
▶ ./run.rb --width=127 --brightness=127 --thinning-factor=0.001 --effective-scattering-length=0.1 --absorption-length=100 --hole-ice-radius-in-dom-radii=40

[2018-07-18 17:59:21] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/FlasherSimulation.tmp
▶ rm -r $SCRATCH/hole-ice-study/results/flasher_arrival_time_large_cylinder/*
▶ mv results/Run-2018-veib3gee $SCRATCH/hole-ice-study/results/flasher_arrival_time_large_cylinder

[fiedl]

There were too few hits in the results file. Start again with more photons, but a smaller scattering length.

[2018-07-18 18:58:39] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/FlasherSimulation.tmp
▶ ./run.rb --width=127 --brightness=127 --thinning-factor=0.1 --effective-scattering-length=1.0 --absorption-length=100 --hole-ice-radius-in-dom-radii=40

[2018-07-18 21:58:20] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/FlasherSimulation.tmp
▶ rm -r $SCRATCH/hole-ice-study/results/flasher_arrival_time_large_cylinder/*
▶ mv results/Run-2018-shai7Cai $SCRATCH/hole-ice-study/results/flasher_arrival_time_large_cylinder

[fiedl]

[2018-07-18 22:24:04] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ python create_time_distribution_histogram.py ../../results/flasher_arrival_time_large_cylinder

[fiedl]

This looks fine. Let's simulate the other side:

[2018-07-18 23:27:06] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/FlasherSimulation.tmp
▶ ./run.rb --width=127 --brightness=127 --thinning-factor=0.1 --effective-scattering-length=1.0 --absorption-length=100 --hole-ice-radius-in-dom-radii=1.0

[2018-07-19 07:13:59] fiedl@wgs16 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/FlasherSimulation.tmp
▶ mkdir $SCRATCH/hole-ice-study/results/flasher_arrival_time
▶ mv results/Run-2018-eis7phoW $SCRATCH/hole-ice-study/results/flasher_arrival_time/

[fiedl]

[2018-07-19 07:26:25] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ python create_time_distribution_histogram.py ../../results/flasher_arrival_time_large_cylinder ../../results/flasher_arrival_time

[fiedl]

• [x] Fix plot-axes labels
• [x] Interpret histogram features on the left-hand side

[fiedl]

[2018-07-19 10:00:27] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ python create_time_distribution_histogram.py ../../results/flasher_arrival_time_large_cylinder ../../results/flasher_arrival_time

[fiedl]

When there is more scattering then there are more steps where the photon chooses a new direction. Thus, there are more possibilities when to arrive at the DOM. Thus, for strong hole ice, i.e. where the scattering length is small compared to the hole-ice radius, the left-hand side of the histogram is less steep.