Open fiedl opened 6 years ago
The simplest geometric approach would be to assume the acceptance A(theta) proportional to the visible fraction of the detector surface.
# python
import matplotlib.pyplot as plt
import numpy as np
theta = np.arange(0, np.pi, 0.01)
costheta = np.cos(theta)
A = 1 - theta / np.pi
f, ax = plt.subplots()
ax.plot(costheta, A)
ax.set_xlabel("cos theta")
ax.set_ylabel("relative sensitivity")
plt.show()
ax.set_yscale("log")
plt.show()
Martin's approach (Rongen, New York, 2016-04-19, slide 17:
# python
A = 16.02 * np.cos(theta) - 4 * np.sin(theta) + 16.51
f, ax = plt.subplots()
ax.plot(costheta, A)
ax.set_xlabel("cos theta")
ax.set_ylabel("relative sensitivity")
plt.show()
ax.set_yscale("log")
plt.show()
[2018-07-31 09:18:08] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptance
▶ rm -r tmp
▶ ./run.rb --cluster --no-hole-ice --distance=1.0 --number-of-photons=1e6 --plane-wave
[2018-07-31 12:42:10] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study
▶ mv scripts/AngularAcceptance/results/current results/angular_acceptance/direct_detection_no_hole_ice_plane_wave
[2018-07-31 12:48:01] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptance
▶ mv tmp /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/results/angular_acceptance/direct_detection_no_hole_ice_plane_wave/
[2018-07-31 12:52:36] 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
▶ ./plot_angular_acceptance.py ../../results/angular_acceptance/direct_detection_no_hole_ice_plane_wave
[2018-08-01 21:09:27] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./plot_angular_acceptance.py --no-hole-ice ../../results/angular_acceptance/direct_detection_no_hole_ice_plane_wave/
[2018-08-01 19:14:44] fiedl@fiedl-mbp ~/hole-ice-study/scripts/FiringRange master ⚡
▶ ./run.rb --distance=1.0 --number-of-photons=1e3 --save-photon-paths --cpu --angle=45 --no-hole-ice --plane-wave
▶ steamshovel tmp/propagated_photons.i3
[2018-07-31 12:48:31] fiedl@kepler00 /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
[2018-07-31 12:57:20] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptance
▶ mv results/current ../../results/angular_acceptance/direct_detection_no_hole_ice_pencil_beam
▶ mv tmp ../../results/angular_acceptance/direct_detection_no_hole_ice_pencil_beam/
[2018-07-31 13:07:45] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./convert_options_txt_to_json.rb ../../results/angular_acceptance/direct_detection_no_hole_ice_pencil_beam
▶ ./plot_angular_acceptance.py ../../results/angular_acceptance/direct_detection_no_hole_ice_pencil_beam
[2018-08-01 20:58:54] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./plot_angular_acceptance.py --no-hole-ice --pencil-beam ../../results/angular_acceptance/direct_detection_no_hole_ice_pencil_beam
[2018-08-01 19:02:49] fiedl@fiedl-mbp ~/hole-ice-study/scripts/FiringRange master ⚡
▶ ./run.rb --distance=1.0 --number-of-photons=1e3 --save-photon-paths --cpu --angle=45 --no-hole-ice
▶ steamshovel tmp/propagated_photons.i3
Regarding the opening angle:
# GeneratePhotonFlashModule.py
pulse = I3CLSimFlasherPulse()
pulse.SetPulseWidth(1. * I3Units.ns)
pulse.SetAngularEmissionSigmaPolar( 0.001 * I3Units.deg )
pulse.SetAngularEmissionSigmaAzimuthal( 0.001 * I3Units.deg )
[2018-08-01 17:20:21] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/clsim
▶ git pull
[2018-08-01 17:22:03] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/software/icecube-simulation-V05-00-07-sl7/build
▶ make -j 4
[2018-08-01 17:28:33] fiedl@kepler16 /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
[2018-08-01 17:47:33] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptance
▶ mv results/current ../../results/angular_acceptance/no_direct_detection_no_hole_ice_plane_waves
▶ mv tmp ../../results/angular_acceptance/no_direct_detection_no_hole_ice_plane_waves/
[2018-08-01 17:52:50] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./convert_options_txt_to_json.rb ../../results/angular_acceptance/no_direct_detection_no_hole_ice_plane_waves
▶ ./plot_angular_acceptance.py --no-hole-ice ../../results/angular_acceptance/no_direct_detection_no_hole_ice_plane_waves
[2018-08-01 20:56:16] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./plot_angular_acceptance.py --no-hole-ice --no-direct-detection ../../results/angular_acceptance/no_direct_detection_no_hole_ice_plane_waves
[2018-08-01 17:52:24] fiedl@kepler16 /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
[2018-08-01 18:03:54] fiedl@kepler00 /afs/ifh.de/group/amanda/scratch/fiedl/hole-ice-study/scripts/AngularAcceptance
▶ mv results/current ../../results/angular_acceptance/no_direct_detection_no_hole_ice_pencil_beam
▶ mv tmp ../../results/angular_acceptance/no_direct_detection_no_hole_ice_pencil_beam/
[2018-08-01 18:25:01] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./convert_options_txt_to_json.rb ../../results/angular_acceptance/no_direct_detection_no_hole_ice_pencil_beam
▶ ./plot_angular_acceptance.py --no-hole-ice ../../results/angular_acceptance/no_direct_detection_no_hole_ice_pencil_beam/
[2018-08-01 19:10:46] fiedl@fiedl-mbp ~/hole-ice-study/scripts/lib master ⚡
▶ ./plot_angular_acceptance.py --no-hole-ice --pencil-beam --no-direct-detection ../../results/angular_acceptance/no_direct_detection_no_hole_ice_pencil_beam
Perform an angular-acceptance scan without hole-ice cylinder. Use direct detection (#32)
Compare to reference curve (#10).