mx/Jungfraujoch
2
1
Fork 0
Code Issues Pull Requests Actions 1 Packages Projects Releases 67 Wiki Activity

BraggIntegrate2D: Lorentz-Polarization correction is not directly applied

Browse Source
This commit is contained in:
leonarski_f
2026-02-17 12:16:58 +01:00
parent 411276415e
commit 175e938616
1 changed files with 1 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files
image_analysis/bragg_integration/BraggIntegrate2D.cpp
+1 -12
View File
@@ -63,8 +63,6 @@ void IntegrateReflection(Reflection &r, const T *image, size_t xpixel, size_t yp
}
}
#include <iostream>
template<class T>
std::vector<Reflection> IntegrateInternal(const DiffractionExperiment &experiment,
const CompressedImage &image,
@@ -88,22 +86,13 @@ std::vector<Reflection> IntegrateInternal(const DiffractionExperiment &experimen
const float r_2_sq = settings.GetR2() * settings.GetR2();
const float r_3_sq = settings.GetR3() * settings.GetR3();
const Coord S0 = geom.GetScatteringVector();
Coord m2 = Coord(0,0,0); // zero vector
if (experiment.GetGoniometer().has_value())
m2 = experiment.GetGoniometer()->GetAxis().Normalize();
for (int i = 0; i < npredicted; i++) {
auto r = predicted.at(i);
IntegrateReflection(r, ptr, image.GetWidth(), image.GetHeight(), special_value, saturation,
r_3, r_1_sq, r_2_sq, r_3_sq);
if (r.observed) {
if (experiment.GetPolarizationFactor())
r.rlp *= geom.CalcAzIntPolarizationCorr(r.predicted_x, r.predicted_y,
experiment.GetPolarizationFactor().value());
r.I *= r.rlp;
r.bkg *= r.rlp;
r.sigma *= r.rlp;
r.rlp *= geom.CalcAzIntPolarizationCorr(r.predicted_x, r.predicted_y, experiment.GetPolarizationFactor().value());
r.image_number = static_cast<float>(image_number);
ret.emplace_back(r);
}