Modern Lab Fixes

statics/physics/modern-phys-lab/counting/run.py

@@ -1,9 +1,10 @@
 import numpy as np
+import matplotlib.pyplot as plt
 
 def read_file(file: str) -> np.typing.ArrayLike:
     with open(f'{file}.csv', 'r') as file:
         lines = file.read().split('Channel,Energy,Counts')[1].strip().split('\n')
-        return np.array([int(line.split(',,')[1]) for line in lines])
+        return np.array([int(line.split(',,')[1]) for line in lines], dtype=np.uint16)
 
 all_data = {
     '200': read_file('200'),
@@ -11,10 +12,44 @@ all_data = {
     '40': sum([read_file(f'40-{i}') for i in range(1, 6)])
 }
 
-for time in all_data:
+fig, axs = plt.subplots(len(all_data), sharex='all')
-    data = all_data[time]
+i = 0
+
+for dwell_time in all_data:
+    print(f'\n==[{dwell_time}ms]==\n')
+    data = all_data[dwell_time]
+
+    N = len(data)
+    print(f'Loaded {sum(data)} events across {N} samples')
 
     mean = np.mean(data)
-    print(f'Found sample mean {mean}')
+    print(f'Found sample mean {mean:0.2f}')
 
-    stdev = np.stdev(mean, ddof=1)
+    stdev = np.std(data, ddof=1)
+    print(f'Found sample standard deviation {stdev:0.2f}')
+
+    sigma = np.sqrt(mean)
+    print(f'Found sigma {sigma:0.2f}')
+
+    P = 0.6745 * sigma
+    print(f'Found P {P:0.2f}')
+    
+    times_dev_more_than_s = (np.abs(data - mean) > sigma).sum()
+    print(f'Found {times_dev_more_than_s} / {N} ({times_dev_more_than_s / N * 100:0.2f}%) samples deviating from the mean more than sigma {sigma:0.2f}')
+    
+    times_dev_more_than_P = (np.abs(data - mean) > P).sum()
+    print(f'Found {times_dev_more_than_P} / {N} ({times_dev_more_than_P / N * 100:0.2f}%) samples deviating from the mean more than P {P:0.2f}')
+
+    ax = axs[i]
+    i += 1
+
+    ax.set_title(f'Cumulative Average ({dwell_time}ms Dwell time)')
+    ax.set_xlabel('After N Runs')
+    ax.set_ylabel('Cumulative Average')
+
+    sample_num = np.arange(1, N + 1)
+    cumulative_average = np.cumsum(data) / sample_num
+    ax.plot(sample_num, cumulative_average, label=f'Final Value: {cumulative_average[-1]:0.2f}')
+    ax.legend()
+
+plt.show()

statics/physics/modern-phys-lab/diffraction.py

@@ -2,17 +2,17 @@ import numpy as np
 import matplotlib.pyplot as plt
 
 # Units: V
-voltages = np.array(np.arange(5000, 2500 - 1, -500)) #-1 to include 2500
+voltages = np.array(np.arange(2500, 5000 + 1, +500)) #+1 to include 5000
 voltages_inv_sqrt = 1 / np.sqrt(voltages)
 
 # New units: m
-diameter_measured_error = 0.02
+diameter_measured_error = 0.002
 
 # Average inner and outer
 # New format: d[0 = small, 1 = large][voltage]
 diameter_measured = np.array([
-    [0.0, 0.024, 0.022, 0.019, 0.02],
+    [0.026, 0.0235, 0.024, 0.022, 0.019, 0.02],
-    [0.0, 0.04, 0.0385, 0.036, 0.035]
+    [0.049, 0.0435, 0.04, 0.0385, 0.036, 0.035]
 ])
 
 diameter_error = 100