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

JFJochHDF5Reader: Read delta phi

Browse Source
This commit is contained in:
leonarski_f
2026-02-08 13:00:56 +01:00
parent 70343ec058
commit 666b52ffc1
1 changed files with 7 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
reader/JFJochHDF5Reader.cpp
View File

@@ -682,7 +682,7 @@ bool JFJochHDF5Reader::LoadImage_i(std::shared_ptr<JFJochReaderDataset> &dataset
auto bkg = source_file->ReadOptVector<float>(image_group_name + "/background_mean");
auto lp = source_file->ReadOptVector<float>(image_group_name + "/lp");
auto partiality = source_file->ReadOptVector<float>(image_group_name + "/partiality");
auto phi = source_file->ReadOptVector<float>(image_group_name + "/delta_phi");
if (h.size() != l.size() || h.size() != k.size() || h.size() != d.size()
|| h.size() != predicted_x.size() || h.size() != predicted_y.size()
|| h.size() != int_sum.size() || h.size() != int_err.size() || h.size() != bkg.size())
@@ -693,14 +693,18 @@ bool JFJochHDF5Reader::LoadImage_i(std::shared_ptr<JFJochReaderDataset> &dataset
if (lp.size() > i && lp[i] != 0.0f)
lp_val = 1.0f / lp[i];
float partiality_val = 0.0f;
if (partiality.size() > i && partiality[i] != 0.0f)
float partiality_val = -1.0f;
if (partiality.size() > i && partiality[i] >= 0.0f)
partiality_val = partiality[i];
float delta_phi_val = NAN;
if (phi.size() > i)
delta_phi_val = phi[i];
Reflection r{
.h = h.at(i),
.k = k.at(i),
.l = l.at(i),
.delta_phi = delta_phi_val,
.predicted_x = predicted_x.at(i),
.predicted_y = predicted_y.at(i),
.d = d.at(i),