Changeset 1264


Ignore:
Timestamp:
Mar 26, 2014 3:08:59 PM (10 years ago)
Author:
vondreele
Message:

Fix SASD plotting when sqPlot is True: log(Q vs Y*Q4) plot

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIplot.py

    r1252 r1264  
    890890                    D = xye[5]-Ymax*G2frame.delOffset
    891891                elif 'SASD' in plottype:
    892                     D = xye[4]
     892                    if G2frame.sqPlot:
     893                        D = xye[4]*X**4
     894                        Z = xye[3]*X**4
     895                    else:
     896                        D = xye[4]
    893897                    Plot.set_yscale("log",nonposy='mask')
    894898                    Plot.set_ylim(bottom=np.min(np.trim_zeros(Y))/2.,top=np.max(Y)*2.)
     
    903907                        Plot.set_yscale("log",nonposy='mask')
    904908                        if G2frame.ErrorBars:
    905                             Plot.errorbar(X,Y,yerr=Sample['Scale'][0]*np.sqrt(1./(Pattern[0]['wtFactor']*xye[2])),
    906                                 ecolor=colors[N%6],picker=3.,clip_on=False)
     909                            if G2frame.sqPlot:
     910                                Plot.errorbar(X,Y,yerr=X**4*Sample['Scale'][0]*np.sqrt(1./(Pattern[0]['wtFactor']*xye[2])),
     911                                    ecolor=colors[N%6],picker=3.,clip_on=False)
     912                            else:
     913                                Plot.errorbar(X,Y,yerr=Sample['Scale'][0]*np.sqrt(1./(Pattern[0]['wtFactor']*xye[2])),
     914                                    ecolor=colors[N%6],picker=3.,clip_on=False)
    907915                        else:
    908916                            Plot.plot(X,Y,colors[N%6]+'+',picker=3.,clip_on=False)
    909                         Plot.plot(X,D,colors[(N+1)%6],picker=False)
     917                        Plot.plot(X,D,colors[(N+2)%6],picker=False)
    910918                        Plot.plot(X,Z,colors[(N+1)%6],picker=False)
    911919                elif G2frame.Weight and 'PWDR' in plottype:
Note: See TracChangeset for help on using the changeset viewer.