Ignore:
Timestamp:
Apr 24, 2017 2:04:29 PM (6 years ago)
Author:
srkline
Message:

Adding utilities to calculate transmission. Added V_Transmission.ipf file. Getting the functionality to work, user interface or incorporation into the data reduction protocol is still to be decided.

File:
1 edited

Legend:

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

    r1026 r1033  
    313313        fname = S_path + fname 
    314314 
    315         Variable nRows = 13 
     315        Variable nRows = 14 
    316316        Redimension/N=(nRows,3) ListWave 
    317317        Redimension/N=(nRows,3) selWave 
     
    349349        listWave[8][1] = "group_id (sample)" 
    350350        listWave[8][2] = num2str(V_getSample_group_ID(fname)) 
    351  
    352         listWave[9][1] = "box_count" 
    353         listWave[9][2] = num2str(V_getBoxCounts(fname)) 
    354          
    355         listWave[10][1] = "box_count_error" 
    356         listWave[10][2] = num2str(V_getBoxCountsError(fname)) 
    357          
    358         listWave[11][1] = "whole_trans" 
    359         listWave[11][2] = num2str(V_getSampleTransWholeDetector(fname)) 
    360          
    361         listWave[12][1] = "whole_trans_error" 
    362         listWave[12][2] = num2str(V_getSampleTransWholeDetErr(fname)) 
    363          
    364          
     351         
     352        listWave[9][1] = "Box Coordinates" 
     353        WAVE boxCoord = V_getBoxCoordinates(fname) 
     354        listWave[9][2] = V_NumWave2List(boxCoord,";") 
     355 
     356        listWave[10][1] = "box_count" 
     357        listWave[10][2] = num2str(V_getBoxCounts(fname)) 
     358         
     359        listWave[11][1] = "box_count_error" 
     360        listWave[11][2] = num2str(V_getBoxCountsError(fname)) 
     361         
     362        listWave[12][1] = "whole_trans" 
     363        listWave[12][2] = num2str(V_getSampleTransWholeDetector(fname)) 
     364         
     365        listWave[13][1] = "whole_trans_error" 
     366        listWave[13][2] = num2str(V_getSampleTransWholeDetErr(fname))    
    365367                 
    366368         
     
    970972        switch(V_Value) // numeric switch 
    971973                case 0: // execute if case matches expression 
    972                         V_WriteHeaderForPatch_0(tempName) 
     974                        V_WriteHeaderForPatch_0(tempName)               //control 
    973975                        break           // exit from switch 
    974976                case 1:  
    975                         V_WriteHeaderForPatch_1(tempName) 
     977                        V_WriteHeaderForPatch_1(tempName)               //reduction 
    976978                        break 
    977979                case 2:  
    978                         V_WriteHeaderForPatch_2(tempName) 
     980                        V_WriteHeaderForPatch_2(tempName)               // sample 
    979981                        break 
    980982                case 3:  
    981                         V_WriteHeaderForPatch_3(tempName) 
     983                        V_WriteHeaderForPatch_3(tempName)               // instrument 
    982984                        break 
    983985                case 4:  
    984                         V_WriteHeaderForPatch_4(tempName) 
     986                        V_WriteHeaderForPatch_4(tempName)               //detectors 
    985987                        break 
    986988                case 5:  
    987                         V_WriteHeaderForPatch_5(tempName) 
     989                        V_WriteHeaderForPatch_5(tempName)               // polSANS 
    988990                        break 
    989991                default:                        // optional default expression executed 
     
    10991101                err = V_writeSample_GroupID(fname,val) 
    11001102        endif    
    1101          
    1102         if ((selWave[9][0] & 2^4) != 0)         //"box_count" 
    1103                 val = str2num(listWave[9][2]) 
     1103 
     1104 
     1105 
     1106        if ((selWave[9][0] & 2^4) != 0)         //"box coordinates" 
     1107                str = listWave[9][2] 
     1108                err = V_writeBoxCoordinates(fname,V_List2NumWave(str,";","inW")) 
     1109        endif    
     1110         
     1111 
     1112         
     1113        if ((selWave[10][0] & 2^4) != 0)                //"box_count" 
     1114                val = str2num(listWave[10][2]) 
    11041115                err = V_writeBoxCounts(fname,val) 
    11051116        endif    
    11061117         
    1107         if ((selWave[10][0] & 2^4) != 0)                //"box_count_error" 
    1108                 val = str2num(listWave[10][2]) 
     1118        if ((selWave[11][0] & 2^4) != 0)                //"box_count_error" 
     1119                val = str2num(listWave[11][2]) 
    11091120                err = V_writeBoxCountsError(fname,val) 
    11101121        endif    
    11111122         
    1112         if ((selWave[11][0] & 2^4) != 0)                //"whole_trans" 
    1113                 val = str2num(listWave[11][2]) 
     1123        if ((selWave[12][0] & 2^4) != 0)                //"whole_trans" 
     1124                val = str2num(listWave[12][2]) 
    11141125                err = V_writeSampleTransWholeDetector(fname,val) 
    11151126        endif    
    11161127         
    1117         if ((selWave[12][0] & 2^4) != 0)                //"whole_trans_error" 
    1118                 val = str2num(listWave[12][2]) 
     1128        if ((selWave[13][0] & 2^4) != 0)                //"whole_trans_error" 
     1129                val = str2num(listWave[13][2]) 
    11191130                err = V_writeSampleTransWholeDetErr(fname,val) 
    11201131        endif    
Note: See TracChangeset for help on using the changeset viewer.