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

Compare angular acceptance from clsim simulation to POCAM ppc simulation #103

Closed fiedl closed 6 years ago

fiedl commented 6 years ago

"Calibrating the hole ice with the POCAM" Rongen, 2017

In the context of feasability studies for POCAM, Martin has created angular-acceptance curves with ppc direct hole ice simulations. I can use those results to check whether the new clsim-medium-propagation algorithm produces the same angular-acceptance curves as ppc.

bildschirmfoto 2018-08-09 um 00 13 25
fiedl commented 6 years ago

Rongen, "Calibrating the hole ice with the POCAM", 2017, slide 11:

bildschirmfoto 2018-08-05 um 20 53 03

Simulation parameters: 1.8 DOMs, 170cm esca HD baseline Standard geometry

fiedl commented 6 years ago

I'm performing the same simulation:

[2018-08-05 20:56:59] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptance
▶ ./run.rb --cluster --hole-ice-scattering-length=0.102 --hole-ice-radius=0.29718 --distance=1.0 --number-of-photons=1e5 --plane-wave --hole-ice-absorption-length=100.0
▶ mv results/current ../../results/angular_acceptance_for_1.8rdom_esca170cm

[2018-08-05 21:13:53] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./convert_options_txt_to_json.rb ../../results/angular_acceptance_for_1.8rdom_esca170cm
▶ ./plot_angular_acceptance.py --no-log-scale --no-reference-curve ../..
/results/angular_acceptance_for_1.8rdom_esca170cm
bildschirmfoto 2018-08-05 um 21 19 27
fiedl commented 6 years ago

Plot overlay:

bildschirmfoto 2018-08-05 um 21 33 08
fiedl commented 6 years ago

Martin's data files:

[2018-08-08 16:56:40] fiedl@fiedl-mbp ~/icecube/pocam-data
▶ mv ~/Downloads/*.p ./
▶ ls
HD06.14_idealbulk_stdgeo.p  HD18.125_idealbulk_stdgeo.p HD18.170_idealbulk_stdgeo.p noHI_idealbulk_stdgeo.p

Format: https://docs.python.org/2/library/pickle.html

fiedl commented 6 years ago

First test plotting POCAM data: POCAM simulation without hole ice vs. H2 a priori angular sensitivity curve

[2018-08-08 17:15:28] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./plot_angular_acceptance.py --pocam-file ~/icecube/pocam-data/noHI_idealbulk_stdgeo.p --pocam-label "POCAM simulation without hole ice" --h0-reference
bildschirmfoto 2018-08-08 um 17 19 40

Data source: Rongen, https://doi.org/10.22323/1.301.0934, 2017

fiedl commented 6 years ago

Compare it to clsim results: direct detection, plane waves, distance 1m

[2018-08-08 17:21:06] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./plot_angular_acceptance.py --pocam-file ~/icecube/pocam-data/noHI_idealbulk_stdgeo.p --pocam-label "POCAM simulation without hole ice" --h0-reference --no-hole-ice ../../results/angular_acceptance/direct_detection_no_hole_ice_plane_wave
bildschirmfoto 2018-08-08 um 17 22 03 
[2018-08-08 20:57:49] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./plot_angular_acceptance.py ../../results/angular_acceptance/direct_detection_no_hole_ice_plane_wave --pocam-file ~/icecube/pocam-data/noHI_idealbulk_stdgeo.p --pocam-label "POCAM simulation without hole ice" --h0-reference --no-hole-ice
bildschirmfoto 2018-08-08 um 20 58 38
fiedl commented 6 years ago

esca = 170cm, r = 1.8 r_dom:

[2018-08-08 17:23:28] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./plot_angular_acceptance.py --pocam-file ~/icecube/pocam-data/HD18.170_idealbulk_stdgeo.p --pocam-label "POCAM simulation with hole ice: esca=170cm, r=1.8r_dom" --h2-reference
bildschirmfoto 2018-08-08 um 17 26 08
fiedl commented 6 years ago 
[2018-08-08 17:39:06] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./plot_angular_acceptance.py --pocam-file ~/icecube/pocam-data/HD18.170_idealbulk_stdgeo.p --pocam-label "POCAM simulation with hole ice: esca=170cm, r=1.8r_dom" --h2-reference --hole-ice ../../results/angular_acceptance_for_1.8rdom_esca170cm
bildschirmfoto 2018-08-08 um 17 39 16
fiedl commented 6 years ago

r = 1.8 DOMs, 14cm = esca

[2018-08-08 17:51:35] fiedl@kepler17 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptance
▶ ./run.rb --cluster --hole-ice-scattering-length=0.0084 --hole-ice-absorption-length=100.0 --hole-ice-radius=0.29718 --distance=1.0 --number-of-photons=1e5 --plane-wave
▶ mv results/current ../../results/angular_acceptance_for_1.8rdom_esca14cm

[2018-08-08 18:30:23] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./convert_options_txt_to_json.rb ../../results/angular_acceptance_for_1.8rdom_esca14cm
▶ ./plot_angular_acceptance.py --pocam-file ~/icecube/pocam-data/HD06.14_idealbulk_stdgeo.p --pocam-label "POCAM simulation with hole ice: esca=14cm, r=0.6r_dom" --h2-reference --hole-ice ../../results/angular_acceptance_for_1.8rdom_esca14cm
bildschirmfoto 2018-08-08 um 18 38 47

I've used the wrong simulation parameter for the hole-ice radius.

fiedl commented 6 years ago

esca = 14cm, r = 0.6 r_dom = 0.6 * 0.16510m = 0.09906m

[2018-08-08 18:40:28] fiedl@kepler17 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptance
▶ ./run.rb --cluster --hole-ice-scattering-length=0.0084 --hole-ice-absorption-length=100.0 --hole-ice-radius=0.09906 --distance=1.0 --number-of-photons=1e6 --plane-wave
▶ mv results/current ../../results/angular_acceptance_for_0.6rdom_esca14cm

[2018-08-08 19:55:26] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./convert_options_txt_to_json.rb ../../results/angular_acceptance_for_0.6rdom_esca14cm
[2018-08-08 19:56:49] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./plot_angular_acceptance.py --pocam-file ~/icecube/pocam-data/HD06.14_idealbulk_stdgeo.p --pocam-label "POCAM simulation with hole ice: esca=14cm, r=0.6r_dom" --h2-reference --hole-ice ../../results/angular_acceptance_for_0.6rdom_esca14cm --dima-reference
bildschirmfoto 2018-08-08 um 19 57 41 
[2018-08-08 19:57:55] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./plot_angular_acceptance.py --pocam-file ~/icecube/pocam-data/HD06.14_idealbulk_stdgeo.p --pocam-label "POCAM simulation with hole ice: esca=14cm, r=0.6r_dom" --h2-reference --hole-ice ../../results/angular_acceptance_for_0.6rdom_esca14cm
bildschirmfoto 2018-08-08 um 20 03 57
fiedl commented 6 years ago

r = 1.8 r_dom = 0.29718m, esca = 125cm, sca = 0.075m

[2018-08-08 19:55:04] fiedl@kepler17 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptance
▶ ./run.rb --cluster --hole-ice-scattering-length=0.075 --hole-ice-absorption-length=100.0 --hole-ice-radius=0.29718 --distance=1.0 --number-of-photons=1e6 --plane-wave
▶ mv results/current ../../results/angular_acceptance_for_1.8rdom_esca125cm

[2018-08-08 20:54:07] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./convert_options_txt_to_json.rb ../../results/angular_acceptance_for_1.8rdom_esca125cm
▶ ./plot_angular_acceptance.py ../../results/angular_acceptance_for_1.8rdom_esca125cm --pocam-file ~/icecube/pocam-data/HD18.125_idealbulk_stdgeo.p --pocam-label "POCAM simulation with hole ice: esca=125cm, r=1.8r_dom" --hole-ice --h2-reference
bildschirmfoto 2018-08-08 um 20 55 09