Changeset 238

Timestamp:

Jan 14, 2008 4:25:16 PM (15 years ago)

Author:

ajj

Message:

Update to USANS slit-smearing on data load.

File:

• sans/Analysis/branches/ajj_23APR07/IGOR_Package_Files/Put in User Procedures/SANS_Models_v3.00/USANS_SlitSmearing.ipf (modified) (10 diffs)

sans/Analysis/branches/ajj_23APR07/IGOR_Package_Files/Put in User Procedures/SANS_Models_v3.00/USANS_SlitSmearing.ipf

@@ -22 +22 @@
         Variable/G USANS_N=numpnts($(basestr+"_q"))
         Variable/G USANS_dQv = dQv
-        String/G USANS_basestr = basestr
+        String/G root:myGlobals:USANS_basestr = basestr
 
         Make/D/O/N=(USANS_N,USANS_N) $(basestr+"_res")
@@ -34 +34 @@
         Variable/G USANS_slope_numpnts = 10
         
+        Print "---- Calculating Weighting Matrix for USANS Data ----"
+        
         EnterSlope(basestr)
-
+        
+        //Deal with broken Pauseforuser
+        SetDataFolder $("root:"+basestr)
+        
         print USANS_m
 
@@ -61 +66 @@
         String baseStr
         
+        NVAR USANS_m,USANS_slope_numpnts
+        
 //      Variable slope=-4
 
@@ -71 +78 @@
 //      return slope
 
-        NewPanel /K=1/W=(1141,543,1548,910)/N=USANS_Slope as "USANS Slope Extrapolation"
+        NewPanel /W=(600,300,1000,700)/N=USANS_Slope as "USANS Slope Extrapolation"
         SetDrawLayer UserBack
         Button button_OK,pos={294,317},size={100,20},title="Accept Slope",font="Geneva"
         Button button_OK,proc=USANS_Slope_ButtonProc
-        Button button_Cancel,pos={233,317},size={55,20},title="Cancel",font="Geneva"
-        Button button_Cancel,proc=USANS_Slope_ButtonProc
+        //Button button_Cancel,pos={233,317},size={55,20},title="Cancel",font="Geneva"
+        //Button button_Cancel,proc=USANS_Slope_ButtonProc
         SetVariable setvar_numpnts, pos={43,350}, size={150,19}, title="# Points",fSize=13
         SetVariable setvar_numpnts, value= USANS_slope_numpnts 
@@ -86 +93 @@
         SetActiveSubwindow ##
         
-        Print GetDataFolder(1)
+        //Print TraceNameList("USANS_Slope#SlopePlot",";",1)
+        
+        USANS_CalculateSlope(basestr,USANS_slope_numpnts)
+        
         PauseForUser  USANS_Slope       
-        Print GetDataFolder(1)
 //      return slope
         
 End
+
+Function USANS_CalculateSlope(basestr, calcpoints)
+        String basestr
+        Variable calcpoints
+        
+        Wave iw = $(basestr+"_i")
+        Wave qw = $(basestr+"_i")
+        Wave sw = $(basestr+"_i")
+
+        Variable num_extr=25
+        NVAR nend = USANS_slope_numpnts
+        // Make extra waves for extrapolation 
+        // Taken from DSM_SetExtrWaves
+                
+        Make/O/D/N=(num_extr) extr_hqq,extr_hqi
+        extr_hqi=1              //default values
+
+        //set the q-range
+        Variable qmax,num
+
+        num=numpnts(qw)
+        qmax=6*qw[num-1]
+        
+        extr_hqq = qw[num-1] + x * (qmax-qw[num-1])/num_extr
+                
+        
+        // Modifed from DSM_DoExtraploate in LakeDesmearing_JB.ipf 
+        // which is part of the USANS Reduction macros
+        
+                
+        //      Wave extr_lqi=extr_lqi
+        //      Wave extr_lqq=extr_lqq
+        //      Wave extr_hqi=extr_hqi
+        //      Wave extr_hqq=extr_hqq
+                Variable/G V_FitMaxIters=300
+                Variable/G V_fitOptions=4               //suppress the iteration window
+                Variable retval
+                num=numpnts(iw)
+        
+                Make/O/D P_coef={0,1,-4}                        //input
+        //      Make/O/T Constr={"K2<0","K2 > -20"}
+                //(set background to zero and hold fixed)
+                CurveFit/H="100" Power kwCWave=P_coef  iw[(num-1-nend),(num-1)] /X=qw /W=sw /D 
+                extr_hqi=P_coef[0]+P_coef[1]*extr_hqq^P_coef[2]
+        
+                Printf "Smeared Power law exponent = %g\r",P_coef[2]
+                Printf "**** For Desmearing, use a Power law exponent of %5.1f\r",P_coef[2]-1
+        
+                retVal = P_coef[2]-1                    
+                return(retVal)
+                
+End
+
 
 Function USANS_Slope_ButtonProc(ctrlName) : ButtonControl
@@ -115 +177 @@
 
         Variable i,j
-        SVAR USANS_basestr
+        SVAR USANS_basestr=root:myGlobals:USANS_basestr 
+        SetDataFolder $("root:"+USANS_basestr)
+
+        
+        
         NVAR dQv = USANS_dQv
         
@@ -138 +204 @@
         Variable i,j
         
-        SVAR USANS_basestr
+        SVAR USANS_basestr=root:myGlobals:USANS_basestr 
+        SetDataFolder $("root:"+USANS_basestr)
+
+
+
         NVAR dQv = USANS_dQv
         
@@ -161 +231 @@
         Variable i
 
-        SVAR USANS_basestr
+        SVAR USANS_basestr=root:myGlobals:USANS_basestr 
+        SetDataFolder $("root:"+USANS_basestr)
+
+
         NVAR m = USANS_m
         NVAR N = USANS_N
@@ -190 +263 @@
         Variable i
         
-        SVAR USANS_basestr
+        SVAR USANS_basestr=root:myGlobals:USANS_basestr 
+        SetDataFolder $("root:"+USANS_basestr)
+
         NVAR m = USANS_m
         NVAR qi = USANS_intQpt
@@ -202 +277 @@
 
 End
-
-Window USANS_Slope() : Panel
-        PauseUpdate; Silent 1           // building window...
-        NewPanel /K=1/W=(1141,543,1548,910) as "USANS Slope Extrapolation"
-        SetDrawLayer UserBack
-        Button button_OK,pos={294,317},size={100,20},title="Accept Slope",font="Geneva"
-        Button button_Cancel,pos={233,317},size={55,20},title="Cancel",font="Geneva"
-        SetVariable setvar_Slope,pos={43,321},size={100,19},title="Slope",fSize=13
-        SetVariable setvar_Slope,limits={-inf,inf,0},value= USANS_m
-        Display/W=(9,6,402,305)/HOST=# 
-        RenameWindow #,G0
-        SetActiveSubwindow ##
-EndMacro