output-skew figure script

72090aae · Anders Wallin · fca5283e · 72090aae
Commit 72090aae authored May 24, 2021 by Anders Wallin
Hide whitespace changes
Inline Side-by-side

Showing with 72 additions and 0 deletions

pda_output_skew.py pda_output_skew.py +72 -0

No files found.
--- a/pda_output_skew.py
+++ b/pda_output_skew.py
+"""
+    measurement data for
+    Pulse distribution amplifier, https://ohwr.org/project/pda-8ch-fda-8ch
+"""
+import numpy
+import matplotlib.pyplot as plt
+sn12_b = [14526, 14542, 14313, 14263, 14313, 14321, 14420, 14422]
+sn13_b = [14399, 14278, 14153, 14081, 14113, 14185, 14241, 14278]
+sn16_b = [13505, 13491, 13302, 13198, 13284, 13322, 13453, 13374]
+sn15_b = [13310, 13230, 13101, 13036, 13057, 13070, 13236, 13210]
+sn14_b = [13347, 13274, 13108, 13029, 13094, 13119, 13219, 13210]
+sn17_b = [17172, 17109, 16922, 16876, 16943, 16928, 17047, 17036 ]
+sn22_b = [12833, 12780, 12585, 12558, 12606, 12593, 12730, 12681]
+sn23_b = [337, 465, 331, 234, 367, 243, 406, 262]
+sn24_b = [16322, 16180, 16051, 16028, 16228, 16017, 16153, 16038]
+sn27l_b = [724, 577, 403, 395, 471, 444, 534, 533 ]
+sn27_b = [677, 723, 576, 524, 630, 514, 712, 666]
+sn28_b = [18505, 18743, 18664, 18557, 18459, 18739, 18545, 18847 ]
+sn30_b = [635, 620, 519, 395, 652, 446, 668, 807]
+sn33_b = [ 19008, 19091, 18869, 18723, 19125, 18865, 19124, 18859 ]
+sn36 = [ 16699, 16570, 16481, 16464, 16705, 16489, 16724, 16638 ]
+sn37 = [ 618, 427, 627, 485, 656, 635, 625, 592 ]
+before = []
+before.append( sn37 )
+before.append( sn36 )
+before.append( sn33_b )
+before.append( sn30_b )
+before.append( sn28_b )
+before.append( sn27_b )
+before.append( sn27l_b )
+before.append( sn24_b )
+before.append( sn23_b )
+before.append( sn22_b )
+before.append( sn17_b )
+before.append( sn14_b )
+before.append( sn15_b )
+before.append( sn16_b )
+before.append( sn13_b )
+before.append( sn12_b )
+plt.figure()
+d = [ [], [], [], [], [], [], [], []]
+for b in before:
+    b = numpy.array(b)
+    plt.plot( b-numpy.average(b) , 'o')
+    for idx, dd in enumerate(b-numpy.average(b)):
+        d[idx].append(dd)
+print d
+mu=[]
+s=[]
+for dd in d:
+    mu.append( numpy.mean(dd) )
+    s.append( numpy.std(dd) )
+plt.errorbar( range(8), mu, yerr=s)
+i=0
+for (m,std) in zip(mu, s):
+    plt.text(i-0.25, 0-i*20, "%d +/- %d" %(m, std))
+    i=i+1
+print(mu)
+print(s)
+plt.title('Pulse distribution amplifier, output-to-output skew\n AW2021-05-07, https://ohwr.org/project/pda-8ch-fda-8ch')
+plt.xlabel('Output #')
+plt.ylabel('Output skew / ps')
+plt.grid()
+plt.show()