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

Visualize flasher scenario with steamshovel #107

Open fiedl opened 6 years ago

fiedl commented 6 years ago

Flasher study: #59

Simulate less photons locally and record photon paths. Visualize in steamshovel.

fiedl commented 6 years ago 
# FlasherSimulation/run.rb
# ...
propagation_options = {
  save_photon_paths: true,
  cpu: true,
  # ...
}
[2018-08-09 00:54:31] fiedl@fiedl-mbp ~/hole-ice-study/scripts/FlasherSimulation master ⚡
▶ ./run.rb --width=127 --brightness=127 --thinning-factor=0.00000001 --effective-scattering-length=0.14 --hole-ice-radius-in-dom-radii=0.6
▶ steamshovel tmp/propagated_photons.i3
bildschirmfoto 2018-08-09 um 01 01 59 bildschirmfoto 2018-08-09 um 01 02 50
fiedl commented 6 years ago

I need more photons to photograph them hitting a receiving DOM.

[2018-08-09 01:06:49] fiedl@fiedl-mbp ~/hole-ice-study/scripts/FlasherSimulation master ⚡
▶ ./run.rb --width=127 --brightness=127 --thinning-factor=0.000001 --effective-scattering-length=0.14 --hole-ice-radius-in-dom-radii=0.6
▶ steamshovel tmp/propagated_photons.i3
bildschirmfoto 2018-08-09 um 01 13 53

This is not easy.

Also, I would have expected photons being caught in a cylinder and then scattered into the DOM.

fiedl commented 6 years ago

I'm removing the cylinder around the sending DOM. Then restarting the simulation with extra-large cylinders.

[2018-08-09 01:17:05] fiedl@fiedl-mbp ~/hole-ice-study/scripts/FlasherSimulation master ⚡
▶ ./run.rb --width=127 --brightness=127 --thinning-factor=0.0000001 --effective-scattering-length=0.14 --hole-ice-radius-in-dom-radii=5.0
▶ steamshovel tmp/propagated_photons.i3
bildschirmfoto 2018-08-09 um 01 25 56

Still no luck.

fiedl commented 6 years ago

Next step: Record only photons that hit a DOM, but save photon paths. Does this work in combination?

fiedl commented 6 years ago

In propagate_photons_with_clsim.py, I'm using I3CLSimMakePhotons as tray segment, but setting StopDetectedPhotons = True and SaveAllPhotons = False, hoping only to visualize the photons that hit a receiving DOM.

[2018-08-09 09:36:26] fiedl@fiedl-mbp ~/hole-ice-study/scripts/FlasherSimulation master ⚡
▶ ./run.rb --width=127 --brightness=127 --thinning-factor=0.00001 --effective-scattering-length=0.14 --hole-ice-radius-in-dom-radii=5.0
▶ steamshovel tmp/propagated_photons.i3

bildschirmfoto 2018-08-09 um 09 45 04

fiedl commented 6 years ago

More photons, but smaller cylinder:

[2018-08-09 09:50:30] fiedl@fiedl-mbp ~/hole-ice-study/scripts/FlasherSimulation master ⚡
▶ ./run.rb --width=127 --brightness=127 --thinning-factor=0.001 --effective-scattering-length=0.14 --hole-ice-radius-in-dom-radii=0.6