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

BraggIntegrate: Add Lorentz factor (1/sin(2*theta)) for rotation measurment

Browse Source
This commit is contained in:
leonarski_f
2025-12-18 13:13:39 +01:00
parent 2971b20b87
commit 25f068e111
1 changed files with 10 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files
image_analysis/bragg_integration/BraggIntegrate2D.cpp
+10 -6
View File
@@ -91,14 +91,18 @@ std::vector<Reflection> IntegrateInternal(const DiffractionExperiment &experimen
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()) {
float pol = geom.CalcAzIntPolarizationCorr(r.predicted_x, r.predicted_y,
float corr = 1.0;
if (experiment.GetPolarizationFactor())
corr /= geom.CalcAzIntPolarizationCorr(r.predicted_x, r.predicted_y,
experiment.GetPolarizationFactor().value());
r.I /= pol;
r.bkg /= pol;
r.sigma /= pol;
}
if (experiment.GetGoniometer().has_value()) {
float two_theta = geom.TwoTheta_rad(r.predicted_x, r.predicted_y);
corr *= std::sin(two_theta);
}
r.I *= corr;
r.bkg *= corr;
r.sigma *= corr;
r.image_number = static_cast<float>(image_number);
ret.emplace_back(r);
}