Tests sp correction

ohmpi/hardware_system.py

@@ -214,7 +214,7 @@ class OhmPiHardware:
     @property
     def last_dev(self):
         if len(self.readings) > 1:
-            v = self.readings[:, 2] != 0  # exclude sample where the is no injection
+            v = self.readings[:, 2] != 0  # exclude sample where there is no injection
             return 100. * np.std(self.readings[v, 2] * (self.readings[v, 4] - self.sp) / self.readings[v, 3]) / self.last_rho
         else:
             return np.nan
@@ -316,7 +316,7 @@ class OhmPiHardware:
 
     def _plot_readings(self, save_fig=False):
         # Plot graphs
-        fig, ax = plt.subplots(nrows=4, sharex=True)
+        fig, ax = plt.subplots(nrows=5, sharex=True)
         ax[0].plot(self.readings[:, 0], self.readings[:, 3], '-r', marker='.', label='iab')
         ax[0].set_ylabel('Iab [mA]')
         ax[1].plot(self.readings[:, 0], self.readings[:, 2] * (self.readings[:, 4] - self.sp) , '-b', marker='.', label='vmn')
@@ -327,6 +327,9 @@ class OhmPiHardware:
         ax[3].plot(self.readings[v, 0], (self.readings[v, 2] * (self.readings[v, 4] - self.sp)) / self.readings[v, 3],
                    '-m', marker='.', label='R [ohm]')
         ax[3].set_ylabel('R [ohm]')
+        ax[4].plot(self.readings[v, 0], (self.readings[v, 2] * (self.readings[v, 4] + self.sp)) / self.readings[v, 3],
+                   '-m', marker='.', label='R [ohm]')
+        ax[4].set_ylabel('R [ohm]')
         fig.legend()
         if save_fig:
             fig.savefig(f'figures/test.png')