# -*- coding: utf-8 -*-
'''Classes in :mod:`~GSASII.exports.G2export_FIT2D` follow:
'''
from __future__ import division, print_function
import os.path
import numpy as np
from .. import GSASIIfiles as G2fil
from .. import GSASIIobj as G2obj
[docs]
class ExportPowderCHI(G2fil.ExportBaseclass):
'''Used to create a CHI 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 = 'Fit2D chi file',
extension='.chi',
longFormatName = 'Export powder data as Fit2D .chi file'
)
self.exporttype = ['powder']
self.multiple = True
def Writer(self,TreeName,filename=None):
self.OpenFile(filename)
histblk = self.Histograms[TreeName]
self.Write(str(TreeName)[5:]) # drop 'PWDR '
self.Write("2-Theta Angle (Degrees)")
self.Write("Intensity")
self.Write(" "+str(len(histblk['Data'][0])))
for X,Y in zip(histblk['Data'][0],histblk['Data'][1]):
line = " %5.7e" % X
line += " %5.7e" % Y
self.Write(line)
self.CloseFile()
[docs]
def Exporter(self,event=None):
'''Export a set of powder data as a Fit2D .qchi 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='single' # get a file name/directory to save in
): return
filenamelist = []
for hist in self.histnam:
if len(self.histnam) == 1:
name = self.filename
else: # multiple files: create a unique name from the histogram
name = self.MakePWDRfilename(hist)
fileroot = os.path.splitext(G2obj.MakeUniqueLabel(name,filenamelist))[0]
self.filename = os.path.join(self.dirname,fileroot + self.extension)
self.Writer(hist)
print('Histogram '+hist+' written to file '+self.fullpath)
[docs]
class ExportPowderQCHI(G2fil.ExportBaseclass):
'''Used to create a q-binned CHI 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 = 'Fit2D q-bin chi file',
extension='.qchi',
longFormatName = 'Export powder data as q-bin Fit2D .qchi file'
)
self.exporttype = ['powder']
self.multiple = True
def Writer(self,TreeName,filename=None):
from .. import GSASIIlattice as G2lat
self.OpenFile(filename)
histblk = self.Histograms[TreeName]
inst = histblk['Instrument Parameters'][0]
self.Write(str(TreeName)[5:]) # drop 'PWDR '
if 'Lam1' in inst:
print('Do you really want to write a multi-wavelength pattern in Q?')
lam = 0.
else:
lam = inst['Lam'][1]
self.Write("Q{:>20.6f}".format(lam))
self.Write("Intensity")
self.Write(" "+str(len(histblk['Data'][0])))
for X,Y in zip(histblk['Data'][0],histblk['Data'][1]):
line = " %5.7e" % (2.*np.pi/G2lat.Pos2dsp(inst,X))
line += " %5.7e" % Y
self.Write(line)
self.CloseFile()
[docs]
def Exporter(self,event=None):
'''Export a set of powder data as a q-bin Fit2D .qchi 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='single' # get a file name/directory to save in
): return
filenamelist = []
for hist in self.histnam:
if len(self.histnam) == 1:
name = self.filename
else: # multiple files: create a unique name from the histogram
name = self.MakePWDRfilename(hist)
fileroot = os.path.splitext(G2obj.MakeUniqueLabel(name,filenamelist))[0]
self.filename = os.path.join(self.dirname,fileroot + self.extension)
self.Writer(hist)
print('Histogram '+hist+' written to file '+self.fullpath)