aboutsummaryrefslogtreecommitdiffstats
path: root/python
diff options
context:
space:
mode:
authorAaron LI <aaronly.me@outlook.com>2016-06-07 22:26:47 +0800
committerAaron LI <aaronly.me@outlook.com>2016-06-07 22:26:47 +0800
commit63ccfae1b5f4ff9c2b43f206c6b350b123376204 (patch)
treef65f28721d26fde65b48ebe1467bf525f1827857 /python
parent34ed04686523ecdafb96c6cdff379df176cf1984 (diff)
downloadatoolbox-63ccfae1b5f4ff9c2b43f206c6b350b123376204.tar.bz2
crosstalk_deprojection.py: explain error calculation approach
Diffstat (limited to 'python')
-rwxr-xr-xpython/crosstalk_deprojection.py14
1 files changed, 9 insertions, 5 deletions
diff --git a/python/crosstalk_deprojection.py b/python/crosstalk_deprojection.py
index d5bab05..b08a66a 100755
--- a/python/crosstalk_deprojection.py
+++ b/python/crosstalk_deprojection.py
@@ -21,9 +21,11 @@
#
# Weitian LI
# Created: 2016-03-26
-# Updated: 2016-04-20
+# Updated: 2016-06-07
#
-# ChangeLog:
+# Change log:
+# 2016-06-07:
+# * Explain the errors/uncertainties calculation approach
# 2016-04-20:
# * Add argument 'add_history' to some methods (to avoid many duplicated
# histories due to Monte Carlo)
@@ -65,8 +67,8 @@
# * Split classes ARF, RMF, Spectrum, and SpectrumSet to a separate module
#
-__version__ = "0.5.2"
-__date__ = "2016-04-20"
+__version__ = "0.5.3"
+__date__ = "2016-06-07"
"""
@@ -1421,7 +1423,9 @@ def calc_median_errors(results):
Calculate the median and errors for the spectral data gathered
through Monte Carlo simulations.
- TODO: investigate the errors calculation approach used here!
+ NOTE:
+ Errors calculation just use the quantiles,
+ i.e., 1sigma ~= 68.3% = Q(84.15%) - Q(15.85%)
"""
results = np.array(results)
# `results' now has shape: (mc_times, num_spec, num_channel)