Changeset 797 for sans/Dev/trunk/NCNR_User_Procedures/Reduction/SANS
- Timestamp:
- Apr 8, 2011 12:58:26 PM (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sans/Dev/trunk/NCNR_User_Procedures/Reduction/SANS/WriteQIS.ipf
r794 r797 817 817 //********************* 818 818 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. 820 825 // 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 // 827 828 // WaveStats/Q sw 828 829 // sw = numtype(sw[p]) == 0 ? sw[p] : V_avg 829 830 // 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 832 834 Duplicate/O data_err sw 833 835 Redimension/N=(pixelsX*pixelsY) sw
Note: See TracChangeset
for help on using the changeset viewer.