Open fiedl opened 6 years ago
New branch: sf/perfect-bulk-ice
Increase bulk-ice scattering length to very large value by replacing calls to getScatteringLength
with 10000.0
in the propagation-kernel files.
# example: hole_ice.c
# local_scattering_lengths[*number_of_medium_changes] =
# getScatteringLength(photonLayerAtTheCylinderBorder, photonDirAndWlen.w);
local_scattering_lengths[*number_of_medium_changes] =
10000.0;
I'm re-running the simulation with plane waves and direct detection (right upper plot).
[2018-08-09 12:09:08] fiedl@kepler17 /afs/ifh.de/group/amanda/scratch/fiedl/clsim
▶ git co sf/perfect-bulk-ice
[2018-08-09 12:10:12] fiedl@kepler17 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptance
▶ ./run.rb --cluster --no-hole-ice --distance=1.0 --number-of-photons=1e6 --plane-wave
▶ mv results/current ../../results/angular_acceptance/direct_detection_no_hole_ice_plane_wave_perfect_bulk_ice
[2018-08-09 12:37:28] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./convert_options_txt_to_json.rb ../../results/angular_acceptance/direct_detection_no_hole_ice_plane_wave_perfect_bulk_ice
▶ ./plot_angular_acceptance.py ../../results/angular_acceptance/direct_detection_no_hole_ice_plane_wave_perfect_bulk_ice --h0-reference --direct-detection
No hits. The simulation crashes.
Check if the simulation crashes after reducing the bulk-ice scattering length to 100m.
[2018-08-09 15:03:43] fiedl@kepler17 /afs/ifh.de/group/amanda/scratch/fiedl/clsim
▶ git pull
[2018-08-09 15:04:14] fiedl@kepler17 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptance
▶ ./run.rb --cluster --no-hole-ice --distance=1.0 --number-of-photons=1e6 --plane-wave --angles=0
▶ cat results/current/data/angle_hits_and_photons_000.txt
000 4028 1000000.0
Ok, there are hits, now.
Rerun the simulation for all angles with geometric bulk-ice scattering-length 100m.
[2018-08-09 15:11:50] fiedl@kepler17 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptance
▶ ./run.rb --cluster --no-hole-ice --distance=1.0 --number-of-photons=1e6 --plane-wave
The simulation for 180° crashed, the others have worked.
Last log
----------------------------------------------------------
$ tail -n 20 tmp/propagation.log tmp/clsim_180.log
/afs/ifh.de/group/amanda/scratch/fiedl/.rbenv/versions/2.5.0/lib/ruby/gems/2.5.0/gems/fiedl-log-0.3.0/lib/fiedl/log/log.rb:62: warning: Insecure world writable dir /afs/ifh.de/group/amanda/scratch/fiedl in PATH, mode 040777
==> tmp/propagation.log <==
$ tail -n 30 tmp/clsim_180.log
/afs/ifh.de/group/amanda/scratch/fiedl/.rbenv/versions/2.5.0/lib/ruby/gems/2.5.0/gems/fiedl-log-0.3.0/lib/fiedl/log/log.rb:62: warning: Insecure world writable dir /afs/ifh.de/group/amanda/scratch/fiedl in PATH, mode 040777
NOTICE (I3Tray): I3Tray finishing... (I3Tray.cxx:473 in void I3Tray::Execute(unsigned int))
FATAL (I3PhotonToMCPEConverter): distance not 1.000000*165.100000=165.100000mm.. it is 196.447357mm (diff=31.3474mm) (OMKey=(1,19) (photon @ pos=(-255.859,-521.349,194.086)m) (DOM @ pos=(-256.023,-521.282,194)m) (I3PhotonToMCPEConverter.cxx:432 in virtual void I3PhotonToMCPEConverter::DAQ(I3FramePtr))
ERROR (I3Module): I3CLSimMakeHits_0000_makeHitsFromPhotons_clsim_make_hits: Exception thrown (I3Module.cxx:116 in void I3Module::Do(void (I3Module::*)()))
ERROR (I3Module): I3CLSimMakeHits_0000_makePhotons_clsim: Exception thrown (I3Module.cxx:116 in void I3Module::Do(void (I3Module::*)()))
***** running unweighted simulation with a photon pre-scaling of 1.0
number of spectra (1x Cherenkov + Nx flasher): 7
Traceback (most recent call last):
File "/afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/lib/propagate_photons_with_clsim.py", line 173, in <module>
PerformPropagationSimulation()
File "/afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/lib/propagate_photons_with_clsim.py", line 154, in PerformPropagationSimulation
tray.Execute()
File "/afs/ifh.de/group/amanda/scratch/fiedl/software/icecube-simulation-V05-00-07-sl7/build/lib/I3Tray.py", line 234, in Execute
super(I3Tray, self).Execute()
RuntimeError: distance not 1.000000*165.100000=165.100000mm.. it is 196.447357mm (diff=31.3474mm) (OMKey=(1,19) (photon @ pos=(-255.859,-521.349,194.086)m) (DOM @ pos=(-256.023,-521.282,194)m) (in virtual void I3PhotonToMCPEConverter::DAQ(I3FramePtr))
../lib/propagate_photons_and_count_hits.rb:113:in `raise_if_log_file_contains': Simulation has failed with RuntimeError (RuntimeError)
from ../lib/propagate_photons_and_count_hits.rb:118:in `block in <main>'
from ../lib/propagate_photons_and_count_hits.rb:76:in `each'
from ../lib/propagate_photons_and_count_hits.rb:76:in `<main>'
Resume reading the other angles out:
[2018-08-09 15:25:58] fiedl@kepler17 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptance
▶ ./run.rb --cluster --no-hole-ice --distance=1.0 --number-of-photons=1e6 --plane-wave --skip-propagation --resume
▶ rm -r ../../results/angular_acceptance/direct_detection_no_hole_ice_plane_wave_perfect_bulk_ice
▶ mv results/current ../../results/angular_acceptance/direct_detection_no_hole_ice_plane_wave_perfect_bulk_ice
[2018-08-09 15:34:05] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./convert_options_txt_to_json.rb ../../results/angular_acceptance/direct_detection_no_hole_ice_plane_wave_perfect_bulk_ice
▶ ./plot_angular_acceptance.py ../../results/angular_acceptance/direct_detection_no_hole_ice_plane_wave_perfect_bulk_ice --h0-reference --direct-detection --no-hole-ice
The photon number for upper angles is rather high when using plane waves as photon sources, also for simulations without hole ice (#98):
Check whether the the high number of photons for upper angles when using plane waves decreases when setting a very high scattering length in the bulk ice. Then there should be a far less chance of photons scattering in from the side.