Ignore:
Timestamp:
Apr 17, 2013 11:36:13 AM (10 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.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • sans/Dev/trunk/NCNR_User_Procedures/Reduction/USANS/BT5_Loader.ipf

    r807 r907  
    215215// now normalizes to the monitor counts 
    216216// 
     217// 17 APR 2013 SRK - coverted the 2 degree "cutoff" to use q-values instead. Then the  
     218//  location is generic and can be used for either NCNR or KUSANS 
     219// 
    217220Function FindTWideCts(type) 
    218221        String type 
     
    226229        Wave TransCts = $(USANSFolder+":"+type+":TransCts") 
    227230        Wave monCts = $(USANSFolder+":"+type+":MonCts") 
     231        Wave Qvals = $(USANSFolder+":"+type+":Qvals") 
    228232         
    229233        num=numpnts(TransCts) 
    230234 
    231  
    232         FindLevel/Q/P angle,2           //use angles greater than 2 deg 
     235//      FindLevel/Q/P angle,2           //use angles greater than 2 deg 
     236//      Print "Using angle, pt = ",V_levelX 
     237         
     238        FindLevel/Q/P Qvals,1e-4                //use angles greater than 2 deg = 2*5.55e-5 = 1e-4 (1/A) 
     239        Print "Using Qval, pt = ",V_levelX 
     240 
    233241        levNotFound=V_Flag              //V_Flag==1 if no pk found 
    234242         
    235243        if(levNotFound) 
    236244                //post a warning, and use just the last 4 points... 
    237                 DoAlert 0,"You don't have "+type+" data past 2 degrees - so Twide may not be reliable" 
     245//              DoAlert 0,"You don't have "+type+" data past 2 degrees - so Twide may not be reliable" 
     246                DoAlert 0,"You don't have "+type+" data past 1e-4 A-1 (~2 degrees) - so Twide may not be reliable" 
    238247                levPt = num-4 
    239248        else 
Note: See TracChangeset for help on using the changeset viewer.