Changeset 1933
- Timestamp:
- Jul 17, 2015 11:37:24 AM (8 years ago)
- Location:
- trunk
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/GSASIIphsGUI.py
r1931 r1933 102 102 reflSets = G2frame.PatternTree.GetItemPyData(G2gd.GetPatternTreeItemId(G2frame,PatternId,'Reflection Lists')) 103 103 reflData = reflSets[phaseName] 104 if 'list' in str(type(reflData)): #patch for old reflection data105 RefData = {'RefList':[],'FF':[]}106 for ref in reflDict:107 RefData['RefList'].append(ref[:11]+[ref[13],])108 RefData['FF'].append(ref[14])109 RefData['RefList'] = np.array(RefData['RefList'])110 reflData = RefData111 104 elif 'HKLF' in reflName: 112 105 PatternId = G2gd.GetPatternTreeItemId(G2frame,G2frame.root, reflName) -
trunk/GSASIIplot.py
r1932 r1933 969 969 plotting available. 970 970 ''' 971 global HKL972 971 global exclLines 973 972 global DifLine … … 1147 1146 if pickIdText in ['Index Peak List','Unit Cells List',]: 1148 1147 indx = -2 1149 if len( HKL):1148 if len(G2frame.HKL): 1150 1149 view = Page.toolbar._views.forward()[0][:2] 1151 1150 wid = view[1]-view[0] 1152 found = HKL[np.where(np.fabs(HKL.T[indx]-xpos) < 0.002*wid)]1151 found = G2frame.HKL[np.where(np.fabs(G2frame.HKL.T[indx]-xpos) < 0.002*wid)] 1153 1152 if len(found): 1154 1153 if len(found[0]) > 6: #SS reflections … … 1341 1340 Pattern[0]['delOffset'] = -ypos/Ymax 1342 1341 G2frame.itemPicked = None 1343 G2frame.PatternTree.SetItemPyData(PickId,data)1344 1342 wx.CallAfter(PlotPatterns,G2frame,plotType=plottype) 1345 1343 return … … 1369 1367 limits[1][0] = min(max(limits[0][0],limits[1][0]),limits[1][1]) 1370 1368 limits[1][1] = max(min(limits[0][1],limits[1][1]),limits[1][0]) 1371 G2frame.PatternTree.SetItemPyData(LimitId,limits)1372 1369 if G2frame.PatternTree.GetItemText(G2frame.PickId) == 'Limits': 1373 1370 G2pdG.UpdateLimitsGrid(G2frame,limits,plottype) … … 1385 1382 peaks['peaks'][lineNo-2][0] = xpos 1386 1383 peaks['sigDict'] = {} #no longer valid 1387 G2frame.PatternTree.SetItemPyData(PeakId,peaks)1388 1384 G2pdG.UpdatePeakGrid(G2frame,peaks) 1389 1385 elif G2frame.PatternTree.GetItemText(PickId) in ['Models',] and xpos: … … 1401 1397 data[1][0] = min(max(data[0][0],data[1][0]),data[1][1]) 1402 1398 data[1][1] = max(min(data[0][1],data[1][1]),data[1][0]) 1403 G2frame.PatternTree.SetItemPyData(LimitId,data)1404 1399 elif (G2frame.PatternTree.GetItemText(PickId) == 'Reflection Lists' or \ 1405 1400 'PWDR' in G2frame.PatternTree.GetItemText(PickId)) and xpos: … … 1413 1408 else: #1st row of refl ticks 1414 1409 data[0]['refOffset'] = event.ydata 1415 G2frame.PatternTree.SetItemPyData(PickId,data)1416 1410 PlotPatterns(G2frame,plotType=plottype) 1417 1411 G2frame.itemPicked = None … … 1522 1516 item, cookie = G2frame.PatternTree.GetNextChild(G2frame.root, cookie) 1523 1517 lenX = 0 1524 if PickId:1525 if G2frame.PatternTree.GetItemText(PickId) in ['Reflection Lists']:1526 Phases = G2frame.PatternTree.GetItemPyData(G2gd.GetPatternTreeItemId(G2frame,PatternId,'Reflection Lists'))1527 HKL = []1528 if Phases:1529 try:1530 for peak in Phases[G2frame.RefList]['RefList']:1531 if len(peak) > 15:1532 HKL.append(peak[:7]) #SS reflection list - need peak[:7]1533 else:1534 HKL.append(peak[:6])1535 except TypeError: #old style patch1536 for peak in Phases[G2frame.RefList]:1537 HKL.append(peak[:6])1538 HKL = np.array(HKL)1539 else:1540 HKL = np.array(G2frame.HKL)1541 1518 Ymax = None 1542 1519 for Pattern in PlotList: … … 1813 1790 Phases = G2frame.PatternTree.GetItemPyData(G2gd.GetPatternTreeItemId(G2frame,PatternId,'Reflection Lists')) 1814 1791 for pId,phase in enumerate(Phases): 1815 try: #patch for old style reflection lists 1816 peaks = Phases[phase]['RefList'] 1817 except TypeError: 1818 peaks = Phases[phase] 1792 peaks = Phases[phase]['RefList'] 1819 1793 if not len(peaks): 1820 1794 continue … … 2426 2400 ''' plotting of powder lines (i.e. no powder pattern) as sticks 2427 2401 ''' 2428 global HKL2429 2402 2430 2403 def OnMotion(event): … … 2435 2408 if G2frame.PickId and G2frame.PatternTree.GetItemText(G2frame.PickId) in ['Index Peak List','Unit Cells List']: 2436 2409 found = [] 2437 if len( HKL):2410 if len(G2frame.HKL): 2438 2411 view = Page.toolbar._views.forward()[0][:2] 2439 2412 wid = view[1]-view[0] 2440 found = HKL[np.where(np.fabs(HKL.T[-1]-xpos) < 0.002*wid)]2413 found = G2frame.HKL[np.where(np.fabs(G2frame.HKL.T[-1]-xpos) < 0.002*wid)] 2441 2414 if len(found): 2442 2415 h,k,l = found[0][:3] … … 2465 2438 for peak in peaks: 2466 2439 Plot.axvline(peak[0],color='b') 2467 HKL = np.array(G2frame.HKL)2468 2440 for hkl in G2frame.HKL: 2469 2441 Plot.axvline(hkl[-1],color='r',dashes=(5,5)) -
trunk/GSASIIpwdGUI.py
r1913 r1933 2299 2299 maxH = ssopt['maxH'] 2300 2300 G2frame.HKL = G2pwd.getHKLMpeak(dmin,Inst,SGData,SSGData,Vec,maxH,A) 2301 G2frame.HKL = np.array(G2frame.HKL) 2301 2302 data[0] = G2indx.IndexSSPeaks(data[0],G2frame.HKL)[1] 2302 2303 else: #select cell from table - no SS … … 2308 2309 for hkl in G2frame.HKL: 2309 2310 hkl.insert(4,G2lat.Dsp2pos(Inst,hkl[3])) 2311 G2frame.HKL = np.array(G2frame.HKL) 2310 2312 data[0] = G2indx.IndexPeaks(data[0],G2frame.HKL)[1] 2311 2313 break … … 2598 2600 M20 = X20 = 0. 2599 2601 G2frame.HKL = G2pwd.getHKLpeak(dmin,SGData,A,Inst) 2602 G2frame.HKL = np.array(G2frame.HKL) 2600 2603 if len(G2frame.HKL): 2601 2604 print ' new M20,X20: %.2f %d fraction found: %.3f'%(M20,X20,float(len(peaks[0]))/len(G2frame.HKL)) … … 2685 2688 peaks = [G2indx.IndexPeaks(peaks[0],G2frame.HKL)[1],peaks[1]] #put peak fit esds back in peaks 2686 2689 Lhkl,M20,X20,Aref,Zero = G2indx.refinePeaksT(peaks[0],difC,ibrav,A,controls[1],controls[0]) 2690 G2frame.HKL = np.array(G2frame.HKL) 2687 2691 controls[1] = Zero 2688 2692 controls[6:12] = G2lat.A2cell(Aref) … … 2708 2712 for hkl in G2frame.HKL: 2709 2713 hkl[ip] = G2lat.Dsp2pos(Inst,hkl[ip-1])+controls[1] 2714 G2frame.HKL = np.array(G2frame.HKL) 2710 2715 if 'PKS' in G2frame.PatternTree.GetItemText(G2frame.PatternId): 2711 2716 G2plt.PlotPowderLines(G2frame) … … 2755 2760 for hkl in G2frame.HKL: 2756 2761 hkl.insert(4,G2lat.Dsp2pos(Inst,hkl[3])+controls[1]) 2762 G2frame.HKL = np.array(G2frame.HKL) 2757 2763 if 'PKS' in G2frame.PatternTree.GetItemText(G2frame.PatternId): 2758 2764 G2plt.PlotPowderLines(G2frame) … … 2782 2788 for hkl in G2frame.HKL: 2783 2789 hkl.insert(4,G2lat.Dsp2pos(Inst,hkl[3])+controls[1]) 2790 G2frame.HKL = np.array(G2frame.HKL) 2784 2791 if 'PKS' in G2frame.PatternTree.GetItemText(G2frame.PatternId): 2785 2792 G2plt.PlotPowderLines(G2frame) … … 3030 3037 for hkl in G2frame.HKL: 3031 3038 hkl.insert(4,G2lat.Dsp2pos(Inst,hkl[3])+controls[1]) 3039 G2frame.HKL = np.array(G2frame.HKL) 3032 3040 table.append(row) 3033 3041 UnitCellsTable = G2G.Table(table,rowLabels=rowLabels,colLabels=colLabels,types=Types) … … 3150 3158 elif 'T' in Inst['Type'][0]: 3151 3159 refs = np.vstack((refList.T[:18+Super],I100)).T 3152 for i in range(len(refs)): rowLabels.append(str(i)) 3160 G2frame.HKL = np.vstack((refList.T[:6+Super])).T #build for plots 3161 rowLabels = [str(i) for i in range(len(refs))] 3153 3162 Types = (4+Super)*[wg.GRID_VALUE_LONG,]+4*[wg.GRID_VALUE_FLOAT+':10,4',]+ \ 3154 3163 2*[wg.GRID_VALUE_FLOAT+':10,2',]+[wg.GRID_VALUE_FLOAT+':10,3',]+ \
Note: See TracChangeset
for help on using the changeset viewer.