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

Location:
sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS
Files:
1 added
9 edited

Legend:

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

    r1052 r1056  
    101101// 
    102102#include "V_Instrument_Resolution" 
     103#include "V_IQ_Annular" 
    103104 
    104105 
     106 
  • sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/V_EventModeProcessing.ipf

    r1050 r1056  
    570570        SetDataFolder root:Packages:NIST:VSANS:Event: 
    571571 
    572 s_tic() 
     572v_tic() 
    573573        if(WaveExists($"root:Packages:NIST:VSANS:Event:OscSortIndex") == 0 ) 
    574574                Duplicate/O rescaledTime OscSortIndex 
     
    581581         
    582582printf "sort time = " 
    583 s_toc() 
     583v_toc() 
    584584 
    585585        Wave index = root:Packages:NIST:VSANS:Event:SavedIndex          //this is the histogram index 
    586586 
    587 s_tic() 
     587v_tic() 
    588588        for(ii=0;ii<nslices;ii+=1) 
    589589                if(ii==0) 
     
    620620        endfor 
    621621printf "histogram time = " 
    622 s_toc() 
     622v_toc() 
    623623 
    624624        Duplicate/O slicedData,root:Packages:NIST:VSANS:Event:dispsliceData,root:Packages:NIST:VSANS:Event:logSlicedData 
     
    986986                Redimension/D xLoc,yLoc,timePt   
    987987                 
    988 //s_tic() 
     988//v_tic() 
    989989//      V_SortAndSplitEvents() 
    990990// 
    991991//Printf "File sort and split time (s) = " 
    992 //s_toc() 
     992//v_toc() 
    993993 
    994994// 
     
    12331233        Open/R refnum as filepathstr 
    12341234         
    1235 s_tic() 
     1235v_tic() 
    12361236 
    12371237        FBinRead refnum, gVSANSStr 
     
    12621262 
    12631263Printf "Time to read file (s) = " 
    1264 s_toc()  
     1264v_toc()  
    12651265 
    12661266 
     
    12741274 
    12751275 
    1276 s_tic() 
     1276v_tic() 
    12771277        WAVE V_Events = V_Events 
    12781278        uint64 val,b1,b2,btime 
     
    13091309 
    13101310Printf "File decode time (s) = " 
    1311 s_toc() 
     1311v_toc() 
    13121312 
    13131313//      KillWaves/Z timePt,xLoc,yLoc 
     
    14191419//      NVAR removeBadEvents = root:Packages:NIST:VSANS:Event:gRemoveBadEvents 
    14201420// 
    1421 //s_tic() 
     1421//v_tic() 
    14221422// 
    14231423////    Wave/Z wave0=wave0 
     
    15001500//      dispStr += tmpStr 
    15011501// 
    1502 //      s_toc() 
     1502//      v_toc() 
    15031503//       
    15041504//      Print "Events removed (XOP) = ",numRemoved 
  • sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/V_HDF5_RW_Utils.ipf

    r1055 r1056  
    435435        endif 
    436436         
    437 //s_tic()               //fast  
     437//v_tic()               //fast  
    438438         
    439439        SVAR tmpStr=root:file_name 
     
    536536        HDF5CloseFile fileID 
    537537         
    538 //s_toc() 
     538//v_toc() 
    539539 
    540540        // save a global string with the file name to be picked up for the status on the display 
  • 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         
  • sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/V_Menu.ipf

    r1051 r1056  
    2525                "Patch Detector Deadtime",V_PatchDetectorDeadtimePanel() 
    2626                "Patch Detector Calibration",V_PatchDetectorCalibrationPanel() 
     27                "-" 
     28                "Annular Binning",Annular_Binning() 
    2729        End 
    2830        Submenu "Work Files" 
  • sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/V_Protocol_Reduction.ipf

    r1055 r1056  
    675675        String match="BLOCKED BEAM" 
    676676 
    677 //s_tic() 
     677//v_tic() 
    678678// get the list from the file catalog (=0.0007s) 
    679679//       
     
    690690        List = SortList(List,";",0) 
    691691        Printf "BGD files = %s\r",list   
    692 //S_toc()        
    693          
    694          
    695 ////S_tic()                     // from grep = 3.3s 
     692//v_toc()        
     693         
     694         
     695////v_tic()                     // from grep = 3.3s 
    696696//      newList = V_GetRawDataFileList() 
    697697//      num=ItemsInList(newList) 
     
    719719//      List = SortList(List,";",0) 
    720720//      Printf "BGD files = %s\r",list 
    721 ////S_toc() 
     721////v_toc() 
    722722         
    723723        return(list) 
     
    26762676                //              TODO 
    26772677                // -- get all of the proper values for the calculation 
     2678                // -- currently the attenuation is incorrect 
     2679                //   such that kappa_err = 1*kappa 
    26782680                // -- verify the calculation (no solid angle needed) 
     2681                 
     2682                DoAlert 0,"This calculation is not reliable - attenuation is not calibrated" 
    26792683                 
    26802684                // get the attenuation factor for the empty beam 
  • sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/V_Utilities_General.ipf

    r1051 r1056  
    12001200        return(newList) 
    12011201End 
     1202 
     1203 
     1204///// 
     1205// @ IgorExchange 
     1206//TicToc 
     1207//Posted April 16th, 2009 by bgallarda 
     1208//      ¥       in Programming 6.10.x 
     1209 
     1210////duplicated here -- for VSANS use 
     1211function v_tic() 
     1212        variable/G tictoc = startMSTimer 
     1213end 
     1214  
     1215function v_toc() 
     1216        NVAR/Z tictoc 
     1217        variable ttTime = stopMSTimer(tictoc) 
     1218        printf "%g seconds\r", (ttTime/1e6) 
     1219        killvariables/Z tictoc 
     1220end 
     1221 
     1222 
  • sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/V_VSANS_Event_Testing.ipf

    r1048 r1056  
    253253        Wave w 
    254254 
    255 s_tic() 
     255v_tic() 
    256256//      WAVE w = eventWave 
    257257        uint64 val,b1,b2,btime 
     
    290290//      endfor 
    291291 
    292 s_toc() 
     292v_toc() 
    293293                 
    294294        return(0) 
     
    383383        filename = S_fileName 
    384384 
    385 s_tic() 
     385v_tic() 
    386386 
    387387        FBinRead refnum, vsansStr 
     
    411411        Duplicate/O $(StringFromList(0,S_waveNames)) V_Events 
    412412        KillWaves/Z $(StringFromList(0,S_waveNames)) 
    413 s_toc()  
     413v_toc()  
    414414         
    415415        Print vsansStr 
  • sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/V_VSANS_Preferences.ipf

    r1051 r1056  
    6262        Variable/G root:Packages:NIST:VSANS:Globals:gBinWidth=val 
    6363         
    64         val = NumVarOrDefault("root:Packages:NIST:VSANS:Globals:gNPhiSteps", 72 ) 
     64        val = NumVarOrDefault("root:Packages:NIST:VSANS:Globals:gNPhiSteps", 90 ) 
    6565        Variable/G root:Packages:NIST:VSANS:Globals:gNPhiSteps=val 
    6666         
Note: See TracChangeset for help on using the changeset viewer.