Ignore:
Timestamp:
Apr 16, 2018 10:50:33 AM (4 years ago)
Author:
srkline
Message:

Added two model functions for white beam smearing.

Many other small changes for processing of the back detector, shuffling of VSANS menu items, and consistent naming of V_ procedures.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/V_IQ_Annular.ipf

    r1073 r1095  
    1818// -- integrate this with a more general "average panel" 
    1919// -- draw the q-center and width on the image (as a contour? - if no, use a draw layer...) 
    20 //   Can I "flag" the pixels that contribute to the annualr average, and overlay them like a  
     20// -- Can I "flag" the pixels that contribute to the annualr average, and overlay them like a  
    2121//    Mask(translucent?), like a "thresholding" operation, but act on the Q-value, not the Z- value. 
     22// 
     23// TODO!!! 
     24// x- none of these procedures are aware of the BACK detector 
     25// 
    2226 
    2327// 
     
    4650        String detStr 
    4751         
    48         if(cmpstr(detGroup,"F") == 0) 
    49                 detStr = "FLRTB" 
    50         else 
    51                 detStr = "MLRTB" 
    52         endif 
     52        strswitch(detGroup) 
     53                case "F": 
     54                        detStr = "FLRTB" 
     55                        break 
     56                case "M": 
     57                        detStr = "MLRTB" 
     58                        break 
     59                case "B": 
     60                        detStr = "B" 
     61                        break 
     62                default: 
     63                        DoAlert 0,"No detGroup match in V_QBinAllPanels_Annular" 
     64                        return(0) 
     65        endswitch 
     66 
    5367         
    5468// right now, use all of the detectors. There is a lot of waste in this and it could be  
     
    7892                RemoveFromGraph/Z iPhiBin_qxqy_FLRTB 
    7993                RemoveFromGraph/Z iPhiBin_qxqy_MLRTB 
     94                RemoveFromGraph/Z iPhiBin_qxqy_B 
    8095        endif 
    8196 
     
    89104                Label bottom "Phi (degrees)" 
    90105                Legend 
    91         else 
     106        endif 
     107        if(cmpstr(detGroup,"M") == 0) 
    92108                AppendToGraph iPhiBin_qxqy_MLRTB vs phiBin_qxqy_MLRTB 
    93109                ModifyGraph mode=4 
     
    98114                Legend 
    99115        endif 
    100  
     116        if(cmpstr(detGroup,"B") == 0) 
     117                AppendToGraph iPhiBin_qxqy_B vs phiBin_qxqy_B 
     118                ModifyGraph mode=4 
     119                ModifyGraph marker=19 
     120                ErrorBars iPhiBin_qxqy_B Y,wave=(ePhiBin_qxqy_B,ePhiBin_qxqy_B) 
     121                Label left "Counts" 
     122                Label bottom "Phi (degrees)" 
     123                Legend 
     124        endif 
    101125         
    102126        SetDataFolder fldrSav0 
     
    247271                        break                                                                    
    248272                                         
     273                case "B": 
     274                        if(isVCALC) 
     275                                WAVE inten = $(folderPath+instPath+"B"+":det_"+"B") 
     276                                WAVE/Z iErr = $("iErr_"+"B")                    // 2D errors -- may not exist, especially for simulation                 
     277                        else 
     278                                Wave inten = V_getDetectorDataW(folderStr,"B") 
     279                                Wave iErr = V_getDetectorDataErrW(folderStr,"B") 
     280 
     281                                Wave/Z mask = $("root:Packages:NIST:VSANS:MSK:entry:instrument:detector_"+"B"+":data") 
     282                                 
     283                                if(WaveExists(mask) == 1) 
     284                                        maskMissing = 0 
     285                                endif 
     286                        endif    
     287//                      NVAR delQ = $(folderPath+instPath+"ML"+":gDelQ_ML") 
     288                         
     289                        Wave qTotal = $(folderPath+instPath+"B"+":qTot_"+"B")                   // 2D q-values   
     290         
     291                        Wave qx = $(folderPath+instPath+"B"+":qx_"+"B")                 // 2D qx-values  
     292 
     293                        Wave qy = $(folderPath+instPath+"B"+":qy_"+"B")                 // 2D qy-values  
     294                                         
     295                        nSets = 1 
     296                        break                                                            
     297 
     298                                         
    249299                default: 
    250300                        nSets = 0                                                       // optional default expression executed 
     
    654704                Wave/Z iw = iPhiBin_qxqy_FLRTB 
    655705                Wave/Z sw = ePhiBin_qxqy_FLRTB 
    656         else 
     706        endif 
     707        if(cmpstr(detGroup,"M") == 0) 
    657708                Wave/Z pw = phiBin_qxqy_MLRTB 
    658709                Wave/Z iw = iPhiBin_qxqy_MLRTB 
    659710                Wave/Z sw = ePhiBin_qxqy_MLRTB 
    660711        endif 
    661  
     712        if(cmpstr(detGroup,"B") == 0) 
     713                Wave/Z pw = phiBin_qxqy_B 
     714                Wave/Z iw = iPhiBin_qxqy_B 
     715                Wave/Z sw = ePhiBin_qxqy_B 
     716        endif 
    662717         
    663718        String dataSetFolderParent,basestr 
Note: See TracChangeset for help on using the changeset viewer.