Changeset 1459 for trunk/GSASII.py


Ignore:
Timestamp:
Aug 8, 2014 2:59:09 PM (8 years ago)
Author:
vondreele
Message:

add instrument parameters (flight path & detector 2-theta) needed for TOF
rework reflection list for TOF
change default diff curve & reflection marker offsets
change weighting to instrument constants calibration to be 1/esd2 from peak fit positions - works a lot better
1st shot at TOF Rietveld refinement with derivatives - need to be checked for correctness (some are wrong)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASII.py

    r1453 r1459  
    979979                return [G2IO.makeInstDict(names,data,codes),{}]
    980980            elif 'T' in DataType:
    981                 names = ['Type','2-theta','difC','difA','difB','Zero','alpha','beta-0','beta-1',
    982                     'beta-q','sig-0','sig-1','sig-q','X','Y','Azimuth']
    983                 codes = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
     981                names = ['Type','fltPath','2-theta','difC','difA', 'difB','Zero','alpha','beta-0','beta-1',
     982                    'beta-q','sig-0','sig-1','sig-q','X', 'Y','Azimuth',]
     983                codes = [0,0,0,0,0, 0,0,0,0,0, 0,0,0,0,0, 0,0,]
    984984                azm = 0.
    985985                if 'INS  1DETAZM' in Iparm:
    986986                    azm = float(Iparm['INS  1DETAZM'])
     987                s = Iparm['INS   FPATH1'].split()
     988                fltPath0 = G2IO.sfloat(s[0])
    987989                s = Iparm['INS  1BNKPAR'].split()
     990                fltPath1 = G2IO.sfloat(s[0])
     991                data.extend([fltPath0+fltPath1,])               #Flight path source-sample-detector
    988992                data.extend([G2IO.sfloat(s[1]),])               #2-theta for bank
    989993                s = Iparm['INS  1 ICONS'].split()
     
    18891893        self.Offset = [0.0,0.0]
    18901894        self.delOffset = .02
    1891         self.refOffset = -100.0
     1895        self.refOffset = -1.0
    18921896        self.refDelt = .01
    18931897        self.Weight = False
     
    27912795                            if name2 == 'Peak List':
    27922796                                peaks = self.PatternTree.GetItemPyData(item2)['peaks']
    2793                                 file.write("%s \n" % (name+' Peak List'))               
     2797                                file.write("%s \n" % (name+' Peak List'))
     2798                                if len(peaks[0]) == 8:
     2799                                    file.write('%10s %12s %10s %10s\n'%('pos','int','sig','gam'))
     2800                                else:
     2801                                    file.write('%10s %12s %10s %10s %10s %10s\n'%('pos','int','alp','bet','sig','gam'))                                   
    27942802                                for peak in peaks:
    2795                                     file.write("%10.5f %12.2f %10.3f %10.3f \n" % \
    2796                                         (peak[0],peak[2],peak[4],peak[6]))
     2803                                    if len(peak) == 8:  #CW
     2804                                        file.write("%10.5f %12.2f %10.3f %10.3f \n" % \
     2805                                            (peak[0],peak[2],peak[4],peak[6]))
     2806                                    else:               #TOF - more cols
     2807                                        file.write("%10.5f %12.2f %10.3f %10.3f %10.3f %10.3f\n" % \
     2808                                            (peak[0],peak[2],peak[4],peak[6],peak[8],peak[10]))
    27972809                            item2, cookie2 = self.PatternTree.GetNextChild(item, cookie2)                           
    27982810                    item, cookie = self.PatternTree.GetNextChild(self.root, cookie)                           
Note: See TracChangeset for help on using the changeset viewer.