1 | # global |
---|
2 | import sys |
---|
3 | import os |
---|
4 | import logging |
---|
5 | import copy |
---|
6 | import webbrowser |
---|
7 | |
---|
8 | from PyQt5 import QtCore |
---|
9 | from PyQt5 import QtGui, QtWidgets |
---|
10 | |
---|
11 | from twisted.internet import threads |
---|
12 | from twisted.internet import reactor |
---|
13 | |
---|
14 | # sas-global |
---|
15 | from sas.sascalc.invariant import invariant |
---|
16 | from sas.qtgui.Plotting.PlotterData import Data1D |
---|
17 | import sas.qtgui.Utilities.GuiUtils as GuiUtils |
---|
18 | |
---|
19 | # import sas.qtgui.Plotting.PlotHelper as PlotHelper |
---|
20 | |
---|
21 | # local |
---|
22 | from .UI.TabbedInvariantUI import Ui_tabbedInvariantUI |
---|
23 | from .InvariantDetails import DetailsDialog |
---|
24 | from .InvariantUtils import WIDGETS |
---|
25 | |
---|
26 | # The minimum q-value to be used when extrapolating |
---|
27 | Q_MINIMUM = 1e-5 |
---|
28 | # The maximum q-value to be used when extrapolating |
---|
29 | Q_MAXIMUM = 10 |
---|
30 | # the ratio of maximum q value/(qmax of data) to plot the theory data |
---|
31 | Q_MAXIMUM_PLOT = 3 |
---|
32 | # Default number of points of interpolation: high and low range |
---|
33 | NPOINTS_Q_INTERP = 10 |
---|
34 | # Default power law for interpolation |
---|
35 | DEFAULT_POWER_LOW = 4 |
---|
36 | |
---|
37 | # Background of line edits if settings OK or wrong |
---|
38 | BG_WHITE = "background-color: rgb(255, 255, 255);" |
---|
39 | BG_RED = "background-color: rgb(244, 170, 164);" |
---|
40 | |
---|
41 | class InvariantWindow(QtWidgets.QDialog, Ui_tabbedInvariantUI): |
---|
42 | # The controller which is responsible for managing signal slots connections |
---|
43 | # for the gui and providing an interface to the data model. |
---|
44 | name = "Invariant" # For displaying in the combo box in DataExplorer |
---|
45 | |
---|
46 | def __init__(self, parent=None): |
---|
47 | super(InvariantWindow, self).__init__() |
---|
48 | self.setupUi(self) |
---|
49 | |
---|
50 | self.setWindowTitle("Invariant Perspective") |
---|
51 | |
---|
52 | # initial input params |
---|
53 | self._background = 0.0 |
---|
54 | self._scale = 1.0 |
---|
55 | self._contrast = 1.0 |
---|
56 | self._porod = None |
---|
57 | |
---|
58 | self.parent = parent |
---|
59 | |
---|
60 | self._manager = parent |
---|
61 | self._reactor = reactor |
---|
62 | self._model_item = QtGui.QStandardItem() |
---|
63 | |
---|
64 | self.detailsDialog = DetailsDialog(self) |
---|
65 | self.detailsDialog.cmdOK.clicked.connect(self.enabling) |
---|
66 | |
---|
67 | self._low_extrapolate = False |
---|
68 | self._low_guinier = True |
---|
69 | self._low_fit = False |
---|
70 | self._low_power_value = False |
---|
71 | self._low_points = NPOINTS_Q_INTERP |
---|
72 | self._low_power_value = DEFAULT_POWER_LOW |
---|
73 | |
---|
74 | self._high_extrapolate = False |
---|
75 | self._high_power_value = False |
---|
76 | self._high_fit = False |
---|
77 | self._high_points = NPOINTS_Q_INTERP |
---|
78 | self._high_power_value = DEFAULT_POWER_LOW |
---|
79 | |
---|
80 | # no reason to have this widget resizable |
---|
81 | self.resize(self.minimumSizeHint()) |
---|
82 | |
---|
83 | self.communicate = self._manager.communicator() |
---|
84 | |
---|
85 | self._data = None |
---|
86 | self._path = "" |
---|
87 | |
---|
88 | self._allow_close = False |
---|
89 | |
---|
90 | # Modify font in order to display Angstrom symbol correctly |
---|
91 | new_font = 'font-family: -apple-system, "Helvetica Neue", "Ubuntu";' |
---|
92 | self.lblTotalQUnits.setStyleSheet(new_font) |
---|
93 | self.lblSpecificSurfaceUnits.setStyleSheet(new_font) |
---|
94 | self.lblInvariantTotalQUnits.setStyleSheet(new_font) |
---|
95 | self.lblContrastUnits.setStyleSheet(new_font) |
---|
96 | self.lblPorodCstUnits.setStyleSheet(new_font) |
---|
97 | self.lblExtrapolQUnits.setStyleSheet(new_font) |
---|
98 | |
---|
99 | # To remove blue square around line edits |
---|
100 | self.txtBackgd.setAttribute(QtCore.Qt.WA_MacShowFocusRect, False) |
---|
101 | self.txtContrast.setAttribute(QtCore.Qt.WA_MacShowFocusRect, False) |
---|
102 | self.txtScale.setAttribute(QtCore.Qt.WA_MacShowFocusRect, False) |
---|
103 | self.txtPorodCst.setAttribute(QtCore.Qt.WA_MacShowFocusRect, False) |
---|
104 | self.txtNptsHighQ.setAttribute(QtCore.Qt.WA_MacShowFocusRect, False) |
---|
105 | self.txtNptsLowQ.setAttribute(QtCore.Qt.WA_MacShowFocusRect, False) |
---|
106 | self.txtPowerLowQ.setAttribute(QtCore.Qt.WA_MacShowFocusRect, False) |
---|
107 | self.txtPowerHighQ.setAttribute(QtCore.Qt.WA_MacShowFocusRect, False) |
---|
108 | |
---|
109 | self.txtExtrapolQMin.setText(str(Q_MINIMUM)) |
---|
110 | self.txtExtrapolQMax.setText(str(Q_MAXIMUM)) |
---|
111 | |
---|
112 | # Let's choose the Standard Item Model. |
---|
113 | self.model = QtGui.QStandardItemModel(self) |
---|
114 | |
---|
115 | # Connect buttons to slots. |
---|
116 | # Needs to be done early so default values propagate properly. |
---|
117 | self.setupSlots() |
---|
118 | |
---|
119 | # Set up the model. |
---|
120 | self.setupModel() |
---|
121 | |
---|
122 | # Set up the mapper |
---|
123 | self.setupMapper() |
---|
124 | |
---|
125 | # Default enablement |
---|
126 | self.cmdCalculate.setEnabled(False) |
---|
127 | |
---|
128 | # validator: double |
---|
129 | self.txtBackgd.setValidator(GuiUtils.DoubleValidator()) |
---|
130 | self.txtContrast.setValidator(GuiUtils.DoubleValidator()) |
---|
131 | self.txtScale.setValidator(GuiUtils.DoubleValidator()) |
---|
132 | self.txtPorodCst.setValidator(GuiUtils.DoubleValidator()) |
---|
133 | |
---|
134 | # validator: integer number |
---|
135 | self.txtNptsLowQ.setValidator(QtGui.QIntValidator()) |
---|
136 | self.txtNptsHighQ.setValidator(QtGui.QIntValidator()) |
---|
137 | self.txtPowerLowQ.setValidator(GuiUtils.DoubleValidator()) |
---|
138 | self.txtPowerHighQ.setValidator(GuiUtils.DoubleValidator()) |
---|
139 | |
---|
140 | def enabling(self): |
---|
141 | """ """ |
---|
142 | self.cmdStatus.setEnabled(True) |
---|
143 | |
---|
144 | def setClosable(self, value=True): |
---|
145 | """ Allow outsiders close this widget """ |
---|
146 | assert isinstance(value, bool) |
---|
147 | |
---|
148 | self._allow_close = value |
---|
149 | |
---|
150 | def closeEvent(self, event): |
---|
151 | """ |
---|
152 | Overwrite QDialog close method to allow for custom widget close |
---|
153 | """ |
---|
154 | if self._allow_close: |
---|
155 | # reset the closability flag |
---|
156 | self.setClosable(value=False) |
---|
157 | # Tell the MdiArea to close the container |
---|
158 | self.parentWidget().close() |
---|
159 | event.accept() |
---|
160 | else: |
---|
161 | event.ignore() |
---|
162 | # Maybe we should just minimize |
---|
163 | self.setWindowState(QtCore.Qt.WindowMinimized) |
---|
164 | |
---|
165 | def updateFromModel(self): |
---|
166 | """ Update the globals based on the data in the model """ |
---|
167 | self._background = float(self.model.item(WIDGETS.W_BACKGROUND).text()) |
---|
168 | self._contrast = float(self.model.item(WIDGETS.W_CONTRAST).text()) |
---|
169 | self._scale = float(self.model.item(WIDGETS.W_SCALE).text()) |
---|
170 | if self.model.item(WIDGETS.W_POROD_CST).text() != 'None' \ |
---|
171 | and self.model.item(WIDGETS.W_POROD_CST).text() != '': |
---|
172 | self._porod = float(self.model.item(WIDGETS.W_POROD_CST).text()) |
---|
173 | |
---|
174 | # Low extrapolating |
---|
175 | self._low_extrapolate = str(self.model.item(WIDGETS.W_ENABLE_LOWQ).text()) == 'true' |
---|
176 | self._low_points = float(self.model.item(WIDGETS.W_NPTS_LOWQ).text()) |
---|
177 | self._low_guinier = str(self.model.item(WIDGETS.W_LOWQ_GUINIER).text()) == 'true' |
---|
178 | self._low_fit = str(self.model.item(WIDGETS.W_LOWQ_FIT).text()) == 'true' |
---|
179 | self._low_power_value = float(self.model.item(WIDGETS.W_LOWQ_POWER_VALUE).text()) |
---|
180 | |
---|
181 | # High extrapolating |
---|
182 | self._high_extrapolate = str(self.model.item(WIDGETS.W_ENABLE_HIGHQ).text()) == 'true' |
---|
183 | self._high_points = float(self.model.item(WIDGETS.W_NPTS_HIGHQ).text()) |
---|
184 | self._high_fit = str(self.model.item(WIDGETS.W_HIGHQ_FIT).text()) == 'true' |
---|
185 | self._high_power_value = float(self.model.item(WIDGETS.W_HIGHQ_POWER_VALUE).text()) |
---|
186 | |
---|
187 | def calculateInvariant(self): |
---|
188 | """ Use twisted to thread the calculations away """ |
---|
189 | # Find out if extrapolation needs to be used. |
---|
190 | extrapolation = None |
---|
191 | if self._low_extrapolate and not self._high_extrapolate: |
---|
192 | extrapolation = "low" |
---|
193 | elif not self._low_extrapolate and self._high_extrapolate: |
---|
194 | extrapolation = "high" |
---|
195 | elif self._low_extrapolate and self._high_extrapolate: |
---|
196 | extrapolation = "both" |
---|
197 | |
---|
198 | # modify the Calculate button to indicate background process |
---|
199 | self.cmdCalculate.setText("Calculating...") |
---|
200 | self.cmdCalculate.setEnabled(False) |
---|
201 | |
---|
202 | # Send the calculations to separate thread. |
---|
203 | d = threads.deferToThread(self.calculateThread, extrapolation) |
---|
204 | |
---|
205 | # Add deferred callback for call return |
---|
206 | d.addCallback(self.deferredPlot) |
---|
207 | d.addErrback(self.calculationFailed) |
---|
208 | |
---|
209 | def calculationFailed(self, reason): |
---|
210 | print("calculation failed: ", reason) |
---|
211 | pass |
---|
212 | |
---|
213 | def deferredPlot(self, model): |
---|
214 | """ |
---|
215 | Run the GUI/model update in the main thread |
---|
216 | """ |
---|
217 | reactor.callFromThread(lambda: self.plotResult(model)) |
---|
218 | |
---|
219 | def plotResult(self, model): |
---|
220 | """ Plot result of calculation """ |
---|
221 | # Set the button back to available |
---|
222 | self.cmdCalculate.setEnabled(True) |
---|
223 | self.cmdCalculate.setText("Calculate") |
---|
224 | self.cmdStatus.setEnabled(True) |
---|
225 | |
---|
226 | self.model = model |
---|
227 | self.mapper.toFirst() |
---|
228 | self._data = GuiUtils.dataFromItem(self._model_item) |
---|
229 | |
---|
230 | # Send the modified model item to DE for keeping in the model |
---|
231 | # Currently -unused |
---|
232 | # self.communicate.updateModelFromPerspectiveSignal.emit(self._model_item) |
---|
233 | |
---|
234 | plot_data = GuiUtils.plotsFromCheckedItems(self._manager.filesWidget.model) |
---|
235 | #self.communicate.plotRequestedSignal.emit([item, plot], self.tab_id) |
---|
236 | |
---|
237 | self._manager.filesWidget.plotData(plot_data) |
---|
238 | |
---|
239 | # Update the details dialog in case it is open |
---|
240 | self.updateDetailsWidget(model) |
---|
241 | |
---|
242 | def updateDetailsWidget(self, model): |
---|
243 | """ |
---|
244 | On demand update of the details widget |
---|
245 | """ |
---|
246 | if self.detailsDialog.isVisible(): |
---|
247 | self.onStatus() |
---|
248 | |
---|
249 | def calculateThread(self, extrapolation): |
---|
250 | """ |
---|
251 | Perform Invariant calculations. |
---|
252 | TODO: Create a dictionary of results to be sent to DE on completion. |
---|
253 | """ |
---|
254 | self.updateFromModel() |
---|
255 | msg = '' |
---|
256 | |
---|
257 | qstar_low = 0.0 |
---|
258 | qstar_low_err = 0.0 |
---|
259 | qstar_high = 0.0 |
---|
260 | qstar_high_err = 0.0 |
---|
261 | qstar_total = 0.0 |
---|
262 | qstar_total_error = 0.0 |
---|
263 | |
---|
264 | temp_data = copy.deepcopy(self._data) |
---|
265 | |
---|
266 | # Prepare the invariant object |
---|
267 | inv = invariant.InvariantCalculator(data=temp_data, |
---|
268 | background=self._background, |
---|
269 | scale=self._scale) |
---|
270 | if self._low_extrapolate: |
---|
271 | |
---|
272 | function_low = "power_law" |
---|
273 | if self._low_guinier: |
---|
274 | function_low = "guinier" |
---|
275 | if self._low_fit: |
---|
276 | self._low_power_value = None |
---|
277 | |
---|
278 | inv.set_extrapolation(range="low", |
---|
279 | npts=int(self._low_points), |
---|
280 | function=function_low, |
---|
281 | power=self._low_power_value) |
---|
282 | |
---|
283 | if self._high_extrapolate: |
---|
284 | function_low = "power_law" |
---|
285 | inv.set_extrapolation(range="high", |
---|
286 | npts=int(self._high_points), |
---|
287 | function=function_low, |
---|
288 | power=self._high_power_value) |
---|
289 | # Compute invariant |
---|
290 | calculation_failed = False |
---|
291 | |
---|
292 | try: |
---|
293 | qstar_total, qstar_total_error = inv.get_qstar_with_error() |
---|
294 | except Exception as ex: |
---|
295 | msg += str(ex) |
---|
296 | calculation_failed = True |
---|
297 | # Display relevant information |
---|
298 | item = QtGui.QStandardItem("ERROR") |
---|
299 | self.model.setItem(WIDGETS.W_INVARIANT, item) |
---|
300 | item = QtGui.QStandardItem("ERROR") |
---|
301 | self.model.setItem(WIDGETS.W_INVARIANT_ERR, item) |
---|
302 | |
---|
303 | try: |
---|
304 | volume_fraction, volume_fraction_error = \ |
---|
305 | inv.get_volume_fraction_with_error(self._contrast) |
---|
306 | |
---|
307 | except Exception as ex: |
---|
308 | calculation_failed = True |
---|
309 | msg += str(ex) |
---|
310 | # Display relevant information |
---|
311 | item = QtGui.QStandardItem("ERROR") |
---|
312 | self.model.setItem(WIDGETS.W_VOLUME_FRACTION, item) |
---|
313 | item = QtGui.QStandardItem("ERROR") |
---|
314 | self.model.setItem(WIDGETS.W_VOLUME_FRACTION_ERR, item) |
---|
315 | |
---|
316 | if self._porod: |
---|
317 | try: |
---|
318 | surface, surface_error = \ |
---|
319 | inv.get_surface_with_error(self._contrast, self._porod) |
---|
320 | except Exception as ex: |
---|
321 | calculation_failed = True |
---|
322 | msg += str(ex) |
---|
323 | # Display relevant information |
---|
324 | item = QtGui.QStandardItem("ERROR") |
---|
325 | self.model.setItem(WIDGETS.W_SPECIFIC_SURFACE, item) |
---|
326 | item = QtGui.QStandardItem("ERROR") |
---|
327 | self.model.setItem(WIDGETS.W_SPECIFIC_SURFACE_ERR, item) |
---|
328 | else: |
---|
329 | surface = None |
---|
330 | |
---|
331 | if (calculation_failed): |
---|
332 | logging.warning('Calculation failed: {}'.format(msg)) |
---|
333 | return self.model |
---|
334 | else: |
---|
335 | |
---|
336 | if self._low_extrapolate: |
---|
337 | # for presentation in InvariantDetails |
---|
338 | qstar_low, qstar_low_err = inv.get_qstar_low() |
---|
339 | extrapolated_data = inv.get_extra_data_low(self._low_points) |
---|
340 | power_low = inv.get_extrapolation_power(range='low') |
---|
341 | |
---|
342 | # Plot the chart |
---|
343 | title = "Low-Q extrapolation" |
---|
344 | |
---|
345 | # Convert the data into plottable |
---|
346 | extrapolated_data = self._manager.createGuiData(extrapolated_data) |
---|
347 | |
---|
348 | extrapolated_data.name = title |
---|
349 | extrapolated_data.title = title |
---|
350 | extrapolated_data.style = "Line" |
---|
351 | extrapolated_data.has_errors = False |
---|
352 | extrapolated_data.plot_role = Data1D.ROLE_DEFAULT |
---|
353 | |
---|
354 | # copy labels and units of axes for plotting |
---|
355 | extrapolated_data._xaxis = temp_data._xaxis |
---|
356 | extrapolated_data._xunit = temp_data._xunit |
---|
357 | extrapolated_data._yaxis = temp_data._yaxis |
---|
358 | extrapolated_data._yunit = temp_data._yunit |
---|
359 | |
---|
360 | # Add the plot to the model item |
---|
361 | # This needs to run in the main thread |
---|
362 | reactor.callFromThread(GuiUtils.updateModelItemWithPlot, |
---|
363 | self._model_item, |
---|
364 | extrapolated_data, |
---|
365 | title) |
---|
366 | |
---|
367 | if self._high_extrapolate: |
---|
368 | # for presentation in InvariantDetails |
---|
369 | qmax_plot = Q_MAXIMUM_PLOT * max(temp_data.x) |
---|
370 | |
---|
371 | if qmax_plot > Q_MAXIMUM: |
---|
372 | qmax_plot = Q_MAXIMUM |
---|
373 | qstar_high, qstar_high_err = inv.get_qstar_high() |
---|
374 | power_high = inv.get_extrapolation_power(range='high') |
---|
375 | high_out_data = inv.get_extra_data_high(q_end=qmax_plot, npts=500) |
---|
376 | |
---|
377 | # Plot the chart |
---|
378 | title = "High-Q extrapolation" |
---|
379 | |
---|
380 | # Convert the data into plottable |
---|
381 | high_out_data = self._manager.createGuiData(high_out_data) |
---|
382 | high_out_data.name = title |
---|
383 | high_out_data.title = title |
---|
384 | high_out_data.style = "Line" |
---|
385 | high_out_data.has_errors = False |
---|
386 | high_out_data.plot_role = Data1D.ROLE_DEFAULT |
---|
387 | |
---|
388 | # copy labels and units of axes for plotting |
---|
389 | high_out_data._xaxis = temp_data._xaxis |
---|
390 | high_out_data._xunit = temp_data._xunit |
---|
391 | high_out_data._yaxis = temp_data._yaxis |
---|
392 | high_out_data._yunit = temp_data._yunit |
---|
393 | |
---|
394 | # Add the plot to the model item |
---|
395 | # This needs to run in the main thread |
---|
396 | reactor.callFromThread(GuiUtils.updateModelItemWithPlot, |
---|
397 | self._model_item, high_out_data, title) |
---|
398 | |
---|
399 | item = QtGui.QStandardItem(str(float('%.3g'% volume_fraction))) |
---|
400 | self.model.setItem(WIDGETS.W_VOLUME_FRACTION, item) |
---|
401 | item = QtGui.QStandardItem(str(float('%.3g'% volume_fraction_error))) |
---|
402 | self.model.setItem(WIDGETS.W_VOLUME_FRACTION_ERR, item) |
---|
403 | if surface: |
---|
404 | item = QtGui.QStandardItem(str(float('%.3g'% surface))) |
---|
405 | self.model.setItem(WIDGETS.W_SPECIFIC_SURFACE, item) |
---|
406 | item = QtGui.QStandardItem(str(float('%.3g'% surface_error))) |
---|
407 | self.model.setItem(WIDGETS.W_SPECIFIC_SURFACE_ERR, item) |
---|
408 | item = QtGui.QStandardItem(str(float('%.3g'% qstar_total))) |
---|
409 | self.model.setItem(WIDGETS.W_INVARIANT, item) |
---|
410 | item = QtGui.QStandardItem(str(float('%.3g'% qstar_total_error))) |
---|
411 | self.model.setItem(WIDGETS.W_INVARIANT_ERR, item) |
---|
412 | |
---|
413 | item = QtGui.QStandardItem(str(float('%.3g'% qstar_low))) |
---|
414 | self.model.setItem(WIDGETS.D_LOW_QSTAR, item) |
---|
415 | item = QtGui.QStandardItem(str(float('%.3g'% qstar_low_err))) |
---|
416 | self.model.setItem(WIDGETS.D_LOW_QSTAR_ERR, item) |
---|
417 | item = QtGui.QStandardItem(str(float('%.3g'% qstar_high))) |
---|
418 | self.model.setItem(WIDGETS.D_HIGH_QSTAR, item) |
---|
419 | item = QtGui.QStandardItem(str(float('%.3g'% qstar_high_err))) |
---|
420 | self.model.setItem(WIDGETS.D_HIGH_QSTAR_ERR, item) |
---|
421 | |
---|
422 | return self.model |
---|
423 | |
---|
424 | def title(self): |
---|
425 | """ Perspective name """ |
---|
426 | return "Invariant panel" |
---|
427 | |
---|
428 | def onStatus(self): |
---|
429 | """ |
---|
430 | Display Invariant Details panel when clicking on Status button |
---|
431 | """ |
---|
432 | self.detailsDialog.setModel(self.model) |
---|
433 | self.detailsDialog.showDialog() |
---|
434 | self.cmdStatus.setEnabled(False) |
---|
435 | |
---|
436 | def onHelp(self): |
---|
437 | """ Display help when clicking on Help button """ |
---|
438 | treeLocation = "/user/qtgui/Perspectives/Invariant/invariant_help.html" |
---|
439 | self.parent.showHelp(treeLocation) |
---|
440 | |
---|
441 | def setupSlots(self): |
---|
442 | """ """ |
---|
443 | self.cmdCalculate.clicked.connect(self.calculateInvariant) |
---|
444 | self.cmdStatus.clicked.connect(self.onStatus) |
---|
445 | self.cmdHelp.clicked.connect(self.onHelp) |
---|
446 | |
---|
447 | self.chkLowQ.stateChanged.connect(self.stateChanged) |
---|
448 | self.chkLowQ.stateChanged.connect(self.checkQExtrapolatedData) |
---|
449 | |
---|
450 | self.chkHighQ.stateChanged.connect(self.stateChanged) |
---|
451 | self.chkHighQ.stateChanged.connect(self.checkQExtrapolatedData) |
---|
452 | |
---|
453 | # slots for the Guinier and PowerLaw radio buttons at low Q |
---|
454 | # since they are not auto-exclusive |
---|
455 | self.rbGuinier.toggled.connect(self.lowGuinierAndPowerToggle) |
---|
456 | |
---|
457 | self.rbPowerLawLowQ.toggled.connect(self.lowGuinierAndPowerToggle) |
---|
458 | |
---|
459 | self.rbFitHighQ.toggled.connect(self.hiFitAndFixToggle) |
---|
460 | |
---|
461 | self.rbFitLowQ.toggled.connect(self.lowFitAndFixToggle) |
---|
462 | |
---|
463 | self.model.itemChanged.connect(self.modelChanged) |
---|
464 | |
---|
465 | # update model from gui editing by users |
---|
466 | self.txtBackgd.textChanged.connect(self.updateFromGui) |
---|
467 | |
---|
468 | self.txtScale.textChanged.connect(self.updateFromGui) |
---|
469 | |
---|
470 | self.txtContrast.textChanged.connect(self.updateFromGui) |
---|
471 | |
---|
472 | self.txtPorodCst.textChanged.connect(self.updateFromGui) |
---|
473 | |
---|
474 | self.txtPowerLowQ.textChanged.connect(self.updateFromGui) |
---|
475 | |
---|
476 | self.txtPowerHighQ.textChanged.connect(self.updateFromGui) |
---|
477 | |
---|
478 | self.txtNptsLowQ.textChanged.connect(self.updateFromGui) |
---|
479 | |
---|
480 | self.txtNptsHighQ.textChanged.connect(self.updateFromGui) |
---|
481 | |
---|
482 | # check values of n_points compared to distribution length |
---|
483 | if self.txtNptsLowQ.isEnabled(): |
---|
484 | self.txtNptsLowQ.textChanged.connect(self.checkLength) |
---|
485 | |
---|
486 | if self.txtNptsHighQ.isEnabled(): |
---|
487 | self.txtNptsHighQ.textChanged.connect(self.checkLength) |
---|
488 | |
---|
489 | def stateChanged(self): |
---|
490 | """ |
---|
491 | Catch modifications from low- and high-Q extrapolation check boxes |
---|
492 | """ |
---|
493 | sender = self.sender() |
---|
494 | |
---|
495 | itemf = QtGui.QStandardItem(str(sender.isChecked()).lower()) |
---|
496 | if sender.text() == 'Enable Low-Q extrapolation': |
---|
497 | self.model.setItem(WIDGETS.W_ENABLE_LOWQ, itemf) |
---|
498 | |
---|
499 | if sender.text() == 'Enable High-Q extrapolation': |
---|
500 | self.model.setItem(WIDGETS.W_ENABLE_HIGHQ, itemf) |
---|
501 | |
---|
502 | def checkLength(self): |
---|
503 | """ |
---|
504 | Validators of number of points for extrapolation. |
---|
505 | Error if it is larger than the distribution length |
---|
506 | """ |
---|
507 | try: |
---|
508 | int_value = int(self.sender().text()) |
---|
509 | except ValueError: |
---|
510 | self.sender().setStyleSheet(BG_RED) |
---|
511 | self.cmdCalculate.setEnabled(False) |
---|
512 | return |
---|
513 | |
---|
514 | if self._data: |
---|
515 | if len(self._data.x) < int_value: |
---|
516 | self.sender().setStyleSheet(BG_RED) |
---|
517 | logging.warning('The number of points must be smaller than {}'.format(len(self._data.x))) |
---|
518 | self.cmdCalculate.setEnabled(False) |
---|
519 | else: |
---|
520 | self.sender().setStyleSheet(BG_WHITE) |
---|
521 | self.cmdCalculate.setEnabled(True) |
---|
522 | else: |
---|
523 | # logging.info('no data is loaded') |
---|
524 | self.cmdCalculate.setEnabled(False) |
---|
525 | |
---|
526 | def modelChanged(self, item): |
---|
527 | """ Update when model changed """ |
---|
528 | if item.row() == WIDGETS.W_ENABLE_LOWQ: |
---|
529 | toggle = (str(item.text()) == 'true') |
---|
530 | self._low_extrapolate = toggle |
---|
531 | self.lowQToggle(toggle) |
---|
532 | elif item.row() == WIDGETS.W_ENABLE_HIGHQ: |
---|
533 | toggle = (str(item.text()) == 'true') |
---|
534 | self._high_extrapolate = toggle |
---|
535 | self.highQToggle(toggle) |
---|
536 | |
---|
537 | def checkQExtrapolatedData(self): |
---|
538 | """ |
---|
539 | Match status of low or high-Q extrapolated data checkbox in |
---|
540 | DataExplorer with low or high-Q extrapolation checkbox in invariant |
---|
541 | panel |
---|
542 | """ |
---|
543 | # name to search in DataExplorer |
---|
544 | if 'Low' in str(self.sender().text()): |
---|
545 | name = "Low-Q extrapolation" |
---|
546 | if 'High' in str(self.sender().text()): |
---|
547 | name = "High-Q extrapolation" |
---|
548 | |
---|
549 | GuiUtils.updateModelItemStatus(self._manager.filesWidget.model, |
---|
550 | self._path, name, |
---|
551 | self.sender().checkState()) |
---|
552 | |
---|
553 | def updateFromGui(self): |
---|
554 | """ Update model when new user inputs """ |
---|
555 | possible_senders = ['txtBackgd', 'txtContrast', 'txtPorodCst', |
---|
556 | 'txtScale', 'txtPowerLowQ', 'txtPowerHighQ', |
---|
557 | 'txtNptsLowQ', 'txtNptsHighQ'] |
---|
558 | |
---|
559 | related_widgets = [WIDGETS.W_BACKGROUND, WIDGETS.W_CONTRAST, |
---|
560 | WIDGETS.W_POROD_CST, WIDGETS.W_SCALE, |
---|
561 | WIDGETS.W_LOWQ_POWER_VALUE, WIDGETS.W_HIGHQ_POWER_VALUE, |
---|
562 | WIDGETS.W_NPTS_LOWQ, WIDGETS.W_NPTS_HIGHQ] |
---|
563 | |
---|
564 | related_internal_values = [self._background, self._contrast, |
---|
565 | self._porod, self._scale, |
---|
566 | self._low_power_value, |
---|
567 | self._high_power_value, |
---|
568 | self._low_points, self._high_points] |
---|
569 | |
---|
570 | item = QtGui.QStandardItem(self.sender().text()) |
---|
571 | |
---|
572 | index_elt = possible_senders.index(self.sender().objectName()) |
---|
573 | |
---|
574 | self.model.setItem(related_widgets[index_elt], item) |
---|
575 | try: |
---|
576 | related_internal_values[index_elt] = float(self.sender().text()) |
---|
577 | self.sender().setStyleSheet(BG_WHITE) |
---|
578 | self.cmdCalculate.setEnabled(True) |
---|
579 | except ValueError: |
---|
580 | # empty field, just skip |
---|
581 | self.sender().setStyleSheet(BG_RED) |
---|
582 | self.cmdCalculate.setEnabled(False) |
---|
583 | |
---|
584 | def lowGuinierAndPowerToggle(self, toggle): |
---|
585 | """ |
---|
586 | Guinier and Power radio buttons cannot be selected at the same time |
---|
587 | If Power is selected, Fit and Fix radio buttons are visible and |
---|
588 | Power line edit can be edited if Fix is selected |
---|
589 | """ |
---|
590 | if self.sender().text() == 'Guinier': |
---|
591 | self._low_guinier = toggle |
---|
592 | |
---|
593 | toggle = not toggle |
---|
594 | self.rbPowerLawLowQ.setChecked(toggle) |
---|
595 | |
---|
596 | self.rbFitLowQ.toggled.connect(self.lowFitAndFixToggle) |
---|
597 | self.rbFitLowQ.setVisible(toggle) |
---|
598 | self.rbFixLowQ.setVisible(toggle) |
---|
599 | |
---|
600 | self.txtPowerLowQ.setEnabled(toggle and (not self._low_fit)) |
---|
601 | |
---|
602 | else: |
---|
603 | self._low_guinier = not toggle |
---|
604 | |
---|
605 | self.rbGuinier.setChecked(not toggle) |
---|
606 | |
---|
607 | self.rbFitLowQ.toggled.connect(self.lowFitAndFixToggle) |
---|
608 | self.rbFitLowQ.setVisible(toggle) |
---|
609 | self.rbFixLowQ.setVisible(toggle) |
---|
610 | |
---|
611 | self.txtPowerLowQ.setEnabled(toggle and (not self._low_fit)) |
---|
612 | |
---|
613 | def lowFitAndFixToggle(self, toggle): |
---|
614 | """ Fit and Fix radiobuttons cannot be selected at the same time """ |
---|
615 | self._low_fit = toggle |
---|
616 | |
---|
617 | toggle = not toggle |
---|
618 | self.txtPowerLowQ.setEnabled(toggle) |
---|
619 | |
---|
620 | def hiFitAndFixToggle(self, toggle): |
---|
621 | """ |
---|
622 | Enable editing of power exponent if Fix for high Q is checked |
---|
623 | Disable otherwise |
---|
624 | """ |
---|
625 | self.txtPowerHighQ.setEnabled(not toggle) |
---|
626 | |
---|
627 | def highQToggle(self, clicked): |
---|
628 | """ Disable/enable High Q extrapolation """ |
---|
629 | self.rbFitHighQ.setEnabled(clicked) |
---|
630 | self.rbFixHighQ.setEnabled(clicked) |
---|
631 | self.txtNptsHighQ.setEnabled(clicked) |
---|
632 | self.txtPowerHighQ.setEnabled(clicked) |
---|
633 | |
---|
634 | def lowQToggle(self, clicked): |
---|
635 | """ Disable / enable Low Q extrapolation """ |
---|
636 | self.rbGuinier.setEnabled(clicked) |
---|
637 | self.rbPowerLawLowQ.setEnabled(clicked) |
---|
638 | self.txtNptsLowQ.setEnabled(clicked) |
---|
639 | # Enable subelements |
---|
640 | self.rbFitLowQ.setVisible(self.rbPowerLawLowQ.isChecked()) |
---|
641 | self.rbFixLowQ.setVisible(self.rbPowerLawLowQ.isChecked()) |
---|
642 | self.rbFitLowQ.setEnabled(clicked) # and not self._low_guinier) |
---|
643 | self.rbFixLowQ.setEnabled(clicked) # and not self._low_guinier) |
---|
644 | |
---|
645 | self.txtPowerLowQ.setEnabled(clicked |
---|
646 | and not self._low_guinier |
---|
647 | and not self._low_fit) |
---|
648 | |
---|
649 | def setupModel(self): |
---|
650 | """ """ |
---|
651 | # filename |
---|
652 | item = QtGui.QStandardItem(self._path) |
---|
653 | self.model.setItem(WIDGETS.W_FILENAME, item) |
---|
654 | |
---|
655 | # add Q parameters to the model |
---|
656 | qmin = 0.0 |
---|
657 | item = QtGui.QStandardItem(str(qmin)) |
---|
658 | self.model.setItem(WIDGETS.W_QMIN, item) |
---|
659 | qmax = 0.0 |
---|
660 | item = QtGui.QStandardItem(str(qmax)) |
---|
661 | self.model.setItem(WIDGETS.W_QMAX, item) |
---|
662 | |
---|
663 | # add custom input params |
---|
664 | item = QtGui.QStandardItem(str(self._background)) |
---|
665 | self.model.setItem(WIDGETS.W_BACKGROUND, item) |
---|
666 | item = QtGui.QStandardItem(str(self._contrast)) |
---|
667 | self.model.setItem(WIDGETS.W_CONTRAST, item) |
---|
668 | item = QtGui.QStandardItem(str(self._scale)) |
---|
669 | self.model.setItem(WIDGETS.W_SCALE, item) |
---|
670 | # leave line edit empty if Porod constant not defined |
---|
671 | if self._porod != None: |
---|
672 | item = QtGui.QStandardItem(str(self._porod)) |
---|
673 | else: |
---|
674 | item = QtGui.QStandardItem(str('')) |
---|
675 | self.model.setItem(WIDGETS.W_POROD_CST, item) |
---|
676 | |
---|
677 | # Dialog elements |
---|
678 | itemf = QtGui.QStandardItem("false") |
---|
679 | self.model.setItem(WIDGETS.W_ENABLE_HIGHQ, itemf) |
---|
680 | itemf = QtGui.QStandardItem("false") |
---|
681 | self.model.setItem(WIDGETS.W_ENABLE_LOWQ, itemf) |
---|
682 | |
---|
683 | item = QtGui.QStandardItem(str(NPOINTS_Q_INTERP)) |
---|
684 | self.model.setItem(WIDGETS.W_NPTS_LOWQ, item) |
---|
685 | item = QtGui.QStandardItem(str(NPOINTS_Q_INTERP)) |
---|
686 | self.model.setItem(WIDGETS.W_NPTS_HIGHQ, item) |
---|
687 | |
---|
688 | itemt = QtGui.QStandardItem("true") |
---|
689 | self.model.setItem(WIDGETS.W_LOWQ_GUINIER, itemt) |
---|
690 | |
---|
691 | itemt = QtGui.QStandardItem("true") |
---|
692 | self.model.setItem(WIDGETS.W_LOWQ_FIT, itemt) |
---|
693 | item = QtGui.QStandardItem(str(DEFAULT_POWER_LOW)) |
---|
694 | self.model.setItem(WIDGETS.W_LOWQ_POWER_VALUE, item) |
---|
695 | |
---|
696 | itemt = QtGui.QStandardItem("true") |
---|
697 | self.model.setItem(WIDGETS.W_HIGHQ_FIT, itemt) |
---|
698 | item = QtGui.QStandardItem(str(DEFAULT_POWER_LOW)) |
---|
699 | self.model.setItem(WIDGETS.W_HIGHQ_POWER_VALUE, item) |
---|
700 | |
---|
701 | def setupMapper(self): |
---|
702 | # Set up the mapper. |
---|
703 | self.mapper = QtWidgets.QDataWidgetMapper(self) |
---|
704 | self.mapper.setOrientation(QtCore.Qt.Vertical) |
---|
705 | self.mapper.setModel(self.model) |
---|
706 | |
---|
707 | # Filename |
---|
708 | self.mapper.addMapping(self.txtName, WIDGETS.W_FILENAME) |
---|
709 | |
---|
710 | # Qmin/Qmax |
---|
711 | self.mapper.addMapping(self.txtTotalQMin, WIDGETS.W_QMIN) |
---|
712 | self.mapper.addMapping(self.txtTotalQMax, WIDGETS.W_QMAX) |
---|
713 | |
---|
714 | # Background |
---|
715 | self.mapper.addMapping(self.txtBackgd, WIDGETS.W_BACKGROUND) |
---|
716 | |
---|
717 | # Scale |
---|
718 | self.mapper.addMapping(self.txtScale, WIDGETS.W_SCALE) |
---|
719 | |
---|
720 | # Contrast |
---|
721 | self.mapper.addMapping(self.txtContrast, WIDGETS.W_CONTRAST) |
---|
722 | |
---|
723 | # Porod constant |
---|
724 | self.mapper.addMapping(self.txtPorodCst, WIDGETS.W_POROD_CST) |
---|
725 | |
---|
726 | # Lowq/Highq items |
---|
727 | self.mapper.addMapping(self.chkLowQ, WIDGETS.W_ENABLE_LOWQ) |
---|
728 | self.mapper.addMapping(self.chkHighQ, WIDGETS.W_ENABLE_HIGHQ) |
---|
729 | |
---|
730 | self.mapper.addMapping(self.txtNptsLowQ, WIDGETS.W_NPTS_LOWQ) |
---|
731 | self.mapper.addMapping(self.rbGuinier, WIDGETS.W_LOWQ_GUINIER) |
---|
732 | self.mapper.addMapping(self.rbFitLowQ, WIDGETS.W_LOWQ_FIT) |
---|
733 | self.mapper.addMapping(self.txtPowerLowQ, WIDGETS.W_LOWQ_POWER_VALUE) |
---|
734 | |
---|
735 | self.mapper.addMapping(self.txtNptsHighQ, WIDGETS.W_NPTS_HIGHQ) |
---|
736 | self.mapper.addMapping(self.rbFitHighQ, WIDGETS.W_HIGHQ_FIT) |
---|
737 | self.mapper.addMapping(self.txtPowerHighQ, WIDGETS.W_HIGHQ_POWER_VALUE) |
---|
738 | |
---|
739 | # Output |
---|
740 | self.mapper.addMapping(self.txtVolFract, WIDGETS.W_VOLUME_FRACTION) |
---|
741 | self.mapper.addMapping(self.txtVolFractErr, WIDGETS.W_VOLUME_FRACTION_ERR) |
---|
742 | self.mapper.addMapping(self.txtSpecSurf, WIDGETS.W_SPECIFIC_SURFACE) |
---|
743 | self.mapper.addMapping(self.txtSpecSurfErr, WIDGETS.W_SPECIFIC_SURFACE_ERR) |
---|
744 | self.mapper.addMapping(self.txtInvariantTot, WIDGETS.W_INVARIANT) |
---|
745 | self.mapper.addMapping(self.txtInvariantTotErr, WIDGETS.W_INVARIANT_ERR) |
---|
746 | |
---|
747 | self.mapper.toFirst() |
---|
748 | |
---|
749 | def setData(self, data_item=None, is_batch=False): |
---|
750 | """ |
---|
751 | Obtain a QStandardItem object and dissect it to get Data1D/2D |
---|
752 | Pass it over to the calculator |
---|
753 | """ |
---|
754 | assert data_item is not None |
---|
755 | |
---|
756 | if self.txtName.text() == data_item[0].text(): |
---|
757 | logging.info('This file is already loaded in Invariant panel.') |
---|
758 | return |
---|
759 | |
---|
760 | if not isinstance(data_item, list): |
---|
761 | msg = "Incorrect type passed to the Invariant Perspective" |
---|
762 | raise AttributeError(msg) |
---|
763 | |
---|
764 | if not isinstance(data_item[0], QtGui.QStandardItem): |
---|
765 | msg = "Incorrect type passed to the Invariant Perspective" |
---|
766 | raise AttributeError(msg) |
---|
767 | |
---|
768 | # only 1 file can be loaded |
---|
769 | self._model_item = data_item[0] |
---|
770 | |
---|
771 | # Extract data on 1st child - this is the Data1D/2D component |
---|
772 | data = GuiUtils.dataFromItem(self._model_item) |
---|
773 | self.model.item(WIDGETS.W_FILENAME).setData(self._model_item.text()) |
---|
774 | # update GUI and model with info from loaded data |
---|
775 | self.updateGuiFromFile(data=data) |
---|
776 | |
---|
777 | def updateGuiFromFile(self, data=None): |
---|
778 | """ |
---|
779 | update display in GUI and plot |
---|
780 | """ |
---|
781 | self._data = data |
---|
782 | |
---|
783 | # plot loaded file |
---|
784 | if not isinstance(self._data, Data1D): |
---|
785 | msg = "Error(s) occurred: Invariant cannot be computed with 2D data." |
---|
786 | raise AttributeError(msg) |
---|
787 | |
---|
788 | try: |
---|
789 | filename = data.filename |
---|
790 | except: |
---|
791 | msg = 'No filename' |
---|
792 | raise ValueError(msg) |
---|
793 | try: |
---|
794 | qmin = min(self._data.x) |
---|
795 | qmax = max(self._data.x) |
---|
796 | except: |
---|
797 | msg = "Unable to find q min/max of \n data named %s" % \ |
---|
798 | data.filename |
---|
799 | raise ValueError(msg) |
---|
800 | |
---|
801 | # update model with input form files: filename, qmin, qmax |
---|
802 | self.model.item(WIDGETS.W_FILENAME).setText(filename) |
---|
803 | self.model.item(WIDGETS.W_QMIN).setText(str(qmin)) |
---|
804 | self.model.item(WIDGETS.W_QMAX).setText(str(qmax)) |
---|
805 | self._path = filename |
---|
806 | |
---|
807 | # Calculate and add to GUI: volume fraction, invariant total, |
---|
808 | # and specific surface if porod checked |
---|
809 | self.calculateInvariant() |
---|
810 | |
---|
811 | def allowBatch(self): |
---|
812 | """ |
---|
813 | Tell the caller that we don't accept multiple data instances |
---|
814 | """ |
---|
815 | return False |
---|