Ignore:
Timestamp:
May 7, 2018 1:14:22 PM (5 years ago)
Author:
srkline
Message:

added Beaucage model to White Beam so that I could model the CFG sample.

cleaned up a few routines to get rid of NaN values in the output I(q) data sets if they were not masked out completely

File:
1 edited

Legend:

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

    r1097 r1099  
    15101510end 
    15111511 
    1512  
    1513  
    1514  
    1515  
     1512// from WM procedure, extended to three waves for I(q) data sets 
     1513// 
     1514// RemoveNaNsQIS(theQWave, theIWave, theSWave) 
     1515//      Removes all points in an XYZ triple if any of the three waves has the value NaN. 
     1516//      A NaN represents a blank or missing value. 
     1517//      Returns the number of points removed. 
     1518Function V_RemoveNaNsQIS(theXWave, theYWave, theZWave) 
     1519        Wave theXWave 
     1520        Wave theYWave 
     1521        Wave theZWave 
     1522 
     1523        Variable p, numPoints, numNaNs 
     1524        Variable xval, yval, zval 
     1525         
     1526        numNaNs = 0 
     1527        p = 0                                                                                   // the loop index 
     1528        numPoints = numpnts(theXWave)                   // number of times to loop 
     1529 
     1530        do 
     1531                xval = theXWave[p] 
     1532                yval = theYWave[p] 
     1533                zval = theZWave[p] 
     1534                if ((numtype(xval)==2) %| (numtype(yval)==2) %| (numtype(zval)==2) )            // either is NaN? 
     1535                        numNaNs += 1 
     1536                else                                                                            // if not an outlier 
     1537                        theXWave[p - numNaNs] = xval            // copy to input wave 
     1538                        theYWave[p - numNaNs] = yval            // copy to input wave 
     1539                        theZWave[p - numNaNs] = zval            // copy to input wave 
     1540                endif 
     1541                p += 1 
     1542        while (p < numPoints) 
     1543         
     1544        // Truncate the wave 
     1545        DeletePoints numPoints-numNaNs, numNaNs, theXWave, theYWave, theZWave 
     1546         
     1547        return(numNaNs) 
     1548End 
     1549 
     1550 
Note: See TracChangeset for help on using the changeset viewer.