source: trunk/exports/G2export_examples.py @ 1902

#!/usr/bin/env python
# -*- coding: utf-8 -*-
########### SVN repository information ###################
# $Date: 2015-06-22 16:42:30 +0000 (Mon, 22 Jun 2015) $
# $Author: vondreele $
# $Revision: 1902 $
# $URL: trunk/exports/G2export_examples.py $
# $Id: G2export_examples.py 1902 2015-06-22 16:42:30Z vondreele $
########### SVN repository information ###################
'''
*Module G2export_examples: Examples*
-------------------------------------------

Code to demonstrate how GSAS-II data export routines are created. The
classes defined here, :class:`ExportPhaseText`, 
:class:`ExportSingleText`, :class:`ExportPowderReflText`, 
and :class:`ExportPowderText` each demonstrate a different type
of export. Also see :class:`G2export_map.ExportMapASCII` for an
example of a map export.

'''
import os.path
import numpy as np
import GSASIIpath
GSASIIpath.SetVersionNumber("$Revision: 1902 $")
import GSASIIIO as G2IO
import GSASIIpy3 as G2py3
#import GSASIIgrid as G2gd
#import GSASIIstrIO as G2stIO
import GSASIImath as G2mth
import GSASIIpwd as G2pwd
#import GSASIIlattice as G2lat
#import GSASIIspc as G2spc
#import GSASIIphsGUI as G2pg
#import GSASIIstrMain as G2stMn

class ExportPhaseText(G2IO.ExportBaseclass):
    '''Used to create a text file for a phase

    :param wx.Frame G2frame: reference to main GSAS-II frame
    '''
    def __init__(self,G2frame):
        super(self.__class__,self).__init__( # fancy way to say <parentclass>.__init__
            G2frame=G2frame,
            formatName = 'Text file',
            extension='.txt',
            longFormatName = 'Export phase as text file'
            )
        self.exporttype = ['phase']
        self.multiple = True # allow multiple phases to be selected

    def Exporter(self,event=None):
        '''Export a phase as a text file
        '''
        # the export process starts here
        self.InitExport(event)
        # load all of the tree into a set of dicts
        self.loadTree()
        # create a dict with refined values and their uncertainties
        self.loadParmDict()
        if self.ExportSelect(): return # set export parameters; prompt for file name
        self.OpenFile()
        # if more than one format is selected, put them into a single file
        for phasenam in self.phasenam:
            phasedict = self.Phases[phasenam] # pointer to current phase info            
            i = self.Phases[phasenam]['pId']
            self.Write('\n'+80*'=')
            self.Write("Phase "+str(phasenam)+" from "+str(self.G2frame.GSASprojectfile))
            self.Write("\nSpace group = "+str(phasedict['General']['SGData']['SpGrp'].strip()))
            # get cell parameters & print them
            cellList,cellSig = self.GetCell(phasenam)
            prevsig = 0
            for lbl,defsig,val,sig in zip(
                ['a','b','c','alpha','beta ','gamma','volume'],
                3*[-0.00001] + 3*[-0.001] + [-0.01], # sign values to use when no sigma
                cellList,cellSig
                ):
                if sig:
                    txt = G2mth.ValEsd(val,sig)
                    prevsig = -sig # use this as the significance for next value
                else:
                    txt = G2mth.ValEsd(val,min(defsig,prevsig),True)
                self.Write(lbl+' = '+txt)
            # get atoms and print them in nice columns
            AtomsList = self.GetAtoms(phasenam)
            fmt = "{:8s} {:4s} {:4s} {:12s} {:12s} {:12s} {:10s} {:10s}"
            self.Write('\nAtoms\n'+80*'-')
            self.Write(fmt.format("label","elem","mult","x","y","z","frac","Uiso"))
            self.Write(80*'-')
            aniso = False
            for lbl,typ,mult,xyz,td in AtomsList:
                vals = [lbl,typ,str(mult)]
                if xyz[3][0] == 0: continue
                for val,sig in xyz:
                    vals.append(G2mth.ValEsd(val,sig))
                if len(td) == 1:
                    vals.append(G2mth.ValEsd(td[0][0],td[0][1]))
                else:
                    vals.append("aniso")
                    aniso = True
                self.Write(fmt.format(*vals))
            # print anisotropic values, if any
            if aniso:
                self.Write('\nAnisotropic parameters')
                self.Write(80*'-')
                fmt = "{:8s} {:4s} {:10s} {:10s} {:10s} {:10s} {:10s} {:10s}"
                self.Write(fmt.format("label","elem",'U11','U22','U33','U12','U13','U23'))
                self.Write(80*'-')
                for lbl,typ,mult,xyz,td in AtomsList:
                    if len(td) == 1: continue
                    if xyz[3][0] == 0: continue
                    vals = [lbl,typ]
                    for val,sig in td:
                        vals.append(G2mth.ValEsd(val,sig))
                    self.Write(fmt.format(*vals))
            print('Phase '+str(phasenam)+' written to file '+str(self.fullpath))
        self.CloseFile()

class ExportPowderText(G2IO.ExportBaseclass):
    '''Used to create a text file for a powder data set

    :param wx.Frame G2frame: reference to main GSAS-II frame
    '''
    def __init__(self,G2frame):
        super(self.__class__,self).__init__( # fancy way to say <parentclass>.__init__
            G2frame=G2frame,
            formatName = 'Text file',
            extension='.txt',
            longFormatName = 'Export powder data as text file'
            )
        self.exporttype = ['powder']
        self.multiple = False # only allow one histogram to be selected

    def Exporter(self,event=None):
        '''Export a set of powder data as a text file
        '''
        # the export process starts here
        self.InitExport(event)
        # load all of the tree into a set of dicts
        self.loadTree()
        if self.ExportSelect( # set export parameters
            AskFile='default' # base name on the GPX file name
            ): return 
        self.OpenFile()
        hist = self.histnam[0] # there should only be one histogram, in any case take the 1st
        histblk = self.Histograms[hist]
        hfmt = 5*"{:12s} "
        digitList = 2*((13,3),) + ((13,5),) + 2*((13,3),)
        
        self.Write(hfmt.format("x","y_obs","weight","y_calc","y_bkg"))
        for vallist in zip(histblk['Data'][0],
                           histblk['Data'][1],
                           histblk['Data'][2],
                           histblk['Data'][3],
                           histblk['Data'][4],
                           #histblk['Data'][5],
                           ):
            strg = ''
            for val,digits in zip(vallist,digitList):
                strg += G2py3.FormatPadValue(val,digits)
            self.Write(strg)
        self.CloseFile()
        print(str(hist)+' written to file '+str(self.fullpath))
        
class ExportPowderReflText(G2IO.ExportBaseclass):
    '''Used to create a text file of reflections from a powder data set

    :param wx.Frame G2frame: reference to main GSAS-II frame
    '''
    def __init__(self,G2frame):
        super(self.__class__,self).__init__( # fancy way to say <parentclass>.__init__
            G2frame=G2frame,
            formatName = 'reflection list as text',
            extension='.txt',
            longFormatName = 'Export powder reflection list as a text file'
            )
        self.exporttype = ['powder']
        self.multiple = False # only allow one histogram to be selected

    def Exporter(self,event=None):
        '''Export a set of powder reflections as a text file
        '''
        self.InitExport(event)
        # load all of the tree into a set of dicts
        self.loadTree()
        if self.ExportSelect( # set export parameters
            AskFile='default' # base name on the GPX file name
            ): return 
        self.OpenFile()
        hist = self.histnam[0] # there should only be one histogram, in any case take the 1st
        self.Write('\nHistogram '+hist)
        histblk = self.Histograms[hist]
        for phasenam in histblk['Reflection Lists']:
            phasDict = histblk['Reflection Lists'][phasenam]
            tname = {'T':'TOF','C':'2-theta'}[phasDict['Type'][2]]
            self.Write('\nPhase '+str(phasenam))
            if phasDict.get('Super',False):
                self.Write(96*'=')
                hklfmt = "{:.0f},{:.0f},{:.0f},{:.0f}"
                hfmt = "{:>10s} {:>8s} {:>12s} {:>12s} {:>7s} {:>6s} {:>8s} {:>8s} {:>8s} {:>8s}"
                if 'T' in phasDict['Type']:
                    fmt = "{:>10s} {:8.3f} {:12.3f} {:12.3f} {:7.2f} {:6.0f} {:8.3f} {:8.3f} {:8.3f} {:8.4f}"
                else:
                    fmt = "{:>10s} {:8.3f} {:12.3f} {:12.3f} {:7.2f} {:6.0f} {:8.5f} {:8.5f} {:8.5f} {:8.4f}"
                self.Write(hfmt.format("h,k,l,m",tname,"F_obs","F_calc","phase","mult","sig","gam","FWHM","Prfo"))
                self.Write(96*'=')
                refList = phasDict['RefList']
                for refItem in refList:
                    if 'T' in phasDict['Type']:
                        h,k,l,m,mult,dsp,pos,sig,gam,Fobs,Fcalc,phase,x,x,x,x,prfo = refItem[:17]
                        FWHM = 2.*G2pwd.getgamFW(gam,sig)
                        self.Write(fmt.format(hklfmt.format(h,k,l,m),pos,Fobs,Fcalc,phase,mult,sig,gam,FWHM,prfo))
                    else:
                        h,k,l,m,mult,dsp,pos,sig,gam,Fobs,Fcalc,phase,x,prfo = refItem[:14]
                        FWHM = 2.*G2pwd.getgamFW(gam,sig)
                        self.Write(fmt.format(hklfmt.format(h,k,l,m),pos,Fobs,Fcalc,phase,mult, \
                            np.sqrt(max(sig,0.0001))/100.,gam/100.,FWHM/100.,prfo))
            else:
                self.Write(94*'=')
                hklfmt = "{:.0f},{:.0f},{:.0f}"
                hfmt = "{:>8s} {:>8s} {:>12s} {:>12s} {:>7s} {:>6s} {:>8s} {:>8s} {:>8s} {:>8s}"
                if 'T' in phasDict['Type']:
                    fmt = "{:>8s} {:8.3f} {:12.3f} {:12.3f} {:7.2f} {:6.0f} {:8.3f} {:8.3f} {:8.3f} {:8.4f}"
                else:
                    fmt = "{:>8s} {:8.3f} {:12.3f} {:12.3f} {:7.2f} {:6.0f} {:8.5f} {:8.5f} {:8.5f} {:8.4f}"
                self.Write(hfmt.format("h,k,l",tname,"F_obs","F_calc","phase","mult","sig","gam","FWHM","Prfo"))
                self.Write(94*'=')
                refList = phasDict['RefList']
                for refItem in refList:
                    if 'T' in phasDict['Type']:
                        h,k,l,mult,dsp,pos,sig,gam,Fobs,Fcalc,phase,x,x,x,x,prfo = refItem[:16]
                        FWHM = 2.*G2pwd.getgamFW(gam,sig)
                        self.Write(fmt.format(hklfmt.format(h,k,l),pos,Fobs,Fcalc,phase,mult,np.sqrt(max(sig,0.0001)),gam,FWHM,prfo))
                    else:
                        h,k,l,mult,dsp,pos,sig,gam,Fobs,Fcalc,phase,x,prfo = refItem[:13]
                        FWHM = 2.*G2pwd.getgamFW(gam,sig)
                        self.Write(fmt.format(hklfmt.format(h,k,l),pos,Fobs,Fcalc,phase,mult,   \
                            np.sqrt(max(sig,0.0001))/100.,gam/100.,FWHM/100.,prfo))
        self.CloseFile()
        print(str(hist)+'reflections written to file '+str(self.fullpath))                        

class ExportSingleText(G2IO.ExportBaseclass):
    '''Used to create a text file with single crystal reflection data

    :param wx.Frame G2frame: reference to main GSAS-II frame
    '''
    def __init__(self,G2frame):
        super(self.__class__,self).__init__( # fancy way to say <parentclass>.__init__
            G2frame=G2frame,
            formatName = 'Text file',
            extension='.txt',
            longFormatName = 'Export reflection list as a text file'
            )
        self.exporttype = ['single']
        self.multiple = False # only allow one histogram to be selected

    def Exporter(self,event=None):
        '''Export a set of single crystal data as a text file
        '''
        # the export process starts here
        self.InitExport(event)
        # load all of the tree into a set of dicts
        self.loadTree()
        if self.ExportSelect( # set export parameters
            AskFile='default' # base name on the GPX file name
            ): return 
        self.OpenFile()
        hist = self.histnam[0] # there should only be one histogram, in any case take the 1st
        histblk = self.Histograms[hist]
        hklfmt = "{:.0f},{:.0f},{:.0f}"
        hfmt = "{:>10s} {:>8s} {:>12s} {:>12s} {:>12s} {:>7s} {:>6s}"
        fmt = "{:>10s} {:8.3f} {:12.2f} {:12.4f} {:12.2f} {:7.2f} {:6.0f}"
        self.Write(80*'=')
        self.Write(hfmt.format("h,k,l","d-space","F_obs","sig(Fobs)","F_calc","phase","mult"))
        self.Write(80*'=')
        for (
            h,k,l,mult,dsp,Fobs,sigFobs,Fcalc,FobsT,FcalcT,phase,Icorr
            ) in histblk['Data']['RefList']:
            self.Write(fmt.format(hklfmt.format(h,k,l),dsp,Fobs,sigFobs,Fcalc,phase,mult))
        self.CloseFile()
        print(str(hist)+' written to file '+str(self.fullpath))