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

Angular-acceptance simulations for three azimuthal angles with displaced DOM #117

Open fiedl opened 5 years ago

fiedl commented 5 years ago

Summer has requested the following set of simulations:

Create angular-acceptance plots for incoming plane waves from three different azimuthal directions. The bubble column with given radius and given geometric scattering length needs to be displaced relative to the DOM ("offset").

offset [cm] radius [cm] scattering [cm]
15  7.5 10  
15  15  10
15  15  70
15  30  70

https://www.youtube.com/watch?v=M_QGNdGG9Ew

fiedl commented 5 years ago

image from ios

summerscenario 002 summerscenario 003 summerscenario 004 summerscenario 005

fiedl commented 5 years ago

What are the polar angles?

The string spacing is delta_x = 20.0 meters. The vertical DOM spacing is 10.0 meters. The string length is at least 300 meters. Also, the DOM above and the DOM below may emit light towards the central DOM.

img_3115

import numpy as np
delta_x = 20.0
delta_z_range = np.arange(-150, 151, 10)
incoming_polar_angles = [0] + [90 + np.arctan(1.0 * delta_z / delta_x) / np.pi * 180 for delta_z in delta_z_range] + [180]
fiedl commented 5 years ago

What is perfect bulk ice?

Or, not use perfect bulk ice, but regular bulk ice.

fiedl commented 5 years ago

Regular bulk ice:

Starting photons from 20m horizontal distance at z=500 relative to the IceCube center, most photons are scattered away before reaching the receiving DOM.

bildschirmfoto 2018-10-18 um 18 46 18

The "ice paper", 2013, "Measurement of South Pole ice transparency with the IceCube LED calibration system", shows the depth dependence of scattering and absorption:

bildschirmfoto 2018-10-18 um 18 54 41

The two translate into each other that the ice surface is at z = 1948.07 m (https://wiki.icecube.wisc.edu/index.php/Coordinate_system). z=500m equals a depth=1448.07.

This is the very top of the detector. The effective scattering length should be about 14m there according to the ice-paper figure. This corresponds to a geometric scattering length of 84cm, which explains the observed behaviour.

At the center of the detector, it's even worse. One would have to be lucky to sit in a more transparent layer. Thus, just keep this ice and start more photons, or start them from a closer distance to the receiving DOM.

fiedl commented 5 years ago

Submitting a test job to the cluster:

[2018-10-18 19:17:13] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study
▶ git pull

[2018-10-18 19:19:08] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptanceFromSeveralAzimuthalAngles
▶ qsub -l gpu -l tmpdir_size=10G -l s_rt=4:0:00 -l h_rss=2G -m ae -t 1-1 batch-job.sh

In this configuration, the simulation takes about 2mins. But this is still for only two polar angles.

fiedl commented 5 years ago

Submit array job to cluster:

[2018-10-18 19:38:33] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptanceFromSeveralAzimuthalAngles
▶ git pull
▶ qsub -l gpu -l tmpdir_size=10G -l s_rt=4:0:00 -l h_rss=2G -m ae -t 1-12 batch-job.sh
fiedl commented 5 years ago

Gauging simulation: Count the number of photons arriving for the same starting distance and plane extend when photons are coming from below and no hole ice is present.

[2018-10-18 19:48:10] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/FiringRange
▶ $ICESIM/env-shell.sh
▶ rm -r tmp
▶ ./run.rb --distance=3.0 --number-of-photons=1e7 --number-of-runs=1 --number-of-parallel-runs=1 --angle=0 --plane-wave --no-hole-ice
fiedl commented 5 years ago 
[2018-10-18 19:56:13] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/FiringRange
▶ git pull
▶ ./run.rb --distance=3.0 --number-of-photons=1e7 --number-of-runs=1 --number-of-parallel-runs=1 --angle=0 --plane-wave
SIMULATION RESULT
[3814]

gauging_factor = 3814 / 1e7 = 0.0003814

fiedl commented 5 years ago

Concerning the cluster run, I just got no space left on device. Start with more space requirements:

[2018-10-18 20:18:58] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptanceFromSeveralAzimuthalAngles
▶ qsub -l gpu -l tmpdir_size=20G -l s_rt=1:0:00 -l h_rss=2G -m ae -t 1-12 batch-job.sh
fiedl commented 5 years ago

Convert options files:

[2018-10-18 21:04:21] fiedl@fiedl-mbp ~/hole-ice-study/scripts master ⚡
▶ lib/convert_options_txt_to_json.rb AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results
fiedl commented 5 years ago

Comparison plots, each plot for different hole-ice parameters, azimuthal angles in one plot.

[2018-10-18 21:59:22] fiedl@fiedl-mbp ~/hole-ice-study/scripts master ⚡
▶ lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_7.5cm_sca_10.0cm_azi* --label "r=7.5cm sca=10cm azi=0" "r=7.5cm sca=10cm azi=120deg" "r=7.5cm sca=10cm azi=240deg"
bildschirmfoto 2018-10-18 um 22 00 33 
[2018-10-18 21:54:10] fiedl@fiedl-mbp ~/hole-ice-study/scripts master ⚡
▶ lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_10.0cm_azi0 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_10.0cm_azi120 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_10.0cm_azi240 --label "r=15cm sca=10cm azi=0" "r=15cm sca=10cm aiz=120deg" "r=15cm sca=10cm aiz=240deg"

[2018-10-18 21:54:36] fiedl@fiedl-mbp ~/hole-ice-study/scripts master ⚡
▶ lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_10.0cm_azi* --label "r=15cm sca=10cm azi=0" "r=15cm sca=10cm aiz=120deg" "r=15cm sca=10cm aiz=240deg"
bildschirmfoto 2018-10-18 um 21 54 16 
[2018-10-18 21:56:03] fiedl@fiedl-mbp ~/hole-ice-study/scripts master ⚡
▶ lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_70.0cm_azi* --label "r=15cm sca=70cm azi=0" "r=15cm sca=70cm azi=120deg" "r=15cm sca=70cm azi=240deg"
bildschirmfoto 2018-10-18 um 21 57 16 
[2018-10-18 21:57:25] fiedl@fiedl-mbp ~/hole-ice-study/scripts master ⚡
▶ lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_30.0cm_sca_70.0cm_azi* --label "r=30cm sca=70cm azi=0" "r=30cm sca=70cm azi=120deg" "r=30cm sca=70cm azi=240deg"
bildschirmfoto 2018-10-18 um 21 59 11
fiedl commented 5 years ago 
[2018-10-18 22:06:38] fiedl@fiedl-mbp ~/hole-ice-study/scripts master ⚡
▶ lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_10.0cm_azi0 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_70.0cm_azi0 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_30.0cm_sca_70.0cm_azi0 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_7.5cm_sca_10.0cm_azi0 --label "r=15cm sca=10cm azi=0" "r=15cm sca=70cm azi=0" "r=30cm sca=70 azi=0" "r=7.5cm sca=10cm azi=0"
bildschirmfoto 2018-10-18 um 22 06 52 
[2018-10-18 22:07:00] fiedl@fiedl-mbp ~/hole-ice-study/scripts master ⚡
▶ lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_10.0cm_azi120 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_70.0cm_azi120 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_30.0cm_sca_70.0cm_azi120 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_7.5cm_sca_10.0cm_azi120 --label "r=15cm sca=10cm azi=120deg" "r=15cm sca=70cm azi=120deg" "r=30cm sca=70 azi=120deg" "r=7.5cm sca=10cm azi=120deg"
bildschirmfoto 2018-10-18 um 22 08 23 
[2018-10-18 22:08:39] fiedl@fiedl-mbp ~/hole-ice-study/scripts master ⚡
▶ lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_10.0cm_azi240 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_70.0cm_azi240 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_30.0cm_sca_70.0cm_azi240 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_7.5cm_sca_10.0cm_azi240 --label "r=15cm sca=10cm azi=240deg" "r=15cm sca=70cm azi=240deg" "r=30cm sca=70 azi=240deg" "r=7.5cm sca=10cm azi=240deg"
bildschirmfoto 2018-10-18 um 22 09 54
fiedl commented 5 years ago

Restart cluster simulation with more (1e7 rather than 1e6) photons per angle.

[2018-10-18 22:53:33] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptanceFromSeveralAzimuthalAngles
▶ git pull
▶ qsub -l gpu -l tmpdir_size=30G -l s_rt=4:0:00 -l h_rss=2G -m ae -t 1-12 batch-job.sh

Request more memory:

[2018-10-18 23:00:27] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptanceFromSeveralAzimuthalAngles
▶ qsub -l gpu -l tmpdir_size=30G -l s_rt=4:0:00 -l h_rss=20G -m ae -t 1-12 batch-job.sh

These jobs ran 2.5h each, but produced no files. From the logs, the first polar angle succeeded, but the second already experienced "RuntimeError: output stream error". Disk space again?

fiedl commented 5 years ago

Local run to record photon paths:

[2018-10-18 23:13:03] fiedl@fiedl-mbp ~/hole-ice-study/scripts/AngularAcceptanceFromSeveralAzimuthalAngles master ⚡
▶ ./run.rb
▶ steamshovel results/*/propagated_photons_from_angle_125.i3
bildschirmfoto 2018-10-18 um 23 40 36 bildschirmfoto 2018-10-18 um 23 40 58
fiedl commented 5 years ago

Steamshovel visualization on youtube: https://www.youtube.com/watch?v=M_QGNdGG9Ew

fiedl commented 5 years ago

Replot without log scale:

[2018-10-19 12:27:03] fiedl@fiedl-mbp ~/hole-ice-study/scripts master ⚡

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_7.5cm_sca_10.0cm_azi* --label "r=7.5cm sca=10cm azi=0" "r=7.5cm sca=10cm azi=120deg" "r=7.5cm sca=10cm azi=240deg" --no-log-scale --outfile "r7.5cm_sca10cm.png"

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_10.0cm_azi* --label "r=15cm sca=10cm azi=0" "r=15cm sca=10cm aiz=120deg" "r=15cm sca=10cm aiz=240deg" --no-log-scale --outfile "r15cm_sca10cm.png"

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_70.0cm_azi* --label "r=15cm sca=70cm azi=0" "r=15cm sca=70cm azi=120deg" "r=15cm sca=70cm azi=240deg" --no-log-scale --outfile "r15cm_sca70cm.png"

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_30.0cm_sca_70.0cm_azi* --label "r=30cm sca=70cm azi=0" "r=30cm sca=70cm azi=120deg" "r=30cm sca=70cm azi=240deg" --no-log-scale --outfile "r30cm_sca70cm.png"

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_10.0cm_azi0 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_70.0cm_azi0 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_30.0cm_sca_70.0cm_azi0 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_7.5cm_sca_10.0cm_azi0 --label "r=15cm sca=10cm azi=0" "r=15cm sca=70cm azi=0" "r=30cm sca=70cm azi=0" "r=7.5cm sca=10cm azi=0" --no-log-scale --outfile "aiz0.png"

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_10.0cm_azi120 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_70.0cm_azi120 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_30.0cm_sca_70.0cm_azi120 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_7.5cm_sca_10.0cm_azi120 --label "r=15cm sca=10cm azi=120deg" "r=15cm sca=70cm azi=120deg" "r=30cm sca=70cm azi=120deg" "r=7.5cm sca=10cm azi=120deg" --no-log-scale --outfile "azi120deg.png"

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_10.0cm_azi240 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_70.0cm_azi240 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_30.0cm_sca_70.0cm_azi240 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_7.5cm_sca_10.0cm_azi240 --label "r=15cm sca=10cm azi=240deg" "r=15cm sca=70cm azi=240deg" "r=30cm sca=70cm azi=240deg" "r=7.5cm sca=10cm azi=240deg" --no-log-scale --outfile "azi240deg.png"

r7 5cm_sca10cm r15cm_sca10cm r15cm_sca70cm r30cm_sca70cm aiz0 azi120deg azi240deg

fiedl commented 5 years ago

Rerun on the cluster with parameters for effective scattering length.

offset [cm] radius [cm] effective scattering length [cm]
15  7.5 10  
15  15  10
15  15  70
15  30  70
[2018-10-19 18:20:22] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptanceFromSeveralAzimuthalAngles
▶ git pull
▶ rm -r cluster-results
▶ qsub -l gpu -l tmpdir_size=50G -l s_rt=10:0:00 -l h_rss=30G -m ae -t 1-12 batch-job.sh
fiedl commented 5 years ago 
[2018-10-19 18:54:55] fiedl@fiedl-mbp ~/hole-ice-study/scripts master ⚡
▶ git pull
▶ lib/convert_options_txt_to_json.rb AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results
▶ 

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_7.5cm_sca_0.6cm_azi* --label "r=7.5cm esca=10cm azi=0" "r=7.5cm esca=10cm azi=120deg" "r=7.5cm esca=10cm azi=240deg" --no-log-scale --outfile "r7.5cm_sca10cm.png"

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_0.6cm_azi* --label "r=15cm esca=10cm azi=0" "r=15cm esca=10cm aiz=120deg" "r=15cm esca=10cm aiz=240deg" --no-log-scale --outfile "r15cm_sca10cm.png"

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_4.199999999999999cm_azi* --label "r=15cm esca=70cm azi=0" "r=15cm esca=70cm azi=120deg" "r=15cm esca=70cm azi=240deg" --no-log-scale --outfile "r15cm_sca70cm.png"

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_30.0cm_sca_4.199999999999999cm_azi* --label "r=30cm esca=70cm azi=0" "r=30cm esca=70cm azi=120deg" "r=30cm esca=70cm azi=240deg" --no-log-scale --outfile "r30cm_sca70cm.png"

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_0.6cm_azi0 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_4.199999999999999cm_azi0 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_30.0cm_sca_4.199999999999999cm_azi0 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_7.5cm_sca_0.6cm_azi0 --label "r=15cm esca=10cm azi=0" "r=15cm esca=70cm azi=0" "r=30cm esca=70cm azi=0" "r=7.5cm esca=10cm azi=0" --no-log-scale --outfile "aiz0.png"

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_0.6cm_azi120 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_4.199999999999999cm_azi120 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_30.0cm_sca_4.199999999999999cm_azi120 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_7.5cm_sca_0.6cm_azi120 --label "r=15cm esca=10cm azi=120deg" "r=15cm esca=70cm azi=120deg" "r=30cm esca=70cm azi=120deg" "r=7.5cm esca=10cm azi=120deg" --no-log-scale --outfile "azi120deg.png"

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_0.6cm_azi240 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_15.0cm_sca_4.199999999999999cm_azi240 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_30.0cm_sca_4.199999999999999cm_azi240 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results/offset_15.0cm_radius_7.5cm_sca_0.6cm_azi240 --label "r=15cm esca=10cm azi=240deg" "r=15cm esca=70cm azi=240deg" "r=30cm esca=70cm azi=240deg" "r=7.5cm esca=10cm azi=240deg" --no-log-scale --outfile "azi240deg.png"

Total photon hit count: 12020
Total photon hit count: 9572
Total photon hit count: 11012
Total photon hit count: 7235
Total photon hit count: 15566
Total photon hit count: 12171
Total photon hit count: 12102

r7 5cm_sca10cm r15cm_sca10cm r15cm_sca70cm r30cm_sca70cm

aiz0 azi120deg azi240deg

fiedl commented 5 years ago

Recreating plots with 1e7 photons per angle.

[2018-11-07 14:36:43] fiedl@fiedl-mbp ~/hole-ice-study/scripts master ⚡
▶ git pull
▶ lib/convert_options_txt_to_json.rb AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results-1e7
▶ 

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results-1e7/offset_15.0cm_radius_7.5cm_sca_0.6cm_azi* --label "r=7.5cm esca=10cm azi=0" "r=7.5cm esca=10cm azi=120deg" "r=7.5cm esca=10cm azi=240deg" --no-log-scale --outfile "r7.5cm_sca10cm.png"
Total photon hit count: 118441

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results-1e7/offset_15.0cm_radius_15.0cm_sca_0.6cm_azi* --label "r=15cm esca=10cm azi=0" "r=15cm esca=10cm aiz=120deg" "r=15cm esca=10cm aiz=240deg" --no-log-scale --outfile "r15cm_sca10cm.png"
Total photon hit count: 94866

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results-1e7/offset_15.0cm_radius_15.0cm_sca_4.199999999999999cm_azi* --label "r=15cm esca=70cm azi=0" "r=15cm esca=70cm azi=120deg" "r=15cm esca=70cm azi=240deg" --no-log-scale --outfile "r15cm_sca70cm.png"
Total photon hit count: 108198

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results-1e7/offset_15.0cm_radius_30.0cm_sca_4.199999999999999cm_azi* --label "r=30cm esca=70cm azi=0" "r=30cm esca=70cm azi=120deg" "r=30cm esca=70cm azi=240deg" --no-log-scale --outfile "r30cm_sca70cm.png"
Total photon hit count: 72595

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results-1e7/offset_15.0cm_radius_15.0cm_sca_0.6cm_azi0 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results-1e7/offset_15.0cm_radius_15.0cm_sca_4.199999999999999cm_azi0 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results-1e7/offset_15.0cm_radius_30.0cm_sca_4.199999999999999cm_azi0 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results-1e7/offset_15.0cm_radius_7.5cm_sca_0.6cm_azi0 --label "r=15cm esca=10cm azi=0" "r=15cm esca=70cm azi=0" "r=30cm esca=70cm azi=0" "r=7.5cm esca=10cm azi=0" --no-log-scale --outfile "aiz0.png"
Total photon hit count: 153886

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results-1e7/offset_15.0cm_radius_15.0cm_sca_0.6cm_azi120 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results-1e7/offset_15.0cm_radius_15.0cm_sca_4.199999999999999cm_azi120 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results-1e7/offset_15.0cm_radius_30.0cm_sca_4.199999999999999cm_azi120 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results-1e7/offset_15.0cm_radius_7.5cm_sca_0.6cm_azi120 --label "r=15cm esca=10cm azi=120deg" "r=15cm esca=70cm azi=120deg" "r=30cm esca=70cm azi=120deg" "r=7.5cm esca=10cm azi=120deg" --no-log-scale --outfile "azi120deg.png"
Total photon hit count: 120035

lib/plot_angular_acceptance.py AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results-1e7/offset_15.0cm_radius_15.0cm_sca_0.6cm_azi240 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results-1e7/offset_15.0cm_radius_15.0cm_sca_4.199999999999999cm_azi240 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results-1e7/offset_15.0cm_radius_30.0cm_sca_4.199999999999999cm_azi240 AngularAcceptanceFromSeveralAzimuthalAngles/cluster-results-1e7/offset_15.0cm_radius_7.5cm_sca_0.6cm_azi240 --label "r=15cm esca=10cm azi=240deg" "r=15cm esca=70cm azi=240deg" "r=30cm esca=70cm azi=240deg" "r=7.5cm esca=10cm azi=240deg" --no-log-scale --outfile "azi240deg.png"
Total photon hit count: 120179

aiz0 azi120deg azi240deg r7 5cm_sca10cm r15cm_sca10cm r15cm_sca70cm r30cm_sca70cm

fiedl commented 5 years ago

Creating visualisation:

[2018-11-07 15:11:32] fiedl@fiedl-mbp ~/hole-ice-study/scripts/AngularAcceptanceFromSeveralAzimuthalAngles master ⚡
▶ ./run.rb

[2018-11-07 15:21:36] fiedl@fiedl-mbp ~/hole-ice-study/scripts/AngularAcceptanceFromSeveralAzimuthalAngles/results/offset_15.0cm_radius_15.0cm_sca_0.6cm_azi120 master ⚡
▶ steamshovel propagated_photons_from_angle_125.i3

video: https://youtu.be/Wiu8CpVQn14