Changeset 907

Timestamp:

Apr 17, 2013 11:36:13 AM (9 years ago)

Author:

srkline

Message:

Changed the calculation of the wide angle transmission to look for q-values greater than 1e-4 (1/A) rather than > 2 degrees so that it would be more generic, and can be used for KUSANS. The 2 degree point is where the transmission counts levels off (empirically) and I simply converted this to Q.

Location:

sans/Dev/trunk/NCNR_User_Procedures/Reduction/USANS

Files:

• BT5_Loader.ipf (modified) (2 diffs)
• COR_Graph.ipf (modified) (3 diffs)

sans/Dev/trunk/NCNR_User_Procedures/Reduction/USANS/BT5_Loader.ipf

@@ -215 +215 @@
 // now normalizes to the monitor counts
 //
+// 17 APR 2013 SRK - coverted the 2 degree "cutoff" to use q-values instead. Then the 
+//  location is generic and can be used for either NCNR or KUSANS
+//
 Function FindTWideCts(type)
         String type
@@ -226 +229 @@
         Wave TransCts = $(USANSFolder+":"+type+":TransCts")
         Wave monCts = $(USANSFolder+":"+type+":MonCts")
+        Wave Qvals = $(USANSFolder+":"+type+":Qvals")
         
         num=numpnts(TransCts)
 
-
-        FindLevel/Q/P angle,2           //use angles greater than 2 deg
+//      FindLevel/Q/P angle,2           //use angles greater than 2 deg
+//      Print "Using angle, pt = ",V_levelX
+        
+        FindLevel/Q/P Qvals,1e-4                //use angles greater than 2 deg = 2*5.55e-5 = 1e-4 (1/A)
+        Print "Using Qval, pt = ",V_levelX
+
         levNotFound=V_Flag              //V_Flag==1 if no pk found
         
         if(levNotFound)
                 //post a warning, and use just the last 4 points...
-                DoAlert 0,"You don't have "+type+" data past 2 degrees - so Twide may not be reliable"
+//              DoAlert 0,"You don't have "+type+" data past 2 degrees - so Twide may not be reliable"
+                DoAlert 0,"You don't have "+type+" data past 1e-4 A-1 (~2 degrees) - so Twide may not be reliable"
                 levPt = num-4
         else

sans/Dev/trunk/NCNR_User_Procedures/Reduction/USANS/COR_Graph.ipf

@@ -111 +111 @@
         if(strlen(list)!=0)
                 //Print "SAM already on graph"
+                SetDataFolder root:
                 return(0)
         Endif
@@ -136 +137 @@
         if(strlen(list)!=0)
         //      Print "EMP already on graph"
+                SetDataFolder root:
                 return(0)
         Endif
@@ -162 +164 @@
         if(strlen(list)!=0)
         //      Print "COR already on graph"
+                SetDataFolder root:
                 return(0)
         Endif