diff --git a/reader/JFJochHDF5Reader.cpp b/reader/JFJochHDF5Reader.cpp
index 985859db..0a28b0de 100644
--- a/reader/JFJochHDF5Reader.cpp
+++ b/reader/JFJochHDF5Reader.cpp
@@ -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),