Index: branches/trigger_burst_research/plot_more_stuff.py
===================================================================
--- branches/trigger_burst_research/plot_more_stuff.py	(revision 18306)
+++ branches/trigger_burst_research/plot_more_stuff.py	(revision 18306)
@@ -0,0 +1,188 @@
+# coding: utf-8
+from pylab import *
+import pandas as pd
+import matplotlib.dates as mdates
+from matplotlib.colors import LogNorm
+
+r = pd.DataFrame.from_csv('burst_ratio.csv', 
+	infer_datetime_format=True, 
+	parse_dates=['fRunStart', 'fRunStop'])
+r = r.dropna()
+
+t = r.fRunStart.values.astype('float64')/1e9
+x = r.fBoardTriggerRateRatioAboveThreshold
+tt = mdates.epoch2num(t)
+
+plt.ion()
+fig, axes = plt.subplots(3,1, sharex=True)
+fig.subplots_adjust(hspace=0.)
+for ax in axes:
+	ax.xaxis.set_major_locator(mdates.MonthLocator())
+	ax.xaxis.set_major_formatter(mdates.DateFormatter('%Y-%m-%d'))
+	ax.xaxis.set_minor_locator( mdates.DayLocator())
+
+plot0 = axes[0].hist2d(tt, x.values, 
+	cmap='Greys', 
+	bins=(np.arange(int(tt[0])-2, int(tt[-1])+2), 
+		np.arange(0, 0.6+1/60., 1./60.)), 
+	norm=LogNorm())
+#plt.colorbar(plot0[-1])
+axes[0].set_ylabel("Ratio")
+
+plot1 = axes[1].hist2d(tt, r.fZenithDistanceMean.values,
+	cmap='Greys', 
+	bins=(np.arange(int(tt[0])-2, int(tt[-1])+2), 
+		np.arange(0, 95, 1.)), 
+	norm=LogNorm())
+#plt.colorbar(plot1[-1])
+axes[1].set_ylabel("Zenith [deg]")
+plot2 = axes[2].hist2d(tt, r.fAzimuthMean.values, 
+	cmap='Greys', 
+	bins=(np.arange(int(tt[0])-2, int(tt[-1])+2), 
+		np.arange(-270, 90 , 4)), 
+	norm=LogNorm())
+#plt.colorbar(plot2[-1])
+axes[2].set_ylabel("Azimut [deg]")
+#ax.xaxis.set_major_formatter(mdates.AutoDateFormatter(mdates.AutoDateLocator()))
+#ax.xaxis.set_major_locator(mdates.AutoDateLocator())
+
+#ax.format_xdata = mdates.DateFormatter('%Y-%m-%d %H:%M')
+fig.autofmt_xdate()
+plt.grid()
+
+plt.xlabel("Time")
+#plt.tight_layout()
+plt.draw()
+
+"""
+H, b = np.histogram(
+	r.fZenithDistanceMean.values, 
+	bins = np.arange(0, 180))
+
+H2, bx, by = np.histogram2d(
+	r.fZenithDistanceMean.values, 
+	x.values,
+	bins=(np.arange(0, 180), np.arange(0, 0.6, 1/60.)),
+	)
+
+plt.figure()
+plt.title("normiert")
+plt.imshow(
+	(H2.T/H).T,
+	origin='upper',
+	#extent=[bx[0],bx[-1],by[0],by[-1]], 
+	#aspect=180,
+	cmap='Greys',
+	)
+
+"""
+H, b = np.histogram(
+	r.fZenithDistanceMean.values, 
+	bins = np.arange(0, 180))
+
+zenit_bin_index = np.digitize(r.fZenithDistanceMean.values, np.arange(0, 180))-1
+w = 1./H[zenit_bin_index]
+
+plt.figure()
+plt.hist2d(
+	r.fZenithDistanceMean.values, 
+	x.values, 
+	weights=w,
+	#cmap='Greys',
+	bins = (np.arange(0, 180), np.arange(0, 0.6, 1./60.)),
+	norm=LogNorm(),
+	)
+plt.xlabel("fZenithDistanceMean")
+plt.xlim(0, 90)
+plt.ylabel("fBoardTriggerRateRatioAboveThreshold")
+plt.title("fZenithDistanceMean Dependence of Ratio")
+plt.colorbar()
+plt.grid()
+
+H, b = np.histogram(
+	r.fAzimuthMean.values, 
+	bins = np.arange(-360, 360))
+
+az_bin_index = np.digitize(r.fAzimuthMean.values, np.arange(-360, 360))-1
+w = 1./H[az_bin_index]
+
+plt.figure()
+plt.hist2d(
+	r.fAzimuthMean.values, 
+	x.values, 
+	weights=w,
+	#cmap='Greys', 
+	bins = (np.arange(-360, 360), np.arange(0, 0.6, 1./60.)),
+	norm=LogNorm()
+	)
+plt.xlabel("fAzimuthMean")
+plt.ylabel("fBoardTriggerRateRatioAboveThreshold")
+plt.title("Azimuth Dependence of Ratio")
+plt.colorbar()
+plt.grid()
+
+nbins = (
+	np.arange(-2*np.pi, 2*np.pi, np.pi/180*5), 
+	np.arange(0, 90)
+	)
+
+H1, bx, by = np.histogram2d(
+	x=r.fAzimuthMean.values*np.pi/180.,
+	y=r.fZenithDistanceMean.values,
+	weights=x.values,
+	bins=nbins)
+
+H2, bx, by = np.histogram2d(
+	x=r.fAzimuthMean.values*np.pi/180.,
+	y=r.fZenithDistanceMean.values,
+	bins=nbins)
+
+normed_H1 = H1/H2
+
+above_thr = normed_H1 > 0.05
+
+runs_above_thr = H2.copy()
+runs_above_thr[~above_thr]=0
+
+plt.figure()
+plt.title("Ratio vs. zenith and azimut (normalized for #runs)")
+im = plt.imshow(normed_H1.T,
+	origin='low',
+	interpolation='none',
+	extent=[bx[0], bx[-1], by[0], by[-1]],
+	aspect='auto',
+	norm=LogNorm())
+#im.set_cmap('Greys')
+cbar = plt.colorbar()
+plt.grid()
+plt.xlim(-1.5*np.pi, np.pi/2.)
+plt.xlabel("Azimut [rad]")
+plt.ylabel("Zenith [deg]")
+
+
+
+plt.figure()
+plt.title("Above 5% Ratio")
+im = plt.imshow(above_thr.T,
+	origin='low',
+	interpolation='none',
+	extent=[bx[0], bx[-1], by[0], by[-1]],
+	aspect='auto',
+	norm=LogNorm())
+#im.set_cmap('Greys')
+plt.colorbar()
+plt.grid()
+plt.xlim(-1.5*np.pi, np.pi/2.)
+
+plt.figure()
+plt.title("# runs above thr")
+im = plt.imshow(runs_above_thr.T,
+	origin='low',
+	interpolation='none',
+	extent=[bx[0], bx[-1], by[0], by[-1]],
+	aspect='auto')
+#im.set_cmap('Greys')
+plt.colorbar()
+plt.grid()
+plt.xlim(-1.5*np.pi, np.pi/2.)
+