Changeset 4571 for trunk/GSASIIimage.py

Timestamp:

Sep 17, 2020 3:50:24 PM (3 years ago)

Author:

vondreele

Message:

Improvements to polarization stuff

File:

• trunk/GSASIIimage.py (modified) (3 diffs)

trunk/GSASIIimage.py

@@ -1217 +1217 @@
     #prepare masked arrays of bins with pixels for interpolation setup
     H2msk = [ma.array(H2[:-1],mask=np.logical_not(nst)) for nst in NST]
-    H0msk = [ma.array(np.divide(h0,nst),mask=np.logical_not(nst)) for nst,h0 in zip(NST,H0)]
+#    H0msk = [ma.array(np.divide(h0,nst),mask=np.logical_not(nst)) for nst,h0 in zip(NST,H0)]
+    H0msk = [ma.array(h0/nst,mask=np.logical_not(nst)) for nst,h0 in zip(NST,H0)]
     #make linear interpolators; outside limits give NaN
     H0 = []
@@ -1512 +1513 @@
     '''
     from scipy.optimize import minimize_scalar
-    Arc = [arcTth,[45.,135.],2.0]
+    print(' Polarization fit for mask at %.1f deg 2-theta'%arcTth)
+    data = copy.deepcopy(imageData)
+    data['IOtth'] = [arcTth-2.,arcTth+2.]
+    data['fullIntegrate'] = True
+    data['LRazimuth'] = [0.,360.]
+    data['outChannels'] = 200
+    data['outAzimuths'] = 1
+    Arc = [arcTth,[85.,95.],2.0]
+    print(' Integration 2-theta test range %.1f - %.1f in 200 steps'%(data['IOtth'][0],data['IOtth'][1]))
     Masks = {'Points':[],'Rings':[],'Arcs':[],'Polygons':[],'Frames':[],
         'Thresholds':imageData['range'],'SpotMask':{'esdMul':3.,'spotMask':None}}
     AMasks = {'Points':[],'Rings':[],'Arcs':[Arc,],'Polygons':[],'Frames':[],
         'Thresholds':imageData['range'],'SpotMask':{'esdMul':3.,'spotMask':None}}
-    data = copy.deepcopy(imageData)
     
     def func(p):
@@ -1524 +1532 @@
         H0 = ImageIntegrate(ImageZ,data,Masks)[0]
         H0A = ImageIntegrate(ImageZ,data,AMasks)[0]
-        M = np.sum(H0-H0A)
-        print('polarization %.4f, fxn: %.1f'%(p,M))
+        M = np.sum(H0)-np.sum(H0A)
+        print(' Polarization %.4f, fxn: %.1f'%(p,M))
         return M**2