Ticket #3193 (closed: fixed)

Opened 9 years ago

Last modified 5 years ago

Anvred corrections for single crystal

Reported by: Vickie Lynch Owned by: Vickie Lynch
Priority: major Milestone: Iteration 29
Component: Mantid Keywords:
Cc: Blocked By:
Blocking: Tester: Michael Reuter

Description

New algorithm, AnvredCorrections, to do spherical absorption corrections. The pixel and wavelength corrections will be optional if the data has not been divided by vanadium. Following A.J.Schultz's anvred, the weight factors should be:

  • sin2(theta) / (lamda4 * spec * eff * trans)
  • where theta = scattering_angle/2
  • lamda = wavelength (in angstroms?)
  • spec = incident spectrum correction
  • eff = pixel efficiency
  • trans = absorption correction

Change History

comment:1 Changed 9 years ago by Vickie Lynch

  • Status changed from new to accepted
  • Owner set to Vickie Lynch

comment:2 Changed 9 years ago by Vickie Lynch

(In [12604]) Refs #3193 Anvred correction are implemented for historgrams, but not finished for events; have not compared results with ISAW

comment:3 Changed 9 years ago by Vickie Lynch

(In [12610]) Refs #3193 Same anvred corrections for histograms and events

comment:4 Changed 9 years ago by Vickie Lynch

(In [12633]) Refs #3193 Power correct for weights; interpolation removed; checkbox renamed to OnlySphericalAbsorption

comment:5 Changed 9 years ago by Vickie Lynch

(In [12634]) Refs #3193 BuildLamdaWeights arguments changed

comment:6 Changed 9 years ago by Vickie Lynch

(In [12649]) Refs #3193 add detector id to EventLists and do corrections in place for events

comment:7 Changed 9 years ago by Vickie Lynch

Matched these corrections from IsawEV with this algorithm. Data was TOPAZ_2511, power=3, Radius=0.17, smu scattering =1.302, amu absorption=1.686 for detector id = 3245368.

max_grid_id = 58 Built APPROXIMATE weights with wl power 3.0 Made Q mapper in 1477.4 ms LOADED 4194304 EVENTS IN 58.4 ms LOADED 8388608 EVENTS IN 46.5 ms lamda = 1.9289504weight = 0.014794302 lamda = 1.5949249weight = 0.02419433 lamda = 3.3021934weight = 0.0039280728 lamda = 2.9713235weight = 0.005017997 LOADED 12582912 EVENTS IN 1131.4 ms lamda = 0.91291904weight = 0.11080247 lamda = 2.0927753weight = 0.011891115 lamda = 1.443294weight = 0.031490322 lamda = 0.72511417weight = 0.21441211 lamda = 3.006553weight = 0.004906128 LOADED 16777216 EVENTS IN 2674.0 ms lamda = 1.0574653weight = 0.07451244 lamda = 0.9779252weight = 0.09283592 lamda = 1.8250223weight = 0.016975451 lamda = 0.9119101weight = 0.110777274 lamda = 2.5400462weight = 0.0073058195 lamda = 1.6446884weight = 0.022292988 LOADED 20971520 EVENTS IN 112.6 ms lamda = 1.8626993weight = 0.016036544 lamda = 2.2585537weight = 0.009883552 lamda = 1.1580445weight = 0.058009826 lamda = 2.3457797weight = 0.008955821 lamda = 0.8593342weight = 0.13432859 lamda = 1.8878386weight = 0.015616583 LOADED 25165824 EVENTS IN 124.3 ms LOADED 25552337 EVENTS IN 11.2 ms lamda = 0.9953789weight = 0.08766594 lamda = 2.80196weight = 0.0057731555 lamda = 2.1811604weight = 0.0106836995 lamda = 2.0104656weight = 0.013128343 lamda = 2.1067512weight = 0.0117581105 lamda = 0.23956504weight = 5.8878136 lamda = 2.667289weight = 0.006539576 lamda = 2.8211741weight = 0.0056747673 lamda = 2.5578222weight = 0.0072482186 lamda = 3.1017652weight = 0.004538702 lamda = 0.81347775weight = 0.15362896 lamda = 2.5024123weight = 0.0075995494 lamda = 3.170335weight = 0.0043077185 lamda = 2.3631904weight = 0.008763356 lamda = 1.0908056weight = 0.06710697

comment:8 Changed 9 years ago by Vickie Lynch

(In [12663]) Refs #3193 Test for AnvredCorrection

comment:9 Changed 9 years ago by Vickie Lynch

  • Status changed from accepted to verify
  • Resolution set to fixed

Agrees with corrections in IsawEV. Test and wiki page added.

comment:10 Changed 9 years ago by Michael Reuter

  • Status changed from verify to verifying
  • Tester set to Michael Reuter

comment:11 Changed 9 years ago by Michael Reuter

  • Status changed from verifying to closed

Everything is within acceptable limits.

comment:12 Changed 5 years ago by Stuart Campbell

This ticket has been transferred to github issue 4040

Note: See TracTickets for help on using tickets.