Ignore:
Timestamp:
Jul 21, 2017 1:37:44 PM (5 years ago)
Author:
srkline
Message:

added Annular averaging routine. still need writer, and link with protocols.

replaced tic/toc with v_tic/v_toc to avoid missing functions

File:
1 edited

Legend:

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

    r1044 r1056  
    230230        Variable xzsum,yzsum,zsum,xctr,yctr 
    231231        Variable left,right,bottom,top,ii,jj,counts 
    232          
     232        Variable x_mm_sum,y_mm_sum,x_mm,y_mm 
    233233 
    234234         
     
    259259                Wave data = V_getDetectorDataW(gCurDispType,detStr) 
    260260                Wave data_err = V_getDetectorDataErrW(gCurDispType,detStr) 
     261                 
     262                // get the real-space information 
     263                String destPath = "root:Packages:NIST:VSANS:"+gCurDispType 
     264                Wave data_realDistX = $(destPath + ":entry:instrument:detector_"+detStr+":data_realDistX") 
     265                Wave data_realDistY = $(destPath + ":entry:instrument:detector_"+detStr+":data_realDistY") 
    261266         
    262267                xzsum = 0 
    263268                yzsum = 0 
    264269                zsum = 0 
     270                x_mm_sum = 0 
     271                y_mm_sum = 0 
     272                 
    265273                // count over rectangular selection, doing each row, L-R, bottom to top 
    266274                ii = bottom -1 
     
    274282                                yzsum += ii*counts 
    275283                                zsum += counts 
     284                                 
     285                                x_mm_sum += data_realDistX[jj][ii]*counts 
     286                                y_mm_sum += data_realDistY[jj][ii]*counts 
    276287                        while(jj<right) 
    277288                while(ii<top) 
     
    280291                yctr = yzsum/zsum 
    281292                 
     293                x_mm = x_mm_sum/zsum 
     294                y_mm = y_mm_sum/zsum 
    282295                // add 1 to each to get to detector coordinates (1,128) 
    283296                // rather than the data array which is [0,127] 
     
    287300                Print "X-center (in array coordinates 0->n-1 ) = ",xctr 
    288301                Print "Y-center (in array coordinates 0->n-1 ) = ",yctr 
     302                 
     303                Print "X-center (mm) = ",x_mm 
     304                Print "Y-center (mm) = ",y_mm 
    289305        endif 
    290306         
Note: See TracChangeset for help on using the changeset viewer.