Ignore:
Timestamp:
Feb 5, 2019 3:13:21 PM (4 years ago)
Author:
srkline
Message:

adding procedures for calculating sector averages of the data

still need to add overlay of sector angles.

File:
1 edited

Legend:

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

    r1111 r1121  
    496496        Wave qb = tmp_qb 
    497497        Wave fs = tmp_fs 
     498         
     499        // zero point waves can be generated from sector averages, exit now 
     500        if(numpnts(qw) == 0) 
     501                SetDataFolder root: 
     502                return(1) 
     503        endif 
    498504         
    499505        Variable q1,q2,tol,ii 
     
    748754        SetDataFolder $("root:Packages:NIST:VSANS:"+folderStr) 
    749755 
    750         Printf "%d points removed from beginning, %d points from the end  of %s \r",nbeg,nend,detStr 
    751756 
    752757// TODO  
     
    768773                Wave/Z qBar = $("qBar_qxqy_"+detStr) 
    769774                Wave/Z fSubS = $("fSubS_qxqy_"+detStr) 
    770                          
     775 
     776// check for the existence of every wave + non-zero length 
     777// zero length/null waves can be generated with sector averages 
     778// exit if something is wrong 
     779                if(waveExists(qw) == 0 || numpnts(qw) == 0) 
     780                        SetDataFolder root: 
     781                        return(1) 
     782                endif 
     783                if(waveExists(iw) == 0 || numpnts(iw) == 0) 
     784                        SetDataFolder root: 
     785                        return(1) 
     786                endif 
     787                if(waveExists(ew) == 0 || numpnts(ew) == 0) 
     788                        SetDataFolder root: 
     789                        return(1) 
     790                endif 
     791                if(waveExists(sigQ) == 0 || numpnts(sigQ) == 0) 
     792                        SetDataFolder root: 
     793                        return(1) 
     794                endif 
     795                if(waveExists(qBar) == 0 || numpnts(qBar) == 0) 
     796                        SetDataFolder root: 
     797                        return(1) 
     798                endif 
     799                if(waveExists(fSubS) == 0 || numpnts(fSubS) == 0) 
     800                        SetDataFolder root: 
     801                        return(1) 
     802                endif 
     803 
     804 
     805        // waves all exist               
    771806                DeletePoints 0,nBeg, qw,iw,ew,sigQ,qBar,fSubS 
    772807 
     
    775810                DeletePoints npt-nEnd,nEnd, qw,iw,ew,sigQ,qBar,fSubS 
    776811         
     812                Printf "%d points removed from beginning, %d points from the end  of %s \r",nbeg,nend,detStr 
     813                SetDataFolder root: 
    777814        return(0) 
    778815End 
Note: See TracChangeset for help on using the changeset viewer.