Load and save POLAR scalers as double

Author: Arthur Glowacki

src/io/file/hdf5_io.h

@@ -1038,9 +1038,9 @@ class DLL_EXPORT HDF5_IO
     }
 
     //-----------------------------------------------------------------------------
-
-    template<typename T_real>
-    bool load_spectra_vol_polar_energy_scan(std::string path, std::string filename, size_t detector_num, data_struct::Spectra_Volume<T_real>* spec_vol, data_struct::Scan_Info<T_real> &scan_info, bool logerr = true)
+    // load spectra volume as float (T_real) but scan_info as double (T_real2)
+    template<typename T_real, typename T_real2>
+    bool load_spectra_vol_polar_energy_scan(std::string path, std::string filename, size_t detector_num, data_struct::Spectra_Volume<T_real>* spec_vol, data_struct::Scan_Info<T_real2> &scan_info, bool logerr = true)
     {
         std::stack<std::pair<hid_t, H5_OBJECTS> > close_map;
         hid_t    file_id, dset_id, space_id, x_dset_id, y_dset_id, x_space_id, y_space_id;
@@ -1256,13 +1256,13 @@ class DLL_EXPORT HDF5_IO
             scan_info.meta_info.y_axis.resize(scan_info.meta_info.requested_rows);
             spec_vol->resize_and_zero(scan_info.meta_info.requested_rows, scan_info.meta_info.requested_cols, dims3[2]);
 
-            struct data_struct::Scaler_Map<T_real> energy_map;
+            struct data_struct::Scaler_Map<T_real2> energy_map;
             energy_map.name = "Energy";
             energy_map.unit = "";
             energy_map.time_normalized = false;
             energy_map.values.resize(scan_info.meta_info.requested_rows, scan_info.meta_info.requested_cols);
 
-            struct data_struct::Scaler_Map<T_real> i0_map;
+            struct data_struct::Scaler_Map<T_real2> i0_map;
             i0_map.name = "I0";
             i0_map.unit = "counts";
             i0_map.time_normalized = false;
@@ -1273,8 +1273,8 @@ class DLL_EXPORT HDF5_IO
             data_struct::ArrayTr<T_real> icr_array(dims3[0]);
             data_struct::ArrayTr<T_real> ocr_array(dims3[0]);
 
-            data_struct::ArrayTr<T_real> i0_array(dims3[0]);
-            data_struct::ArrayTr<T_real> energy_array(dims3[0]);
+            data_struct::ArrayTr<T_real2> i0_array(dims3[0]);
+            data_struct::ArrayTr<T_real2> energy_array(dims3[0]);
 
             err = _read_h5d<T_real>(elt_dset_id, H5S_ALL, elt_space_id, H5P_DEFAULT, elt_array.data());
             if(err < 0)
@@ -1301,13 +1301,13 @@ class DLL_EXPORT HDF5_IO
                 ocr_array.setOnes();
             }
 
-            err = _read_h5d<T_real>(energy_dset_id, H5S_ALL, energy_space_id, H5P_DEFAULT, energy_array.data());
+            err = _read_h5d<T_real2>(energy_dset_id, H5S_ALL, energy_space_id, H5P_DEFAULT, energy_array.data());
             if(err < 0)
             {
                 logE<< "Could not read energy data, setting it all to 1's\n";
                 energy_array.setOnes();
             }
-            err = _read_h5d<T_real>(i0_dset_id, H5S_ALL, i0_space_id, H5P_DEFAULT, i0_array.data());
+            err = _read_h5d<T_real2>(i0_dset_id, H5S_ALL, i0_space_id, H5P_DEFAULT, i0_array.data());
             if(err < 0)
             {
                 logE<< "Could not read i0 data, setting it all to 1's\n";
@@ -6980,6 +6980,10 @@ class DLL_EXPORT HDF5_IO
 
             _save_scalers(maps_grp_id, &(scan_info->scaler_maps), params_override->us_amp_sens_num, params_override->us_amp_sens_unit, params_override->ds_amp_sens_num, params_override->ds_amp_sens_unit);
         }
+        else
+        {
+            _save_scalers(maps_grp_id, &(scan_info->scaler_maps), (T_real)0.0, "0.0", (T_real)0.0, "0.0");
+        }
         _close_h5_objects(_global_close_map);
 
         end = std::chrono::system_clock::now();

src/io/file/hl_file_io.h

@@ -848,13 +848,14 @@ DLL_EXPORT bool load_spectra_volume(std::string dataset_directory,
     if (dataset_file[dlen - 4] == '.' && dataset_file[dlen - 3] == 'h' && dataset_file[dlen - 2] == 'd' && dataset_file[dlen - 1] == 'f')
     {
         // try to load polar hdf master file
-        data_struct::Scan_Info<T_real> scan_info;
+        data_struct::Scan_Info<double> scan_info;
         if(true == io::file::HDF5_IO::inst()->load_spectra_vol_polar_energy_scan(dataset_directory, dataset_file, detector_num, spectra_volume, scan_info))
         {
             scan_type = STR_SCAN_TYPE_POLAR_XANES;
             if(io::file::HDF5_IO::inst()->start_save_seq(true))
             {
-                io::file::HDF5_IO::inst()->save_scan_scalers(detector_num, &scan_info, params_override);
+                data_struct::Params_Override<double>* null_po = nullptr;
+                io::file::HDF5_IO::inst()->save_scan_scalers(detector_num, &scan_info, null_po);
                 return true;
             }
             else