source: sans/Dev/trunk/NCNR_User_Procedures/Common/NIST_XML_v40.ipf @ 481

Last change on this file since 481 was 481, checked in by ajj, 14 years ago

Sundry changes...

File size: 10.9 KB
RevLine 
[403]1#pragma rtGlobals=1             // Use modern global access method.
2#pragma version=1.00
3#pragma IgorVersion=6.0
4
5#include "cansasXML_v11", version >= 1.10
6
7
[461]8// The function is called "LoadNISTXMLData" but is actually more generic
9// and will load any canSAS XML v1 dataset
10// Dec 2008 :
11//                       Caveats - Assumes Q in /A and I in /cm
[403]12function LoadNISTXMLData(filestr,doPlot)
13        String filestr
14        Variable doPlot
15       
16       
17        Variable rr,gg,bb
18        NVAR dQv = root:Packages:NIST:USANS_dQv
19
20               
21        Print "Trying to load canSAS XML format data"
22        Variable result = CS_XMLReader(filestr)
23       
24        String xmlReaderFolder = "root:Packages:CS_XMLreader:"
25       
26        if (result == 0)
27                        SetDataFolder xmlReaderFolder
28                                               
[461]29                        Variable i,j,numDataSets
30                        Variable np
[403]31                       
[461]32                        String w0,w1,w2,basestr,fileName
33                        String xmlDataFolder,xmlDataSetFolder
34                       
35                       
36                       
[403]37                        for (i = 0; i < CountObjects(xmlReaderFolder,4); i+=1)
38                                                               
[461]39                                xmlDataFolder = xmlReaderFolder+GetIndexedObjName(xmlReaderFolder,4,i)+":"
40                                numDataSets = CountObjects(xmlDataFolder,4)
41                                if (numDataSets > 0)
42                                        //Multiple SASData sets in this SASEntry
43                                        for (j = 0; j < numDataSets; j+=1)
44                                               
45                                                xmlDataSetFolder = xmlDataFolder+GetIndexedObjName(xmlDataFolder,4,j)+":"
46                                               
47                                                SetDataFolder xmlDataSetFolder
48                                       
49                                                basestr = CleanupName(getXMLDataSetTitle(xmlDataSetFolder,j),0)
50                                                //String basestr = ParseFilePath(3, ParseFilePath(5,filestr,":",0,0),":",0,0)                           
51                                                fileName =  ParseFilePath(0,ParseFilePath(5,filestr,":",0,0),":",1,0)
52                                                       
53                                                //print "In NIST XML Loader"
54                                                //print "fileStr: ",fileStr
55                                                //print "basestr: ",basestr
56                                                //print "fileName: ",fileName
57                                                //remove the semicolon AND period from files from the VAX
58                                                w0 = basestr + "_q"
59                                                w1 = basestr + "_i"
60                                                w2 = basestr + "_s"
61                                               
62                                                if(DataFolderExists("root:"+baseStr))
63                                                                DoAlert 1,"The data set " + basestr + " from file "+fileName+" has already been loaded. Do you want to load the new data file, overwriting the data in memory?"
64                                                                if(V_flag==2)   //user selected No, don't load the data
65                                                                        SetDataFolder root:
66                                                                        if(DataFolderExists("root:Packages:NIST"))
67                                                                                String/G root:Packages:NIST:gLastFileName = filename
68                                                                        endif           //set the last file loaded to the one NOT loaded
69                                                                        return  0       //quits the macro
70                                                                endif
71                                                                SetDataFolder $("root:"+baseStr)
72                                                else
73                                                        NewDataFolder/S $("root:"+baseStr)
74                                                endif
[403]75                       
[461]76                                                Duplicate/O $(xmlDataSetFolder+"Qsas") $w0
77                                                Duplicate/O $(xmlDataSetFolder+"Isas") $w1
78                                                Duplicate/O $(xmlDataSetFolder+"Idev") $w2
79               
80               
81                                                       
82                                                if (exists(xmlDataSetFolder+"Qdev"))
83                                                        Wave Qsas = $(xmlDataSetFolder+"Qsas")
84                                                        Wave Qdev = $(xmlDataSetFolder+"Qdev")
85                                               
86                                                // make a resolution matrix for SANS data
87                                                         np=numpnts($w0)
88                                                        Make/D/O/N=(np,4) $(baseStr+"_res")
89                                                        Wave reswave =  $(baseStr+"_res")
90                                                       
91                                                        reswave[][0] = Qdev[p]          //sigQ
92                                                        reswave[][3] = Qsas[p]  //Qvalues
93                                                        if(exists(xmlDataSetFolder+"Qmean"))
94                                                                Wave Qmean = $(xmlDataSetFolder+"Qmean")
95                                                                reswave[][1] = Qmean[p]         //qBar
96                                                        endif
97                                                        if(exists(xmlDataSetFolder+"Shadowfactor"))
98                                                                Wave Shadowfactor = $(xmlDataSetFolder+"Shadowfactor")
99                                                                reswave[][2] = Shadowfactor[p]          //fShad
100                                                        endif
101                                                elseif(exists(xmlDataSetFolder+"dQl"))
102                                                        //USAS Data
103                                                        Wave dQl = $(xmlDataSetFolder+"dQl")
104                                                        dQv = dQl[0]
105                                               
106                                                        USANS_CalcWeights(baseStr,dQv)
107                                                else
108                                                        //No resolution data
109                                                endif
110                                                        //get rid of the resolution waves that are in the matrix
[403]111                                       
[461]112                                                        SetScale d,0,0,"1/A",$w0
113                                                        SetScale d,0,0,"1/cm",$w1
114                                               
115                                                       
[403]116                               
[461]117                                                //////
118                                                if(DataFolderExists("root:Packages:NIST"))
119                                                        String/G root:Packages:NIST:gLastFileName = filename
[403]120                                                endif
[461]121                                       
122                                               
123                                                //plot if desired
124                                                if(doPlot)
125                                                        // assign colors randomly
126                                                        rr = abs(trunc(enoise(65535)))
127                                                        gg = abs(trunc(enoise(65535)))
128                                                        bb = abs(trunc(enoise(65535)))
129                                                       
130                                                        // if target window is a graph, and user wants to append, do so
131                                                   DoWindow/B Plot_Manager
132                                                        if(WinType("") == 1)
133                                                                DoAlert 1,"Do you want to append this data to the current graph?"
134                                                                if(V_Flag == 1)
135                                                                        AppendToGraph $w1 vs $w0
136                                                                        ModifyGraph mode($w1)=3,marker($w1)=19,msize($w1)=2,rgb($w1) =(rr,gg,bb),tickUnit=1
[477]137                                                                        ErrorBars/T=0 $w1 Y,wave=($w2,$w2)
[461]138                                                                        ModifyGraph tickUnit(left)=1
139                                                                else
140                                                                //new graph
141                                                                        Display $w1 vs $w0
142                                                                        ModifyGraph log=1,mode($w1)=3,marker($w1)=19,msize($w1)=2,rgb($w1)=(rr,gg,bb),tickUnit=1
143                                                                        ModifyGraph grid=1,mirror=2,standoff=0
[477]144                                                                        ErrorBars/T=0 $w1 Y,wave=($w2,$w2)
[461]145                                                                        ModifyGraph tickUnit(left)=1
146                                                                        Legend
147                                                                endif
148                                                        else
149                                                        // graph window was not target, make new one
150                                                                Display $w1 vs $w0
151                                                                ModifyGraph log=1,mode($w1)=3,marker($w1)=19,msize($w1)=2,rgb($w1)=(rr,gg,bb),tickUnit=1
152                                                                ModifyGraph grid=1,mirror=2,standoff=0
[477]153                                                                ErrorBars/T=0 $w1 Y,wave=($w2,$w2)
[461]154                                                                ModifyGraph tickUnit(left)=1
155                                                                Legend
156                                                        endif
157                                                endif
158                                       
159                                        endfor
160                                       
161                                       
[403]162                                else
[461]163                                        //No multiple SASData sets for this SASEntry
164                                        SetDataFolder xmlDataFolder
165                                       
166                                        basestr = CleanupName(getXMLDataSetTitle(xmlDataFolder,0),0)
167                                        //String basestr = ParseFilePath(3, ParseFilePath(5,filestr,":",0,0),":",0,0)                           
168                                        fileName =  ParseFilePath(0,ParseFilePath(5,filestr,":",0,0),":",1,0)
169                                               
170                                        //print "In NIST XML Loader"
171                                        //print "fileStr: ",fileStr
172                                        //print "basestr: ",basestr
173                                        //print "fileName: ",fileName
174                                        w0 = basestr + "_q"
175                                        w1 = basestr + "_i"
176                                        w2 = basestr + "_s"
177                                       
178                                        if(DataFolderExists("root:"+baseStr))
179                                                        DoAlert 1,"The data set " + basestr + " from file "+fileName+" has already been loaded. Do you want to load the new data file, overwriting the data in memory?"
180                                                        if(V_flag==2)   //user selected No, don't load the data
181                                                                SetDataFolder root:
182                                                                if(DataFolderExists("root:Packages:NIST"))
183                                                                        String/G root:Packages:NIST:gLastFileName = filename
184                                                                endif           //set the last file loaded to the one NOT loaded
185                                                                return  0       //quits the macro
186                                                        endif
187                                                        SetDataFolder $("root:"+baseStr)
188                                        else
189                                                NewDataFolder/S $("root:"+baseStr)
190                                        endif
191               
192                                        Duplicate/O $(xmlDataFolder+"Qsas") $w0
193                                        Duplicate/O $(xmlDataFolder+"Isas") $w1
194                                        Duplicate/O $(xmlDataFolder+"Idev") $w2
[403]195       
[461]196       
197                                               
198                                        if (exists(xmlDataFolder+"Qdev"))
199                                                Wave Qsas = $(xmlDataFolder+"Qsas")
200                                                Wave Qdev = $(xmlDataFolder+"Qdev")
[403]201                                       
[461]202                                        // make a resolution matrix for SANS data
203                                                np=numpnts($w0)
204                                                Make/D/O/N=(np,4) $(baseStr+"_res")
205                                                Wave reswave =  $(baseStr+"_res")
206                                               
207                                                reswave[][0] = Qdev[p]          //sigQ
208                                                reswave[][3] = Qsas[p]  //Qvalues
209                                                if(exists(xmlDataFolder+"Qmean"))
210                                                        Wave Qmean = $(xmlDataFolder+"Qmean")
211                                                        reswave[][1] = Qmean[p]         //qBar
212                                                endif
213                                                if(exists(xmlDataFolder+"Shadowfactor"))
214                                                        Wave Shadowfactor = $(xmlDataFolder+"Shadowfactor")
215                                                        reswave[][2] = Shadowfactor[p]          //fShad
216                                                endif
217                                        elseif(exists(xmlDataFolder+"dQl"))
218                                                //USAS Data
219                                                Wave dQl = $(xmlDataFolder+"dQl")
220                                                dQv = dQl[0]
[403]221                                       
[461]222                                                USANS_CalcWeights(baseStr,dQv)
223                                        else
224                                                //No resolution data
[403]225                                        endif
[461]226                                                //get rid of the resolution waves that are in the matrix
[403]227                               
[461]228                                                SetScale d,0,0,"1/A",$w0
229                                                SetScale d,0,0,"1/cm",$w1
[403]230                                       
[461]231                                               
[403]232                       
[461]233                                        //////
234                                        if(DataFolderExists("root:Packages:NIST"))
235                                                String/G root:Packages:NIST:gLastFileName = filename
236                                        endif
[403]237                               
238                                       
[461]239                                        //plot if desired
240                                        if(doPlot)
241                                                // assign colors randomly
242                                                rr = abs(trunc(enoise(65535)))
243                                                gg = abs(trunc(enoise(65535)))
244                                                bb = abs(trunc(enoise(65535)))
245                                               
246                                                // if target window is a graph, and user wants to append, do so
247                                           DoWindow/B Plot_Manager
248                                                if(WinType("") == 1)
249                                                        DoAlert 1,"Do you want to append this data to the current graph?"
250                                                        if(V_Flag == 1)
251                                                                AppendToGraph $w1 vs $w0
252                                                                ModifyGraph mode($w1)=3,marker($w1)=19,msize($w1)=2,rgb($w1) =(rr,gg,bb),tickUnit=1
[477]253                                                                ErrorBars/T=0 $w1 Y,wave=($w2,$w2)
[461]254                                                                ModifyGraph tickUnit(left)=1
255                                                        else
256                                                        //new graph
257                                                                Display $w1 vs $w0
258                                                                ModifyGraph log=1,mode($w1)=3,marker($w1)=19,msize($w1)=2,rgb($w1)=(rr,gg,bb),tickUnit=1
259                                                                ModifyGraph grid=1,mirror=2,standoff=0
[477]260                                                                ErrorBars/T=0 $w1 Y,wave=($w2,$w2)
[461]261                                                                ModifyGraph tickUnit(left)=1
262                                                                Legend
263                                                        endif
[403]264                                                else
[461]265                                                // graph window was not target, make new one
[403]266                                                        Display $w1 vs $w0
267                                                        ModifyGraph log=1,mode($w1)=3,marker($w1)=19,msize($w1)=2,rgb($w1)=(rr,gg,bb),tickUnit=1
268                                                        ModifyGraph grid=1,mirror=2,standoff=0
[477]269                                                        ErrorBars/T=0 $w1 Y,wave=($w2,$w2)
[403]270                                                        ModifyGraph tickUnit(left)=1
271                                                        Legend
272                                                endif
273                                        endif
[461]274                               
[403]275                                endif
276                        endfor
277        endif
278
279        //go back to the root folder and clean up before leaving
280        SetDataFolder root:
281        //KillDataFolder xmlReaderFolder
282       
283
284end
285
286
[461]287function/S getXMLDataSetTitle(xmlDF,dsNum)
[403]288        String xmlDF
[461]289        Variable dsNum
[403]290
[481]291        SVAR title = root:Packages:NIST:gXMLLoader_Title
[461]292
293
[403]294        String mdstring = xmlDF+"metadata"
295
296        Wave/T meta = $mdstring
[461]297        //Check for value when there are multiple datasets
298        FindValue/TEXT="Title"/TXOP=4/Z meta
[481]299        title = TrimWS(meta[V_Value][1])
300        print meta[V_Value][1]
301        print title
[461]302         //Check for Run value
303         FindValue/TEXT="Run_"+num2str(dsNum)/TXOP=4/Z meta
304        if (V_Value >= 0)
[481]305                title = title+" "+TrimWS(meta[V_Value][1])
306                print title
[461]307        else
308                FindValue/TEXT="Run"/TXOP=4/Z meta
[481]309                title = title+" "+TrimWS(meta[V_Value][1])     
310                print title
[461]311        endif
312
313        if (strlen(title) > 28)
314                //Prompt title, "Set New Sample Name"
315                //DoPrompt "Sample Name Is Too Long", title
316                do
317                        Execute "getXMLShorterTitle()"
318                while (strlen(title) > 28)                     
319        endif
[403]320       
321        return title
322end
323
324function isXML(filestr)
325        String filestr
326       
327        String line
328        Variable fileref
329       
330        Open/R fileref as filestr
331        FReadLine fileref,  line
332        Close fileref
333       
334        //Hopefully this will distinguish between other formats and the XML
335        return stringmatch(line, "*xml*")       
336
[461]337end
338
339
340Proc getXMLShorterTitle()
341       
342         //NVAR title = root:myGlobals:gXMLLoader_Title
343       
344        DoWindow/K getNewTitle
345        getNewTitle()
346
347        PauseforUser getNewTitle
348end
349
350Window getNewTitle()
351
352        PauseUpdate; Silent 1           // building window...
353        NewPanel /W=(166,90,666,230) as "Sample Title Too Long!"
354        SetDrawLayer UserBack
355        DrawText 11,22,"The sample title is too long."
356        DrawText 11,42,"Please enter a new one with a maximum length 28 characters"
357        DrawText 11,72,"Current Sample Title:"
358        GroupBox group0 pos={8,55},size={484,50}
[481]359        TitleBox tb_CurrentTitle,pos={150,57}, variable=root:Packages:NIST:gXMLLoader_Title,fSize=12,frame=0   
[461]360        SetVariable sv_NewTitle,pos={11,77},size={476,18},title="New Sample Title"
[481]361        SetVariable sv_NewTitle,fSize=12,value=root:Packages:NIST:gXMLLoader_Title
[461]362        Button btn_SetNewTitle title="Set New Title",pos={150,110},size={200,20}
363        Button btn_SetNewTitle proc=SetNewTitleButtonProc
364
365EndMacro
366
367
368Proc SetNewTitleButtonProc(ctrlName) : ButtonControl
369        String ctrlName
370
371        DoWindow/K getNewTitle
372
373End
Note: See TracBrowser for help on using the repository browser.