Changeset 6060a86 in sasview for src/sas

Timestamp:

Oct 26, 2017 5:11:00 AM (7 years ago)

Author:

Adam Washington <adam.washington@…>

Branches:

ESS_GUI, ESS_GUI_Docs, ESS_GUI_batch_fitting, ESS_GUI_bumps_abstraction, ESS_GUI_iss1116, ESS_GUI_iss879, ESS_GUI_iss959, ESS_GUI_opencl, ESS_GUI_ordering, ESS_GUI_sync_sascalc

Children:

3e20577

Parents:

60bf0db

Message:

Plot propper values

File:

• src/sas/qtgui/Perspectives/PrInversion/dmax.py (modified) (4 diffs)

src/sas/qtgui/Perspectives/PrInversion/dmax.py

@@ -10 +10 @@
 import sys
 import os
+import logging
+import numpy as np
 from PyQt4 import QtCore
 from PyQt4 import QtGui
@@ -25 +27 @@
 # from InvariantDetails import DetailsDialog
 # from InvariantUtils import WIDGETS
+
+logger = logging.getLogger(__name__)
 
 def enum(*sequential, **named):
@@ -56 +60 @@
         # Let's choose the Standard Item Model.
         self.model = QtGui.QStandardItemModel(self)
+        self.mapper = None
 
         # # Connect buttons to slots.
@@ -94 +99 @@
 
     def modelChanged(self, item):
-        pass
+        if not self.mapper:
+            return
+
+        iq0 = []
+        rg = []
+        pos = []
+        pos_err = []
+        osc = []
+        xs = np.linspace(float(self.model.item(W.DMIN).text()),
+                         float(self.model.item(W.DMAX).text()),
+                         float(self.model.item(W.NPTS).text()))
+
+        for x in xs:
+            self.pr_state.d_max = x
+            try:
+                out, cov = self.pr_state.invert(self.nfunc)
+
+                iq0.append(self.pr_state.iq0(out))
+                rg.append(self.pr_state.rg(out))
+                pos.append(self.pr_state.get_positive(out))
+                pos_err.append(self.pr_state.get_pos_err(out, cov))
+                osc.append(self.pr_state.oscillations(out))
+            except:
+                # This inversion failed, skip this D_max value
+                msg = "ExploreDialog: inversion failed "
+                msg += "for D_max=%s\n%s" % (str(x), sys.exc_value)
+                print(msg)
+                # logger.error(msg)
+        self.plot.plot(data=Data1D(xs, rg), marker="-")