aboutsummaryrefslogtreecommitdiff
path: root/utils/plot-dc
diff options
context:
space:
mode:
Diffstat (limited to 'utils/plot-dc')
-rwxr-xr-xutils/plot-dc7
1 files changed, 6 insertions, 1 deletions
diff --git a/utils/plot-dc b/utils/plot-dc
index d2e461d..e20d5cc 100755
--- a/utils/plot-dc
+++ b/utils/plot-dc
@@ -28,7 +28,7 @@ from scipy.stats import iqr, poisson
from matplotlib.lines import Line2D
from scipy.stats import iqr, norm, beta
from scipy.special import spence
-from itertools import izip_longest
+from itertools import chain
particle_id = {20: 'e', 22: r'\mu'}
@@ -52,6 +52,11 @@ def plot_hist2(df, muons=False):
bins = np.logspace(np.log10(20),np.log10(10e3),21)
plt.hist(df_id.ke.values, bins=bins, histtype='step')
plt.gca().set_xscale("log")
+ major = np.array([10,100,1000,10000])
+ minor = np.unique(list(chain(*list(range(i,i*10,i) for i in major[:-1]))))
+ minor = np.setdiff1d(minor,major)
+ plt.gca().set_xticks(major)
+ plt.gca().set_xticks(minor,minor=True)
plt.xlabel("Energy (MeV)")
plt.title('$' + ''.join([particle_id[int(''.join(x))] for x in grouper(str(id),2)]) + '$')