From 57e59ea4d8f22c1c28e7d59bdaeecbdfba68f12d Mon Sep 17 00:00:00 2001
From: su530201 <olivier.kaufmann@umons.ac.be>
Date: Tue, 26 Sep 2023 08:30:54 +0200
Subject: [PATCH] Tests sp correction

---
 ohmpi/hardware_system.py | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

diff --git a/ohmpi/hardware_system.py b/ohmpi/hardware_system.py
index f7934e54..b5ab0d85 100644
--- a/ohmpi/hardware_system.py
+++ b/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')
-- 
GitLab