Changeset 798 for trunk/GSASII.py


Ignore:
Timestamp:
Nov 14, 2012 4:37:26 PM (10 years ago)
Author:
vondreele
Message:

more peak fitting work including TOF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASII.py

    r797 r798  
    658658                return [G2IO.makeInstDict(names,data,codes),{}]
    659659            elif 'T' in DataType:
    660                 names = ['Type','2-theta','difC','difA','Zero','alpha','beta-0','beta-1','var-inst','X','Y','Azimuth']
    661                 codes = [0,0,0,0,0,0,0,0,0,0,0,0]
     660                names = ['Type','2-theta','difC','difA','Zero','alpha','beta-0','beta-1','sig-0','sig-1','X','Y','Azimuth']
     661                codes = [0,0,0,0,0,0,0,0,0,0,0,0,0]
    662662                azm = 0.
    663663                if 'INS  1DETAZM' in Iparm:
     
    673673                    data.extend([G2IO.sfloat(s[1]),G2IO.sfloat(s[2]),G2IO.sfloat(s[3])])
    674674                    s = Iparm['INS  1PRCF12'].split()
    675                     data.extend([G2IO.sfloat(s[1]),0.0,0.0,azm])
     675                    data.extend([0.0,G2IO.sfloat(s[1]),0.0,0.0,azm])
    676676                elif abs(pfType) in [3,4,5]:
    677677                    data.extend([G2IO.sfloat(s[0]),G2IO.sfloat(s[1]),G2IO.sfloat(s[2])])
    678678                    if abs(pfType) == 4:
    679                         data.extend([G2IO.sfloat(s[3]),0.0,0.0,azm])
     679                        data.extend([0.0,G2IO.sfloat(s[3]),0.0,0.0,azm])
    680680                    else:
    681681                        s = Iparm['INS  1PRCF12'].split()
    682                         data.extend([G2IO.sfloat(s[0]),0.0,0.0,azm])                       
     682                        data.extend([0.0,G2IO.sfloat(s[0]),0.0,0.0,azm])                       
    683683                Inst1 = G2IO.makeInstDict(names,data,codes)
    684684                Inst2 = {}
     
    692692                        Inst2['Pdabc'].append([float(t) for t in s])
    693693                    Inst2['Pdabc'] = np.array(Inst2['Pdabc'])
    694                     Inst2['Pdabc'][3] += Inst2['Pdabc'][0]*Inst1['difC'][0] #turn 3rd col into TOF
     694                    Inst2['Pdabc'].T[3] += Inst2['Pdabc'].T[0]*Inst1['difC'][0] #turn 3rd col into TOF
    695695                if 'INS  1I ITYP' in Iparm:
    696696                    s = Iparm['INS  1I ITYP'].split()
     
    866866            if 'T' in Iparm1['Type'][0]:
    867867                if not rd.clockWd and rd.GSAS:
    868                     rd.powderdata[0] *= 100.
     868                    rd.powderdata[0] *= 100.        #put back the CW centideg correction
    869869                cw = np.diff(rd.powderdata[0])
    870870                rd.powderdata[0] = rd.powderdata[0][:-1]+cw/2.
    871871                rd.powderdata[1] = rd.powderdata[1][:-1]/cw
    872                 rd.powderdata[2] = rd.powderdata[2][:-1]/cw**2
     872                rd.powderdata[2] = rd.powderdata[2][:-1]*cw**2  #1/var=w at this point
    873873                if 'Itype' in Iparm2:
    874874                    Ibeg = np.searchsorted(rd.powderdata[0],Iparm2['Tminmax'][0])
     
    878878                    rd.powderdata[1] = rd.powderdata[1][Ibeg:Ifin]/YI
    879879                    var = 1./rd.powderdata[2][Ibeg:Ifin]
    880                     var += rd.powderdata[1]**2+WYI
     880                    var += WYI*rd.powderdata[1]**2
    881881                    var /= YI**2
    882882                    rd.powderdata[2] = 1./var
Note: See TracChangeset for help on using the changeset viewer.