- Timestamp:
- Nov 25, 2009 2:40:15 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sans/Dev/trunk/NCNR_User_Procedures/Reduction/SANS/WorkFileUtils.ipf
r588 r605 484 484 //rad is the distance of pixel ij from the sample 485 485 //domega is the ratio of the solid angle of pixel ij versus center pixel 486 // product xy = 1 for a detector with a linear spatial response (modern Ordela) 487 // solid angle calculated, dW^3 >=1, so multiply data to raise measured values to correct values. 486 488 rad = sqrt(dtdis2 + xd^2 + yd^2) 487 489 domega = rad/dtdist 488 490 ratio = domega^3 489 // solidAngle[ii][jj] = ratio //testing only490 491 xy = xx[ii]*yy[jj] 491 492 data[ii][jj] *= xy*ratio 493 // solidAngle[ii][jj] = xy*ratio //testing only 494 492 495 493 // multiplicative correctionfor detector efficiency JBG memo det_eff_cor2.doc 3/20/07496 // correction factor for detector efficiency JBG memo det_eff_cor2.doc 3/20/07 494 497 // correction inserted 11/2007 SRK 495 // large angle detector efficiency correction is >= 1 and will "bump up"the highest angles498 // large angle detector efficiency is >= 1 and will "bump up" the measured value at the highest angles 496 499 // so divide here to get the correct answer (5/22/08 SRK) 497 500 if(doEfficiency) … … 501 504 #else 502 505 data[ii][jj] /= DetEffCorr(lambda,dtdist,xd,yd) 503 // solidAngle[ii][jj]= DetEffCorr(lambda,dtdist,xd,yd) //testing only506 // solidAngle[ii][jj] /= DetEffCorr(lambda,dtdist,xd,yd) //testing only 504 507 #endif 505 508 endif 506 509 507 // large angle transmission correction is <= 1 and will "bump up" the highest angles 510 // large angle transmission calculation is <= 1 and will "bump down" the measured value at the highest angles 511 // so divide here to get the correct answer 508 512 if(doTrans) 509 513 … … 520 524 521 525 data[ii][jj] /= LargeAngleTransmissionCorr(trans,dtdist,xd,yd) //moved from 1D avg SRK 11/2007 522 solidAngle[ii][jj]= LargeAngleTransmissionCorr(trans,dtdist,xd,yd) //testing only526 // solidAngle[ii][jj] /= LargeAngleTransmissionCorr(trans,dtdist,xd,yd) //testing only 523 527 endif 524 528
Note: See TracChangeset
for help on using the changeset viewer.