Changeset 55db501 in sasview

Timestamp:

Nov 20, 2016 5:36:43 AM (7 years ago)

Author:

jhbakker

Branches:

master, 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, costrafo411, magnetic_scatt, release-4.1.1, release-4.1.2, release-4.2.2, ticket-1009, ticket-1094-headless, ticket-1242-2d-resolution, ticket-1243, ticket-1249, ticket885, unittest-saveload

Children:

e4078ed

Parents:

a01af35

Message:

test branch for SESANS class in sesans.py

Location:

src/sas

Files:

• sascalc/data_util/qsmearing.py (modified) (4 diffs)
• sasgui/perspectives/fitting/basepage.py (modified) (3 diffs)
• sasgui/perspectives/fitting/fitpage.py (modified) (2 diffs)

src/sas/sascalc/data_util/qsmearing.py

@@ -15 +15 @@
 import time
 from sasmodels import sesans
-
 import numpy as np  # type: ignore
 from numpy import pi, exp  # type: ignore
@@ -61 +60 @@
     # TODO: change other sanity checks to check for file loader instead of data structure?
     _found_sesans = False
-    if data.dx is not None and data.meta_data['loader']=='SESANS':
+    #if data.dx is not None and data.meta_data['loader']=='SESANS':
+    if data.dx is not None and data.isSesans:
         if data.dx[0] > 0.0:
             _found_sesans = True
@@ -67 +67 @@
     if _found_sesans == True:
         #Pre-compute the Hankel matrix (H)
-        H0,H, q_calc = sesans.Hankelconstructor(data)
+        qmax, qunits = data.sample.zacceptance
+        Hankelinst=sesans.SesansTransform()
+        sesans.SesansTransform.set_transform(Hankelinst,
+        SE = Converter(data._xunit)(data.x, "A"),
+        zaccept = Converter(qunits)(qmax, "1/A"),
+        Rmax = 1000000)
+        H=sesans.SesansTransform._H
+        H0=sesans.SesansTransform._H0
+        q=sesans.SesansTransform.q
         # Then return the actual transform, as if it were a smearing function
-        return PySmear(SESANS1D(data, H0, H, q_calc), model)
+        # applying evalDistribution to a model, with a q-space as param, returns the I(q) values that go with the q-values
+
+        return PySmear(SESANS1D(data, H0, H, q), model)
 
     _found_resolution = False
@@ -175 +185 @@
     width = data.dx if data.dx is not None else 0
     return PySmear(Pinhole1D(q, width), model)
-
-def sesans_smear(data, model=None):
-    #This should be calculated characteristic length scale
-    #Probably not a data prameter either
-    #Need function to calculate this based on model
-    #Here assume a number
-    Rmax = 1000000
-    q_calc = sesans.make_q(data.sample.zacceptance, Rmax)
-    SElength=Converter(data._xunit)(data.x, "A")
-    #return sesans.HankelTransform(q_calc, SElength)
-    #Old return statement, running through the smearer
-    #return PySmear(SESANS1D(data,q_calc),model)

src/sas/sasgui/perspectives/fitting/basepage.py

@@ -431 +431 @@
         self.sizer2 = wx.BoxSizer(wx.VERTICAL)
         self.sizer3 = wx.BoxSizer(wx.VERTICAL)
-        self.sizerTrafo = wx.BoxSizer(wx.VERTICAL)
         self.sizer4 = wx.BoxSizer(wx.VERTICAL)
         self.sizer5 = wx.BoxSizer(wx.VERTICAL)
@@ -440 +439 @@
         self.sizer2.SetMinSize((PANEL_WIDTH, -1))
         self.sizer3.SetMinSize((PANEL_WIDTH, -1))
-        self.sizerTrafo.SetMinSize((PANEL_WIDTH, -1))
         self.sizer4.SetMinSize((PANEL_WIDTH, -1))
         self.sizer5.SetMinSize((PANEL_WIDTH, -1))
@@ -449 +447 @@
         self.vbox.Add(self.sizer2)
         self.vbox.Add(self.sizer3)
-        self.vbox.Add(self.sizerTrafo)
         self.vbox.Add(self.sizer4)
         self.vbox.Add(self.sizer5)

src/sas/sasgui/perspectives/fitting/fitpage.py

@@ -55 +55 @@
         self.weightbt_string = None
         self.m_name = None
-        # transform implementation
-        self._fill_Trafo_sizer()
-       # self.Trafobt_string()
         # get smear info from data
         self._get_smear_info()
@@ -95 +92 @@
         self.parent.on_set_focus(event)
         self.on_tap_focus()
-
-    def onTrafo(self, event):
-        """
-        On Trafo radio button event, sets the Trafobt_string
-        """
-        self.Trafobt_string = event.GetEventObject().GetLabelText()
-        self._set_Trafo()
-
-    def _fill_Trafo_sizer(self):
-
-        title = "Transform"
-        box_description_trafo = wx.StaticBox(self, wx.ID_ANY, str(title))
-        box_description_trafo.SetForegroundColour(wx.BLUE)
-        #boxsizer_trafo = wx.StaticBoxSizer(box_description_trafo, wx.VERTICAL)
-        boxsizer_trafo = wx.StaticBoxSizer(box_description_trafo, wx.HORIZONTAL)
-        #sizer_trafo = wx.StaticBoxSizer(box_description_trafo, wx.HORIZONTAL)
-        #weighting_set_box = wx.StaticBox(self, wx.ID_ANY,
-         #                              'Select the type of SESANS analysis')
-
-        #sizer_weighting = wx.BoxSizer(wx.HORIZONTAL)
-          #      weighting_box.SetMinSize((_DATA_BOX_WIDTH, 60))
-
-        #For every radio button (each statement x3):
-        self.no_transform = wx.RadioButton(self, wx.ID_ANY,
-                                                  'None', style=wx.RB_GROUP)
-
-        #self.Bind(wx.EVT_RADIOBUTTON, self.onTrafo,
-         #                 id=self.no_transform.GetId())
-        self.hankel = wx.RadioButton(self, wx.ID_ANY,
-                                                  'Hankel')
-        #self.Bind(wx.EVT_RADIOBUTTON, self.onTrafo,
-        #                  id=self.hankel.GetId())
-        self.cosine = wx.RadioButton(self, wx.ID_ANY,
-                                                  'Cosine')
-        #self.Bind(wx.EVT_RADIOBUTTON, self.onTrafo,
-         #                 id=self.cosine.GetId())
-
-        #Not sure about this (only once though)
-        self.no_transform.SetValue(True)
-
-        #For every radio button (each statement x3):
-        boxsizer_trafo.Add(self.no_transform, 0, wx.LEFT, 10)
-        boxsizer_trafo.Add((14, 10))
-        boxsizer_trafo.Add(self.hankel)
-        boxsizer_trafo.Add((14, 10))
-        boxsizer_trafo.Add(self.cosine)
-        boxsizer_trafo.Add((14, 10))
-            #Default for weighting is False, but these need to be on by default!
-        self.no_transform.Enable(True)
-
-        #Not sure about this (only once though)
-        #weighting_box.Add(sizer_trafo)
-
-        self.sizerTrafo.Clear(True)
-        self.sizerTrafo.Add(boxsizer_trafo, 0, wx.EXPAND | wx.ALL, 10)
-        #self.sizerTrafo.Add(sizer_trafo, 0, wx.EXPAND | wx.ALL, 10)
-        self.sizerTrafo.Layout()
 
     def _fill_data_sizer(self):