Ignore:
Timestamp:
Apr 8, 2011 12:58:26 PM (12 years ago)
Author:
srkline
Message:

Removed Smear_Model_N_AAO() and reverted back to the non-AAO version. The AAO version is not cursor-aware, and then smeares incorrectly when using cursors. The AAO version can be corrected in the future if necessary. Using the Debye Sphere functions as fit functions need AAO smearing to be useable at all.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • sans/Dev/trunk/NCNR_User_Procedures/Reduction/SANS/WriteQIS.ipf

    r794 r797  
    817817//********************* 
    818818 
    819 //      // generate my own error wave for I(qx,qy) 
     819        // generate my own error wave for I(qx,qy) 
     820        //      //      sw = 0.05*sw            // uniform 5% error? tends to favor the low intensity too strongly 
     821        // after teh SQRT(),  get rid of the "bad" errors by replacing the NaN, Inf, and zero with V_avg 
     822        // THIS IS EXTREMEMLY IMPORTANT - if this is not done, there are some "bad" values in the  
     823        // error wave (things that are not numbers) - and this wrecks the smeared model fitting. 
     824        // It appears to have no effect on fitting with the unsmeared model. 
    820825//      Duplicate/O z_val sw 
    821 //      sw = sqrt(z_val)                //assumes Poisson statistics for each cell (counter) 
    822 //      //      sw = 0.05*sw            // uniform 5% error? tends to favor the low intensity too strongly 
    823 //      // get rid of the "bad" errorsby replacing the NaN, Inf, and zero with V_avg 
    824 //      // THIS IS EXTREMEMLY IMPORTANT - if this is not done, there are some "bad" values in the  
    825 //      // error wave (things that are not numbers) - and this wrecks the smeared model fitting. 
    826 //      // It appears to have no effect on the unsmeared model. 
     826//      sw = sqrt(z_val)                //assumes Poisson statistics for each cell (counter) - which is wrong, since this is corrected data 
     827// 
    827828//      WaveStats/Q sw 
    828829//      sw = numtype(sw[p]) == 0 ? sw[p] : V_avg 
    829830//      sw = sw[p] != 0 ? sw[p] : V_avg 
    830          
    831         // now use the properly propagated 2D error 
     831 
     832 
     833//      // April 2011 -- now use the properly propagated 2D error 
    832834        Duplicate/O data_err sw 
    833835        Redimension/N=(pixelsX*pixelsY) sw 
Note: See TracChangeset for help on using the changeset viewer.