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 scan with cable rather than hole ice #101

Open fiedl opened 6 years ago

fiedl commented 6 years ago

What if the cable is the dominant factor?

bildschirmfoto 2018-08-11 um 02 19 22

Image source: https://gallery.icecube.wisc.edu/internal/v/GraphicRe/graphics/arraygraphics2011/sketchup/DOMCloseUp.jpg.html, Jamie Yang

fiedl commented 6 years ago

From this image (source: https://icecube.wisc.edu/gallery/view/153)

I'm trying to model the cable as a first step, by using a cylinder part besides the DOM and one along the string axis, both instantly absorbing.

[2018-08-05 00:16:46] fiedl@fiedl-mbp ~/hole-ice-study/scripts/AngularAcceptance master ⚡
▶ ./run.rb --cable --no-hole-ice
▶ steamshovel tmp/gcd_with_hole_ice.i3
bildschirmfoto 2018-08-05 um 00 20 57
fiedl commented 6 years ago
[2018-08-05 00:24:28] fiedl@kepler16 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptance
▶ ./run.rb --cluster --plane-wave --number-of-photons=1e5 --distance=1.0 --no-hole-ice --cable
▶ mv results/current results/cable_test_ephahM3e
fiedl commented 6 years ago
[2018-08-05 00:30:54] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./convert_options_txt_to_json.rb ../AngularAcceptance/results/cable_test_ephahM3e
▶ ./plot_angular_acceptance.py --no-hole-ice ../AngularAcceptance/results/cable_test_ephahM3e
bildschirmfoto 2018-08-05 um 00 33 31
fiedl commented 6 years ago
[2018-08-05 00:34:41] fiedl@fiedl-mbp ~/hole-ice-study/scripts/AngularAcceptance master ⚡
▶ ./run.rb --cluster --plane-wave --number-of-photons=100 --distance=1.0 --no-hole-ice --cable --cpu
▶ steamshovel tmp/propagated_photons_from_angle_*i3
bildschirmfoto 2018-08-05 um 00 48 12

This confirms that the cable is simulated at the correct position for absorbing photons coming from the side.

fiedl commented 6 years ago
bildschirmfoto 2018-08-05 um 00 58 04
fiedl commented 6 years ago
[2018-08-05 00:59:56] fiedl@kepler16 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptance
▶ ./run.rb --cluster --plane-wave --number-of-photons=1e5 --distance=1.0 --no-hole-ice --cable
▶ mv results/current results/cable_test_uwoh9oaL
fiedl commented 6 years ago
[2018-08-05 01:05:27] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./convert_options_txt_to_json.rb ../AngularAcceptance/results/cable_test_uwoh9oaL
▶ ./plot_angular_acceptance.py --no-hole-ice ../AngularAcceptance/results/cable_test_uwoh9oaL
bildschirmfoto 2018-08-05 um 01 06 26
fiedl commented 6 years ago

Compare this to the H2 model (sca = 50cm, r = 30cm): https://github.com/fiedl/hole-ice-study/issues/80#issuecomment-410475494

[2018-08-05 01:33:26] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./plot_angular_acceptance.py --no-hole-ice --no-reference-curve --no-hole-ice ../AngularAcceptance/results/cable_test_* ../../results/angular_acceptance_karle_h2
bildschirmfoto 2018-08-05 um 01 35 45
fiedl commented 6 years ago

I am rerunning the cable simulation with a cable made of three cylinder sections, this time with 1e6 photons.

[2018-08-11 01:46:13] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptance
▶ ./run.rb --cluster --no-hole-ice --number-of-photons=1e6 --plane-wave --cable
▶ mv results/current ../../results/angular_acceptance_cable_only_Thohgai0

[2018-08-11 02:42:39] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./convert_options_txt_to_json.rb ../../results/angular_acceptance_cable_only_Thohgai0
▶ ./plot_angular_acceptance.py ../../results/angular_acceptance_cable_only_Thohgai0 --no-hole-ice
bildschirmfoto 2018-08-11 um 02 44 45
fiedl commented 6 years ago
[2018-08-11 02:47:34] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./plot_angular_acceptance.py ../../results/angular_acceptance_cable_only_Thohgai0 ../../results/angular_acceptance_karle_h2_1e6_photons ../../results/angular_acceptance/direct_detection_no_hole_ice_plane_wave --no-hole-ice --label "Simulation with cable only" "Simulation with hole ice H2 sca=50cm r=30cm" "Simulation without hole ice"
bildschirmfoto 2018-08-11 um 02 48 49