Changeset 4269


Ignore:
Timestamp:
Jan 29, 2020 7:51:54 AM (21 months ago)
Author:
vondreele
Message:

clear out old bond files before next run of RMCProfile
fix conversion from size/mustrain to X & Y for x-ray data
add experimental G(R) to partial G(R) plot scaled & shifted to allow easy comparison

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIphsGUI.py

    r4268 r4269  
    51595159            if os.path.isfile(pName+'.triplets'):
    51605160                os.remove(pName+'.triplets')
     5161            i = 1
     5162            while True:
     5163                if os.path.isfile(pName+'.bondodf_%d'%i):
     5164                    os.remove(pName+'.bondodf_%d'%i)
     5165                    os.remove(pName+'_bondplot_%d.ppm'%i)
     5166                    i += 1
     5167                else:
     5168                    break               
    51615169            G2frame.OnFileSave(event)
    51625170            print (' GSAS-II project saved')
     
    52155223                '_FQ1.csv':[r'$\mathsf{Q,\AA^{-1}}$','F(Q)','RMCP F(Q) for '],
    52165224                '_FT_XFQ1.csv':[r'$\mathsf{R,\AA}$','G(R)','RMCP x-ray G(R) for '],
    5217                 '_XFQ1.csv':[r'$\mathsf{Q,\AA^{-1}}$','F(Q)','RMCP x-ray F(Q) for '],
     5225#                '_XFQ1.csv':[r'$\mathsf{Q,\AA^{-1}}$','F(Q)','RMCP x-ray F(Q) for '],
    52185226                '_bragg.csv':[r'$\mathsf{TOF,\mu s}$','Normalized Intensity','RMCP bragg for ']}
     5227            Ysave = []
    52195228            for label in Labels:
    52205229                X = []
     
    52305239                    Yobs = np.array([X,Yobs])
    52315240                    Ycalc = np.array([X,Ycalc])
     5241                    if 'G(R)' in Labels[label][1]:
     5242                        Ysave.append(Yobs)
     5243                        Ymin = Ysave[0][1][0]
    52325244                    if 'bragg' in label:
    52335245                        Ydiff = np.array([X,(Yobs-Ycalc)[1]])
     
    52615273                        Partials.append([float(item) for item in items[1:]])
    52625274                    X = np.array(X)
     5275                    DX = X[1]-X[0]  #should be 0.02
    52635276                    Partials = np.array(Partials).T
    52645277                    if 'Q' in label:
     
    52665279                        Names = [name for name in Names if 'Va' not in name]
    52675280                    else:
    5268                         XY = [[X.T,Y.T] for Y in Partials]
    5269                     G2plt.PlotXY(G2frame,XY,labelX=Labels[label][0],
    5270                         labelY=Labels[label][1],newPlot=True,Title=Labels[label][2]+pName,
    5271                         lines=True,names=Names[1:])
     5281                        XY = [[X.T,(DX*Y.T)+Ymin] for Y in Partials]
     5282                    if 'G(R)' in Labels[label][1]:
     5283                        Xmax = np.searchsorted(Ysave[0][0],XY[0][0][-1])
     5284                        G2plt.PlotXY(G2frame,XY2=XY,XY=[Ysave[0][:,0:Xmax],],labelX=Labels[label][0],
     5285                            labelY=Labels[label][1],newPlot=True,Title=Labels[label][2]+pName,
     5286                            lines=False,names=[r'   $G(R)_{obs}$',]+Names[1:])
     5287                    else:                       
     5288                        G2plt.PlotXY(G2frame,XY,labelX=Labels[label][0],
     5289                            labelY=Labels[label][1],newPlot=True,Title=Labels[label][2]+pName,
     5290                            lines=True,names=Names[1:])
    52725291#chi**2 plot
    52735292            X = []
  • trunk/GSASIIpwd.py

    r4268 r4269  
    21292129            wave = G2mth.getWave(inst)
    21302130            if 'ellipsoidal' not in Size[0]:    #take the isotropic term only
    2131                 Xsb = 1.8*wave/(np.pi*Size[1][0])/2.
     2131                Xsb = 1.8*wave/(np.pi*Size[1][0])
    21322132        if useSamBrd[1]:
    21332133            if 'generalized' not in Mustrain[0]:    #take the isotropic term only
    2134                 Ysb = 0.018*np.pi*Mustrain[1][0]/2.
     2134                Ysb = 0.0180*Mustrain[1][0]/np.pi
    21352135        prms = ['Bank',
    21362136                'Lam','Zero','Polariz.',
Note: See TracChangeset for help on using the changeset viewer.