def plot_afstelling_pmt():
"""Plot voor afstellen spanning"""
multiplot = MultiPlot(1, 3, width=r'0.4\linewidth')
x = pylab.linspace(1e-10, 11, 200)
signal = pylab.zeros(x.shape)
for N in range(1, 15):
scale = 100. / N ** 3
pdf = 80 * scale * pylab.normpdf(x, N, pylab.sqrt(N) * 0.35)
signal += pdf
gammas = 1e2 * x ** -3
signal += gammas
p = multiplot.get_subplot_at(0, 0)
v = x * 35
p.plot(v, pylab.where(v >= 30, signal, 1), mark=None)
p.draw_vertical_line(200, linestyle='gray')
p.set_label(r"$V_\mathrm{PMT}$ te laag")
p = multiplot.get_subplot_at(0, 1)
v = x * 200
p.plot(v, pylab.where(v >= 30, signal, 1), mark=None)
p.draw_vertical_line(200, linestyle='gray')
p.set_label(r"$V_\mathrm{PMT}$ correct")
p = multiplot.get_subplot_at(0, 2)
v = x * 400
p.plot(v, pylab.where(v >= 30, signal, 1), mark=None)
p.draw_vertical_line(200, linestyle='gray')
p.set_label(r"$V_\mathrm{PMT}$ te hoog")
multiplot.set_xlabel(r"Pulseheight [\si{\milli\volt}]")
multiplot.set_ylabel("Counts")
multiplot.set_xlimits_for_all(min=0, max=1000)
multiplot.set_ylimits_for_all(min=1)
multiplot.show_xticklabels_for_all()
multiplot.set_xticklabels_position(0, 1, 'top')
multiplot.set_yticks_for_all(ticks=None)
multiplot.save('afstelling_pmt')
def multiplot():
# data series
x = [0, 40, 60, 69, 80, 90, 100]
y = [0, 0, 0.5, 2.96, 2, 1, .5]
multiplot = MultiPlot(1, 2)
# make graph
graph1 = multiplot.get_subplot_at(0, 0)
graph2 = multiplot.get_subplot_at(0, 1)
# make Plot
graph1.plot(x, y, mark=None, linestyle='smooth,very thick')
graph2.plot(x, y, mark=None, linestyle='smooth,very thick')
multiplot.show_xticklabels_for_all()
multiplot.show_yticklabels(0, 0)
multiplot.set_xticklabels_position(0, 1, 'top')
multiplot.set_yticks_for_all(ticks=range(6))
graph1.set_xlimits(0, 100)
graph2.set_xlimits(60, 80)
multiplot.set_ylimits_for_all(min=0, max=5)
# set scale: 1cm equals 10 or 5 units along the x-axis
graph1.set_xscale(cm=10)
graph2.set_xscale(cm=5)
# set scale: 1cm equals 1 unit along the y-axis
graph1.set_yscale(cm=1)
graph2.set_yscale(cm=1)
# set graph paper
graph1.use_graph_paper()
graph2.use_graph_paper()
# set labels
multiplot.set_xlabel(r"$f [\si{\mega\hertz}]$")
multiplot.set_ylabel("signal strength")
# save graph to file
multiplot.save('mm_paper_multiplot')
# Plot voor afstellen spanning
multiplot = MultiPlot(1, 3, width=r'.4\linewidth')
p = multiplot.get_subplot_at(0, 1)
V = x * 200
p.plot(V, where(V >= 30, signal, 1), mark=None)
p.draw_vertical_line(200, linestyle='gray')
p.set_label(r"$V_\mathrm{PMT}$ correct")
p = multiplot.get_subplot_at(0, 2)
V = x * 400
p.plot(V, where(V >= 30, signal, 1), mark=None)
p.draw_vertical_line(200, linestyle='gray')
p.set_label(r"$V_\mathrm{PMT}$ te hoog")
p = multiplot.get_subplot_at(0, 0)
V = x * 35
p.plot(V, where(V >= 30, signal, 1), mark=None)
p.draw_vertical_line(200, linestyle='gray')
p.set_label(r"$V_\mathrm{PMT}$ te laag")
multiplot.set_xlabel(r"Pulseheight [\si{\milli\volt}]")
multiplot.set_ylabel("Counts")
multiplot.set_xlimits_for_all(min=0, max=1000)
multiplot.set_ylimits_for_all(min=1)
multiplot.show_xticklabels_for_all()
multiplot.set_xticklabels_position(0, 1, 'top')
multiplot.set_yticks_for_all(ticks=None)
multiplot.save('afstelling_pmt')
