| 1 | import os, sys |
|---|
| 2 | import wx |
|---|
| 3 | from dataFitting import Data1D, Theory1D |
|---|
| 4 | from danse.common.plottools.plottables import Data2D |
|---|
| 5 | |
|---|
| 6 | from DataLoader.loader import Loader |
|---|
| 7 | |
|---|
| 8 | def choose_data_file(parent, location=None): |
|---|
| 9 | path = None |
|---|
| 10 | if location==None: |
|---|
| 11 | location = os.getcwd() |
|---|
| 12 | |
|---|
| 13 | l = Loader() |
|---|
| 14 | cards = l.get_wildcards() |
|---|
| 15 | wlist = '|'.join(cards) |
|---|
| 16 | |
|---|
| 17 | dlg = wx.FileDialog(parent, "Choose a file", location, "", wlist, wx.OPEN) |
|---|
| 18 | if dlg.ShowModal() == wx.ID_OK: |
|---|
| 19 | path = dlg.GetPath() |
|---|
| 20 | mypath = os.path.basename(path) |
|---|
| 21 | dlg.Destroy() |
|---|
| 22 | |
|---|
| 23 | return path |
|---|
| 24 | |
|---|
| 25 | |
|---|
| 26 | def load_ascii_1D(path): |
|---|
| 27 | """ |
|---|
| 28 | Load a 1D ascii file, with errors |
|---|
| 29 | """ |
|---|
| 30 | import numpy |
|---|
| 31 | if path and os.path.isfile(path): |
|---|
| 32 | |
|---|
| 33 | file_x = numpy.zeros(0) |
|---|
| 34 | file_y = numpy.zeros(0) |
|---|
| 35 | file_dy = numpy.zeros(0) |
|---|
| 36 | |
|---|
| 37 | input_f = open(path,'r') |
|---|
| 38 | buff = input_f.read() |
|---|
| 39 | lines = buff.split('\n') |
|---|
| 40 | |
|---|
| 41 | has_dy = False |
|---|
| 42 | |
|---|
| 43 | for line in lines: |
|---|
| 44 | try: |
|---|
| 45 | toks = line.split() |
|---|
| 46 | x = float(toks[0]) |
|---|
| 47 | y = float(toks[1]) |
|---|
| 48 | if len(toks)==3: |
|---|
| 49 | has_dy = True |
|---|
| 50 | err = float(toks[2]) |
|---|
| 51 | else: |
|---|
| 52 | err = 0.0 |
|---|
| 53 | file_x = numpy.append(file_x, x) |
|---|
| 54 | file_y = numpy.append(file_y, y) |
|---|
| 55 | file_dy = numpy.append(file_dy, err) |
|---|
| 56 | except: |
|---|
| 57 | print "READ ERROR", line |
|---|
| 58 | |
|---|
| 59 | if has_dy==False: |
|---|
| 60 | file_dy = None |
|---|
| 61 | |
|---|
| 62 | return file_x, file_y, file_dy |
|---|
| 63 | return None, None, None |
|---|
| 64 | |
|---|
| 65 | def plot_data(parent, path, name="Loaded Data"): |
|---|
| 66 | from sans.guicomm.events import NewPlotEvent, StatusEvent |
|---|
| 67 | |
|---|
| 68 | from DataLoader.loader import Loader |
|---|
| 69 | import numpy |
|---|
| 70 | #Instantiate a loader |
|---|
| 71 | L = Loader() |
|---|
| 72 | |
|---|
| 73 | #Recieves data |
|---|
| 74 | try: |
|---|
| 75 | output=L.load(path) |
|---|
| 76 | |
|---|
| 77 | except: |
|---|
| 78 | wx.PostEvent(parent, StatusEvent(status="Problem loading file: %s" % sys.exc_value)) |
|---|
| 79 | return |
|---|
| 80 | filename = os.path.basename(path) |
|---|
| 81 | if not output.__class__.__name__=="list": |
|---|
| 82 | try: |
|---|
| 83 | dxl=output.dxl |
|---|
| 84 | dxw=output.dxw |
|---|
| 85 | except: |
|---|
| 86 | dxl=None |
|---|
| 87 | dxw=None |
|---|
| 88 | if hasattr(output,'data'): |
|---|
| 89 | new_plot = Data2D(image=output.data,err_image=output.err_data, |
|---|
| 90 | xmin=output.xmin,xmax=output.xmax, |
|---|
| 91 | ymin=output.ymin,ymax=output.ymax) |
|---|
| 92 | new_plot.x_bins=output.x_bins |
|---|
| 93 | new_plot.y_bins=output.y_bins |
|---|
| 94 | #print "data_loader",output |
|---|
| 95 | else: |
|---|
| 96 | if not hasattr(output,"dy"): |
|---|
| 97 | new_plot = Theory1D(output.x,output.y, dxl, dxw) |
|---|
| 98 | else: |
|---|
| 99 | new_plot = Data1D(x=output.x,y=output.y,dy=output.dy, dxl=dxl, dxw=dxw) |
|---|
| 100 | #print "dataloader",output[0],output[1] |
|---|
| 101 | |
|---|
| 102 | new_plot.source=output.source |
|---|
| 103 | new_plot.name = filename |
|---|
| 104 | new_plot.interactive = True |
|---|
| 105 | new_plot.info= output |
|---|
| 106 | if hasattr(output, "dxl"): |
|---|
| 107 | new_plot.dxl = output.dxl |
|---|
| 108 | if hasattr(output, "dxw"): |
|---|
| 109 | new_plot.dxw = output.dxw |
|---|
| 110 | #print "loader output.detector",output.source |
|---|
| 111 | new_plot.detector =output.detector |
|---|
| 112 | |
|---|
| 113 | # If the data file does not tell us what the axes are, just assume... |
|---|
| 114 | new_plot.xaxis(output._xaxis,output._xunit) |
|---|
| 115 | new_plot.yaxis(output._yaxis,output._yunit) |
|---|
| 116 | new_plot.group_id = filename |
|---|
| 117 | wx.PostEvent(parent, NewPlotEvent(plot=new_plot, title=filename)) |
|---|
| 118 | else: |
|---|
| 119 | i="res" |
|---|
| 120 | for item in output: |
|---|
| 121 | try: |
|---|
| 122 | dxl=item.dxl |
|---|
| 123 | dxw=item.dxw |
|---|
| 124 | except: |
|---|
| 125 | dxl=None |
|---|
| 126 | dxw=None |
|---|
| 127 | if not hasattr(item,"dy"): |
|---|
| 128 | new_plot = Theory1D(item.x,item.y,dxl,dxw) |
|---|
| 129 | else: |
|---|
| 130 | new_plot = Data1D(x=item.x,y=item.y,dy=item.dy,dxl=dxl,dxw=dxw) |
|---|
| 131 | |
|---|
| 132 | new_plot.source=item.source |
|---|
| 133 | new_plot.info=output |
|---|
| 134 | new_plot.name = filename+" "+ str(i) |
|---|
| 135 | new_plot.interactive = True |
|---|
| 136 | |
|---|
| 137 | #print "loader output.detector",output.source |
|---|
| 138 | new_plot.detector =item.detector |
|---|
| 139 | # If the data file does not tell us what the axes are, just assume... |
|---|
| 140 | new_plot.xaxis(item._xaxis,item._xunit) |
|---|
| 141 | new_plot.yaxis(item._yaxis,item._yunit) |
|---|
| 142 | new_plot.group_id = filename |
|---|
| 143 | wx.PostEvent(parent, NewPlotEvent(plot=new_plot, title=filename)) |
|---|
| 144 | i="slit" |
|---|
| 145 | |
|---|
| 146 | |
|---|
