Changeset 1038 for trunk/exports


Ignore:
Timestamp:
Aug 16, 2013 12:30:54 PM (8 years ago)
Author:
vondreele
Message:

more work on residuals & cif files

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/exports/G2cif.py

    r1037 r1038  
    142142            R = '%.5f'%(self.OverallParms['Covariance']['Rvals']['Rwp']/100.)
    143143            WriteCIFitem('\n# OVERALL WEIGHTED R-FACTOR')
    144             if len(self.powderDict) > 1:
    145                 WriteCIFitem('_pd_proc_ls_prof_wR_factor',R)
    146                 #WriteCIFitem('_pd_proc_ls_prof_R_factor',TEXT(11:20)) # who cares!
    147             if len(self.xtalDict) > 1:
    148                 WriteCIFitem('_refine_ls_wR_factor_all',R)
     144            WriteCIFitem('_refine_ls_wR_factor_obs',R)
    149145                # _refine_ls_R_factor_all
    150146                # _refine_ls_R_factor_obs               
     
    827823                WriteCIFitem('_reflns_limit_l_max', str(int(hklmax[2])))
    828824            if hklmin is not None:
    829                 WriteCIFitem('_reflns_d_resolution_low', G2mth.ValEsd(dmax,-0.0009))
    830                 WriteCIFitem('_reflns_d_resolution_high', G2mth.ValEsd(dmin,-0.009))
     825                WriteCIFitem('_reflns_d_resolution_low  ', G2mth.ValEsd(dmax,-0.009))
     826                WriteCIFitem('_reflns_d_resolution_high ', G2mth.ValEsd(dmin,-0.009))
    831827
    832828        def WritePowderData(histlbl):
     
    893889                            )
    894890                    WriteCIFitem('loop_' +
    895                                  '\n\t_pd_phase_id' +
    896                                  '\n\t_refine_ls_R_F_factor' +
    897                                  '\n\t_refine_ls_R_Fsqd_factor')
     891                                 '\n\t_pd_proc_ls_R_F_factor' +
     892                                 '\n\t_pd_proc_ls_R_Fsqd_factor')
    898893                    for phasenam in phasebyhistDict.get(histlbl):
    899894                        pfx = str(self.Phases[phasenam]['pId'])+':'+str(hId)+':'
    900895                        WriteCIFitem(
    901896                            '  '+
    902                             str(self.Phases[phasenam]['pId']) +
    903897                            '  '+G2mth.ValEsd(histblk[pfx+'Rf']/100.,-.00009) +
    904898                            '  '+G2mth.ValEsd(histblk[pfx+'Rf^2']/100.,-.00009)
     
    906900            else:
    907901                pfx = '0:'+str(hId)+':'
    908                 WriteCIFitem('_refine_ls_R_F_factor','%.5f'%(histblk[pfx+'Rf']/100.))
    909                 WriteCIFitem('_refine_ls_R_Fsqd_factor','%.5f'%(histblk[pfx+'Rf^2']/100.))
     902                WriteCIFitem('_pd_proc_ls_R_F_factor      ','%.5f'%(histblk[pfx+'Rf']/100.))
     903                WriteCIFitem('_pd_proc_ls_R_Fsqd_factor   ','%.5f'%(histblk[pfx+'Rf^2']/100.))
    910904               
    911             # WriteCIFitem('_pd_proc_ls_prof_R_factor','?')
    912             WriteCIFitem('_pd_proc_ls_prof_wR_factor','%.5f'%(histblk['wR']/100.))
    913             # WriteCIFitem('_pd_proc_ls_prof_wR_expected','?')
     905            WriteCIFitem('_pd_proc_ls_prof_R_factor   ','%.5f'%(histblk['R']/100.))
     906            WriteCIFitem('_pd_proc_ls_prof_wR_factor  ','%.5f'%(histblk['wR']/100.))
     907            WriteCIFitem('_pd_proc_ls_prof_R_B_factor ','%.5f'%(histblk['Rb']/100.))
     908            WriteCIFitem('_pd_proc_ls_prof_wR_B_factor','%.5f'%(histblk['wRb']/100.))
     909            WriteCIFitem('_pd_proc_ls_prof_wR_expected','%.5f'%(histblk['wRmin']/100.))
    914910
    915911            if histblk['Instrument Parameters'][0]['Type'][1][1] == 'X':
     
    11301126                WriteCIFitem(s)
    11311127            WriteReflStat(refcount,hklmin,hklmax,dmin,dmax)
     1128            hId = histblk['hId']
     1129            pfx = '0:'+str(hId)+':'
     1130            WriteCIFitem('_reflns_wR_factor_obs    ','%.4f'%(histblk['wR']/100.))
     1131            WriteCIFitem('_reflns_R_F_factor_obs   ','%.4f'%(histblk[pfx+'Rf']/100.))
     1132            WriteCIFitem('_reflns_R_Fsqd_factor_obs','%.4f'%(histblk[pfx+'Rf^2']/100.))
    11321133
    11331134        #============================================================
Note: See TracChangeset for help on using the changeset viewer.